UNIVERZITET U BEOGRADU HEMIJSKI FAKULTET Marko R. Perić TEORIJSKO PROUČAVANJE STRUKTURNIH I ELEKTRONSKIH UTICAJA NA MAGNETNE OSOBINE KOMPLEKSA PRELAZNIH METALA Doktorska Disertacija Beograd, 2012 UNIVERSITY OF BELGRADE FACULTY OF CHEMISTRY Marko R. Perić THEORETICAL INVESTIGATION OF STRUCTURAL AND ELECTRONIC INFLUENCES ON THE MAGNETIC PROPERTIES OF TRANSITION METAL COMPLEXES Ph.D. Thesis Belgrade, 2012 Mentor __________________________________ Docent dr Maja Gruden-Pavlović, Univerzitet u Beogradu Hemijski fakultet Mentor __________________________________ Viši naučni saradnik dr Sonja Grubišić, Univerzitet u Beogradu IHTM Članovi komisije __________________________________ Redovni profesor dr Katarina Anđelković, Univerzitet u Beogradu Hemijski fakultet __________________________________ Naučni saradnik dr Matija Zlatar, Univerzitet u Beogradu IHTM Datum odbrane ______________2012. Zahvalnica Doktorska disertacija je urađena na Katedri za Opštu i neorgansku hemiju Hemijskog fakulteta Univerziteta u Beogradu. Deo istraživanja urađen je u saradnji sa istraživačkom grupom iz Kantabrije u Španiji i ovom prilikom želim da se zahvalim dr Pablo García Fernández-u za pomoć i saradnju. Zahvaljujem se profesoru Claude Daul-u i Pawel Bednarek-u sa Univerziteta u Friburu, Švajcarska, za korišćenje kompjuterskih resursa i programa. Iskreno se zahvaljujem profesoru Svetozaru Niketiću mom prvom mentoru, kod koga sam uradio diplomski rad. Iskreno se zahvaljujem dr Maji Gruden-Pavlović, dr Sonji Grubišić i dr Matiji Zlatar, koji su predložili temu ove disertacije i rukovodili njenom izradom. Njihova nesebična pomoć i strpljenje bili su presudni za tok i izradu ove disertacije. Veliku zahvalnost dugujem dr Katarini Anđelković, redovnom profesoru, na pomoći i korisnim sugestijama tokom izrade i pisanja disertacije. Posebno se zahvaljujem mojim dragim kolegama Ljubici Anđelković i dr Aleksandru S. Nikoliću za bezrezervnu pomoć i razumevanje. Najveću zahvalnost dugujem mojoj porodici i najdražim prijateljima, za razumevanje i podršku, koju su mi pružali proteklih godina tokom izrade ove disertacije. TEORIJSKO PROUČAVANJE STRUKTURNIH I ELEKTRONSKIH UTICAJA NA MAGNETNE OSOBINE KOMPLEKSA PRELAZNIH METALA REZIME U okviru ove doktorske disertacije analizirane su magnetne osobine mononuklearnih i binuklearnih kompleksa prve serije prelaznih metala LF-DFT i DFT metodama. Poznavanje magnetnih osobina kompleksnih jedinjenja omogućuje predviđanje i dizajniranje molekula sa traženim svojstvima koja se mogu iskoristiti za skladištenje informacija na kvantnom nivou. Molekulski magneti, tj. mononuklearni sistemi sa velikom vrednošću magnetne anizotropije nalaze se u žiži interesovanja s obzirom na njihovu potencijalnu primenu u kvantnom računarstvu. Analizom faktora koji određuju magnetna svojstva u mononuklearnim kompleksima pronašli smo već sintetisani kompleksni jon [NiCl3(Hdabco)2]+ (dabco=1,4-diazabicyclo[2.2.2]-octane) sa veoma velikom vrednosti aksijalnog parametra cepanja, D. Osnovni uzrok ovako velikih vrednosti predstavlja interakcija između dx2-y2 i dxy Magnetna kuplovanja između paramagnetnih centara u binuklearnim kompleksima analizirana su BS (Broken Symmetry) DFT metodom. U svim slučajevima jačina magnetnih interakcija zavisi od elektronske strukture metala, kao i od prirode i orijentacije mostnih liganada. orbitala. Efekti Jahn-Teller-ove distorzije, koji deluju nasuprot spin-orbitalnom kuplovanju, umanjeni su prisustvom voluminoznih dabco liganada, koji se nalaze u aksijalnom položaju. Usled sternog zaklanjanja aksijalnih liganada atomi hlora u ekvatorijalnoj ravni, slabo se pomeraju jedan u odnosu na drugi. Distorzije u regiji mosta dovode do znatnog slabljenja antiferomagnetnih interakcija. Za analizu uticaja elektronske konfiguracije na tip kuplovanja izabrani su binuklearni kompleksi prve serije prelaznih metala sa planarnim i rigidnim oksalatnim mostovima. Najveće vrednosti konstanti kuplovanja izračunate su za komplekse nikla (II) i bakra (II), usled σ interakcija između d orbitala metala i p orbitala mostnih liganada. Nadalje, porast broja nesparenih elektrona po metalu dovodi do pojačanja antiferomagnetnih interakcija, ali istovremo i do pojačanja feromagnetnih interakcija između ortogonalnih orbitala. Antagonistički efekti se međusobno poništavaju, što kao rezultat daje male vrednosti konstanti kuplovanja. Analizom magnetnih interakcija koje se prenose posredstvom vodoničnih veza, dokazano je da jačina kuplovanja između monomernih jedinica zavisi od usmerenosti i prostornog rasporeda vodoničnih mostova. U većini ispitivanih slučajeva Jahn-Teller-ova distorzija predstavlja osnovni uzrok odsustva kuplovanja. Ključne reči: Magnetna anizotropija, Magnetna kuplovanja, DFT, LF-DFT, Binuklearni kompleksi Naučna oblast: Hemija Uža naučna oblast: Neorganska hemija UDK broj: 546.733 THEORETICAL INVESTIGATION OF STRUCTURAL AND ELECTRONIC INFLUENCES ON THE MAGNETIC PROPERTIES OF TRANSITION METAL COMPLEXES SUMMARY Magnetic properties of mononuclear and binuclear transition metal complexes were analyzed by the means of LF-DFT and DFT methods. The understanding of magnetic properties of metal complexes allow us to predict and design molecules with required characteristics, which can be used to store the informations on quantum level. Molecular magnets, i.e. mononuclear systems with large magnetic anisotropy value are in the focus of interest due to their potential applications in quantum computing. Through the analysis of the factors that determine the magnetic properties of mononuclear complexes, we found complex ion [NiCl3(Hdabco)2]+ (dabco=1,4-diazabicyclo[2.2.2]- octane) with very high values of axial splitting parameter, D. The main cause for such a large value lies in interactions between dx2-y2 and dxy Magnetic couplings between paramagnetic centers in binuclear complexes were analyzed using Broken Symmetry (BS) DFT method. In all investigated cases the strength of magnetic interactions depends upon the electronic structure and the structure of bridging ligands. orbitals. The effects of Jahn-Teller distortion, which act opposed to spin-orbit coupling, are suppressed by the presence of voluminous dabco ligands, which are placed in axial position. Due to the steric hindrance of axial ligands, the chlorine atoms in equatorial plane can not move easily. The distortions in the bridge region lead to a significant weakening of antiferromagnetic interactions. In order to analyze the influence of electronic factors on exchange couplings, binuclear complexes with planar and rigid oxalate bridges were chosen. The highest values of coupling constants were calculated for complexes of nickel (II) and copper (II) due to the σ interactions between d metal orbitals and p bridging ligand orbitals. Furthermore, the rising number of unpaired electrons per metal center leads to the increase of antiferomagnetic and feromagnetic interactions at the same time. The antagonistic effects cancel each other, causing very small values of coupling constants. Magnetic interactions that are transmitted through the hydrogen bonds strongly depends on the orientation and arrangement of hydrogen bridges. In most examined cases Jahn-Teller distortion is the main cause of the absence of coupling . Keywords: Magnetic anisotropy, Magnetic couplings, DFT, LF-DFT, Binuclear complexes Area of science: Chemistry Sub-area of science: Inorganic Chemistry UDC number: 546.733 Sadržaj 1 Uvod ........................................................................................................................... 1 2 Teorija Funkcionala Gustine (Density Functional Theory-DFT) ........................ 2 2.1 Kohn-Sham DFT ................................................................................................. 5 2.2 LDA funkcional ................................................................................................... 7 2.3 “Generalized Gradient” Aproksimacija (GGA funkcionali) ............................... 9 2.4 Hibridni funkcionali; B3LYP funkcional .......................................................... 10 3 Uvod u magnetizam ................................................................................................ 12 4 Analiza spinskog hamiltonijana ............................................................................ 15 4.1 Magnetna anizotropija mononuklearnih kompleksa prelaznih metala; D vrednost ........................................................................................................................... 15 4.1.1 Spin dipolarni efekat ..................................................................................... 23 4.1.2 Uticaj Jahn-Teller-ovog efekta na magnetnu anizotropiju ........................... 24 4.1.3 Uticaj elektronske konfiguracije na magnetnu anizotropiju trigonalno bipiramidalnih sistema ................................................................................................. 28 4.2 Izračunavanje D vrednosti u mononuklearnim kompleksima prelaznih metala LF- DFT metodom .................................................................................................... 30 4.3 Magnetna kuplovanja u binuklearnim kompleksima prelaznih metala ............. 32 4.4 DFT proračuni magnetnih kuplovanja u binuklearnim kompleksima prelaznih metala ................................................................................................................. 42 5 Magnetna anizotropija u mononuklearnim kompleksima .................................. 46 5.1 Računski metod ................................................................................................. 47 5.2 Rezultati i diskusija ........................................................................................... 48 5.2.1 Analiza magnetne anizotropije u model sistemima ...................................... 48 5.2.2 Analiza magnetne anizotropije za [NiCl3(Hdabco)2]+ kompleksni jon ....... 51 5.3 Zaključak ........................................................................................................... 53 6 DFT analiza magnetnih interakcija u binuklearnim kompleksima mangana .. 54 6.1 Računski metod ................................................................................................. 55 6.2 Rezultati i diskusija ........................................................................................... 56 6.3 Zaključak ........................................................................................................... 60 7 DFT analiza magnetnih interakcija u binuklearnim oksalatnim kompleksima prve serije prelaznih metala ........................................................................................... 61 7.1 Računski metod ................................................................................................. 62 7.2 Rezultati i diskusija ........................................................................................... 63 7.2.1 Određivanje linearne zavisnosti između konstante kuplovanja i (ε2-ε1)2 i S2 parametara ................................................................................................................... 71 7.3 Zaključak ........................................................................................................... 72 8 Prenošenje magnetnih kuplovanja posredstvom ne-kovalentnih interakcija ... 73 8.1 Računski metod ................................................................................................. 74 8.2 Rezultati i diskusija ........................................................................................... 74 8.2.1 Kuplovanja posredstvom F-H…. F vodoničnih veza ...................................... 75 8.2.2 Kuplovanja posredstvom O-H….F and O-H... ..................O vodoničnih veza 80 8.2.3 Kuplovanja posredstvom Br....H-A-H....Br (A=O,N) i Br.... Br veza ................ 84 8.3 Zaključak ........................................................................................................... 89 9 ................................................................................................................. Zaključak 90 10 Korišćeni programi ................................................................................................. 92 11 Literatura...................................................................................................................93 12 Prilog ...................................................................................................................... 101 1 1 Uvod Nagli razvoj računarskih tehnologija, u proteklih nekoliko godina, omogućio je rešavanje veoma složenih problema primenom kompjuterskih simulacija, ne samo u hemiji već i u drugim oblastima kao što su biologija, medicina i fizika. Upotreba računarskih metoda u hemiji, neophodna je u cilju dobijanja jasnijeg uvida u svojstva ispitivanih sistema, kao i radi predviđanja osobina nepoznatih sistema. Analiza porekla magnetnih interakcija u mononuklearnim i binuklearnim kompleksima d- i f-metala, postala je predmet opsežnog istraživanja u proteklih par godina, a teorijski aspekti su detaljno opisani u mnogim knjigama i radovima1,2,3,4,5,6,7,8,9,10,11 . Bez obzira na postojeći napredak eksperimentalnih tehnika u hemiji, ispitivanje magnetnih osobina kompleksa d- i f-metala, u današnje vreme je značajno olakšano primenom različitih kvantno mehaničkih metoda, a posebno KS-DFT-a (Kohn-Sham Density Functional Theory, videti odeljak 2.1), pomoću kojeg se u razumnom vremenskom roku mogu dobiti rezultati zadovoljavajuće preciznosti. Prvi deo ove doktorske disertacije obuhvatao je ispitivanje magnetnih svojstava mononuklearnih kompleksa prelaznih metala, dok se drugi deo bazirao na magnetnim osobinama binuklearnih kompleksa, kao i međusobnim magnetnim interakcijama između mononuklearnih gradivnih jedinica, korišćenjem DFT i LF-DFT (Ligand Field-DFT) teorijskih metoda. Kako molekulski magnetizam povezuje nauku o materijalima sa biološkim procesima u kojima učestvuju metali kao aktivni centri enzima ovde su se, između ostalog, razmatrale i magnetne interakcije u metaloenzimima bakra. Osnovni cilj ove teze, bio je da se predvide magnetna svojstva materijala potrebna za rad odgovarajućih elektronskih uređaja, kao što su prekidači i uređaji za skladištenje informacija, analizom na elektronskom, mikroskopskom nivou, kao i dizajn novih molekulskih sistema sa željenim magnetnim osobinama. Dalje, kao krajnji cilj bio je dizajn molekulskih magneta i molekulskih sistema sa kontrolisanim magnetnim osobinama, tj. sa mogućnošću ‘’uključivanja /isključivanja’’ magnetnog polja po određenoj potrebi. 2 2 Teorija Funkcionala Gustine (Density Functional Theory- DFT) Pored velikog broja kvantno mehaničkih metoda koji za opisivanje sistema koriste talasnu funkciju, teorija funkcionala gustine, DFT, kao jedan od najviše korišćenih metoda kvantne hemije, koristi raspodelu elektronske gustine za izračunavanje merljivih veličina. Jedna je od osnovnih prednosti DFT-a u odnosu na ostale kvanto mehaničke metode, što se rezultati dobijaju u kratkom vremenskom intervalu uz zadovoljavajuću tačnost i preciznost. Teorijski aspekti DFT-a detaljno su opisani u nekoliko knjiga i revijalnih radova, dok će ovde biti date samo neke osnovne napomene12,13,14,15,16,17,18,19 Razmatranje je neophodno započeti sa vremenski nezavisnom Schrödinger-ovom jednačinom za sistem od M elektrona, koji se kreću pod privlačnim potencijalom N jezgara (Jednačina 2.1). . ( ) ( )MM rrrErrrH  ,....,,....,ˆ 2121 Ψ=Ψ Jednačina 2.1 Hamiltonijan za dati sistem sastoji se od tri člana: kinetičke energije, privlačnih interakcija jezgra i elektrona i međuelektronskih interakcija (Jednačine 2.2-2.5). ∑ ∑ ∑∑ ∑ < = == ∇= ++= M ji ij ee M i M i N A iA A extne M i i eene r V r ZvV T VVTH i   1ˆ ˆˆ ˆ ˆˆˆˆ 2 Jednačine 2.2-2.5 3 Oznake Tˆ , neVˆ i eeVˆ predstavljaju kvantno mehaničke operatore za kinetičku energiju, potencijalnu energiju interakcija jezgra i elektrona i potencijalnu energiju međuelektronskih interakcija. Operatori Tˆ i eeVˆ su univerzalni, budući da su isti za bilo koji sistem. Bilo da je ispitivani sistem atom, molekul, makromolekul ili kristal, Šredinger-ova jednačina zavisi samo od. neVˆ . U DFT termonologiji extvˆ predstavlja eksterni potencijal, odnosno polje koje generišu jezgra. Oznaka ZA iArˆ, predstavlja naelektrisanje jezgra, dok i ijrˆ označavaju rastojanje između jezgra i elektrona i međuelektronsko rastojanje. Po prvoj Hohenberg-Kohn-ovoj teoremi20,21 neVˆ predstavlja jedinstveni funkcional elektronske gustine do aditivne konstante, odnosno elektronska gustina određuje eksterni potencijal, a samim tim jedinstveno određuje hamiltonijan, odnosno energiju ispitivanog sistema. Nadalje, Šredinger-ova jednačina se može prikazati u drugačijem obliku (Jednačina 2.6.). ( ) k kee k i iextne VTEvV Ψ Ψ+ −== ∑ ˆˆˆ Jednačina 2.6 Izraz 2.6. pokazuje da bilo koje svojstveno stanje Ψk određuje operator potencijalne energije, do aditivne konstante tj. do odgovarajuće svojstvene vrednosti Ek Prva Hohenberg-Kohn-ova teorema se može sumirati na sledeći način: . Kao posledica navedenog, elektronska gustina određuje kako talasnu funkciju osnovnog stanja, tako i talasne funkcije ekscitovanih stanja. ( ) ( ) ( )rrvr ext  Ψ↔↔ρ Jednačina 2.7 Talasna funkcija Ψ predstavlja funkcional eksternog potencijala i elektronske gustine. Ovakva formulacija predstavlja osnov Hohenberg-Kohn-ove teoreme, odnosno sva svojstva 4 sistema predstavljaju funkcional elektronske gustine osnovnog stanja. Očekivana vrednost bilo koje obzervable, O, je funkcional ρ(r): [ ] [ ] [ ]ρρρ ΨΨ== OOO ˆ Jednačina 2.8 Kinetička energija, T[ρ] i energija međuelektronskih interakcija, Vee [ ] [ ] [ ] [ ] [ ]ρρρρρ Ψ++Ψ=ΨΨ= eenev VVTHE ˆˆˆˆ [ρ], takođe su funkcionali gustine. Totalna energija sistema, koja je takođe funkcional gustine, može se izraziti na sledeći način: Jednačina 2.9 Drugom Hohenberg-Kohn-ovom teoremom20 [ ] 0' EEv ≥ρ uspostavlja se veza varijacionog principa i elektronske gustine: Jednačina 2.10 , gde ρ' predstavlja neku inicijalno pretpostavljenu gustinu. Energija osnovnog stanja i gustina odgovaraju minimumu nekog funkcionala, Ev Kako su kinetička energija i međuelektronske interakcije opisani univerzalnim operatorom, totalna energija se može napisati u drugačijem obliku. [ρ], pri čemu elektronska gustina mora sadržati odgovarajući broj elektrona. Ovakva formulacija sledi iz obe Hohenberg- Kohn-ove teoreme. [ ] [ ] [ ] [ ] [ ] [ ] [ ] ( ) rdrvFVFVVTE extneeene  ∫+=+=++= ρρρρρρρρ ˆ Jednačina 2.11 Oznaka F[ρ] predstavlja univerzalni funkcional nezavisan od vext . Ukoliko bi F[ρ] funkcional bio poznat, primenom variacionog principa, bilo bi moguće dobiti najtačniju moguću elektronsku gustinu. 5 2.1 Kohn-Sham DFT Nepoznavanje univerzalnog funkcionala, F[ρ], predstavljalo je jedan od osnovnih problema različitih DFT metoda. Većina modela, poput Thomas-Fermi-jevog22, nisu bili u stanju da reprodukuju kinetičku energiju, T[ρ], koja prema Virijalnoj teoremi predstavlja polovinu totalne energije sistema. Tek 1965 godine Kohn i Sham su predložili indirektno ali veoma originalno rešenje ovog problema, kojim je omogućeno izračunavanje ovako velikog doprinosa kinetičke energije u funkcionalu F[ρ]23 svˆ . Njihova ideja pretvorila je DFT u veoma praktično sredstvo za komplikovane i veoma rigorozne proračune, tako da se u današnje vreme većina DFT pristupa u kvantnoj hemiji zasniva na Kohn-Sham-ovoj jednoelektronskoj formulaciji (KS-DFT). U okviru svog modela Kohn i Sham su uveli orbitale (spin-orbitale), na takav način da kinetička energija može biti veoma lako izračunata. Osnovu Kohn-Sham-ovog pristupa predstavlja fiktivni sistem od M elektrona koji međusobno ne interaguju, kreću se u eksternom potencijalu i imaju istu gustinu kao realni sistem. Ovakav fiktivni sistem može se opisati talasnom funkcijom, koja se prikazuje jednom determinantom. Kohn-Sham-ov Hamiltonijan predstavlja sumu jednoelektronskih Hamiltonijana. ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) [ ] ∑ ∑∑ ∑∑ ∇−=       +∇−== = =Ψ == = N i iis i isi i ss iiis Ns exact N i s s T rvihH eh N rrr φφρ φφ φφφ ρφρ 2 2 21 2 1 2 1 2 1ˆˆ 111ˆ ...21   Jednačine 2.12-2.15 Oznaka s u indeksu označava sistem čestica koje međusobno ne interaguju. Ako se uzme u obzir da klasične Coulomb-ove interakcije elektrona, J, čine najveći deo elektron-elektron interakcija, energija realnog sistema se može prikazati jednačinom 2.16: 6 [ ] [ ] [ ] [ ] [ ] [ ] ( ) ( ) 21 21 21 2 1 rdrd rr rrJ EJVTE xcnes    ∫ −= +++= ρρ ρ ρρρρρ Jednačina 2.16 Oznaka Exc[ρ] označava exchange-correlation (XC) funkcional. Poređenjem izraza 2.11 i 2.16 dobija se jednačina 2.17 na osnovu koje se može videti koji je zaista smisao EXC iz KS teorije. [ ] [ ] [ ]( ) [ ] [ ]( )ρρρρρ JVTTE eesxc −+−= Jednačina 2.17 Potrebno je napomenuti da Exc Pored energije izmena i korelacije, ovaj funkcional takođe sadrži korekciju za kinetičku energiju. Uzimajući u obzir prethodne izraze, i primenom varijacione teorije Kohn-Sham- ove jednačine se mogu napisati na sledeći način: energija predstavlja funkcional u kojem je sadržano sve ono što nije poznato, odnosno “mesto gde se odlaže sve ono sa čime ne znamo da rukujemo”. ( ) ( ) iiixc ervrdrr r φφρ =    + − +∇− ∫    ' ' ' 2 1 2 Jednačina 2.18 Oznaka vxc predstavlja potencijal s obzirom na EXC energiju i izražava se kao prvi izvod EXC energije po elektronskoj gustini (dEXC KS-DFT pristup je u principu tačan, ali E /dρ). Rešavanjem ovakvih jednačina dobijaju se Kohn-Sham-ove (KS) orbitale. Za egzaktni funkcional, a samim tim i egzaktni eksterni potencijal, KS orbitale daju tačnu elektronsku gustinu i energiju osnovnog stanja. xc član ostaje eksplicitno nedefinisan. Iako je taj član u izrazu za energiju (Jednačina 2.16) relativno mali on se ne može zanemariti, pa je iz tog razloga neophodno uvesti aproksimacije o kojima će više reči biti u narednim odeljcima. 7 Na kraju ovog odeljka, a pre nego što se uđe u detaljniju diskusiju o različitim aspektima Kohn-Sham-ovog pristupa neophodno je sumirati osnovne korake ove procedure: 1. Na početku se definiše referentni system od M elektrona koje međusobno ne interaguju i čije se egzaktno osnovno stanje može predstaviti jednom Slater-ovom determinantom čija elektronska gustina ρs 2. Orbitale koje formiraju Slater-ovu determinantu predstavljaju rešenja M jedno- elektronskih jednačina (2.12-2.16). Ovim je omogućeno određivanje kinetičke energije neinteragujućih čestica T izjednačena sa gustinom realnog sistema interagujućih čestica. s 3. Energija sistema interagujućih elektrona rastavljena je na kinetičku energiju neinteragujućeg sistema, T . s, energiju interakcija elektrona sa jezgrima, Vne, klasičnu elektrostatičku energiju odbijanja elektrona, J, i ostatak, EXC, koji se satoji od kvantno mehaničkih doprinosa potencijalnoj energiji (self-interaction korekcija, izmene i korelacije) i deo prave kinetičke energije koji nije obuhvaćen sa Ts 4. Rešavanjem Kohn-Sham-ovih jednoelektronskih jednačina (Jednačina 2.18) metodom samo-usaglašenog polja (Self Consistent Field-SCF) . 24 5. Orbitale daju elektronsku gustinu, pa ubacivanjem date gustine u izraz za energiju dobija se egzaktna gustina osnovnog stanja, kao i egzaktna energija osnovnog stanja. Kako u svim realnim aplikacijama E dobijaju se Kohn- Sham-ove (KS) orbitale. Ovim se definiše neinteragujući sistem koji deli istu elektronsku gustinu kao ispitivani realni sistem. XC član nije poznat, potrebno je uvesti aproksimacije za XC potencijal, vXC . 2.2 LDA funkcional Kako je napomenuto u prethodnom odeljku, za praktičnu primenu DFT-a neophodne su aproksimacije za XC energije. Naravno, izbor funkcionala primarno zavisi od hemijske prirode sistema koji se ispituje. LDA funkcional (Local Density 8 Approximation) 25 predstavlja klasu aproksimacija za Exc koja zavisi jedino od vrednosti elektronske gustine u određenoj tački prostora. LDA se zasniva na principu homogenog elektronskog gasa koji ima istu gustinu kao sistem koji se posmatra. Ovakav pristup postiže se smeštanjem M elektrona u prostor zapremine V sa pozitivno naelektrisanim okruženjem. Za homogeni elektronski gas, Exc [ ] ( ) ( )drvrE LDAxcLDAxc ρρρ ∫= se može izraziti jednačinom 2.19. Jednačina 2.19 Energija Exc cxxc EEE += se može razložiti na dva člana: Jednačina 2.20 Prvi član u jednačini poznat je analitički26 i može se izraziti u funkciji od elektronske gustine27 [ ] ( )∫    −= drrE LDAx 3/4 3/13 4 3 ρ π ρ . Jednačina 2.21 Analitički izraz za korelacionu energiju, drugi član u izrazu 2.20, nije poznat sem za granične slučajeve beskonačno slabih i beskonačno jakih korelacija. U slučaju homogenog elektronskog gasa gustine ρ, beskonačno jake korelacije mogu se izraziti jednačinom 2.22. ( ) ( )( )DrCrBrA sssc +++= lnlnε Jednačina 2.22 Korelaciona energija u slučaju beskonačno slabih interakcija izražava se jednačinom 2.23.       ++= ... 2 1 2/3 10 ss c r g r g ε Jednačina 2.23 9 Oznaka rs ρ π 1 3 4 3 =sr predstavlja Wigner-Seitz-ov radius, koji je sa elektronskom gustinom povezan sledećim izrazom. Jednačina 2.24 Veoma precizne Monte Carlo simulacije za energiju homogenog elektronskog gasa izvršene su za nekoliko prelaznih vrednosti gustine, čime su se dobile veoma tačne vrednosti energija korelacija. Najpopularnije lokalne aproksimacije interpoliraju ove precizne vrednosti, dobijene simulacijama. Različiti pristupi, koristeći raličite analitičke obrasce za εc, generisali su nekoliko lokalnih aproksimacija za korelacioni funkcional od kojih je najpoznatija VWN (Vosko-Wilk-Nusair)25 Generalno govoreći LDA daje kvalitativno dobre rezultate pri optimizaciji geometrije iako su dužine veza nešto kraće u odnosu na eksperimentalne vrednosti. Takođe se dobijaju dobri rezultati za vibracione frekvence i gustine naelektrisanja. Međutim, LDA funkcional nije dobra aproksimacija za izračunavanja u sistemima sa slabim vezama ili za predviđanje termodinamičkih i magnetnih osobina. Da bi se dobili precizniji rezultati potrebno je koristiti druge tipove funkcionala. aproksimacija. 2.3 “Generalized Gradient” Aproksimacija (GGA funkcionali) S obzirom na činjenicu da KS elektronska gustina nije lokalizovana, čini se da čisto lokalne aproksimacije kao LDA nisu dovoljne. GGA funkcionali predstavljaju drugi tip aproksimacija, u kojima potencijali zavise, kako od vrednosti elektronske gustine u određenoj tački, tako i od vrednosti gradijenta elektronske gustine. GGA aproksimacije obično koriguju rezultate dobijene na LDA nivou, kao što su geometrije ispitivanih sistema, vibracione frekvence i gustine naelektrisanja. Druga karakteristika GGA funkcionala je da daju dobre rezultate za sisteme kod kojih postoje vodonične veze. Međutim, u nekim slučajevima izbor određenog GGA funkcionala ne dovodi do poboljšanja rezultata u 10 odnosu na rezultate dobijene LDA funkcionalom, pa je neophodno koristiti hibridne funkcionale. 2.4 Hibridni funkcionali; B3LYP funkcional Hibridni funkcionali predstavljaju klasu aproksimacija koje se dobijaju linearnom kombinacijom egzaktnih interakcija izmene iz Hartree-Fock teorije i Ex i Ec ( ) ( ) ( ) ( ) 2122 12 1 * 1 , * 1 2 1 drdrrr r rrE jij ji i HF x φφφφ∑∫∫= iz standardnih funkcionala. Interakcije koje se uzimaju iz Hartree-Fock-a izražene su na sledeći način (Jednačina 2.25). Jednačina 2.25 Za izračunavanje magnetnih osobina kompleksnih jedinjenja najviše se koristi B3LYP (Becke, 3-parametarski, Lee-Yang-Parr) hibridni funkcional (Jednačina 2.26 )28,29 . ( ) ( ) ( )LDAcGGAccLDAxGGAxxLDAxHFxLDAxcLYPBxc EEaEEaEEaEE −+−+−+= 03 Jednačina 2.26 Energije ExGGA i EcGGA predstavljaju Becke 88 funkcional28 i Lee-Yang-Parr29 korelacioni funkcional, EcLDA označava Vosko-Wilk-Nusair (VWN) aproksimaciju25 lokalne gustine korelacionom funkcionalu, dok parametri a0, ax i ac imaju vrednosti a0=0.20, ax=0.72, ac=0.81. B3LYP funkcional predstavlja najbolji izbor za izračunavanje magnetnih osobina koordinacionih jedinjenja, iako u nekim slučajevima LDA i GGA funkcionali daju sasvim zadovoljavajuće rezultate. Osnovni nedostatak hibridnih funkcionala predstavlja dužina proračuna. Odnos vremena i tačnosti proračuna od velike je važnosti za rešavanje složenih problema, pa se u današnje vreme sve više radi na razvoju novih funkcionala, koji su u mogućnosti da sistematski poboljšaju performanse DFT-a (Slika 2.1) 11 Slika 2.1 Šematski prikaz “Jakovovih lestvi” za DFT koje pokazuju odnos brzine i tačnosti proračuna sa različitim funkcionalima 12 3 Uvod u magnetizam Pod pojmom magnetizma podrazumeva se svojstvo materijala da interaguju sa spoljašnjim magnetnim poljem. Uopšteno govoreći, magnetno polje je posrednik uzajamnog delovanja magnetnih sila. Magnetna svojstva materijala zavise od njihove strukture, posebno od njihove elektronske konfiguracije, kao i od spoljašnjih faktora (pritisak, temperatura, …). Na osnovu pomenutog sve supstance se mogu podeliti u pet kategorija: 1. Paramagnetici su supstance sa nesparenim spinovima elektrona. Atomi ovih supstanci poseduju magnetne momente kada se nalaze van spoljašnjeg magnetnog polja, ali su pomenuti momenti haotično raspoređeni pa je ukupni magnetni momenat jednak nuli. Kada se takve supstance nađu u spoljašnjem magnetnom polju, magnetni momenti atoma se orijentišu u smeru polja, zbog čega supstanca biva “privučena” od strane spoljašnjeg magnetnog polja. Po prestanku delovanja polja supstanca nije namagnetisana i magnetni momenti su ponovo haotično raspoređeni. 2. Dijamagnetici su supstance sa potpuno popunjenim elektronskim ljuskama. Kada se supstance sa dijamagnetnim osobinama nalaze van magnetnog polja njihovi atomi ne poseduju magnetne momente. To znači da takvi marterijali ne stvaraju magnetno polje oko sebe. Kada se ovakve supstance nađu u spoljašnjem magnetnom polju, dolazi do indukovanja magnetnih momenata u atomima, ali su oni usmereni suprotno u odnosu na spoljašnje magnetno polje, što dovodi do njegovog slabljenja. Rezultat prethodno opisanog efekta je da supstanca biva “odbijena” od strane magnetnog polja. 3. Feromagnetici slično paramagneticima imaju nesparene elektrone. Međutim za razliku od paramagnetika, feromagnetici formiraju domene, unutar kojih su magnetni momenti paralelni (Slika 3.1). Ukupni magnetni momenti domena su haotično raspoređeni pa supstanca ne stvara magnetno polje oko sebe. Kada se feromagnetni materijal nađe u spoljašnjem magnetnom polju, rastu domeni čiji magnetni moment ima isti smer kao i 13 spoljašnje polje, a oni koji ne zadovoljavaju ovaj uslov, smanjuju se (Slika 3.1). Po prestanku delovanja polja supstanca zadržava magnetne osobine, odnosno magnetni momenti ostaju usmereni. Klasični primeri feromagnetika su gvožđe, kobalt, nikl zatim legure kao što su CuMn, MnBi, MnAs, MnSb i pojedini lantanidi poput gadolinijuma i disprozijuma . Slika 3.1 Šematski prikaz feromagnetnog poretka 4. Antiferomagnetici su supstance kod kojih postoji tendencija da se magnetni momenti susednih elektrona raspoređuju antiparalelno. Uslov da neka supstanca bude antiferomagnetik je da svaki atom u odnosu na susedni ima antiparalelnu orijentaciju magnetnog momenta (Slika 3.2) Antiferomagnetici imaju rezultujući magnetni momenat jednak nuli, što znači da ne stvaraju magnetno polje oko sebe. Slika 3.2 Šematski prikaz antiferomagnetnog poretka. Magnetni momenti su međusobno antiparalelno orijentisani, dajući rezultujući momenat jednak nuli 5. Ferimagnetici su supstance koje poput feromagnetika zadržavaju magnetizaciju uzorka po prestanku delovanja spoljašnjeg polja. Ferimagnetizam se javlja kod materijala koji su sastavljeni iz različitih pod-rešetki ili jona (magnetit Fe3O4), kao i u molekulskim magnetima30,31. Međutim slično antiferomagneticima magnetni momenti atoma različitih pod-rešetki nalaze se u antiparalelnom položaju. Pri optimalnom geometrijskom rasporedu ima više magnetnih momenata elektrona koji su usmereni u istom smeru, nego sa 14 antiparalelnim rasporedom (Slika 3.3). Takođe, suprotni magnetni momenti različitih pod- rešetki nisu jednaki, pa se magnetizacija spontano zadržava. Ferimagnetici pokazuju veoma izraženu magnetnu anizotropiju. Slika 3.3 Šematski prikaz ferimagnetnog poretka. Različite boje označavaju paralelne i antiparalelne orijentacije magnetnih momenata za različite pod-rešetke U cilju boljeg razumevanja magnetnog ponašanja određenih supstanci kao i njihove primene u tehnologiji i elektronici, ali i drugim oblastima kao što su biologija i medicina, potrebno je detaljnije objasniti pojam molekulskog magnetizma koji se odnosi na magnetne osobine izolovanih molekula i molekulskih sistema. Takvi molekuli mogu da sadrže jedan ili više magnetnih centara. Sistemi molekula se najčešće javljaju u molekulskim kristalima sa veoma slabim interakcijama između molekulskih jedinica. Takođe, mogu se naći i u sistemima izgrađenim od molekulskih prekursora, na takav način da se povećavaju interakcije između gradivnih jedinica što dovodi do znatnih promena u magnetnim osobinama na makroskopskom nivou. Ispitivanje molekulskih magneta daje uvid u osobine i svojstva potpuno različite u odnosu na uobičajene neorganske magnete. Molekulski magnetizam ima značajnu ulogu na polju molekularne elektronike, odnosno u primeni molekulskih sistema za izradu elektronskih uređaja i kola32,33,34. Pojava spinskih prelaza, karakterističnih za neke komplekse prelaznih metala, predstavlja značajan primer molekulske bistabilnosti. U sistemima sa velikim vrednostima magnetne anizotropije osnovno stanje je bi-stabilno u odnosu na paralelnu i antiparalelnu orijentaciju spinskih stanja, na relativno visokim temperaturama, što ima znatnu primenu na polju skladištenja informacija i kvantnog računarstva.35,36,37 Detaljnim objašnjenjem magnetnih osobina na kvantnom nivou može se steći uvid o osobinama na makroskopskom nivou. Pored toga, analiza porekla magnetne anizotropije na molekulskom nivou predstavlja put ka dizajniranju materijala sa pomenutim karakteristikama. 15 4 Analiza spinskog hamiltonijana Ispitivanje magnetnih osobina na kvantnom nivou zahteva poznavanje veličina koje su satavni deo spinskog hamiltonijana. Spinski hamiltonijan predstavlja kvantno-mehanički operator koji zavisi od spina, a ne zavisi od bilo koje prostorne koordinate ili vremena. Koncept spinskog hamiltonijana se bazira na radovima Abragam-a i Pryce-a38,39 i može se predstaviti jednačinom 4.1. QPHFNZEZEXZF kk k kkik i k ik T kn k n i ii T BBAii i i HHHHHH IPIIASIBg SgBSSJSDSH ˆˆˆˆˆˆ ˆˆˆˆˆ ˆˆˆˆˆˆ 11 1 , 1 11 +++++= ++− +−= ∑∑∑∑ ∑∑ == == ==   µ µ Jednačina 4.1 Prvi član u jednačini označava cepanje spinskih multipleta u odsustvu spoljašnjeg magnetnog polja (“zero field splitting-zfs” u daljem tekstu) za molekule sa vrednošću spina S>1/2. Drugi član se odnosi na magnetne interakcije između dva ili više metalnih centara. Treći i četvrti član predstavljaju interakcije spina elektrona i spina jezgra sa spoljašnjim magnetnim poljem (Zeeman-ove interakcije), dok peti i šesti član označavaju interakcije spin aktivnog jezgra sa nesparenim elektronima i kvadrupolne interakcije jezgra. Nadalje, detaljno će biti opisana samo prva dva člana, relevantna za problematiku ovog rada. 4.1 Magnetna anizotropija mononuklearnih kompleksa prelaznih metala; D vrednost Kao što je napomenuto zfs predstavlja uklanjanje degeneracije spinskih stanja u sistemima sa S>1/2 u odsustvu spoljašnjeg magnetnog polja (Slika 4.1). Degeneracija se uklanja kao posledica elektronske strukture i raspodele spinske gustine. Cepanje spinskih multipleta uzrokuje magnetnu anizotropiju u sistemima sa nižom-aksijalnom lokalnom 16 simetrijom. U sistemima sa visokom lokalnom simetrijom kao što je oktaedarska i tetraedarska zfs doprinosi podjednako svim nivoima, pa ne dolazi do pojave magnetne anizotropije. Oznaka D u jednačini 4.1 predstavlja zfs tenzor, dok Ŝ predstavlja operator spinskog ugaonog momenta. Energije spinskih mikrostanja usled cepanja se mogu predstaviti jednačinom 4.2. ( ) ( )      −+      +−= 222 3 1 YXZ SSE SSSDE Jednačina 4.2 Skalarne veličine D i E označavaju aksijalni i rombični parametar cepanja, dok S predstavlja kvantni broj ukupnog spinskog ugaonog momenta, a Sx, Sy i Sz ( ) ( )YYXX YYXXZZ DDE DDDD −= +−= 2 1 2 1 projekcije spinskog ugaonog momenta. Aksijalni i rombični parametri cepanja se mogu izraziti preko komponenata D tenzora na sledeći način (Jednačina 4.3): Jednačina 4.3 U slučaju kada je rombični parametar cepanja jednak nuli, a D<0, energije mikrostanja mogu se predstaviti šematski na sledeći način (Slika 4.1): Slika 4.1 Energije mikrostanja za sisteme sa S=1, S=3/2, S=2 i S=5/2 u slučaju kada je E=0, a D<0 17 Kada je energetski najstabilnije stanje sa minimalnom ms vrednosti, D parametar je veći od nule. Aksijalni parametar ujedno predstavlja meru magnetne anizotropije, odnosno što je veća apsolutna vrednost aksijalnog parametra, veća je i magnetna anizotropija sistema. Za tehnološku aplikaciju od primarnog je značaja dizajnirati i sintetisati takve molekule ili sisteme molekula kod kojih je aksijalni parametar izrazito negativan. Ukoliko je visina energetske barijere koja razdvaja stanja sa pozitivnim i negativnim vrednostima magnetnih momenata dovoljno velika, magnetizacija uzorka se može zadržati i po prestanku delovanja spoljašnjeg magnetnog polja, što se koristi za skladištenje informacija na kvantnom nivou (Slika 4.2). a. b. Slika 4.2 Šematski prikaz magnetne anizotropije sistema. U odsustvu spoljašnjeg magnetnog polja (Slika 4.2 a.) obe orijentacije spina su ekvivalentne, usled čega je populacija ravnomerna. U prisustvu spoljašnjeg magnetnog polja (Slika 4.2 b.) jedna od orijentacija je favorizovana. Ukoliko je barijera mala može doći do prelaska iz jedne orijentacije u drugu termalnim popunjavanjem viših stanja. Ako je barijera dovoljno visoka sistem zadržava favorizovanu orijentaciju i na višim temperaturama. Osnovni doprinos zfs-u, pored spin-dipolarnog efekta čiji je doprinos znatno manji, predstavljaju interakcije orbitalnog i spinskog ugaonog momenta. Postoje dva osnovna načina na koje se orbitalni ugaoni momenat manifestuje: ukoliko je osnovno stanje 18 degenerisano, ili spin-orbitalnim mešanjem određenih ekscitovanih stanja u osnovno stanje. U oba slučaja elektron doprinosi orbitalnom angularnom momentu samo ukoliko se orbitala koja sadrži elektron prostom rotacijom oko neke od tri ose koordinatnog sistema prevodi u drugu, po obliku identičnu orbitalu. Ovim se implicira da degeneracija osnovnog stanja nije sama po sebi uzrok pojave orbitalnog ugaonog momenta ako oblik degenerisanih orbitala nije isti (t2g orbitale kod oktaedarskih kompleksa imaju isti oblik, dok eg orbitale nemaju isti oblik). Osnovno selekciono pravilo pri ovakvim interkonverzijama je Δml=0,±1. Da bi se pravilno shvatila veza između orbitalnog ugaonog momenta i spinskog ugaonog momenta, potrebno je razmotriti orbitalno kretanje elektrona po kružnoj putanji u odsustvu spoljašnjeg magnetnog polja (Slika 4.3). Slika 4.3 Interakcije spinskog magnetnog momenta sa magnetnim poljem nastalim usled orbitalnog kretanja elektrona Usled obrtanja u smeru kazaljke na satu elektron poseduje orbitalni ugaoni momenat usmeren ka dole. Kružno kretanje elektrona dovodi do pojave magnetnog polja, usmerenog suprotno od smera orbitalnog ugaonog momenta tj. na gore. Usled pojave usmerenog magnetnog polja dolazi do orijentisanja spinskih magnetnih momenata elektrona u odnosu na polje i do uklanjanja degeneracije spinskih multipleta. Totalni hamiltonijan ispitivanog sistema se na osnovu toga može prikazati jednačinom 4.4. SOel HHH ˆˆˆ += Jednačina 4.4 gde Ĥel predstavlja uobičajeni elektronski hamiltonijan, dok ĤSO označava efektivni jednoelektronski operator spin-orbitalnog kuplovanja. Jednoelektronski operator spin- 19 orbitalnog kuplovanja, najpogodnije je predstaviti u funkciji rastojanja jezgra i elektrona (Jednačina 4.5): ( ) ( ) i i iii i i eff SO srslr Z H   ⋅=⋅= ∑∑ ξα λ λ λ χ 3 2 2 1ˆ Jednačina 4.5 gde α označava konstantu fine strukture spektra, Z(λ)eff efektivno naleketrisanje jezgra, ξ spin-orbitalnu konstantu metala, a riλ rastojanje između jezgra i elektrona. Nadalje, korišćenjem perturbacione teorije (Prilog) može se pokazati veza između energije sistema i spin-orbitalnog kuplovanja (Jednačina 4.6). ∆ ΨΨ −= 2 '0 0 e Z SO Z S S H S S EE Jednačina 4.6 Uopšteno govoreći numerator je proporcionalan spin-orbitalnoj konstanti metala, a energija Δ predstavlja eneriju ekscitovanja elektrona iz stanja Ψ0 u stanje Ψe. Molekulski sistemi sa velikom konstantom spin-orbitalnog kuplovanja i malom razlikom energije Δ pokazuju velike vrednosti magnetne anizotropije. U slučaju kada Δ→0, odnosno kada su stanja Ψ0 i Ψe Za bolje razumevanje problema, kao najjednostavniji primer je potrebno uzeti elongovani i kompresovani oktaedar u slučaju kada je S=2 (Slika 4.4). simetrijski i energetski degenerisana izraz 4.5 prestaje da važi, ali se ne narušava kvalitativnost rezultata. 20 Slika 4.4 Šematski prikaz elongovanog i kompresovanog oktaedra za model system sa S=2 U slučaju elongovanog oktaedra, korišćenjem jednačine 4.6 moguće je izraziti komponente D tenzora na sledeći način: ∆ −≅ − ΨΨΨΨ −= ∆ −≅ − ΨΨΨΨ −= ∆ −≅ − ΨΨΨΨ −= − −− − −− − −− 416 1616 1616 22 22 22 22 2222 22 2222 22 2222 eff xyyx yx M zxyxy M zyxeff ZZ eff xzyx yx M yxzxz M yyxeff YY eff yzyx yx M xyzyz M xyxeff XX ll D ll D ll D ζ εε ζ ζ εε ζ ζ εε ζ Jednačine 4.7-4.9 gde lx, ly i lz predstavljaju efektivni jednoelektronski operator spin-orbitalnog kuplovanja, koji se u datoj grupi simetrije transformiše kao rotacija oko ose koordinatnog sistema, a ζeff ∆ −= 16 3 2effD ζ efektivno naelektrisanje jezgra. Ako se D vrednost izrazi na osnovu jednačine 4.3 i jednačina 4.6-4.8 kao krajnji rezultat za elongovani oktaedar dobija se . Kako DZZ term dominira u odnosu na DXX i DYY izvodi se zaključak da spin-orbitalno kuplovanje 21 oko z ose favorizuje stanja sa maksimalnom ms vrednosti. U slučaju kompresovanog oktaedra izrazi za komponente D tenzora mogu se prikazati jednačinama 4.10-4.12. 0 16 3 16 16 3 16 22 22 2 22 2 22 = ∆ −≅ − ΨΨΨΨ −= ∆ −≅ − ΨΨΨΨ −= ZZ eff xzz z M yxzxz M yzeff YY eff yzz z M xyzyz M xzeff XX D ll D ll D ζ εε ζ ζ εε ζ Jednačine 4.10-4.12 Imajući u vidu da je DZZ=0 t.j da ne postoji popunjena orbitala koja bi se rotacijom oko z ose transformisala u dz2 orbitalu, aksijalni parametar cepanja je pozitivan. Na osnovu analiziranih izraza osnovno selekciono pravilo za kuplovanje preko lz operatora je Δml=0, dok je za kuplovanje preko lx i ly operatora Δm l Ovakvo tumačenje ipak je nepotpuno ako se u obzir uzmu samo prostorne funkcije, a izostave spinske funkcije od kojih zavisi spinski hamiltonijan. Talasna funkcija elektrona koji se ekscituje, uzimajući u obzir spinske funkcije elektrona može se izraziti jednačinom 4.13. =±1. β α ,;2/1,2/1, 2/1 2/1 2/1 2/1 ,;2/1,2/1, 2/1 2/1 2/1 2/1 aSSD S S S S S S aSSD S S S S S S Z ZZ a Z Z ZZ a Z +−      −+ − =Ψ +−−      − − =Ψ Jednačina 4.13 Simboli u zagradama predstavljaju Clebsch-Gordan-ove koeficijente, a označava prostornu funkciju elektrona, dok α i β označavaju spinske funkcije elektrona. Interakcije osnovnog i ekscitovanog stanja usled spin orbitalnog kuplovanja mogu se prikazati na sledeći način (Jednačina 4.14): 22 ββ αβ βα αα ,;2/1,2/1,,;2/1,2/1, 2/12/1 2/12/1 2/12/1 2/12/1 ,;2/1,2/1,,;2/1,2/1, 2/12/1 2/12/1 2/12/1 2/12/1 ,;2/1,2/1,,;2/1,2/1, 2/12/1 2/12/1 2/12/1 2/12/1 ,;2/1,2/1,,;2/1,2/1, 2/12/1 2/12/1 2/12/1 2/12/1 bdSSDHaSSD dS S dS S S S S S bdSSDHaSSD dS S dS S S S S S bdSSDHaSSD dS S dS S S S S S bdSSDHaSSD dS S dS S S S S S dS S H S S Z c SOZ ZZZZ Z c SOZ ZZZZ Z c SOZ ZZZZ Z c SOZ ZZZZ b Z c SOa Z ++−+−      −++ − +       −+ − + +−+−+−      −+ − +       −+ − + +++−−−      −++ − +       − − + +−+−−−      −+ − +       − − = =Ψ + Ψ Jednačina 4.14 Kombinovanjem jednačina 4.5 i 4.14 kao krajnji rezultat dobija se jednačina 4.15. αξβδ βξαδ βξβαξαδ ,, 2/12/1 2/12/1 12/12/1 2/12/1 ,, 2/12/1 2/12/1 12/12/1 2/12/1 ,, 2/12/1 2/12/1 ,, 2/12/1 2/12/1 1, 1, 22 0, bsa S S S S S S S S bsa S S S S S S S S bsa S S S S bsa S S S S dS S H S S cc ZZZZ d cc ZZZZ d cc ZZ cc ZZ d b Z c SOa Z       + − +       −+ − + +      −− − −       − − + +               −+ − +      − − = =Ψ + Ψ + − Jednačina 4.15 Analizom izraza 4.15 dolazi se do zaključka da je član term različit od nule samo ako je c=z, dok su druga dva člana različita od nule kada je c=x,y. Spin orbitalno kuplovanje oko z ose moguće je samo ako se pri tome ne menja spin elektrona (interakcije između stanja sa istom ms vrednosti), dok je kuplovanje oko x i y ose moguće samo uz promenu spina (interakcije između stanja sa različitom ms vrednosti) (Slika 4.5). Ako se razmotri prvi term u jednačini 4.15, jedan od Clebsch-Gordan-ovih koeficijenata uvek je jednak jedinici dok je drugi jednak nuli, ukoliko se međusobno kupluju stanja sa maksimalnom ms vrednosti. 23 Slika 4.5 Spin-orbitalno kuplovanje preko lz operatora (gore), moguće samo između stanja sa istom ms vrednosti i kuplovanje preko lx i ly operatora (dole) između stanja sa različitim ms vrednostima Ovakvo tumačenje je sasvim opravdano ako se uzme u obzir da se maksimalna Sz Kada se kupluju stanja sa nižim m vrednost može prikazati jednom Slater-ovom determinantom. s Na osnovu dosadašnjih teorijskih razmatranja jasno je da će kompleksi kod kojih je moguće kuplovanje stanja preko H vrednostima Clebsch-Gordan-ovi koeficijenti su manji od jedan pa je i doprinos takvih interakcija znatno manji (Slika 4.5). SO z operatora pokazivati veću tendenciju ka negativnim D vrednostima. Interakcije preko HSOz operatora, imajući u vidu selekciono pravilo Δml=0, moguće su između dx2-y2 i dxy orbitala kao između dxz i dyz 4.1.1 Spin dipolarni efekat orbitala. Ukoliko je energetska razlika između pomenutih orbitala manja, utoliko je D parametar negativniji (Jednačina 3.5). Kod kompleksnih jedinjenja kod kojih je efekat spin orbitalnog kuplovanja zanemarljiv i kod kojih je osnovno stanje nedegenerisano (visokospinski kompleksi d5 konfiguracije), potrebno je uzeti u obzir uticaj dipolarnog efekta na magnetnu anizotropiju. Vrednosti magnetne anizotropije su u ovakvim slučajevima veoma male. Dipolarne interakcije potiču od uzajamnih delovanja spinskih magnetnih momenata elektrona, odnosno usled spinskog ugaonog momenta elektroni stvaraju magnetno polje sa kojim 24 interaguju ostali elektroni. Matematički oblik dipolarnih interakcija može se izraziti jednačinom 4.16. ( )( )         −= ∑ ≠ 23 2 31 2 ij jijiij ji ji ij B dip r mrmr mm r E  µ Jednačina 4.16 Oznake mi i mj predstavljaju magnetne momente u tačkama i i j, μB označava Bohr-ov magneton, a rij ( )ij ji ij jiB dip r mm E θµ 23 2 cos31 2 −= ∑ ≠ predstavlja rastojanje između dva položaja. Kako je energetski povoljno da magnetni momenti budu paralelni, izraz se može napisati u redukovanom obliku (Jednačina 4.17): Jednačina 4.17 gde θij predstavlja ugao između magnetizacije i rastojanja rij . 4.1.2 Uticaj Jahn-Teller-ovog efekta na magnetnu anizotropiju Najveći doprinos spin-orbitalnog efekta prisutan je u sistemima sa degenerisanim osnovnim stanjima, usled čega nema utroška energije potrebne za ekscitovanje elektrona. U slučaju kada su dx2-y2 i dxy orbitale degenerisane, mogu se očekivati visoke vrednosti magnetne anizotropije. Nelinearni molekuli sa prostorno degenerisanim osnovnim stanjem energetski su nestabilni i podložni geometrijskoj distorziji, usled koje se snižava simetrija i energija sistema uz otklanjanje degeneracije osnovnog stanja. Ovaj elektronski efekat, nazvan Jahn-Teller-ov (JT) efekat40,41,42,43,44, deluje nasuprot spin-orbitalnom kuplovanju, usled čega se može očekivati smanjenje magnetne anizotropije sistema. Ukoliko su spin- orbitalne interakcije dovoljno jake, Jahn-Teller-ova distorzija je slabije izražena. Na osnovu pomenutih razmatranja jasno je da se radi o dva antagonistička efekta od kojih svaki može biti kontrolisan izborom metala i liganada. Kao pogodan model za objašnjenje uzajamnog 25 dejstva spin-orbitalnih interakcija i Jahn-Teller-ovog efekta potrebno je uzeti penta- koordinovani sistem d8 elektronske konfiguracije, sa D3h simetrijom okruženja (geometrija trigonalne bipiramide) (Slika 4.6). Najpre je potrebno razmotriti slučaj kada je distorzija molekula dominantna u odnosu na spin-orbitalne interakcije, a zatim obrnut slučaj, kada je dominantno spin-orbitalno kuplovanje. Slika 4.6 Jahn-Teller-ova distorzija penta-koordinovanih sistema d8 elektronske konfiguracije Osnovno stanje sistema, 3E sadrži šest mikrostanja 3θ-1, 3θ0, 3θ1, 3ε-1, 3ε0 i 3ε1. Distorzijom duž e vibracije (E⊗e problem) snižava se simetrija od D3h do C2v, a osnovni term 3E razlaže se na 3A1 i 3B2 (Slika 4.6). Usled Jahn-Teller distorzije javlja se energetska razlika između dx2-y2 i dxy Ukoliko bi spin-orbitalno kuplovanje delovalo nezavisno od distorzije, degeneracija šest mikrostanja ( orbitala, što umanjuje efekat spin-orbitalnih interakcija i smanjuje vrednost magnetne anizotropije. 3θ-1, 3θ0, 3θ1, 3ε-1, 3ε0 i 3ε1) bila bi uklonjena u skladu sa prethodnim teorijskim postavkama (Slika 4.7). 26 Slika 4.7 Uklanjanje degeneracije mikrostanja spin-orbitalnim interakcijama u penta-koordinovanim sistemima d8 elektronske konfiguracije i D3h simetrije, ne uzimajući u obzir Jahn-Teller-ov efekat Energetski su najstabilnija stanja 3θ-1 i 3θ1 sa ms=±1 (D<0). Kako je Jahn-Teller elektronski efekat i ne zavisi od spinskih funkcija, moguća je interakcija samo onih stanja sa istim ms vrednostima, odnosno 3θ-1 stanje se može kuplovati samo sa ekscitovanim stanjem 3ε-1, 3θ1 sa stanjem 3ε1 dok se 3θ0 kupluje sa stanjem 3ε0 Uzimajući u obzir JT i spin-orbitalne interakcije, hamiltonijan za dati sistem može se izraziti jednačinom 4.18. . Kuplovanje osnovnog sa ekscitovanim stanjima predstavlja pseudo Jahn-Teller-ov efekat (PJT), koji je znatno slabiji u odnosu na Jahn-Teller-ov efekat prvog reda i direktno zavisi od energetske razlike kuplovanih stanja. ( )                     +− +− +− +− +− +− ++= θε θε θε εθ εθ εθ εθ ξ ξ ξ ξ ξ ξ VQVQi VQVQic VQVQi VQiVQ VQicVQ VQiVQ QQKH 0000 0000 0000 0000 0000 0000 2 1ˆ 22 Jednačina 4.18 Oznaka ξ predstavlja konstantu spin-orbitalnog kuplovanja za Sz=1→Sz=1 i Sz=-1→Sz=-1 stanja, a cξ redukovanu konstantu za stanja Sz=0→Sz=0, dok su V i K Konstanta JT kuplovanja i konstanta sile. Dijagonalizacijom matrice dobijaju se energije stanja (Jednačina 4.19). 27 ( ) ( ) ( ) ( ) ( )2222220 222222 1 2 1 2 1 cfQQVQQKE fQQVQQKE ++±+= ++±+=± εθεθ εθεθ Jednačina 4.19 Energetska razlika između E±1 i E0 nivoa daje vrednost aksijalnog parametra cepanja (Jednačina 4.20). ( ) ( ) ( ) ( ) 0 22 22222222 10 VQ fcffQQVcfQQVEE −≈++−++=− ± εθεθ Jednačina 4.20 Grafički prikaz dva moguća slučaja, kada je D>0 (dominantan Jahn-Teller-ov efekat) i kada je D<0 (dominantno spin-orbitalno kuplovanje) dat je na Slici 4.8. Slika 4.8 Uklanjanje degeneracije spinskih mikrostanja u slučaju kada je dominantan Jahn-Teller-ov efekat (levo) i kada je dominantno spin-orbitalno kuplovanje (desno) U D3h simetriji (maksimumi na krivama potencijalne energije) stanja sa maksimalnom Sz vrednosti (ms=±1) su u oba slučaja energetski stabilnija, bez obzira na intenzitet dva suprotna efekta. Ukoliko je JT efekat dominantan (Slika 4.8 levo), distorzija sistema dovodi do interkonverzije stanja, pa je u minimumu najstabilnije stanje sa minimalnom Sz vrednosti (ms=0). Nasuprot tome jake spin-orbitalne interakcije (Slika 4.8 desno), dovode do veće energetske razlike između stanja sa različitim ms vrednostima t.j do većih D 28 vrednosti, pa su samim tim PJT i JT slabije izraženi. Kao rezultat, interkonverzija stanja nije moguća, pa je samim tim D<0. Međutim treba istaći da JT efekat nije sam po sebi uzrok pozitivnih D vrednosti, kao što ni spin-orbitalno kuplovanje nije uvek uzrok negativnih D vrednosti, već zavisi od elektronske konfiguracije i tipa metala, kao i simetrije okruženja. 4.1.3 Uticaj elektronske konfiguracije na magnetnu anizotropiju trigonalno bipiramidalnih sistema Kuplovanjem dx2-y2 i dxy orbitala, kao i dxz i dyz orbitala mogu se postići visoke vrednosti magnetne anizotropije. U cilju povećanja magnetne anizotropije kompleksa prelaznih metala, a na osnovu dosadašnjih razmatranja, potrebno je izabrati penta- koordinovane molekule, D3h simetrije, kod kojih je osnovno stanje degenerisano (E’ ili E“). Efekti Jahn-Teller-ove distorzije mogu se umanjiti adekvatnim izborom liganada. Pored model kompleksa d8 elektronske konfiguracije, potrebno je analizirati sisteme sa visokospinskom d3 i d6 elektronskom konfiguracijom (Slika 4.9). Slika 4.9 Raspored elektrona u orbitalama za d3 i d6 elektronsku konfiguraciju kod penta- koordinovanih sistema D3h simetrije U slučaju d3 konfiguracije osnovno stanje je 4E’, pa je kuplovanje preko HSOz operatora moguće između dx2-y2 i dxy orbitala. Osnovno stanje sadrži osam mikrostanja: 4θ-3/2, 4θ-1/2, 4θ1/2, 4θ3/2 ,4ε-3/2, 4ε-1/2, 4ε1/2 i 4ε3/2 . Delovanjem spin-orbitalnog kuplovanja uklanja se degeneracija mikrostanja kao što je prikazano na slici 4.10. 29 Slika 4.10 Uklanjanje degeneracije 4E' mikrostanja spin-orbitalnim interakcijama u penta- koordinovanim sistemima d3 elektronske konfiguracije i D3h simetrije, ne uzimajući u obzir Jahn- Teller-ov efekat Strukture sa d3 (S=3/2) i d8 (S=1) elektronskom konfiguracijom predstavljaju sisteme sa jednostepenom barijerom (3/2→1/2 i 1→0). Kod struktura sa d6 konfiguracijom osnovno stanje je 5E’’, koje sadrži deset mikrostanja: 5θ-2, 5θ-1, 5θ0, 5θ1 , 5θ2, 5ε-2, 5ε-1, 5ε0, 5ε1, 5ε2. Kuplovanje preko HSOz operatora moguće je samo između dxz i dyz orbitala. Uklanjanje degeneracije mikrostanja spin-orbitalnim kuplovanjem prikazano je na slici 4.11. Slika 4.11 Uklanjanje degeneracije 5E'' mikrostanja spin-orbitalnim interakcijama u penta- koordinovanim sistemima d6 elektronske konfiguracije i D3h simetrije, ne uzimajući u obzir Jahn- Teller-ov efekat 30 Centralni metalni joni sa d6 Interkonverzija između d konfiguracijom predstavljaju sisteme sa dvostepenom barijerom (2→1→0). xz i dyz orbitala daje znatno manje vrednosti magnetne anizotropije u odnosu na interakcije između dx2-y2 i dxy orbitala. Kretanje elektrona u xy ravni pri transformaciji dx2-y2→ dxy dovodi do indukovanja magnetnog polja usmerenog duž z ose, u odnosu na koje se orijentišu spinski magnetni momenti elektrona (Slika 4.3). Ukoliko se elektron nalazi u dz2 orbitali (kreće se duž z ose), kuplovanje preko HSOz operatora nije moguće u skladu sa selekcionim pravilima, pa ne dolazi do indukovanja magnetnog polja. Kako je elektronska gustina u dxz i dyz orbitalama delimično raspoređena duž z ose, transformacija dxz → dyz indukuje znatno slabije polje. 4.2 Izračunavanje D vrednosti u mononuklearnim kompleksima prelaznih metala LF-DFT metodom Iz dosadašnjih razmatranja jasno se vidi da se prilikom izračunavanja magnetne anizotropije, moraju uzeti u obzir ekscitovana stanja, odnosno interakcije osnovnog stanja sa ekscitovanim stanjima. Za egzaktno rešenje problema potrebno je uzeti u obzir linearnu kombinaciju različitih Slater-ovih determinanti. Kako je DFT jedno-determinantni pristup, tačne vrednosti D parametra ne mogu se dobiti konvencionalnim DFT metodama. Kombinacijom DFT-a i teorije ligandnog polja (Ligand Field Theory-LFT)45, mogu se rešiti postojeći problemi i izračunati D vrednosti46,47,48,49,50,51,52,53,54. Osnovni nedostatak ovakvog pristupa predstavlja izbor sistema, odnosno kompleksi sa izrazito kovalentnim karakterom metal-ligand veze ne daju dobre rezultate. Teorijske postavke teorije ligandnog polja opisane su u literaturi45, tako da će ovde samo ukratko biti napomenut postupak za izračunavanje D vrednosti. Prvi deo procedure zasniva se na DFT proračunima i ekstrakciji neophodnih podataka. Drugi deo obuhvata izračunavanje Racah parametara, LF matrice i orbitalnog redukcionog faktora (orf) na osnovu podataka dobijenih DFT proračunima, dok se u trećem koraku na osnovu Racah parametara B i C, LF matrice, orbitalnog redukcionog faktora (orf) i konstante spin-orbitalnog kuplovanja izračunavaju energije odgovarajućih 31 multipleta, na osnovu kojih se može dobiti D vrednost. DFT proračuni se mogu podeliti u tri koraka. Prvi korak obuhvata geometrijsku optimizaciju ispitivanog sistema i analizu molekulskih orbitala. Drugi korak obuhvata AOC proračune (Average of configuration), na geometrijama dobijenim iz prvog koraka, gde se ukupan broj d elektrona ravnomerno raspoređuje po svim KS MO sa dominantnim karakterom d orbitala metala t.j za sistem sa d8 konfiguracijom osam elektrona se smešta u pet orbitala tako da je svaka orbitala na metalu okupirana sa 8/5=1.6 elektrona. Gustina dobijena u AOC proračunima najbolje je prilagođena modelu kojim se tretiraju međuelektronske odbojne interakcije u okviru teorije ligandnog polja. Energije svih Slater-ovih determinanti (SD), koje potiču od dn Energije Slater-ovih determinanti se dalje parametrizuju u okviru LFDFT dela u cilju dobijanja LF matrice dimenzija 5x5, konfiguracije ispitivanog sistema, računaju se u trećem koraku DFT proračuna, koristeći KS elektronsku gustinu dobijenu iz AOC proračuna. ii dAd ˆ , i Racah parametara B i C, koji se uz orbitalni redukcioni faktor (dobijen takođe u LFDFT koraku), unose u LF program kojim se izračunavaju energije multipleta. Slaganje SD energija dobijenih DFT i LFDFT proračunima dobija se takođe kao jedan od izlaznih podataka na kraju LFDFT koraka (Slika 4.12). 0 5 10 15 20 25 30 35 40 45 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0 x 10 4 Slika 4.12 Poređenje između SD energija dobijenih DFT proračunima (*) i LFDFT proračunima (°). Energija od 0eV odgovara referentnoj vrednosti AOC proračuna 32 Aksijalni parametar cepanja, D, moguće je dobiti iz razlike energija multipleta, ako se u LF delu procedure, pored pomenutih parametara, unese konstanta spin-orbitalnog kuplovanja metala. Navedena LF-DFT procedura se grafički može prikazati slikom 4.13. Slika 4.13 Šematski prikaz LF-DFT procedure. DFT korak podrazumeva tri stupnja na osnovu kojih se dobijaju podaci neohodni za LFDFT korak. LFDFT proračuni daju vrednosti Racah parametara i LF matricu, koji se koriste u LF programu za dobijanje D vrednosti 4.3 Magnetna kuplovanja u binuklearnim kompleksima prelaznih metala Dosadašnja razmatranja bila su usmerena na međusobne magnetne interakcije spinskih i orbitalnih magnetnih momenata, odnosno na magnetne osobine mono-nuklearnih kompleksa prelaznih metala. Kako se u čvrstom stanju metalni centri mogu naći na udaljenosti dovoljnoj za interakciju magnetnih momenata, potrebno je proširiti teorijska razmatranja i detaljnije objasniti magnetne osobine binuklearnih kompleksa. Teorijska Geometrijska optimizacija AOC proračuni SD proračuni LFDFT program LF program Energije Slater-ovih determinanti Racah parametri i LF matrica DFT LFDFT LF MO sa dominantnim d karakterom 33 razmatranja potrebno je započeti sa sistemom od dva metalna jona od kojih svaki ima d1 elektronsku konfiguraciju (Slika 4.14). Slika 4.14 Šematski prikaz interakcije dva metalna jona i oznake odgovarajućih termova Ukoliko su dva katjona na udaljenosti na kojoj mogu da interaguju totalni hamiltonijan se može prikazati jednačinom 4.21. 'ˆˆˆˆ HHHH BA ++= Jednačina 4.21 Hamiltonijani ĤA i ĤB R Z rrr ZH AB 2 1221 111'ˆ ++      +−= predstavljaju jednoelektronske operatore za svaki metalni centar, a Ĥ’ označava perturbacioni operator usled interakcije dva centra (Jednačina 4.22): Jednačina 4.22 Prvi član u perturbacionom operatoru predstavlja privlačne interakcije elektrona jednog atoma i jezgra drugog atoma, dok druga dva člana odgovaraju odbojnim interakcijama između dva elektrona i između dva jezgra. Oznake r1B i r2A predstavljaju rastojanje između jezgra i elektrona, r12 Kada je rastojanje između centara veliko sistem je dvostruko degenerisan. Degeneracija se uklanja perturbacijom preko operatora Ĥ’ pa se talasne funkcije mogu napisati kao simetrična i antisimetrična kombinacija (Jednačina 4.23): označava rastojanje između dva elektrona, dok R označava udaljenost između jezgara. 34 ( ) ( ) ( ) ( ) ( )[ ]122112 1 2 BABAS ΦΦ±ΦΦ ± Jednačina 4.23 Ukoliko se prostornim funkcijama pridodaju spinske funkcije, simetrična funkcija odgovara singletnom stanju, dok antisimetrična funkcija odgovara tripletnom stanju (Jednačina 4.24). ( ) ( ) ( ) ( )[ ] ( ) ( ) ( ) ( )[ ] ( ) ( ) ( ) ( )[ ] ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )           +ΦΦ−ΦΦ − =Ψ −ΦΦ+ΦΦ + =Ψ 21 1221 21 2 11221 )1(2 1 1221 2 11221 )1(2 1 2 3 2 1 ββ βαβα αα βαβα BABA BABA S S Jednačina 4.24 Imajući u vidu jednačine 4.21 i 4.22, kao i oblik talasnih funkcija, energije singletnog i tripletnog stanja mogu se izraziti preko odgovarajućih integrala (Jednačina 4.25): ( ) ( ) 203 2 01 1 2 1 2 S KJEE S KJEE ABAB ABAB − − +=Ψ + + +=Ψ Jednačina 4.25 Oznaka 2E0 predstavlja jedinstvenu vrednost HA + HB, dok JAB i KAB ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )12121 21121 12 12 BABAAB BABAAB r K r J ΦΦΦΦ= ΦΦΦΦ= označavaju Coulomb-ov integral i integral izmene (Jednačina 4.26). Jednačina 4.26 Energetska razlika između singletnog i tripletnog stanja definiše jačinu magnetnih interakcija između metalnih centara i može se izraziti jednačinom 4.27. 35 ( ) ( ) ( )4 2 31 1 2 S SJKEE ABAB − − =Ψ−Ψ Jednačina 4.27 Energetska razlika može se prikazati posredstvom efektivnog hamiltonijana koji je deo ukupnog spinskog hamiltonijana (Jednačina 4.1-drugi član)55,56. Oznake ŜA i ŜB BA SSS ˆˆˆ  += predstavljaju efektivne operatore spinskog ugaonog momenta, a J označava konstantu kuplovanja između paramagnetnih centara A i B. Totalni operator spinskog ugaonog momenta može se izraziti kao zbir pojedinačnih operatora (Jednačina 4.28). Jednačina 4.28 Na osnovu izraza 4.28 efektivni hamiltonijan se može reformulisati u drugačijem obliku (Jednačina 4.29).      −−−= 222 ˆˆˆ 2 ˆ BAEX SSS JH  Jednačina 4.29 Svojstvene vrednosti operatora se mogu izraziti u opštoj formi na sledeći način (Jednačina 4.30): ( ) ( ) ( ) ( )[ ]111 2 12 +−+−+−=Ψ+ BBAA S SSSSSSJE Jednačina 4.30 Oznake S, SA i SB predstavljaju kvantne brojeve povezane sa spinskim operatorima. Ukoliko se za tripletno stanje unesu vrednosti S=1, SA=1/2 i SB=1/2 kao krajnji rezultat dobija se E(3Ψ)=-J/4. Zamenom vrednosti S=0, SA=1/2 i SB=1/2 u izraz 4.30 dobija se energija singletnog stanja E(1Ψ)=3J/4. Energetska razlika jednaka je konstanti kuplovanja J. Ako je konstanta negativna, singletno stanje je stabilnije pa je kuplovanje antiferomagnetno. Kada je konstanta pozitivna stabilnije je tripletno stanje, sa paralelnom orijentacijom spinova, pa su kuplovanja feromagnetna. Za izradu prekidača neophodno je postići antiferomagnetno kuplovanje, jer se promenom spina na jednom metalnom centru 36 automatski menja spin na drugom centru. Ukupna konstanta kuplovanja se može izraziti kao suma feromagnetnih i antiferomagnetnih doprinosa (Jednačina 4.31). AFF JJJ += Jednačina 4.31 Za objašnjenje mehanizma kuplovanja najjednostavnije je zadržati se na modelu sa dva metalna centra elektronske konfiguracije d1 (Slika 4.15)57. Uticaji mostnih liganada će radi jednostavnosti biti zanemareni, mada će kasnije biti detaljnije opisani. Prvo je potrebno razmotriti delokalizaciju koja dovodi do stabilizacije singletnog stanja u odnosu na tripletno. Talasna funkcija neperturbovanog singletnog stanja može se prikazati kako je to naznačeno u jednačini 4.24. Usled kinetičke energije elektrona i privlačnih potencijala jezgra B, elektron iz orbitale ΦA može se delokalizovati u orbitalu ΦB u kojoj se nalazi nespareni elektron (Slika 4.15 (a)). Takođe je moguć i obrnut slučaj delokalizacije elektrona iz ΦB orbitale u ΦA orbitalu. Slika 4.15 Mehanizam delokalizacije elektrona sa jednog metalnog centra na drugi bez posredstva mostnih liganada. Stabilizacija singleta u odnosu na triplet (a) dešava se usled delokalizacije elektrona iz ΦA orbitale u ΦB Φ orbitalu, dok se stabilizacija tripleta (b) dešava usled delokalizacije iz A orbitale u φΒ orbitalu na centru B 37 Delokalizacija dovodi do nastajanja jonskog singletnog stanja (charge transfer singlet-CT) kod koga se oba elektrona nalaze na jednom metalnom centru. Talasne funkcije ekscitovanih jonskih singleta mogu se prikazati jednačinom 4.32. ( ) ( ) ( ) ( )[ ] ( ) ( ) ( ) ( )[ ] ( ) ( ) ( ) ( )[ ] ( ) ( ) ( ) ( )[ ]1221 2 12121 )1(2 1 1221 2 12121 )1(2 1 2 1 2 1 βαβα βαβα +ΦΦ−ΦΦ − =Ψ −ΦΦ+ΦΦ + =Ψ BBAA CT u BBAA CT g S S Jednačina 4.32 Energetska razlika između kovalentnog i jonskog singleta, U, predstavlja razliku u međuelektronskim odbijanjima kada se elektroni nalaze na istom centru i kada se nalaze na različitim centrima (Jednačine 4.33-4.35).       ΦΦΦΦ=       ΦΦΦΦ==      ΦΦΦΦ= −= BBAAAB BBBBBBAAAAAA ABAA r U r U r U UUU 12 1212 1 11 Jednačine 4.33-4.35 Usled delokalizacije elektrona dolazi do mešanja ekscitovanog singleta 1Ψg u osnovno stanje. Ovakvo mešanje dovodi do interakcija osnovnog singleta sa ekscitovanim singletom 1Ψg pri čemu dolazi do povišenja energije jonskog singleta i sniženja energije kovalentnog singleta u odnosu na tripletno stanje. Interakcije su jače ukoliko je energetska razlika U manja. Međusobne interakcije singleta mogu se izraziti preko sledeće matrice (Jednačina 4.36): ( ) ( ) ( ) ( ) ΦΦ==ΨΦ→ΦΨ==ΨΦ→ΦΨ tSVSV ABCTABACTB 20ˆ0ˆ Jednačina 4.36 38 Oznaka tΦΦ predstavlja integral prenosa naelektrisanja, koji je na osnovu perturbacione teorije (Prilog) direkto povezan sa energetskom razlikom između osnovnog singleta i tripleta (Jednačina 4.37). ( ) ( )[ ] U tEE 2 13 4 ΦΦ=Ψ−Ψ Jednačina 4.37 Do stabilizacije tripleta u odnosu na singlet može doći usled transfera elektrona iz ΦA orbitale u φΒ orbitalu na centru B (Slika 4.15 (b)). Talasne funkcije ekscitovanog singleta i tripleta mogu se prikazati jednačinom 4.38. ( ) ( ) ( ) ( )[ ] ( ) ( ) ( ) ( )[ ] ( ) ( ) ( ) ( )[ ] ( ) ( ) ( ) ( )[ ] ( ) ( ) ( ) ( )[ ] ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )[ ] ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )           +Φ−Φ − =Ψ           +Φ−Φ − =Ψ −Φ+Φ + =Ψ −Φ+Φ + =Ψ 21 1221 21 2 11221 )1(2 1 21 1221 21 2 11221 )1(2 1 1221 2 11221 )1(2 1 1221 2 11221 )1(2 1 2 3 2 3 2 1 2 1 ββ βαβα αα φφ ββ βαβα αα φφ βαβαφφ βαβαφφ AAAA CT A BBBB CT B AAAA CT A BBBB CT B S S S S Jednačina 4.38 Energetska razlika između ekscitovanog i osnovnog singleta iznosi U’+Δ+KAB, gde U’ predstavlja razliku u međuelektronskim odbijanjima kada se elektroni nalaze na istom centru ali u različitim orbitalama i kada se nalaze na različitim centrima. Oznaka Δ predstavlja razliku u energijama između Φ i φ orbitala, dok KAB označava integral izmene. Osnovni triplet je odvojen od ekcitovanog tripleta energetskom barijerom U’+Δ-KAB. Kako je 3Ψ CT stabilniji u odnosu na 1Ψ CT, interakcije između osnovnog i ekscitovanog tripleta znatno su izraženije u odnosu na interakcije osnovnog i ekscitovanog singleta. Kao 39 posledica toga dolazi do većeg sniženja energije osnovnog tripleta u odnosu na singlet. Energetska razlika se može izraziti jednačinom 4.39. ( ) ( )[ ]       +∆+ − −∆+ =Ψ−Ψ Φ ABAB KUKU tEE ' 1 ' 12 213 φ Jednačina 4.39 Delokalizacija elektrona takođe se može ostvariti posredstvom dijamagnetnih mostnih liganada (Slika 4.16). Slika 4.16 Mehanizam delokalizacije elektrona posredstvom orbitala mostnih liganada Mehanizam delokalizacije obuhvata prenos naelektrisanja sa liganda na jedan od metalnih centara, praćen prenosom naelektrisanja sa drugog metala na ligand. Odgovarajući doprinos kuplovanju je antiferomagnetni i može biti dominantan ukoliko prethodno opisani mehanizmi delokalizacije nisu jako izraženi (Jednačina 4.40). UE ttJ CT LBLA 2 224 ∆ −← Jednačina 4.40 Delokalizacija elektrona sa liganda na metal, može biti praćena dodatnom delokalizacijom sa orbitala liganda na drugi metal, pri čemu nastaje dvostruka jonska struktura (Slika 4.17). 40 Slika 4.17 Mehanizam delokalizacije elektrona posredstvom mostnih liganada uz formiranje dvostruke jonske strukture Ovakav doprinos je uvek antiferomagnetni, iako izrazito mali, pošto je energetska razlika između osnovnog i dvostruko-jonskog singleta ΔE2CT CTCT LBLA EE ttJ 2 2 228 ∆∆ −← znatno veća u odnosu na energiju U (Jednačina 4.41). Jednačina 4.41 Dvostruki prenos naleketrisanja sa liganda na metal poznat je u fizici čvrstog stanja kao Goodenough-ov mehanizam. Potrebno je napomenuti da delokalizacija preko mostnih liganada ne mora uvek dovesti do antiferomagnetnih interakcija. Ukoliko su orbitale metala između kojih postoje interakcije ortogonalne kuplovanje je feromagnetno (Slika 4.18). Slika 4.18 Mehanizam antiferomagnetnog kuplovanja delokalizacijom preko orbital liganada (levo) i mehanizam feromagnetnog kuplovanja između ortogonalnih orbital metala (desno) Sa porastom broja nesparenih elektrona po metalnom centru povećava se broj kako antiferomagnetnih interakcija, tako i broj feromagnetnih interakcija između ortogonalnih orbitala, pa rezultujuća kuplovanja mogu biti veoma mala. Čeona preklapanja d orbitala metala sa p orbitalama liganada, odnosno σ interakcije dovode do intenzivnijih antiferomagnetnih interakcija u odnosu na bočna π preklapanja. Imajući u vidu jednačinu 4.31 i do sada navedene mehanizme totalna konstanta kuplovanja može se reformulisati na sledeći način (Jednačina 4.42). 41 U TKJ ABAB 242 −= Jednačina 4.42 Prvi član predstavlja feromagnetne doprinose, dok drugi član favorizuje antiferomagnetne interakcije. Pojedinačni doprinosi mogu se dobiti i u slučaju kada se u obzir uzmu molekulske orbitale nastale linearnom kombinacijom d orbitala sa dva metalna centra. Radi jednostavnosti potrebno je pretpostaviti da su d orbitale metala prostorno ortogonalne (Jednačina 4.43). ( ) ( )BAu BAg Φ−Φ= Φ+Φ= 2 1 2 1 σ σ Jednačina 4.43 Raspoređivanjem dva elektrona u dve molekulske orbitale moguće je konstruisati četiri determinante, od kojih tri odgovaraju singletnim stanjima, (σg)2, (σu)2, (σgα σuβ), a jedna odgovara tripletnom stanju (σgα σu guug guguug uuu ggg JE KJE JE JE ++= −++= += += εε εε ε ε 4 3 2 1 2 2 α). Energije datih determinanti mogu se izraziti na sledeći način (Jednačine 4.44-4.47). Jednačine 4.44-4.47 Parametri KAB, U i TAB iz jednačine 4.42 mogu se izraziti preko energija Coulomb-ovih integrala (Jgg, Juu i Jgu), integrala izmene (Kgu) i energija ε(g) i ε(u) (Jednačine 4.48-4.50). 42 ( ) ( ) ( ) ( )12 34 421 4 1 22 2 4 12 4 1 EET EEKU EEEJJJK AB gu guuuggAB −= −== −+=−+= Jednačine 4.48-4.50 Kao krajnji rezultat moguće je ustanoviti linearnu zavisnost između konstante kuplovanja J i razlike u energijama ε(g) i ε(u) ( Jednačina 4.51). ( ) ABAA gu AB JJ KJ − − −= 2 2 εε Jednačina 4.51 Ovakva formulacija obuhvata osnovu Hoffman-ove teorije, koja će biti detaljnije objašnjena kroz praktične primere u daljem tekstu (videti odeljak 7.2.1). 4.4 DFT proračuni magnetnih kuplovanja u binuklearnim kompleksima prelaznih metala Kao i u slučaju izračunavanja magnetne anizotropije, za izračunavanje magnetnih osobina binuklearnih kompleksa, potrebno je uključiti doprinose ekscitovanih stanja. Da bi se rešio ovakav problem potrebno je koristiti modifikovane DFT metode, budući da je standardni DFT jedno-determinantni pristup. Među najviše korišćenim metodama svakako je Noodleman-ov pristup (Broken-Symmetry-BS approach)58,59,60. Ovakav pristup ima za cilj da poveže energiju BS determinante, koja ne predstavlja jedinstveno stanje spinskog ugaonog momenta (BS determinanta je jedinstveno stanje Ms vrednosti), sa energijom tačno određenog spinskog stanja. Da bi se pravilno razumeo BS pristup potrebno je prvo objasniti kako se konstruišu BS orbitale. Kao model sistem mogu poslužiti dve monomerne jedinice [CuBr4]2- postavljene na udaljenost, dovoljnu za međusobne magnetne interakcije (Slika 4.19). 43 Slika 4.19 Šematski prikaz orbitala lokalizovanih na pojedinačnim centrima čijom kombinacijom se dobijaju molekulske orbitale dimera Lokalizovane magnetne orbitale za svaki monomer prikazane na slici mogu se predstaviti u matematičkom obliku (Jednačina 4.52). ( ) ( ) ( )AAzAxA CucdBrbpBrap ++=Φ ( ) ( ) ( )BBzBxB CucdBrbpBrap ++=Φ Jednačina 4.52 Ovakve magnetne orbitale su međusobno neortogonalne, a njihovom linearnom kombinacijom dobijaju se delokalizovane, međusobno ortogonalne, molekulske orbitale dimera kao simetrična i antisimetrična kombinacija (Jednačina 4.53). ( ) ( )[ ] ( ) ( )[ ] ( ) ( )[ ]{ }BABzAzBxAx CudCudcBrpBrpbBrpBrpa S +++++ + =Ψ+ 222 1 ( ) ( )[ ] ( ) ( )[ ] ( ) ( )[ ]{ }BABzAzBxAx CudCudcBrpBrpbBrpBrpa S −+−+− − =Ψ− 222 1 Jednačina 4.53 Iz ovakvog seta molekulskih orbitala moguće je konstruisati par, delimično delokalizovanih, ortogonalnih magnetnih orbitala (OMO) (Jednačina 4.54). ( ) ( )[ ] ( ) ( )[ ] ( ) ( )[ ]{ } ( ) ( )[ ] ( ) ( )[ ] ( ) ( )[ ]{ }               −+−+− − ++++++ +=Φ BABzAzBxAx BABzAzBxAx MO A CudCudcBrpBrpbBrpBrpa S CudCudcBrpBrpbBrpBrpa S 2 2 1 1 1 1 2 1 44 ( ) ( )[ ] ( ) ( )[ ] ( ) ( )[ ]{ } ( ) ( )[ ] ( ) ( )[ ] ( ) ( )[ ]{ }               −+−+− − −+++++ +=Φ BABzAzBxAx BABzAzBxAx MO B CudCudcBrpBrpbBrpBrpa S CudCudcBrpBrpbBrpBrpa S 2 2 1 1 1 1 2 1 Jednačina 4.54 Kombinovanjem ovako dobijenih orbitala konstruišu se BS orbitale koje su parcijalno delokalizovane i neortogonalne (Jednačina 4.55). MO B MO A BS B MO B MO A BS A xy yx Φ+Φ=Φ Φ+Φ=Φ Jednačina 4.55 Šematski prikaz BS orbitala dat je na slici 4.20. Slika 4.20 Šematski prikaz BS orbitala Elektroni smešteni u različitim BS orbitalama imaju različite spinove, odnosno α i β spinske gustine lokalizovane su na različitim paramagnetnim centrima, sa delimičnom delokalizacijom na drugi centar. Drugim rečima dve različite orbitale u koje se smeštaju elektroni različitih spinova, lokalizovane najvećim delom na jednom centru, sa maksimalnom lokalnom ms vrednosti, interaguju međusobno tako da daju minimalnu rezultujuću Mstot vrednost. Ovakav pristup omogućava uvođenje „repova“, odnosno jonskog singleta, kao i interakciju kovalentnog i jonskog singleta što predstavlja osnov 45 antiferomagnetnih doprinosa ukupnom kuplovanju. Talasna funkcija sistema može se predstaviti preko BS determinante (Jednačina 4.56). ( ) ( )[ ]ββββαααα nABnBABS bbbyxaaayxMN ......1 ! 1 2121 Φ+ΦΦ+Φ=Ψ Jednačina 4.56 Kako se stanja sa maksimalnom SZ vrednosti mogu prikazati jednom determinantom, izračunavanje energije stanja najvišeg multipliciteta ne predstavlja problem u okviru DFT- a. Kada se izračunaju energije stanja sa najvišim multiplicitetom (High Spin Energy-EHS) i energije BS stanja (EBS), konstanta kuplovanja može se dobiti na osnovu Yamaguchi-jevog pristupa61 ( ) BSHS BSHS SS EEJ 22 − − = (Jednačina 4.57). Jednačina 4.57 Prednosti BS metode su u tome što se rezultati mogu dobiti u kratkom vremenskom roku i sa zadovoljavajućom preciznošću. Osnovni nedostatak je nemogućnost da se izračunaju doprinosi dvostrukog transfera elektrona sa liganda na metale, koji iako mali mogu igrati odlučujuću ulogu kada su ostali doprinosi zanemarljivi62. 46 5 Magnetna anizotropija u mononuklearnim kompleksima Nagli razvoj molekulskog magnetizma započeo je sa otkrićem da pojedini polinukelarni kompleksi prelaznih metala pokazuju veoma dugo vreme za relaksaciju magnetizacije. Ovakve osobine su iskorišćene za skladištenje informacija na molekulskom nivou63, kao i za izradu kvantnih računara64 2 maxDSUeff = . Osnovni uslov koji je neophodan za postojanje dugog vremena relaksacije magnetizacije jeste postojanje energetske barijere koja razdvaja stanja sa pozitivnim i negativnim vrednostima magnetnih momenata u odnosu na izabranu osu magnetizacije. Energija barijere najlakše se može predstaviti jednačinom 5.1: Jednačina 5.1 , gde D predstavlja aksijalni parametar cepanja, a Smax Nadalje, potrebno je dizajnirati i sintetisati molekule i sisteme molekula koji poseduju tražene osobine. Upotrebom različitih metoda kvantne hemije, moguće je predvideti i objasniti uticaje geometrije i elektronske konfiguracije na magnetnu anizotropiju sistema totalni spin molekula. 53,54. U tom cilju vrednosti magnetne anizotropije su izračunate LF- DFT metodom za seriju model-kompleksa [CrF5]2-, [CrCl5]2-, [CrBr5]2-, [CrI5]2-, [CrCl2F3]2-,[VF5]3-, [VCl5]3-, [VBr5]3-, [FeF5]3-, [FeCl5]3-, [FeBr5]3-, [NiF5]3-, [NiCl5]3-, [NiBr5]3-,[NiCl2F3]2-, [NiBr2F3]2- i [NiI2F3]2-, kao i na prethodno sintetisanom realnom jedinjenju [NiCl3(Hdabco)2]+ (dabco=1,4-diazabicyclo[2.2.2]-octane)65,66 i upoređene sa podacima dobijenim veoma preciznim „ab-initio“ metodama kao što su CASSCF (complete active space self consistent field) i MRCI (multi reference configuration interactions) (Slike 5.1 i 5.2). Slika 5.1 Struktura ispitivanih model-kompleksa tipa MA2E3 (M=V(II), Cr(III), Fe(II), Ni(II); A=F-, Cl- , Br-, I-; E=F-, Cl-, Br-, I-) 47 Slika 5.2 Struktura ispitivanog kompleksa [NiCl3(Hdabco)2]+ sa jasno naznačenim geometrijskim parametrima i rasporedom orbitala (dabco=1,4-diazabicyclo[2.2.2]-octane) 5.1 Računski metod Geometrije model-kompleksa u D3h simetriji optimizovane su LDA funkcionalom25. Na optimizovanim geometrijama urađeni su AOC i SD proračuni OPBE67,68 funkcionalom. Geometrije [NiCl3(Hdabco)2]+ kompleksa u D3 i C2 simetrijama optimizovane su LDA, OPBE i B3LYP28,29 funkcionalima. Na pomenutim geometrijama urađeni su AOC i SD proračuni LDA, OPBE i B3LYP nivoom teorije. Za atome metala i atome nemetala korišćen je TZP (triple zeta bazis set sa polarizacionim funkcijama). U LF- DFT proračunima na model-sistemima korišćena je eksperimentalno dobijena konstanta spin-orbitalnog kuplovanja69 umanjena za orbitalni redukcioni faktor (orf<1). Za ispitivani realni kompleksni jon [NiCl3(Hdabco)2]+ korišćena je konstanta spin orbitalnog kuplovanja dobijena na osnovu MRCI proračuna uz orf=1. Odgovarajući Racah i LF parametri dobijeni u LFDFT koraku dati su u Prilogu. 48 5.2 Rezultati i diskusija Da bi se pravilno razumeli faktori koji određuju intenzitet magnetnih interakcija u mononuklearnim kompleksima, u okviru ovog dela prvo su analizirane magnetne osobine model-sistema, a zatim je urađena detaljna analiza na realnom sistemu [NiCl3(Hdabco)2]+ 5.2.1 Analiza magnetne anizotropije u model sistemima Izračunate vrednosti energija prelaza (Slike 4.6, 4.9 i 4.10) LF-DFT metodom, kao i poređenje sa dostupnim CASSCF rezultatima dobijenim u saradnji sa istraživačkom grupom iz Kantabrije (Španija) za ispitivane model-komplekse u D3h simetriji date su u tabeli 5.1. Kako LF-DFT daje sasvim zadovoljavajuće rezultate, a CASSCF iako veoma precizan metod zahteva veliko kompjutersko vreme, pojedini model-kompleksi koji su ovde uzeti u obzir nisu analizirani CASSCF metodom (prazna mesta u tabeli 5.1). Tabela 5.1 Izračunate energije prelaza za seriju model-kompleksa (cm-1 Struktura ) i poređenje sa CASSCF rezultatima dobijenim u saradnji sa istraživačkom grupom iz Kantabrije (Španija) LF-DFT CASSCF D a 3h, [CrF5] 138,00/344,00/495,00 2- 148,40/302,20/461,50 D3h, [CrCl5] 134,00/338,00/484,00 2- 141,30/288,30/441,00 D3h, [CrBr5] 129,00/327,00/467,00 2- 102,80/209,00/318,90 D3h, [CrI5] 143,00/295,00/453,00 2- / D3h, [CrCl2F3] 153,00/312,00/477,00 2- / D3h, [VF5] 70,00/254,00/330,00 3- / D3h, [VCl5] 80,00/268,00/357,00 3- / D3h, [VBr5] 76,00/258,00/343,00 3- / D3h, [FeF5] 78,00 3- 72,90 D3h, [FeCl5] 73,00 3- 74,20 D3h, [FeBr5] 71,00 3- 74,69 D3h, [NiF5] 471,00/1,047,00 3- 500,20/1181,00 D3h, [NiCl5] 430,00/948,00 3- 408,90/1052,50 D3h, [NiBr5] 407,00/895,00 3- / D3h, [NiCl2F3] 465,00/1030,00 2- / D3h, [NiBr2F3] 454,00/1018,00 2- / D3h, [NiI2F3] 449,00/1005,00 2- / aPablo García-Fernández, Departamento de Ciencias de la Tierra y Física de la Materia Condensada, Universidad de Cantabria, Avenida de los Castros, Santander, Spain 49 Analizirani model-kompleksi poseduju degenerisano osnovno stanje u D3h simetriji, tj. 4E' u slučaju kompleksa vanadijuma(II) i hroma(III), 5E" u slučaju kompleksa gvožđa(II) i 3E' za komplekse nikla(II) (videti teorijski deo odeljci 4.12. i 4.1.3). Degenerisana osnovna stanja indukuju Jahn-Teller-ovu distorziju usled čega se otklanja degeneracija osnovnog stanja i snižava simetrija iz D3h u C2v, uz energetsku stabilizaciju sistema. Na slici 5.3 šematski je prikazan efekat spin-orbitalnog kuplovanja i Jahn-Teller-ove distorzije na šest degenerisanih stanja 3E' multipleta. Slika 5.3 Šematski prikaz delovanja spin-orbitalnog kuplovanja i Jahn-Teller-ove distorzije na šest degenerisanih stanja 3 E' multipleta Kada se analiziraju efekti spin-orbitalnog kuplovanja, nezavisno od distorzije molekula (slučaj dole-levo), može se zaključiti da je energija Ueff koja razdvaja najniža stanja |SZ|=1 i SZ=0, znatno veća nego u slučaju kada je prisutna i distorzija molekula (slučaj dole-desno). Nadalje, spin-orbitalno kuplovanje uklanja degeneraciju nivoa sa |SZ|=1. Kako je Jahn- Teller elektronski efekat i ne zavisi od spinskih funkcija, moguća je interakcija samo onih stanja sa istim ms vrednostima. Na taj način mogu interagovati samo osnovno i ekscitovano stanje sa ms=1, kao i osnovno i ekscitovano stanje sa ms=-1, što predstavlja pseudo Jahn- Teller-ov efekat. Interakcije između stanja sa ms=0 predstavljaju pravi Jahn-Teller-ov efekat jer se degeneracija ovih nivoa ne otklanja spin-orbitalnim kuplovanjem. Veoma važna razlika između pomenutih efekata je u tome da Jahn-Teller uvek indukuje distorziju molekula, dok pseudo Jahn-Teller dovodi do distorzije samo ukoliko je dovoljno jak. Iz tog razloga može se zaključiti da spin-orbitalno kuplovanje sprečava distorziju molekula 40. 50 Kako su sve energije prelaza izračunate za Jahn-Teller aktivna jedinjenja u geometriji visoke simetrije moguće je ispitati efekte spin-orbitalnog kuplovanja nezavisno od distorzije (Tabela 5.1), (Slika 5.3 dole-levo). Dakle, razlike u izračunatim prelazima potiču samo od razlika u elektronskoj konfiguraciji i u tipu metala. Za sve pomenute komplekse u D3h simetriji, poredak molekulskih orbitala sa dominantnim d karakterom je isti (Slika 5.4). Slika 5.4 Poredak i orijentacija orbitala za ispitivane model sisteme Interakcije između dx2-y2 i dxy orbitala koje dovode do visokih vrednosti prelaza moguće su kod kompleksa vanadijuma(II), hroma(II) i nikla(II). Kod kompleksa gvožđa(II) dolazi do kuplovanja između dxz i dyz orbitala pa su vrednosti energija prelaza znatno manje iz razloga koji su objašnjeni u teorijskom delu (videti teorijski deo 4.1.3). Kompleksi vanadijuma(II) i hroma(III) imaju istu, d3 elektronsku konfiguraciju, pa se primetne razlike mogu objasniti samo na osnovu jačeg spin-orbitalnog kuplovanja kod hroma(III) (ξCr(III)=273 cm-1, ξV(II)=167 cm-1)69. Najveće vrednosti energija prelaza dobijene su za model komplekse nikla(II), zbog izuzetno izraženog spin-orbitalnog kuplovanja (ζNi(II)=649 cm-1). Energija prvog prelaza u slučaju kompleksa nikla(II) odgovara D vrednosti, što je u skladu sa razmatranjima u teorijskom delu 4.1 (videti slike 4.1 i 4.7). 51 5.2.2 Analiza magnetne anizotropije za [NiCl3(Hdabco)2]+ Na osnovu detaljne analize model-sistema ustanovili smo da distorzija molekula deluje nasuprot spin orbitalnom kuplovanju, pa je bilo neophodno odabrati sisteme sa krutim helatnim ligandima ili ligandima koji zbog svoje voluminoznosti otežavaju distorziju. Detaljnom pretragom kristalografske banke podataka kompleksni jon 70, pronašli smo mogući odgovarajući sistem-[NiCl3(Hdabco)2]+. Izračunate vrednosti energija prelaza i relevantni geometrijski parametri za [NiCl3(Hdabco)2]+, u D3 i C2 geomerijama na različitim nivoima teorije prezentovani su u tabeli 5.271 Tabela 5.2 Geometrijski parametri za ispitivani kompleks [NiCl . 3(Hdabco)2]+, izračunate vrednosti Jahn-Teller stabilizacione energije (cm-1) i D vrednosti dobijene LF-DFT i MRCI metodom (cm-1 Metod ) CF R(NiCl1 R(NiCl) 2 R(NiN) ) Cl1NiCl E2 UJT ULF-DFT MRCI D 2,276 3 2,276 2,078 120,0 / 554 645 LDA C2 2,271 (min) 2,268 2,086 114,3 412 159 314 C2 2,270 (TS) 2,275 2,085 125,1 355 179 370 D 2,307 3 2,307 2,295 120,0 / 542 641 OPBE C2 2,331 (min) 2,308 2,220 114,9 254 182 358 C2 2,309 (TS) 2,350 2,218 124,6 215 204 372 D 2,337 3 2,337 2,281 120,0 / 530 643 B3LYP C2 2,354 (min) 2,334 2,266 116,0 253 243 411 C2 2,338 (TS) 2,365 2,269 123,7 217 370 426 Exp min 2,309 2,281 2,227 118,9 / / / Dužine veza su reprodukovane sa zadovoljavajućom tačnošću na sva tri nivoa teorije, sa maksimalnom devijacijom od 0,04 Å (LDA), 0,022 Å (OPBE) i 0,045 Å (B3LYP) za R(Ni-Cl1); 0,013 Å (LDA), 0,027 Å (OPBE) i 0,053 Å (B3LYP) za R(Ni-Cl2); 0,141 Å (LDA), 0,007 Å (OPBE) i 0,039 Å (B3LYP) za R(Ni-N). Najveće odstupanje od eksperimentalnih vrednosti primećeno je za R(Ni-N) u geometrijama optimizovanim LDA nivoom teorije. U svim slučajevima R(Ni-Cl1)>R(Ni-Cl2), što je u saglasnosti sa eksperimentalno određenim vrednostima. Vrednosti izračunatih uglova slažu se sa podacima dobijenim difrakcijom x zraka sa greškom od 4.6º (LDA), 4.0º (OPBE) i 2.9º 52 (B3LYP). Izračunate vrednosti energija prelaza u geometriji minimuma su velike i u saglasnosti sa MRCI proračunima. Veoma je važno napomenuti da, iako je magnetna anizotropija u minimumu izražena, iznosi samo 40% od vrednosti u geometriji visoke simetrije gde nije došlo do distorzije. Osnovno stanje [NiCl3(Hdabco)2]+ jona (d8 elektronska konfiguracija) je degenerisano (3E), pa je ova struktura podložna Jahn-Teller-ovoj distorziji. Usled distorzije dolazi do sniženja simetrije iz D3 u C2. Aktivni vibracioni modovi su kvalitativno slični modovima za triangularne molekule kao što je CuF3 72. Međutim, distorzije ugla za takve molekule, znatno su izraženije u odnosu na distorziju ugla u [NiCl3(Hdabco)2]+ Nadalje, poznato je da spin-orbitalno kuplovanje uklanja degeneraciju nivoa sa S kompleksu. U geometriji visoke simetrije svaki hlorov atom, okružen je sa četiri vodonikova atoma aksijalnih liganada. Usled jakih sternih zaklanjanja distorzija ugla nije energetski povoljna. z = ±1, usled čega je distorzija molekula u osnovnom stanju 3E±1 indukovana pseudo Jahn-Teller- ovim efektom. Kako je pseudo Jahn-Teller-ov efekat slabije izražen u odnosu na Jahn- Teller-ov efekat, mogu se očekivati mala odstupanja od geometrije visoke D3 U cilju detaljnije analize MRCI proračuni su takođe odrađeni za odabrane geometrije duž puta najstrmijeg spusta, IDP-a (Intrinsic Distortion Path) simetrije. 42. (Slika 5.5). Osnovni cilj bio je da se utvrdi da li minimum za osnovno stanje 3E±1 odgovara visokoj simetriji ili distorgovanoj C2 geometriji, uzimajući u obzir spin-orbitalno kuplovanje nikla(II). Slika 5.5 Površina potencijalne energije dobijena na osnovu MRCI proračuna za različite tačke duž IDP-a idući od visoke D3 simetrije ka C2 minimumu i prelaznom stanju 53 MRCI proračuni duž IDP-a su pokazali da nema odstupanja od konfiguracije visoke simetrije (α Cl1NiCl2 Na osnovu svega pomenutog može se reći da je distorzija u osnovnom stanju u potpunosti „blokirana“ spin-orbitalnim kuplovanjem. =120º). Imajući u vidu da je eksperimentalna vrednost ugla od 118.9º određena u čvrstoj fazi, gde se molekuli pakuju sa okruženjem niske simetrije, oba rezultata (eksperimentalni i teorijski) u potpunoj su saglasnosti. Uzimajući u obzir gore navedeno, kao i eksperimentalne vrednosti uglova, prelaze ne treba računati u geometriji DFT minimuma, već u strukturi visoke simetrije gde oni imaju veoma visoke vrednosti. Temperatura za prelazak barijere iznosi 925K, što je znatno više od sobne temperature i tri puta više od prethodno zabeleženih vrednosti za druge komplekse (323K)73 . Prisustvo tako visoke barijere može predstavljati problem za aplikaciju, odnosno promena spina može biti otežana. Rešenje ovog problema može se naći u primeni električnog polja, koje favorizuje Jahn-Teller-ovu distorziju, što za posledicu ima smanjenje barijere. Nadalje, za razliku od oktaedarskih kompleksa gde su d-d prelazi simetrijski zabranjeni, to nije slučaj sa trigonalnobipiramidalnim kompleksima. Kombinovanim dejstvom magnetnog polja i svetlosti, može doći do foto-reakcije, odnosno prelaska u ekscitovano stanje, a zatim do relaksacije do osnovnog stanja sa suprotnim spinom. 5.3 Zaključak U okviru ovog dela prezentovana je detaljna anliza magnetnih svojstava trigonalno bipiramidalnih model-kompleksa i sintetisanog [NiCl3(Hdabco)2]+. Korišćenjem LF-DFT pristupa određene su energije prelaza za seriju model-kompleksa, kao i D vrednost u [NiCl3(Hdabco)2]+ molekulu. Uzrok velike D vrednosti u pomenutom kompleksu nikla, koja je veoma blizu konstanti spin-orbitalnog kuplovanja za nikal predstavlja degeneracija parcijalno popunjenih orbitala (dx2-y2 i dxy). Prisustvo voluminoznih aksijalnih liganada kao i spin-orbitalno kuplovanje umanjuju efekte Jahn-Teller-ove distorzije. U cilju daljeg istraživanja potrebno eksperimentalno potvrditi ovako visoku vrednost magnetne anizotropije na ispitivano sistemu [NiCl3(Hdabco)2]+. 54 6 DFT analiza magnetnih interakcija u binuklearnim kompleksima mangana Magnetna svojstva komleksa prelaznih metala koji sadrže više od jednog paramagnetnog centra međusobno premošćena višeatomskim ligandima razlukuju se u odnosu na mononuklearne komplekse i veoma često se ne mogu predstaviti sumom magnetnih osobina pojedinačnih jedinica. Utvrđeno je da magnetne osobine binuklearnih kompleksa zavise od: udaljenosti između metalnih centara, uglova veza između mostnih atoma, torzionih uglova između ravni u kojima se nalaze metalni centri, dužine veza i koordinacionog okruženja oko metala. Sposobnost dicijanamidnih liganada (dicyd2-) da prenose magnetne interakcije, nagoveštava da i fenilcijanamidni derivati (pcyd) takođe mogu imati slična svojstva74,75,76. Na osnovu Extended-Hückel MO proračuna predloženo je da µ1,3-NCN pcyd ligandi mogu imati svojstva koja se nalaze na prelazu između azidnih i dicijanamidnih medijatora magnetnih interakcija75. U svrhu daljih ispitivanja i radi potvrde Extended-Hückel MO rezultata, Escuer et al. su sintetisali binuklearne komplekse mangana (II) sa pcyd mostnim ligandima i 1,10-fenantrolinskim blokirajućim ligandima77 Da bi se detaljnije shvatili uticaji konformacije mostnih liganada na magnetne osobine, primenjen je DFT-BS pristup (videti odeljak 4.4) na seriji binuklearnih kompleksa mangana sa fenilcijanamidnim ligandima: [{Mn(3-Fpcyd)(MeOH)(phen)} . Ovakvi molekuli formiraju supramolekulske jednodimenzionalne mreže kod kojih su binuklearne monomerne jedinice međusobno povezane vodoničnim vezama. 2(µ-3-Fpcyd)2], 1 [{Mn(3-Fpcyd)(EtOH)(phen)}2(µ-3-Fpcyd)2], 2 [{Mn(3- Clpcyd)(MeOH)(phen)}2(µ-3- Clpcyd)2], 3 [{Mn(4-lpcyd)(EtOH)(phen)}2(µ-4-Clpcyd)2] 4 (Slika 6.1) 77,78. 55 Slika 6.1 Strukture ispitivanih binuklearnih kompleksa: [{Mn(3-Fpcyd)(MeOH)(phen)}2(µ-3-Fpcyd)2], 1 [{Mn(3-Fpcyd)(EtOH)(phen)}2(µ-3-Fpcyd)2], 2 [{Mn(3- Clpcyd)(MeOH)(phen)}2(µ-3-Clpcyd)2], 3 [{Mn(4-lpcyd)(EtOH)(phen)}2(µ-4-Clpcyd)2 6.1 Računski metod ] 4 Da bi se izbegli problemi koji postoje pri optimizaciji binuklearnih kompleksa sa singletnim osnovnim stanjem (videti teorijski deo 4.4), strukture ispitivanih kompleksa 1-4 optimizovane su molekulskom mehanikom. Pri radu je korišćena 2010/PC verzija CFF (Consistent Force Field) konformacionog programa (videti spisak korišćenih programa). Konformaciona energija definisana je na sledeći način (Jednačina 6.1): ( ) ( ) ( ) ∑∑ ∑∑∑ +               −      + +++−+−= ∗∗ ij ij ji ij r r total Dr ee r r r r nkkrrkE 612 2020 3 2 cos1 2 1 2 1 2 1 εε φθθ φ φ θ θ Jedanačina 6.1 Prvi član jednačine predstavlja sumu svih energetskih doprinosa koji potiču od deformacije dužine veza. Korišćena je harmonijska funkcija za opis energije istezanja veza (r) (sa 56 konstantama sile kr i ravnotežnim dužinama veza r0). Valencioni potencijal koji potiče od deformacije uglova (θ) (sa konstantama sile kθ i optimalnim uglovima θ0) takođe se opisuje harmonijskom (kvadratnom) funkcijom. Treći član predstavlja sumu svih torzionih potencijala za svaku vezu, gde je kφ rotaciona barijera, a n je parametar koji definiše periodičnost 2π/n, tj. broj minimuma pri zaokretu za 360º oko ose veze. Sledeći član jednačine predstavlja doprinos nevezivnih Van der Waals-ovih interakcija ukupnoj energiji molekula, koji je računat primenom Lennard-Jones-ove ’12-6’ funkcije. Parametar ε* predstavlja dubinu Lennard-Jones-ove potencijalne jame, kojoj odgovara ravnotežno rastojanje r*. Dalje, elektrostatičke interakcije modelirane su Coulomb-ovom energijom interakcija između dva centara naelektrisanja qi i qj, koja se nalaze na rastojanju rij Polje sila u medijumu dialektrične konstante ε. Parcijalne šarže na atomima dobijene su iz DFT proračuna. 79,80,81,82,83,84,85,86,87,88,89 parametrizovano je na osnovu tri različita tipa ugljenikovih atoma (sp hibridizovani Csp atom cijanamidnih mostnih liganada, uglenikov atom aromatičnih prestenova Cp i sp3 hibridizovani Csp3 atom), jednog tipa kiseonikovog atoma, jednog tipa vodonikovog atoma, dva tipa azotovih atoma (sp2 hibridizovani N atom i atom azota iz aromatičnog fenantrolinskog prstena Nph Konstante kuplovanja izračunate su B3LYP ) i centralnog metalnog atoma (Mn(II)). Lista parametara polja sila data je u Prilogu. Optimizacija geometrije izvedena je kombinovanjem tri gradijentne minimizacione metode: metode najstrmijeg spusta (steepest- descent), Davidon-Fletcher-Powell i Newton-Raphson metoda. 28,29 funkcionalom na kristalnim i MM optimizovanim strukturama. Za atome metala korićen je TZVPP (triple zeta valence basis set sa polarizacionim funkcijama) bazis set, dok je za atome nemetala upotrebljen SVP (split valence basis set) bazis set90,91 . 6.2 Rezultati i diskusija Minimizacija energije i optimizacija geometrije kao rezultat daju jedan stabilan konformer za svaku strukturu ispitivanih kompleksa 1-4. Strukturni parametri za 57 minimizovane geometrije u saglasnosti su sa podacima dobijenim difrakcijom x zraka (Tabela 6.1). Tabela 6.1 Strukturni parametri; dužine veza (Å), uglovi (°) i torzioni uglovi (°) za ispitivane komplekse i poređenje sa eksperimentalnim vrednostima M-N M-O N-C NMN Nph CNM MO CNph CNM ph [{Mn( MN 3-Fpcyd)(MeOH)(phen)}2(µ-3-Fpcyd)2 X-ray ] (1) 2,147 2,228 1,162 103,00 82,25 157,7 114,2 12,90 MM 2,137 2,225 1,156 102,70 85,00 157,6 113,1 4,00 [{Mn(3-Fpcyd)(EtOH)(phen)}2(µ-3-Fpcyd)2 X-ray ] (2) 2,143 2,250 1,172 100,60 84,35 140,80 113,30 39,00 MM 2,142 2,260 1,166 99,40 86,65 142,90 111,00 46,80 [{Mn(3-Clpcyd)(MeOH)(phen)}2(µ-3-Clpcyd)2 X-ray ] (3) 2,236 2,220 1,222 95,00 85,10 146,60 122,00 28,50 MM 2,232 2,219 1,223 93,90 86,35 152,40 119,40 36,20 [{Mn(4-Clpcyd)(EtOH)(phen)}2(µ-4-Clpcyd)2 X-ray ] (4) 2,129 2,238 1,172 101,40 83,90 148,40 110,50 26,30 MM 2,132 2,230 1,171 100,60 86,70 152,70 113,00 31,90 Konstante kuplovanja su izračunate na kristalnim strukturama i MM optimizovanim geometrijama. Sve izračunate vrednosti su u saglasnosti sa eksperimentalno dobijenim podacima 77 (Tabela 6.2). Tabela 6.2 Izračunate konstante kuplovanja (cm-1 Struktura ) za seriju kompleksa 1-4 na kristalnim strukturama i optimizovanim geometrijama i poređenje sa eksperimentalnim podacima Jcalc / cm J-1 exp /cm J-1 MM/cm [{Mn(3-Fpcyd)(MeOH)(phen)} - 1 2(l-3-Fpcyd)2 -4,61 ] -3,34 -3,24 [{Mn(3-Fpcyd)(EtOH)(phen)}2(l-3-Fpcyd)2 -2,23 ] -0,41 -0,30 [{Mn(3-Clpcyd)(MeOH)(phen)}2(l-3-Clpcyd)2 -3,52 ] -2,89 -1,30 [{Mn(4-Clpcyd)(EtOH)(phen)}2(l-4-Clpcyd)2 -3,65 ] -2,60 -1,55 Proračunima je pokazano da su interakcije u kompleksima slabog antiferomagnetnog karaktera. Ovakav rezultat posledica je elektronske konfiguracije paramagnetnih centara, njihove udaljenosti i strukture mosta. Svaki metalni centar ima pet nesparenih elektrona 58 usled čega su izražene antiferomagnetne interakcije između neortogonalnih orbitala, ali i feromagnetne interakcije između ortogonalnih orbitala. Različiti doprinosi se međusobno poništavaju, pri čemu antiferomagnetni doprinos preovlađuje, dajući kao krajnji rezultat slabo antiferomagnetno kuplovanje preko fenilcijanamidnog mosta. Udaljenost između metala je ≈6Ǻ što dodatno otežava delokalizaciju elektrona. Faktori koji favorizuju antiferomagnetno kuplovanje povezani su sa distorzijom Mn–(NCN)2–Mn regiona, odnosno veće antiferomagnetno kuplovanje se očekuje kod sistema kod kojih je most manje distorgovan. Distorzija Mn–(NCN)2–Mn dela potiče od deformacije C-Nph -5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0 0 5 10 15 20 25 30 35 40 45 50 Torzioni ugao/stepeni J (c m -1 ) 1 2 3 -Mn-N torzionog ugla. Kompleks 1 pokazuje najmanju distorziju centralnog prstena u odnosu na ostale komplekse, dok je kod kompleksa 2 primećeno najveće odstupanje centralnog prstena od planarnosti. Trend izračunatih vrednosti u korelaciji je sa distorzijom mosta, odnosno kompleks 2 pokazuje najslabije antiferomagnetne interakcije (Slika 6.2). Slika 6.2 Zavisnost konstante kuplovanja od distorzije prstena: 1) eksperimentalne vrednosti 2) vrednosti izračunate na kristalnoj strukturi 3) vrednosti izračunate na MM optimizovanim geometrijama. Simbol (▲) označava kompleks 1, simbol (■) označava kompleks 4, simbol (●) označava kompleks 3, simbol (♦) označava kompleks 2 Eksperimentalne vrednosti pokazuju naglo povećanje antiferomagnetnih interakcija, u slučaju kompleksa 3 što nije u korelaciji sa distorzijom C-Nph-Mn-N torzionog ugla. Konstante izračunate na kristalnim strukturama pokazuju bolje slaganje sa 59 eksperimentalnim vrednostima u slučaju kompleksa 3 i 4, dok se za komplekse 1 i 2 bolji rezultati dobijaju na MM optimizovanim strukturama. Da bi se pravilno objasnio uticaj distorzije prstena na jačinu magnetnih interakcija, potrebno je analizirati magnetne orbitale BS stanja ispitivanih kompleksa (Slika 6.3). Slika 6.3 Magnetne orbitale u BS stanju za analizirane komplekse. Za svaku od četiri ispitivane strukture 1-4 prikazane su samo orbitale a)-d) kod kojih je preklapanje sa orbitalama liganada najizraženije Orbitale ispitivanih kompleksa u BS stanju su lokalizovane na pojedinačnim centrima, bez doprinosa orbitala sa drugog centra. U izvesnoj meri elektronska gustina je delokalizovana sa metala na orbitale mostnih liganada koje su prilagođene za preklapanje. Antiferomagnetne interakcije se dakle ispoljavaju usled delokalizacije preko orbitala mostnih liganada. Preklapanje d orbitala metala sa orbitalama liganada kao i međusobno 60 preklapanje orbitala liganada najjače je izraženo kada se orbitale nalaze u ravni prstena. Distorzija prstena dovodi do narušavanja delokalizacije, a samim tim do smanjenja jačine kuplovanja. 6.3 Zaključak U okviru ovog dela analizirana su magnetna kuplovanja u binuklearnim kompleksima mangana(II) sa fenilcijanamidnim mostnim ligandima DFT-BS pristupom. Da bi se izbegli problemi pri optimizaciji binuklearnih kompleksa sa singletnim osnovnim stanjem, strukture analiziranih kompleksa optimizovane su metodom molekulske mehanike (MM). U svim ispitivanim kompleksima prisutno je slabo antiferomagnetno kuplovanje, što je u saglasnosti sa eksperimentalno dobijenim konstantama kuplovanja. Jačina interakcija opada sa distorzijom mosta, odnosno sa povećanjem CNphMN torzionog ugla. Delokalizacija elektrona na mostne atome najjača je kada se orbitale metala i liganada nalaze u istoj ravni. Distorzijom mosta dolazi do izlaženja orbitala liganada iz ravni maksimalnog preklapanja, usled čega slabe magnetne interakcije između paramagnetnih centara. 61 7 DFT analiza magnetnih interakcija u binuklearnim oksalatnim kompleksima prve serije prelaznih metala Veliki broj naučnih radova u oblasti magnetizma u proteklih par godina, posvećen je analizi oksalatnih mostova kao medijatorima magnetnih kuplovanja između metala prve prelazne serije u homo i heteropolinuklearnim kompleksnim jedinjenjima92,93,94,95 [19-22]. Primećeno je da oksalatni anjoni mogu da formiraju homo- i heterometalne 2-D96 i 3-D97 mreže koje imaju primenu kao molekulski magneti i magnetni materijali. Iz tog razloga u ovom delu će detaljno biti analizirane magnetne osobine DFT-BS (videti odeljak 4.4) pristupom za seriju binuklearnih kompleksa sa oksalatnim mostovima: pet izomera [(VO)2(ox)(SCN)6]4-, trans-(equatorial,equatorial) 1, cis-(equatorial,equatorial) 2, trans- (axial,axial) 3, cis-(axial,axial) 4, (axial,equatorial) 5, kao i kompleksi [Cr2(ox)(SCN)8]4- 6, [Fe2(ox)(SCN)8]4- 7, [CrFe (ox)(SCN)8]4- 8, [Fe2(ox)5]4- 9, [Cr2(ox)5]4- 10, [Ni2(ox)5]6- 11, [Cu2(ox)(C12H8N2)2]2+ 12 (Slika 7.1)98. Različiti izomeri [(VO)2(ox)(SCN)6]4- ispitivani su u cilju određivanja relacije između konstante kuplovanja J i orijentacije dvostruke M=O veze99,100. Takođe su ispitani uticaji različitih funkcionala na geometrije i magnetne osobine ispitivanih binuklearnih kompleksa. Pored navedenog analizirani su i uticaji elektronske konfiguracije i tipa metala na magnetna kuplovanja. 62 Slika 7.1 Strukture ispitivanih kompleksa[(VO)2(ox)(SCN)6]4-, trans-(equatorial,equatorial) 1, cis- (equatorial,equatorial) 2, trans-(axial,axial) 3, cis-(axial,axial) 4, (axial,equatorial) 5, kao i kompleksi [Cr2(ox)(SCN)8]4- 6, [Fe2(ox)(SCN)8]4- 7, [CrFe (ox)(SCN)8]4- 8, [Fe2(ox)5]4- 9, [Cr2(ox)5]4- 10, [Ni2(ox)5]6- 11, [Cu2(ox)(C12H8N2)2]2+ 12 7.1 Računski metod Geometrije ispitivanih kompleksa optimizovane su LDA funkcionalom. Dobijene geometrije su dodatno optimizovane OPBE i PW91101 funkcionalima u cilju poređenja rezultata. Za atome metala korićen je TZVPP (triple zeta valence basis set sa polarizacionim funkcijama) bazis set, dok je za atome nemetala upotrebljen SVP (split valence basis set) bazis set. Konstante kuplovanja su izračunate na dostupnim kristalnim i 63 optimizovanim strukturama, B3LYP hibridnim funkcionalom. U cilju poređenja rezultata konstante su takođe izračunate LDA, OPBE, i PW91 nivoima teorije. 7.2 Rezultati i diskusija Strukture ispitivanih kompleksa optimizovane su BS metodom u cilju dobijanja geometrija antiferomagnetnog stanja. Sve optimizovane strukture predstavljaju minimume na površini potencijalne energije, što je potvrđeno odsustvom vibracionih modova sa negativnim konstantama sila. Strukturni parametri minimizovanih geometrija kao i dostupni kristalografski podaci predstavljeni su u tabeli 7.1. Kod svih molekula LDA funkcional daje kraće vrednosti za dužine veza u odnosu na OPBE i PW91 funkcionale sa izuzetkom M=O veze u kompleksima vanadijuma 1, 2, 3, 4, i 5. Bolje slaganje u M-Oox i M=O dužinama veza između optimizovane i kristalne strukture kompleksa 1, kao i u slučaju M-Oox veze u kompleksu 6, postignuto je LDA funkcionalom. U binuklearnim kompleksima 1 i 2 M-Oox veza koja se nalazi u trans položaju u odnosu na M=O vezu, duža je u odnosu na ostale M-Oox veze usled izraženog trans efekta. Posmatranjem optimizovanih struktura 3 i 4 može se ustanoviti da su M-Oox veze u tim kompleksima, kraće u odnosu na trans M-Oox veze i neznatno duže od cis M- Oox veza u kompleksima 1, 2, i 5. Duže veze u binuklearnim kompleksima 7, 9, i 11 mogu se očekivati usled smeštanja dodatnih elektrona u antivezivne M-L orbitale (σ* Konstante kuplovanja su izračunate sa B3LYP funkcionalom na LDA, OPBE i PW91 optimizovanim geometrijama, kao i na dostupnim kristalnim strukturama binuklearnih kompleksa 1, 6, i 12 ). Najkraće veze su dobijene u slučaju strukture 12 kao posledica izražene kovalentnosti metal-ligand veze. 102,103,104. Sve izračunate konstante su u dobroj saglasnosti sa eksperimentalno dobijenim rezultatima (Tabela 7.2)102,103,105,106. 64 Tabela 7.1 Strukturni parametri za ispitivane komplekse; dužine veza (Å), uglovi (°) i poređenje sa eksperimentalnim vrednostima Struktura M-O M-O ox M-N OoxMOox transNMOox transOMO [(VO) ox 2(ox)(SCN)6]4- (1) X-ray 2,053/2,237 1,582 a 2,054 76,20 164,50 171,60 LDA 2,023/2,237 1,586 a 2,004 75,50 165,30 170,20 OPBE 2,077/2,368 1,579 a 2,085 73,20 164,70 168,70 PW91 2,075/2,308 1,599 a 2,067 74,50 164,80 169,60 [(VO)2(ox)(SCN)6]4- (2) LDA 2,026/2,222 1,587 a 2,005 75,60 166,50 169,40 OPBE 2,109/2,313 1,580 a 2,085 73,60 167,00 166,85 PW91 2,092/2,278 1,600 a 2,067 74,60 166,35 168,25 [(VO)2(ox)(SCN)6]4- (3) LDA 2,066 1,588 2,054 79,10 167,50 / OPBE 2,129 1,578 2,158 77,50 168,00 / PW91 2,115 1,600 2,124 78,30 167,80 / [(VO)2(ox)(SCN)6]4- (4) LDA 2,069 1,590 2,052 79,00 168,10 / OPBE 2,134 1,581 2,152 77,40 168,92 / PW91 2,119 1,602 2,119 78,20 168,40 / [(VO)2(ox)(SCN)6]4- (5) LDA 2,051/2,242 1,587 a 2,031 77,25 167,00 169,50 OPBE 2,110/2,375 1,579 a 2,124 75,40 167,05 167,30 PW91 2,101/2,312 1,599 a 2,097 76,35 166,75 168,60 [Cr2(ox)(SCN)8]4- (6) X-ray 2,017 / 2,007 82,30 173,40 / LDA 2,016 / 1,925 80,80 173,90 / OPBE 2,084 / 1,989 79,20 173,10 / PW91 2,069 / 1,983 79,90 173,30 / [Fe2(ox)(SCN)8] (7) 4- LDA 2,105 / 2,005 77,20 169,80 / OPBE 2,172 / 2,080 75,50 168,70 / PW91 2,159 / 2,061 76,30 169,20 / [CrFe (ox)(SCN)8]4- (8) X-ray 2,177/2,056 / b 2,070/1,992 / b 168,00/173,05 / b LDA 2,135/1,992 / b 2,000/1,928 76,80/80,90b 169,40/174,00b / b OPBE 2,205/2,059 / b 2,069/1,996 74,90/79,60b 168,20/173,50b / b PW91 2,190/2,045 / b 2,053/1,989 75,80/80,30b 168,70/173,60b / b 65 Tabela 7.1 Nastavak [Fe2(ox)5] (9) 4- X-ray 2,092c / /1,977 / 80,80/168,70 / / LDA 2,101c / /1,989 / 79,00/166,60 / / OPBE 2,195c / /2,047 / 77,00/164,70 / / PW91 2,165c / /2,037 / 78,10/165,50 / / [Cr2(ox)5]4- LDA (10) 2,019c / /1,931 / 81,60/172,70 / / OPBE 2,124c / /1,990 / 77,43/171,23 / / PW91 2,089c / /1,980 / 80,80/171,90 / / [Ni2(ox)5] (11) 6- LDA 2,118c / /2,054 / 78,30/169,70 / / OPBE 2,304c / /2,169 / 74,03/164,90 / / PW91 2,218c / /2,135 / 76,30/167,40 / / [Cu2(ox)(C12H8N2)2]2+ (12) X-ray 1,973 / 1,985 85,20 168,00 / LDA 1,943 / 1,933 85,90 179,80 / OPBE 2,009 / 1,999 83,40 179,80 / PW91 1,995 / 1,989 84,50 179,80 / aM-Oox veza usmerena trans u odnosu na M=O vezu; bCr-Oox, Cr-N dužine veza i odgovarajući uglovi; cM-Oox veze u oblasti mosta 66 Tabela 7.2 Izračunate vrednosti konstanti kuplovanja ( cm-1 J ) B3LYP funkcionalom na kristalnoj strukturi i geometrijama optimizovanim LDA, OPBE i PW91 funkcionalima Jexp JX-ray/B3LYP JLDA/B3LYP JOPBE/B3LYP [(VO) PW91/B3LYP 2(ox)(SCN)6]4- -3,85 (1) -2,71 -8,34 -3,95 -4,67 [(VO)2(ox)(SCN)6]4- / (2) / -55,47 -28,07 -31,47 [(VO)2(ox)(SCN)6]4- / (3) / -0,96 -1,40 -1,21 [(VO)2(ox)(SCN)6]4- / (4) / 0,37 -0,68 -0,44 [(VO)2(ox)(SCN)6]4- / (5) / -0,95 -0,74 -0,95 [Cr2(ox)(SCN)8]4- -3,23 (6) -3,94 -5,68 -3,86 -4,14 [Fe2(ox)(SCN)8]4- -3,84 (7) / -5,47 -4,80 -4,86 [CrFe (ox)(SCN)8]4- 1,10 (8) / 0,89 1,45 1,40 [Fe2(ox)5]4- -6,6 (9) / -4,50 -3,60 -3,84 [Cr2(ox)5]4- -6,2 (10) / -6,82 -3,39 -4,15 [Ni2(ox)5]6- -22,8 (11) / -17,45 -13,10 -10,55 [Cu2(ox)(C12H8N2)2]2+ / (12) -351,98 -348,79 -337,97 -330,92 Strukture 2, 3, 4 i 5 predstavljaju izomere koji nisu izolovani, tako da nema dostupnih eksperimentalnih podataka za poređenje. Kao što je napomenuto cilj nam je bio da proučimo vezu između pojedinačnih izomera i ukupne konstante kuplovanja. U cilju dobijanja bolje slike i verodostojnosti razultata konstante su izračunate i LDA, OPBE, i PW91 nehibridnim funkcionalima (Tabela 7.3). 67 Tabela 7.3 Kontante kuplovanja( cm-1 J ) izračunate nehibridnim funkcionalima na geometrijama optimizovanim odgovarajućim funkcionalima Jexp JLDA JOPBE [(VO) PW91 2(ox)(SCN)6]4- -3,85 (1) -11,97 -4,19 -6,51 [(VO)2(ox)(SCN)6]4- / (2) -229,30 -111,92 -138,59 [(VO)2(ox)(SCN)6]4- / (3) -7,16 -6,77 -4,18 [(VO)2(ox)(SCN)6]4- / (4) -1,54 -1,97 0,85 [(VO)2(ox)(SCN)6]4- / (5) 3,85 1,07 1,42 [Cr2(ox)(SCN)8]4- -3,23 (6) -15,78 -11,34 -11,53 [Fe2(ox)(SCN)8]4- -3,84 (7) -16,34 -10,73 -13,37 [CrFe (ox)(SCN)8]4- 1,10 (8) 3,63 6,92 3,37 [Fe2(ox)5]4- -6,6 (9) -10,04 -7,08 -9,24 [Cr2(ox)5]4- -6,2 (10) -14,97 -8,40 -10,51 [Ni2(ox)5]6- -22,8 (11) -129,02 -56,78 -59,86 [Cu2(ox)(C12H8N2)2]2+ / (12) -1747,66 -1267,36 -1253,85 Sa izuzetkom kompleksa 8 gde su primećene feromagnetne interakcije, u svim ostalim strukturama dominantno je antiferomagnetno kuplovanje (Tabela 7.2). Veoma slabe feromagnetne interakcije primećene su za LDA optimizovanu strukturu 4, dok su slabe antiferomagnetne interakcije dobijene za OPBE i PW91 optimizovane geometrije pomenute strukture. Najbolje slaganje sa eksperimentalnim rezultatima za komplekse 8, 9, 10, i 11 postignuto je na LDA geometrijama, dok je za ostale komplekse bolje slaganje dobijeno na OPBE i PW91 geometrijama. Konstante izračunate sa nehibridnim funkcionalima jako su negativne u slučaju kompleksa 2, 11 i 12, dok su za komplekse 1, 8, 9 i 10 slaganja sa eksperimentom znatno bolja. 68 Da bi se stekao bolji uvid u mehanizam kuplovanja kod kompleksnih jedinjenja 1-5, potrebno je detaljnije analizirati molekulske orbitale za pomenute sisteme. Šematski prikaz molekulskih orbitala u tripletnom stanju za strukture 1-5 predstavljen je na slici 7.2. Slika 7.2 Magnetne orbitale u tripletnom stanju za različite izomere kompleksa vanadijuma Slabljenje antiferomagnetnih interakcija u jedinjenjima 3 i 4 u poređenju sa izomernim strukturama 1 i 2 može se smatrati posledicom pomeranja ekvatorijalne M=O veze u aksijalni položaj. U slučaju ekvatorijalne orijentacije M=O veze u kompleksima 1 i 2 d orbitale metala su orijentisane normalno na ravan oksalatnog mosta. Ovakva orijentacija omogućava povoljna preklapanja sa pz orbitalama mostnih atoma. Doprinos pz orbitale sa oksalatnog ugljenikovog atoma, znatno je veći u strukturi 2 u poređenju sa ostalim kompleksima. Magnetne interakcije dva metalna centra u strukturi 2 prostiru se preko pet atoma, pa je samim tim kuplovanje znatno olakšano u odnosu na kompleks 1 kod koga je put kuplovanja duži. U slučaju aksijalnog položaja M=O veze d orbitale metala se nalaze u ravni oksalatnog mosta. Preklapanje sa orbitalama u ravni je izrazito slabo, usled čega 69 dolazi do znatnih redukcija antiferomagnetnog kuplovanja. Za ekvatorijalno-aksijalni položaj M=O veza u kompleksu 5, karakterističa je prostorna ortogonalnost d orbitala, zbog čega se mogu očekivati feromagnetna kuplovanja. Međutim, konstante izračunate B3LYP funkcionalom neznatno su negativne, dok se primetna feromagnetna kuplovanja dobijaju nehibridnim funkcionalima. Iako se B3LYP nivo teorije smatra metodom izbora za izračunavanje magnetnih osobina1 Analizom spinskih gustina može se steći uvid o delokalizaciji elektrona (Tabela 7.4). 07, u ovom slučaju se bolji trendovi dobijaju nehibridnim funkcionalima. Tabela 7.4 Izračunate spinske gustine za komplekse 1-4 na mostnim atomima ugljenika i atomima vanadijuma za različite izomere kompleksa 1-4 Struktura LDA/B3LYP OPBE/B3LYP PW91/B3LYP ρ/C ρ/V ox ρ/C ρ/V ox ρ/C ρ/V ox [(VO)2(ox)(SCN)6]4- (1) 0,008137 0,008138 1,095527 0,004948 0,004948 1,136328 0,005717 0,005717 1,136716 [(VO)2(ox)(SCN)6]4- (2) 0,020930 1,098099 0,012729 1,135605 0,014374 1,140828 [(VO)2(ox)(SCN)6]4- (3) 0,002921 0,002770 1,121096 0,002389 0,002174 1,149023 0,002583 0,002560 1,159604 [(VO)2(ox)(SCN)6] (4) 4- 0,002592 0,002850 1,120923 0,002061 0,002235 1,149926 0,002263 0,002475 1,159615 Generalno govoreći veća delokalizacija molekulskih orbitala na atome liganada povezana je sa manjom spinskom gustinom na vanadijumovim atomima. U slučaju kompleksa 1 i 2 vrednosti spinskih gustina na ugljenikovim mostnim atomima veće su nego za iste atome kod struktura 3 i 4. Nadalje, manje vrednosti spinskih gustina na atomima vanadijuma za ekvatorijalni položaj M=O veze, ukazuju na intenzivniju delokalizaciju. Najmanje vrednosti spinskih gustina na mostnim ugljenikovim atomima dobijene su za kompleks 4, koji ujedno pokazuje i najslabija antiferomagnetna kuplovanja. Izračunate i eksperimentalno dobijene konstante kuplovanja za binuklearne komplekse hroma su negativne, što ukazuje na dominantni antiferomagnetni doprinos. I u ovom slučaju hibridni funkcional daje bolje rezultate. U poređenju sa kompleksima vanadijuma koji imaju jedan nespareni elektron po metalnom centru, dodatak dva elektrona po metalu u ispitivanim kompleksima ne dovodi do značajnijih promena u konstantama kuplovanja. 70 Ova pojava se lako može objasniti činjenicom da porast broja nesparenih elektrona pojačava kako antiferomagnetne, tako i feromagnetne interakcije. Kuplovanje neortogonalnih orbitala dxz→dxz, dyz→dyz, i dxy→dxy vodi ka antiferomagnetnim interakcijama, dok feromagnetne interakcije potiču od kuplovanja ortogonalnih orbitala dxz→dyz, dxz→dxy, i dxy→dyz. Gvožđe (III) jon u kompleksima 7, 8, i 9 ima visoko spinsku konfiguraciju, pa se dva dodatna elektrona nalaze u dx2-y2 i dz2 orbitalama. Dodatni antiferomagnetni (dx2-y2→dx2-y2, dz2→dz2) i feromagnetni (dx2-y2→dz2, dx2-y2→d(xz,yz,xy), dz2→d(xz,yz,xy)) doprinosi vode do umereno negativne vrednosti konstante kuplovanja u kompleksima 7 i 9, slično kao kod kompleksa hroma (III). Sa izuzetkom kompleksa 2, koji još uvek nije sintetisan, najjače antiferomagnetno kuplovanje javlja se u kompleksima 11 i 12. U heterobinuklearnom jedinjenju 8 dodatni feromagnetni doprinos potiče od polupopunjenih dx2-y2 i dz2 orbitala gvožđa (III) i praznih dx2-y2 i dz2 orbitala hrom (III) jona. Dakle, feromagnetno kuplovanje preko oksalatnog mosta se javlja kao dominantni efekat. Nespareni elektroni u jedinjenju 12 nalaze se u dx2-y2 orbitalama sa oba metalna centra, pa je jedini mogući put kuplovanja dx2-y2→dx2-y2. Veoma jake interakcije u binuklearnim kompleksima bakra mogu se objasniti preko σ preklapanja dx2-y2 orbitale metala sa simetrijski adaptitanim p orbitalama mostnog liganda. Odgovarajuće magnetne orbitale BS stanja u kompleksu 12 prikazane su na slici 7.3. Slika 7.3 Magnetne orbitale u BS stanju za kompleks 12 Orbitale metala usmerene su direktno ka orbitalama liganada, što doprinosi znatno jačem preklapanju. U slučaju niklovih kompleksa, jake antiferomagnetne interakcije potiču od σ preklapanja dz2 i dx2-y2 orbitala metala sa p orbitalama liganada. Feromagnetne interakcije između dz2 i dx2-y2 orbitala dovode do redukcije antiferomagnetnih kuplovanja u poređenju sa kompleksima bakra. 71 7.2.1 Određivanje linearne zavisnosti između konstante kuplovanja i (ε2-ε1)2 i S2 Kao dodatak jednostvanoj MO analizi i analizi spinskih gustina, magnetna kuplovanja u binuklearnim kopleksima sa jednim nesparenim elektronom po metalnom centru mogu se analizirati parametara Hoffmann-ovim108 [30] i Kahn-ovim11 [31] pristupom. Primena ova dva pristupa prikazana je za kompleks 1 kao ilustrativni primer. Kao što je napomenuto u teorijskom delu ukupna konstanta kuplovanja se može predstaviti kao suma doprinosa (Jednačina 4.51). Kod Hoffmann-ovog pristupa dvoelektronski integrali su aproksimativno konstantni za male promene u strukturnim parametrima za homobinuklearne komplekse. Na taj način moguće je ustanoviti linearnu zavisnost konstante kuplovanja i parametra (ε2-ε1)2 iz jedna čine 4.51. Vrednosti konstante kuplovanja izračunate na kristalnoj strukturi i geometrijama dobijenim sa različitim funkcionalima za binuklearni kompleks 1 prikazane su na slici 7.4. 0 1 2 3 4 5 6 7 8 9 1.E+05 2.E+05 3.E+05 4.E+05 5.E+05 6.E+05 (ε 1-ε 2)2 / cm-2 │J │/ cm-1 Slika 7.4 Grafička zavisnost konstante kuplovanja od parametra (ε2-ε1)2 izračunata za kompleks 1 na kristalnoj strukturi i geometrijama dobijenim optimizacijom sa LDA, OPBE, PW91, BP86 , BLYP 109,110,111 28,29 , i PWP101 funkcionalima Kod Kahn-ovog modela postoji linearna zavisnost između konstante kuplovanja i kvadrata integrala preklapanja magnetnih orbitala u BS stanju (Slika 7.5). 72 0 1 2 3 4 5 6 7 8 9 0.E+00 1.E-04 2.E-04 3.E-04 4.E-04 5.E-04 6.E-04 S 2 │J │/ cm-1 Slika 7.5 Grafička zavisnost konstante kuplovanja od parametra S2 , BLYP , i PWP funkcionalima izračunata za kompleks 1 na kristalnoj strukturi i geometrijama dobijenim optimizacijom sa LDA, OPBE, PW91, BP86 Ovakva zavisnost pokazuje da se S2 parametar može koristiti za detaljniju analizu magnetnih interakcija u binuklearnim kompleksima. 7.3 Zaključak Tendencija oksalatnih jona ka formiranju 2-D i 3-D mreža sa prelaznim metalima poslužila je u sintezi molekulskih magneta. U cilju dizajniranja molekula sa traženim osobinama, neophodno je prethodno ispitati uticaj različitih metala, odnosno njihove elektronske konfiguracije, na jačinu i tip kuplovanja. Korišćenjem DFT-BS pristupa za seriju binuklearnih kompleksa sa oksalatnim mostom, ustanovljeno je da povećanje broja nesparenih elektrona vodi ka pojačanju kako antiferomagnetnih, tako i feromagnetnih kuplovanja, pa totalne konstante kuplovanja imaju male vrednosti. Kuplovanja između različitih metala skoro uvek imaju slab feromagnetni karakter, zbog dodatnog feromagnetnog doprinosa usled interakcija popunjenih i nepopunjenih d orbitala metala. Nadalje, najjača kuplovanja registrovana su kod kompleksa kod kojih se delokalizacija elektrona ostvaruje usled σ interakcija d orbitala metala sa p orbitalama liganada. 73 8 Prenošenje magnetnih kuplovanja posredstvom ne- kovalentnih interakcija Do sada su naša istraživanja bila fokusirana na magnetnim kuplovanjima preko mostnog ligandnog sistema. Izračunavanje magnetnih kuplovanja koja se prenose posredstvom nekovalentnih interakcija predstavlja pravi izazov u naučnom smislu i pored razvoja teorije o magnetnim osobinama binuklearnih kompleksa112,113,114,115,116,117. U proteklih par godina postignut je veliki napredak u dizajniranju i sintezi koordinacionih polimera kod kojih su gradivne monomerne jedinice međusobno povezane vodoničnim vezama118,119,120. Pored činjenice da se vodonične veze mogu koristiti kao strukturni elementi, one takođe imaju ulogu kao medijatori magnetnih interakcija. Određena usmerenost vodoničnih veza kontroliše raspored susednih paramagnetnih centara, a samim tim orijentaciju njihovih magnetnih momenata. Postojanje stabilnih strukturnih jedinica, kao što je HF2- pokazalo se kao neophodno u dizajniranju materijala sa primenom u elektronici118,119. Pored navedenih aplikacija, prenos naelektrisanja posredstvom vodoničnih veza ima značajnu ulogu u biološkim procesima, kao što su oksidativna fosforilacija i fotosinteza. U ovom delu teksta detaljno će biti analizirani mehanizmi kuplovanja posredstvom nekovalentnih interakcija121, na kristalnim strukturama za seriju ispitivanih kompleksa: [Mn2(F)10HF2]5- (1), [Ni2(HF4)(pyz)8]+ (2a),(2b)(pyz=pyrazine), [Cu2(HF2)(pyz)8]3- (3), [Cu2F(HF)(HF2)(pyz)8]2- (4a) i (4b), [Cu2F4(H2O)2] (5), [Cu2(php)2(H2O)2]2+ (6), [Cu2(L1)2]2+ (7) (HL1=1-(2-hydroxyethylaminomethyl)-3- hydroxyethyl-1,3-diazacyclohexane), (php=6-(pyridine-2-ylhydrazonomethyl)phenol), aktivni centri Escherichia coli (multi-copper oxidase) (8) i peptidylglycine α-hydroxylating monooxygenase (9) i [Cu2Br8]4- (10), [Fe2Br8]2- (10’) (Slika 5.1)118,119,122,123,124,125,126,127,128,129. Monomerne jedinice u svim ispitivanim strukturama međusobno su povezane nekovalentnim interakcijama. 74 Slika 5.1 [Mn2(F)10HF2]5- (1), [Ni2(HF4)(pyz)8]+ (2a),(2b)(pyz=pyrazine), [Cu2(HF2)(pyz)8]3- (3), [Cu2F(HF)(HF2)(pyz)8]2- (4a) i (4b), [Cu2F4(H2O)2] (5), [Cu2(php)2(H2O)2]2+ (6), [Cu2(L1)2]2+ (7) (HL1=1-(2-hydroxyethylaminomethyl)-3-hydroxyethyl-1,3-diazacyclohexane), (php=6-(pyridine-2- ylhydrazonomethyl)phenol), aktivni centri Escherichia coli (multi-copper oxidase) (8) i peptidylglycine α-hydroxylating monooxygenase (9) i [Cu2Br8]4- (10), [Fe2Br8]2- 8.1 Računski metod (10’). Isprekidane linije predstavljaju vodonične veze Konstante kuplovanja za seriju navedenih kompleksa izračunate su na odgovarajućim kristalnim strukturama korišćenjem B3LYP hibridnog funkcionala. U cilju poređenja konstante kuplovanja su takođe izračunate i uz LDA funkcional. Za atome prelaznih metala upotrebljen je TZV (triple-zeta valence) bazis set, dok je za atome nemetala korišćen VDZ (valence double zeta) bazis set90,91 8.2 Rezultati i diskusija . Kod analiziranih kompleksnih jona magnetne interakcije se prenose posredstvom F-H….F, O-H….F, O-H...O, i Br....H-A-H....Br (A=O,N) vodoničnih veza, pa će shodno tome detaljno biti razjašnjen mehanizam kuplovanja za svaku pomenutu grupu. 75 8.2.1 Kuplovanja posredstvom F-H…. F vodoničnih veza Jednu od najjačih poznatih vodoničnih veza može da ostvari HF2- anjon, pa je određivanje mehanizma kuplovanja posredstvom ovakve veze od velikog značaja. U ovom delu detaljno su analizirane magnetne osobine kompleksa sa pomenutim mostnim ligandom t.j strukture 1, 2a, 2b, 3, 4a, 4b (Slika 5.1). Kristalne strukture za komplekse 2a i 2b razlikuju se po prirodi kontra-jona, što uzrokuje razlike u geometriji unutrašnje sfere, dok se kompleksi 4a i 4b razlikuju u strukturnim fragmentima. Svi izračunati podaci sumirani su i upoređeni sa dostupnim eksperimentalnim vrednostima (Tabela 8.1) Tabela 8.1 Konstante kuplovanja( cm-1 ) izračunate na kristalnim strukturama kompleksa 1-4 LDA i B3LYP funkcionalima Struktura J JLDA JB3LYP [Mn exp/calc 2(F)10HF2]5- 0,0 (1) 0,0 0,0 [Ni2(HF4)(pyz)8]+ -11,3 (2a) -2,9 -2,2/-9,2 [Ni2(HF4)(pyz)8]+ -17,7 (2b) -4,1 -4,1 /-14,3 [Cu2(HF2)(pyz)8]3- -1,3 (3) -0,0 -0,1 [Cu2F(HF)(HF2)(pyz)8]2- -0,5 (4a) -0,3 -0,3 [Cu2F(HF)(HF2)(pyz)8]2- 135,3 (4b) 0,8 / Generalno govoreći, rezultati dobijeni B3LYP funkcionalom su u boljoj saglasnosti sa dostupnim eksperimentalnim vrednostima. U slučaju kompleksa 1 i 4a LDA funkcional daje sasvim zadovoljavajuće rezultate. Kompleksi Struktura kompleksa 1 se najlakše može predstaviti kao 2-D mreža izgrađena od kvadratno planarnih [MnF 2a, 2b i 4a pokazuju slabe antiferomagnetne interakcije, dok su za kompleks 4b dobijene slabe feromagnetne interakcije. Odsustvo kuplovanja je primećeno kod jedinjenja 1 i 3. 4]- jedinica međusobno povezanih HF2- vodoničnim mostovima. Visoko spinski kompleksi mangana (III) u Oh simetriji podložu Jahn-Teller-ovoj distorziji usled degeneracije osnovnog stanja (5Eg). Simetrija se snižava do D4h, pri čemu elongovana 76 struktura predstavlja globalni minimum. Kao posledica Jahn-Teller-ovog efekta Mn-F dužine veza u monomernim jedinicama, jako su elongovane u pravcu formiranja vodoničnih mostova (Mn-Felong U cilju objašnjenja odsustva kuplovanja, bilo je potrebno izvršiti detaljnu analizu molekulskih orbitala ispitivanog dimera u visoko spinskom stanju (Slika 8.2). 2.190 Å, Mn-F 1.836 Å). Slika 8.2 Magnetne orbitale u visoko spinskom stanju za kompleks 1 Sa povećanjem broja nesparenih elektrona po metalnom centru, istovremeno se povećava intenzitet kako feromagnetnih, tako i antiferomagnetnih interakcija. Međutim, delokalizacija elektrona posredstvom π interakcija metala i liganada (usled bočnog preklapanja d orbitala metala i p orbitala liganada) nije moguća jer p orbitale mostnih liganada ne učestvuju u preklapanju (Slika 8.2(a)-8.2(c)), što je zasigurno posledica elongacije. Iz tog razloga, jedini mogući put delokalizacije predstavljaju σ interakcije metala i liganada (čeono preklapanje d orbitala metala i p orbitala liganada) (Slika 8.2(d)). Kako vrednost Mn-F….H ugla iznosi 112°, orbitale liganada su smeštene u približno paralelnu poziciju zbog čega dolazi do prekida povoljnih vezivnih σ interakcija između orbitala liganada, a samim tim i prekida antiferomagnetnih interakcija. Iako su vodonikovi atomi delovi mostnih jedinica, oni ne učestvuju u preklapanju i prenošenju magnetnih 77 interakcija. Grafička zavisnost konstante kuplovanja od Mn-F dužine veze prikazana je na slici 8.3. Slika 8.3 Grafička zavisnost konstante kuplovanja od dužine Mn-F veze Kompresovanje Mn-F veze u pravcu formiranja vodoničnih mostova, dovodi do značajnog povećanja antiferomagnetnih interakcija. Doprinosi p orbitala mostnih liganada u bočnom preklapanju postaju sve značajniji skraćivanjem Mn-F veze, usled čega je moguća delokalizacija posredstvom π interakcija metala i liganda. Na dužini veze od 1.30 Å spinska gustina je značajno delokalizovana sa metala na mostne ligande. Promena Mn-F....H ugla ne utiče na bočno preklapanje metala i mostnih liganada. Kako je frontalno preklapanje orbitala metala i mostnih liganada ključno za antiferomagnetne interakcije, povećanje Mn- F....H ugla može dovesti do značajnog pojačanja delokalizacije putem σ interakcija (Slika 8.4). 78 Slika 8.4 Grafička zavisnost konstante kuplovanja od vrednosti Mn-F-H ugla Povoljne vezivne σ interakcije između p orbitala fluora drastično su redukovane za male vrednosti Mn-F....H ugla. Interakcije sa mostnim vodonikovim atomima su nevezivne, što takođe utiče na jačinu magnetnih kuplovanja. Sa povećanjem Mn-F....H ugla doprinosi p orbitala fluora se smanjuju , dok se doprinosi s orbitala povećavaju. Blizu vrednosti ugla od 180° preklapanja orbitala su u potpunosti čeona sa značajnim doprinosom s orbitale sa vodonika. Odsustvo kuplovanja je očigledno posledica kako strukturnih (vrednost Mn-F.... Za razliku od prethodno navedene strukture kompleksa 1, [Ni H ugla), tako i elektronskih faktora (Jahn-Teller-ova distorzija). 2(HF4)(pyz)8]+ kompleksi 2a ( Ni-F….H ~160°) i 2b (Ni-F….H =180°) sa istom strukturom mosta, ali različitim kontra jonima pokazuju veoma slabe antiferomagnetne interakcije. U literaturi su do sada opisani DFT+U proračuni za 2a i 2b, koji nisu bili u saglasnosti sa eksperimentalnim i našim DFT rezultatima (Tabela 8.1)119. Magnetna kuplovanja putem vodoničnih mostova značajno su slabija u poređenju sa tipičnim binuklearnim kompleksima nikla (II) (J ≈ -30cm-1)130, kod kojih se antiferomagnetne interakcije ostvaruju usled σ preklapanja dz2 i dx2-y2 orbitala metala sa p orbitalama mostnih liganada. Pa ipak, određena orijentacija magnetnih orbitala može dovesti do redukcije ili čak odsustva kuplovanja131. U slučaju ispitivanih kompleksa samo jedna d-orbitala metala ima pravilnu orijentaciju, koja omogućava σ interakcije sa p 79 orbitalama mostnih liganada (Slika 8.5), što objašnjava slabe antiferomagnetne interakcije kod ispitivanih kompleksnih jedinjenja. Slika 8.5 Magnetne orbitale u visoko spinskom stanju za kompleks 2 Osnovnu strukturnu razliku u kompleksima 1, 2a i 2b predstavlja vrednost M-F….H ugla. Za komplekse 2a i 2b ovaj ugao je znatno veći (160° u 2a i 180° u 2b). Položaj p orbitala mostnih liganada je skoro u potpunosti frontalni (Slika 8.5(b)), što je važan preduslov za prenos magnetnih kuplovanja putem σ interakcija. U većini binuklearnih kompleksa bakra (II) prisutne su veoma jake magnetne interakcije (J ≈ -350cm -1), usled σ preklapanja dz2 ili dx2-y2 orbitala metala sa simetrijski adaptiranim p orbitalama mostnih liganada132. Struktura kompleksa 3 je jako distorgovana duž Cu-F veze, zbog izraženog Jahn-Teller-ovog efekta, karakterističnog za d9 konfiguraciju (Cu-F 2,280 Å, Cu-N 2,041 Å). Monte Carlo simulacije pokazale su veoma slabo antiferomagnetno kuplovanje122, kao i naši DFT proračuni sa LDA funkcionalom. Proračuni sa B3LYP nivoom teorije pokazuju primetno odsustvo kuplovanja. Kako Cu-F….H ugao iznosi 180°, odsustvo kuplovanja mora biti posledica elongacije. Usled elongacije duž Cu-F veze u kompleksu 3, d orbitale metala su raspoređene tako da nema povoljnih σ preklapanja sa p orbitalama mostnih liganada (Slika 8.6). 80 Slika 8.6 Magnetne orbitale u visoko spinskom stanju za kompleks 3 U slučaju struktura 4a i 4b most je produžen (F-H-F-H-F u poređenju sa F-H-F u kompleksima 1, 2 i 3). Za razliku od strukture 3 u kompleksima 4a i 4b, primećena je elongacija duž Cu-N veze (Cu-Nelong 2,37 Å, Cu-N 2,04 Å, Cu-F 2,01 Å). Odgovarajuće d orbitale metala smeštene su u ravni Cu-F veza, pa su s toga moguće magnetne interakcije usled σ preklapanja (Slika 8.7). Slika 8.7 Magnetne orbitale u visoko spinskom stanju za kompleks 4a Doprinosi p orbitala liganada su veoma mali što dovodi do slabih antiferomagnetnih interakcija u kompleksu 4a i slabih feromagnetnih interakcija u 4b. Sa porastom dužine mosta raste i rastojanje između metalnih centara (Cu-Cu 7,79 Å, Mn-Mn 5,64 Å Ni-Ni 6,30 Å Cu-Cu 6,85 8.2.2 Kuplovanja posredstvom O-H Å) što takođe dovodi do smanjenja magnetnih interakcija. ….F i O-H... O vodoničnih veza Konstante kuplovanja su izračunate na kristalnim strukturama ispitivanih kompleksa 5, 6, 7, 8 i 9 124-127, bez prethodne geometrijske optimizacije. Koordinate za aktivne centre metaloenzima 8 i 9 izvučene su iz odgovarajućih kristalnih struktura dostupnih u okviru proteinske banke podataka (protein data bank PDB) pod kodovima 1VK7 and 1PHM. 81 Izračunate vrednosti konstanti, kao i dostupne eksperimentalne vrednosti predstavljene su u tabeli 8.2. Antiferomagnetno kuplovanje dominira skoro u svim kompleksima, izuzev u strukturama 8 i 9, gde je prisutno feromagnetno kuplovanje. Deprotonovane forme kompleksa 5, 6 i 7 pokazuju znatno jače antiferomagnetne interakcije u poređenju sa protonovanim formama. Tabela 8.2 Konstante kuplovanja( cm-1 ) izračunate LDA i B3LYP funkcionalima za komplekse 5-9 Struktura J JLDA JB3LYP [Cu exp 2F4(H2O)2 -49,4 ] (5) -9,4 -3,9 (5) deprotonovan -127,4 -39,6 / [Cu2(php)2(H2O)2]2+ -71,9 -12,6 -19,8 (6) deprotonovan -304,1 -44,4 / [Cu2(L1)2]2+ -28,5 (7) -14,3 -9,4 (7) deprotonovan -204,3 -76,8 / (8) 517,5 108,0 / (9) 102,2 52,8 / Za razliku od monomernih jedinica povezanih [HF2]- mostovima, kod struktura sa O-H….F i O-H...O vodoničnim mostovima primećene su znatno jače magnetne interakcije između paramagnetnih metalnih centara. Kako F-H….F predstavlja najjaču vodoničnu vezu, razlog ovakvih pojava mora se potražiti u prostornom rasporedu gradivnih jedinica, a ne u jačini vodonične veze. Prostorni raspored monomernih jedinica određuje orijentaciju magnetnih orbitala, pa je bila neophodna detaljna analiza molekulskih orbitala (Slika 8.8) kao i analiza spinskih gustina (Tabela 8.3). 82 Slika 8.8 Magnetne orbitale u visoko spinskom stanju za kompleks 5 (a) i njegovu deprotonovanu strukturu (b) Na slici 8.8 prikazane su magnetne orbitale u tripletnom stanju za strukturu 5, kao i za odgovarajuću deprotonovanu formu. Magnetne orbitale predstavljaju vezivnu i antivezivnu kombinaciju orbitala koje su lokalizovane na monomerima. Odgovarajuće d orbitale metala smeštene su u ravni sa orbitalama liganada, što dovodi do jakih antiferomagnetnih interakcija. Značajni doprinosi p orbitala liganada, ukazuju na intenzivnu delokalizaciju. Doprinosi s orbitala vodonika su zanemarljivi, zbog čega se izvodi zaključak da je uloga vodonične veze da postigne razumnu udaljenost između gradivnih jedinica, koje su neophodne za magnetne interakcije. Drugim rečima, mostni atomi se moraju naći na odgovarajućim rastojanjima, na kojima njihove nekovalentne interakcije dovode do primetnih magnetnih kuplovanja između metalnih centara monomernih jedinica. Uklanjanje protona sa mostnih kiseonikovih atoma, dovodi do veće delokalizacije, a samim tim do jačih antiferomagnetnih interakcija (Table 8.2). Doprinosi p orbitala liganada znatno su veći u poređenju sa protonovanom strukturom (Slika 8.8). Sve ove činjenice ukazuju da vodonikovi atomi ne učestvuju u magnetnim kuplovanjima. Spinske gustine u strukturama 5, 6 i 7 najvećim delom su lokalizovane na metalima, sa izvesnom delokalizacijom na atome liganada (Tabela 8.3). Pojačana delokalizacija spinskih gustina prisutna je u deprotonovanim strukturama uz istovremeno smanjenje spinske gustine na atomima bakra. Aktivni centri metaloenzima bakra 8 i 9 pokazuju veoma jake feromagnetne interakcije (Tabela 8.2). Međusobne orijentacije atoma bakra u analiziranim proteinima direktno 83 zavise od tipa i položaja amino kiselina u blizini aktivnog centra. Svaki metalni centar poseduje geometriju distorgovanog tetraedra, što može da prouzrokuje prostornu ortogonalnost između magnetnih orbitala. Neophodno je napomenuti da se magnetne interakcije u strukturi 8 ne prenose sistemom vodoničnih veza, već preko jednog molekula vode ili OH- jona u zavisnosti od pH sredine. Kako su dužine Cu-O veza 2.43 Å i 2.30 Å redom, ovakva kuplovanja se ipak mogu smatrati magnetnim interakcijama na daljinu. Magnetne orbitale za strukturu 8 lokalizovane su na pojedinačnim metalnim centrima, sa značajnom delokalizacijom na kiseonikov atom, i međusobno su ortogonalne (Slika 8.9). Slika 8.9 Magnetne orbitale u visoko spinskom stanju za kompleks 8 Jako feromagnetno kuplovanje u strukturi 9 prenosi se preko lanca izgrađenog od molekula vode, koji su umetnuti između dva atoma bakra i koji održavaju veliku Cu-Cu udaljenost (9.87 Å). Magnetne orbitale su međusobno ortogonalne, pri čemu je jedna orbitala lokalizovana na kiseonikovom atomu, dok je druga lokalizovana na atomu bakra sa izvesnom delokalizacijom na atome liganada (Slika 8.10). Slika 8.10 Magnetne orbitale u visoko spinskom stanju za kompleks 9 84 Na osnovu toga se može zaključiti da su magnetne interakcije prisutne između atoma bakra sa jedne strane i kiseonikovog atoma sa drugog metalnog centra. Spinske gustine izračunate LDA funkcionalom najvećim delom su lokalizovane na jednom atomu bakra i kiseonikovom atomu vezanom za drugi atom bakra. Mali iznos spinske gustine prisutan je na drugom metalnom centru usled delokalizacije sa kiseonikovog atoma. Međutim, analiza spinskig gustina na B3LYP nivou teorije pokazuje veliki iznos negativne spinske gustine na drugom atomu bakra usled spinske polarizacije. Tabela 8.3 Spinske gustine ispitivanih kompleksa 5-9 Struktura LDA B3LYP ρO/ρF ρCu ρO /ρF ρCu [Cu2F4(H2O)2 0,079/0,114 ] (5) 0,549 0,061/0,083 0,656 (5) deprotonovan 0,350/0,150 0,406 0,320/0,120 0,484 [Cu2(php)2(H2O)2]2+ 0,047 (6) 0,478 0,042 0,600 (6) deprotonovan 0,300 0,400 0,280 0,520 [Cu2(L1)2]2+ 0,062 (7) 0,433 0,053 0,537 (7) deprotonovan 0,272 0,338 0,246 0,412 (8) 0,443 0,400 0,248 0,600 (9) 0,307 0,456 0,695 0,682 8.2.3 Kuplovanja posredstvom Br....H-A-H....Br (A=O,N) i Br.... Magnetne interakcije koje se prenose preko Br Br veza ....H-A-H....Br i Br....Br veza, prisutne su u molekulskim 2-D antiferomagnetima izgrađenim od [CuBr4]2- i [FeBr4]- kompleksnih jona128,129 Magnetna kuplovanja zavise od međumolekulskih geometrijskih parametara (geometrijski parametri između dve [CuBr . 4]2- jedinice), odnosno Br....Br rastojanja i Cu- 85 Br....Br-Cu torzionih uglova. Izračunate vrednosti konstanti za različite položaje (a-g) [CuBr4]2- gradivnih jedinica za strukturu 10, kao i relevantni međumolekulski geometrijski parametri predstavljeni su u tabeli 8.4. Tabela 8.4 Relevantni geometrijski parametri; intermolekulska rastojanja(Å) i torzioni uglovi (°) i konstante kuplovanja( cm-1 Orijentacija ) izračunate LDA i B3LYP funkcionalima za kompleks 10 Br.... CuBr .... Cu-BrCu .... JBr-Cu JLDA JB3LYP 10(a) exp 4,22 7,87 16,00 -15,2 -9,2 -4,4 10(b) 6,49 9,75 -139,5 -0,1 0,0 0,0 10(c) 2,33/3,67 4,51 115,10/44,90 -482,3 -260,7 / 10(d) 5,36 5,39 -7,30 -0,6 0,2 / 10(e) 2,16/2,58 5,33 -32,5/19,1 -729,6 -284,3 / 10 (f) 3,89 8,09 -55,90 -49,8 -27,0 -11,6 10 (g) 4,40/4,52 6,99 -62,40/63,40 -3,3 -1,4 -5,5 Kao i u ostalim slučajevima B3LYP funkcional pokazuje bolja slaganja sa eksperimentalnim vrednostima. Antiferomagnetne interakcije prisutne su za sve orijentacije, sa izuzetkom orijentacija 10 (b) i 10 (d) (Tabela 8.4). Osnovno stanje tetraedarskog [CuBr4]2- jona je 2T2, pa su ovakve strukture podložne Jahn- Teller-ovoj distorziji (T2⊗(e+τ2) JT problem). Distorzija duž e vibracije snižava simetriju iz Td u D2d, pa se 2T2 osnovni term razlaže na 2E i 2B2 stanje, gde 2B2 predstavlja globalni minimum. Distorzija [CuBr4]2- monomernih jedinica veoma je važna za prostorni raspored gradivnih blokova i formiranje 2-D slojeva. Nadalje, pozicija bromovih atoma u kristalu povoljna je za magnetna kuplovanja. Da bi se objasnili uticaji strukturnih efekata na magnetne interakcije, potrebno je analizirati molekulske orbitale dimera u stanju visokog spina za različite orijentacije [CuBr4]2- monomernih jedinica (10a-10g) (Slika 8.11). 86 Slika 8.11 Magnetne orbitale u tripletnom stanju za različite orijentacije strukture 10 Položaji 10a-10g razlikuju se u Br....Br nekovalentnom rastojanju i Cu-Br....Br-Cu torzionim uglovima. Molekulske orbitale predstavljaju simetrične i antisimetrične kombinacije orbitala monomera. Doprinosi p orbitala broma su znatno veći u odnosu na doprinose d orbitale sa bakra, pa jačina kuplovanja direktno zavisi od Br....Br kontakta. Orbitale sa mostnih molekula vode ne doprinose (položaj 10(a)), pa se kuplovanje odvija preko Br....Br interakcija. Najjače antiferomagnetno kuplovanje prisutno je kod orijentacije 10(e), gde su p orbitale broma u čeonom položaju i na dovoljnoj udaljenosti. Povoljna čeona preklapanja drastično se smanjuju promenom Cu-Br....Br-Cu torzionog ugla. Odsustvo kuplovanja ili feromagnetna kuplovanja prisutna su u slučaju kada nema značajnog kontakta između susednih bromovih atoma ili kada su p orbitale u paralelnom položaju (10(d)). Sa porastom Br....Br i Cu....Cu rastojanja jačina antiferomagnetnih interakcija opada. Na udaljenosti od 6.49 Å primetno je odsustvo bilo kakvih interakcija (10(b)). Grafička zavisnost konstante kuplovanja od Br.... Br rastojanja prikazana je na slici 8.12. 87 Slika 8.12 Grafička zavisnost konstante kuplovanja od Br.... Br rastojanja za različite orijentacije monomernih jedinica Absolutna vrednost konstante kuplovanja opada približno eksponencijalno sa porastom rastojanja, što je posledica promene u preklapanju orbitala monomera. Kako se strukture 10(a)-10(g) razlikuju kako u Br....Br rastojanju tako i u Cu-Br....Br-Cu torzionom uglu, u cilju ispitivanja uticaja Br....Br rastojanja, izabrana je 10(f) orijentacija i izračunate su konstante za tačno određene Br....Br udaljenosti (Slika 8.13). Slika 8.13 Grafička zavisnost konstante kuplovanja od Br.... Br rastojanja za 10 (f) orijentaciju monomernih jedinica U visoko spinskom tripletnom stanju jedan elektron se smešta u vezivnu orbitalu, dok se drugi elektron nalazi u antivezivnoj orbitali. Nepovoljne antivezivne interakcije postaju sve jače kako se rastojanje između monomernih jedinica smanjuje. Na udaljenosti od 2.00 Å energija tripletnog stanja naglo raste, dok je porast energije BS stanja manje intenzivan. Samim tim dolazi do porasta antiferomagnetnih interakcija. Udaljavanjem monomernih 88 jedinica antivezivne interakcije u tripletnom stanju slabe, do tačke u kojoj su energije tripletnog i BS stanja izjednačene (Slika 8.14). Slika 8.14 Grafička zavisnost energije tripletnog i BS stanja od udaljenosti između monomernih jedinica Antiferomagnetna kuplovanja između [FeBr4]- monomernih jedinica značajno su slabija u poređenju sa odgovarajućim [CuBr4]2- strukturama (Tabela 8.6). Tabela 8.6 Relevantni geometrijski parametri; intermolekulska rastojanja(Å) i torzioni uglovi (°) i konstante kuplovanja( cm-1 Orijentacija ) izračunate LDA i B3LYP funkcionalima za kompleks 10’ Br.... FeBr .... Fe-BrFe .... JBr-Fe JLDA JB3LYP 10’(a) exp 3,81 6,97 111,41 -2,30 -0,81 / 10’(b) 4,02 6,11 62,86 -1,70 -0,57 / 10’(c) 3,83 7,76 -77,17 -2,94 -1,45 / Gvožđe (III) jon u [FeBr4]- poseduje visoko spinsku konfiguraciju sa pet nesparenih elektrona. Usled simetrične raspodele elektronske gustine, nije moguća pojava Jahn-Teller- ove distotzije. Antiferomagnetni doprinosi potiču od dxz→dxz, dyz→dyz i dxy→dxy interakcija (σ- preklapanja sa orbitalama liganada), kao i dx2-y2→dx2-y2, dz2→dz2 interakcija (π-preklapanja sa orbitalama liganada). Feromagnetni doprinosi prisutni su usled dxz→dyz, 89 dxz→dxy, dxy→dyz, dx2-y2→dz2, dx2-y2→d(xz,yz,xy), dz2→d(xz,yz,xy) 8.3 Zaključak interakcija. Kao posledica dva antagonistička efekta, totalno kuplovanje je veoma slabo. Grafička zavisnost konstante kuplovanja od geometrijskih parametara ne može biti ustanovljena zbog prisustva π- interakcija. Na osnovu analiziranih podataka može se zaključiti, da u većini ispitivanih slučajeva vodonične veze služe kao strukturne jedinice, čija usmerenost u prostoru određuje jačinu kuplovanja između metalnih centara monomernih jedinica. Prostorni raspored vodoničnih veza određuje prostornu orijentaciju magnetnih momenata susednih magnetnih centara. Odsustvo kuplovanja uglavnom je registrovano kod struktura koje su podložne Jahn-Teller-ovoj distorziji. Kao posledica elongacije d orbitale metala se nalaze u položaju u kojem nema povoljnih preklapanja sa orbitalama liganada. Magnetne interakcije koje se prenose preko O-H….F i O-H...O vodoničnih veza značajno su jače u odnosu na kuplovanja preko F-H….F vodoničnih veza. Ovakvi podaci nagoveštavaju da intenzitet kuplovanja ne zavisi uvek od jačine vodoničnih veza. Nadalje, jačina antiferomagnetnih kuplovanja se drastično smanjuje sa povećanjem udaljenosti između monomernih jedinica. Ova zavisnost se posebno uočava u molekulskim 2-D antiferomagnetima izgrađenim od [CuBr4]2-. 90 9 Zaključak U okviru ove doktorske disertacije prezentovana je detaljna analiza magnetnih osobina mononuklearnih i binuklearnih kompleksa prve serije prelaznih metala primenom DFT i LF-DFT metoda. U prvom delu analizirana su magnetna svojstva u trigonalno bipiramidalnim mononuklearnim kompleksima kao i u prethodno sintetisanom [NiCl3(Hdabco)2]+ Magnetna anizotropija mononuklearnih sistema analizirana je LF-DFT metodom, s obzirom da standardne DFT tehnike ne daju dobre rezultate. Izrazito velike vrednosti energija prelaza dobijene su za trigonalno bipiramidalne komplekse nikla (II) usled degeneracije d kompleksu, dok su u drugom delu ispitivane magnetne osobine binuklearnih kompleksa, kao i magnetne interakcije između mononuklearnih gradivnih jedinica međusobno povezanih nekovalentnim interakcijama. Osnovni cilj ovog istraživanja bilo je ispitivanje različitih elektronskih i strukturnih uticaja, na magnetne interakcije u kompleksima prve serije prelaznih metala. Analiza magnetnih svojstava ispitivanih jedinjenja može poslužiti za predviđanje i dizajniranje struktura sa osobinama neophodnim za tehnološku aplikaciju. x 2 -y 2 i dxy orbitala, kao i izraženog spin-orbitalnog kuplovanja karakterističnog za nikal(II) (ζNi(II)=649 cm -1). Detaljnom analizom na model-sistemima utvrđeno je da Jahn-Teller-ova distorzija “gasi” efekte spin-orbitalnog kuplovanja usled čega se smanjuje magnetna anizotropija sistema. Da bi se smanjili efekti Jahn-Teller-ove distorzije, neophodno je bilo izabrati sistem sa nefleksibilnim helatnim ligandima u ekvatorijalnoj ravni ili voluminoznim ligandima u aksijalnoj ravni. Iz pomenutih razloga izabran je realni sistem [NiCl3(Hdabco)2]+, sa voluminoznim aksijalnim ligandima, i na njemu su izračunate energije prelaza za visoku D3 i snižene C2 simetrije, LFDFT i MRCI metodama. Nadalje, ustanovljeno je da spin-orbitalno kuplovanje umanjuje efekte distorzije transformisanjem Jahn-Teller-ovog efekta u pseudo Jahn-Teller-ov efekat, usled čega nema odstupanja od geometrije visoke simetrije, pa se prelazi mogu računati u strukturi visoke simetrije gde oni imaju veoma visoke vrednosti. Problemi visokih energetskih barijera mogu se prevazići ili primenom električnog polja ili kombinovanim delovanjem svetlosti i magnetnog polja. 91 Povezivanje mononuklearnih gradivnih jedinica nekovalentnim interakcijama ili preko kovalentnih mostova dovodi do formiranja 2-D i 3-D mreža. Usled ovakvih povezivanja javljaju se magnetne interakcije između dva ili više magnetnih centara, pa je njihova analiza i predviđanje od velikog značaja za dizajniranje i kasniju sintezu magneta na molekulskom nivou. Na magnetna kuplovanja utiču tip metala i njegova elektronska konfiguracija, struktura mosta, dužine veza i uglovi između veza, udaljenost između paramagnetnih centara, kao i prostorni raspored orbitala, pogodan za kuplovanje. Delokalizacija elektrona sa metala na mostne ligande najizraženija je usled σ preklapanja d orbitala sa orbitalama mostnih liganada. Povećanje broja nesparenih elektrona vodi ka porastu antiferomagnetnih interakcija, pri čemu se istovremeno pojačavaju i feromagnetne interakcije kuplovanjem između ortogonalnih orbitala. Usled poništavanja antagonističkih efekata, totalne konstante kuplovanja su malih vrednosti. Dakle, najjače interakcije se mogu očekivati kod kompleksa nikla (II) i bakra (II). Odsustvo kuplovanja prisutno je usled Jahn- Teller-ove distorzije ili strukturne distorzije mosta. Raspored i pakovanje monomernih jedinica u kristalu mogu značajno uticati na intenzitet kuplovanja, dirigujući krajnji magnetni efekat sistema, kao i sam dizajn. Dakle, elektronski i strukturni faktori predstavljaju glavne činioce pri dizajniranju magneta na molekulskom nivou. Vodeći računa o strukturi tj. geometriji, kao i pri pažljivom izboru prelaznog metala, moguće je konstruisati sisteme koji zadovoljavaju potrebe savremenih tehnologija. 92 10 Korišćeni programi Za optimizacije mononuklearnih model sistema i [NiCl3(Hdabco)2]+ kompleksa korišćen je ADF (Amsterdam Density Functional) softverski paket133. Za izračunavanje Racah parametara korišćen je modifikovani DFT metod, označen kao LF-DFT46. Geometrije binuklearnih kompleksa mangana sa fenilcijanamidnim ligandima optimizovane su CFF programom za molekulsku mehaniku134. Magnetna kuplovanja u binuklearnim kompleksima prelaznih metala izračunata su BS metodom, ORCA programskim paketom135. Molekulske orbitale analizirane su gOpenMol programom136 . Sva izračunavanja urađena su na klasteru Univerziteta u Friburu, koji je obezbedio licence za odgovarajuće programe. 93 11 Literatura: 1 P.W. Anderson, Phys.Rev. 1959, 115, 2. 2 P. J. Hay, J. C. Thibeault, R. Hoffmann, J.Am.Chem.Soc. 1975, 97, 4884. 3 G. Rajaraman, F. Totti, A. Bencini, A. Caneschi, R. Sessoli, D. Gatteschi, Dalton. Tr ans. 2009, 45, 3153. 4 A. Bencini, M. Casarin, D. Forrer, L. Franco, F. Garau, N. Masciocchi, L. Pandolfo, C. Pettinari, M. Ruzzi, A.Vittadini, Inorg. Chem. 2009, 48, 4044. 5 A. Bencini, F. Totti, J. Chem. Theory. Comput. 2009, 5, 44. 6 A. Bencini, F. Totti, Int. J. Quantum Chem. 2005, 101, 819. 7 C. Duboc, D. Ganyushin, K. Sivalingam, M. N. Collomb, F. Neese, J. Phys. Chem. 2010, 114, 10750. 8 D. A. Pantazis, V. Kremwald, M. Orio, F. Neese, Dalton Trans. 2010, 39, 4959. 9 M. Orio, D. A. Pantazis, T. Petrenko, F. Neese, Inorg. Chem. 2009, 48, 7251. 10 E. Ruiz, A. Rodríguez-Fortea, J. Cano, S. Alvarez, P. Alemany, J. Comput. Chem. 2003, 24, 982. 11 O. Kahn (1993) Molecular Magnetism, VCH Publishers, New York. 12 R. M. Dreizler, E. K. U. Gross (1990) Density Functional Theory, An Approach to Quantum Many-Body Problem, Springer-Verlag, Berlin, Heidelberg, New York, London, Paris, Tokyo, Hong Kong, Barcelona. 13 W. Koch, M. C. Holthausen (2001) A Chemist Guide to Density Functional Theory, Willey-VCH. 14 T. Ziegler, Chem. Rev. 1991, 91, 651. 15 E. J. Baerends, O. V. Gritsenko, J. Phys. Chem. 1997, 101, 651. 16 R. Stowasser, R. Hoffman, J. Am. Chem. Soc. 1999, 121, 3414. 17 G. te Velde, F. M. Bickelhaupt, S. J. A. van Gisbergen, C. Fonseca Guerra, E. J. Baerends, J. G. Snijders, T. Ziegler, J. Comput. Chem. 2001, 22, 931. 18 P. Geerlings, F. de Profit, W. Langenaeker, Chem. Rev. 2003, 103, 1793. 94 19 R. G. Parr, W. Yang (1994) Density-Functional Theory of Atoms and Molecules, University Press, Oxford. 20 P. Hohenberg, W. Kohn, Phys. Rev. 1964, 136, 864. 21 M. Levy, Phys. Rev. A, 1982, 26, 1200. 22 L. H. Thomas, Proc. Camb. Phil. Soc. 1926, 23, 542. 23 W. Kohn, L. J. Sham, Phys. Rev. 1965, 140, 1133. 24 C. J. Cramer (2002) Essentials of Computational Chemistry, Theories and Models, John Wiley and Sons LTD 25 S. Vosko, L. Wilk, M. Nusair, Can. J. Phys. 1980, 58, 1200. 26 J. C. Slater, Phys. ReV. 1951, 81, 385. 27 P. A. M. Dirac, Proc. Cambridge Phil. Roy. Soc. 1930, 26, 376. 28 A. D. Becke, Phys. Rev. 1988, 38, 3098. 29 C. Lee, W. Yang, R. G. Parr, Phys. Rev. 1988, 37, 785. 30 R. Sessoli, D. Gatteschi, A. Caneschi, M. A. Novak, Nature 1993, 365, 141. 31 A. Caneschi, D. Gatteschi, R. Sessoli, A. L. Barra, L. C. Brunel, M. Guillot, J. Am. Chem. Soc. 1991, 113, 5873. 32 R. L. White, J. Magn. Magn. Mater. 2000, 209, 1. 33 P. Gambardella, S. Rusponi, M. Veronese, S. S. Dhesi, C. Grazioli, A. Dallmeyer, I. Cabria, R. Zeller, P. H. Dederichs, K. Kern, C. Carbone, H. Brune, Science 2003, 300, 1130. 34 C. F. Hirjibehedin, C.-Y. Lin, A. F. Otte, M. Ternes, C. P. Lutz, B. A. Jones, A. J. Heinrich, Science 2007, 317, 1199. 35 M. N. Leuenberger, D. Loss, Nature 2001, 410, 789. 36 S. Hill, R. S. Edwards, N. Aliaga-Alcalde, G. Christov, Science 2003, 302, 1015. 37 F. Troiani, A. Ghirri, M. Affronte, S. Carretta, P. Santini, G. Amoretti, S. Piligkos, G. Timco, R. E. P. Winpennyet, Phys. Rev. Lett. 2005, 94, 207208. 38 M. H. L. Pryce, Phys. Rev. 1950, 80, 1107. 39 A. Abragam, M. H. L. Pryce, Proc. Roy. Soc. 1951, 205, 135. 40 I. B. Bersuker (2006) The Jahn-Teller Effect, Cambridge University Press, Cambridge. 95 41 M. Zlatar, C. W. Schlaepfer, E. P. Fowe, C. Daul, Pure Appl. Chem. 2009, 81, 1397. 42 M. Zlatar, M. Gruden-Pavlović, C. W. Schlaepfer, C. Daul, J. Mol. Struct. Theochem, 2010, 954, 86. 43 M. Gruden-Pavlović, M. Zlatar, C. W. Schlaepfer, C. Daul, J. Mol. Struct. Theochem, 2010, 954, 80. 44 M. Gruden-Pavlović, P. Garcia-Fernandez, Lj. Andjelković, C. Daul, M. Zlatar, J. Phys. Chem. A, 2011, 115, 10801. 45 H. L. Schläfer, G. Gliemann (1969) Basic Principles of Ligand Field Theory, Wiley Interscience, New York. 46 M. Atanasov, C. Rauzy, P. Baetting, C. Daul, Quant. Chem. 2005, 102, 119. 47 M. Atanasov, C. Daul, C. Rauzy, Struct. Bond. 2004, 106, 97. 48 M. Atanasov. C. A. Daul, Chem. Phys. Lett. 2003, 367, 209. 49 F. Senn, L. Helm, A. Borel, C. Daul, Chimie, 2012, 15, 250. 50 M. Atanasov, C. Daul, M. M. Rohmer, T. Venkatachalam, Chem. Phys. Lett. 2006, 427, 449. 51 L. Petit, A. Borel, C. Daul, P. Maldivi, C. Adamo, Inorg. Chem. 2007, 45, 7382. 52 M. Atanasov, P. Comba, C. Daul, Inorg. Chem. 2008, 47, 2449. 53 A. Trueba, P. Garcia-Fernandez, F. Senn, C. Daul, J. A. Aramburu, M. T. Barriuso, M. Moreno, Phys. Rew. B, 2010, 81, 075107. 54 P. Garcia-Fernandez, F. Senn, C. Daul, J. A. Aramburu, M. T. Barriuso, M. Moreno, Phys. Chem. Chem. Phys. 2009, 11, 7545. 55 W. Heisenberg, Z. Phys.1928, 49, 619. 56 P.A.M. Dirac (1947), The principles of quantum mechanics, Clarendon, Oxford. 57 A. Palli, B. Tsukerblat, S. Klokishner, K. R. Dunbar, J. M. Clemente-Juan, E. Coronado, Chem. Soc. Rev. 2011, 40, 3130. 58 L. Noodleman, C.Y. Peng, D.A. Case, J.M. Mouesca, Coord. Chem. Rev. 1995, 144 199. 59 L. Noodleman, D. A. Case, A. Aizman, J. Am. Chem. Soc. 1988, 110 1001. 96 60 L. Noodleman, J. G. Norman, J. H. Osborne, A. Aizman, D. A. Case, J. Am. Chem. Soc. 1985, 107, 3418. 61 T. Soda, Y. Kitagawa, T. Onishi, Y. Takano, Y. Shigeta, H. Nagao, Y. Yoshika, K. Yamaguchi, Chem. Phys. Lett. 2000, 319, 223. 62 C. J. Calzado, J. Cabrero, J. P. Malrieu, R. Caballol, J. Chem. Phys. 2002, 116, 2728 63 D. C. A. Sessoli, R. Gatteschi, M. A. Novak, Nature, 1993, 365, 64 P. C. E. Stamp, A. Gaita-Ariño, J. Mat. Chem. 2009, 19, 65 M. M. A. Pritchard, R. G. Ali, A. Uddin, Acta Cryst. C, 2006, 62, 66 W. J. Rozell, J. S. Wood, Inorg. Chem. 1977, 16, 67 N. C. Handy, A. J. Cohen, Mol. Phys. 2001, 99, 403. 68 J. P. Perdew, K. Burke, M. Ernzerhof, Phys. Rev. Lett. 1996, 77, 3865 69 J. S. Griffith (1964) The Theory of Transition-Metal Ions, Cambridge University Press, Cambridge. 70 F. H. Allen, The Cambridge Structural Database: a quarter of a million crystal structures and rising, Acta Cryst. B, 2002, 58, 380. 71 M. Gruden-Pavlović, M. Perić, M. Zlatar, P. García-Fernández, rad u pripremi. 72 P. Garcia-Fernandez, I. B. Bersuker, Phys. Rev. Lett. 2011, 108, 73 B. W. S. H. L. W. Z. M. Jiang, S. D. Wang, S. Gao, J. Am. Chem. Soc. 2011, 133, 74 R. J. Crutchley, Coord. Chem. Rew. 2001, 125, 219. 75 A. Escuer, N. Sanz, R. Vicente, F. A. Mautner, Dalton Trans. 2003, 11, 2121. 76 A. Escuer, N. Sanz, R. Vicente, F. A. Mautner Inorg. Chem. 2003, 42, 541. 77 A. Escuer, F. A. Mautner, N. Sanz, R. Vicente, Polihedron 2004, 23, 1409. 78 M. Perić, S. R. Niketić, M. Zlatar, M. Gruden-Pavlović, S. Grubišić, Monatsh. Chem. 2011, 142, 585. 79 M. Gruden-Pavlović, S. Grubišić, A. G. Coutsolelos, S. R. Niketić, J. Mol. Struct. 2001, 595, 209. 80 S. Grubišić, M. Gruden-Pavlović, S. R. Niketić, N. Sakagami-Yoshida, S. Kaizaki, J. Mol. Struct. 2002, 609, 1. 97 81 S. Grubišić, M. Gruden-Pavlović, S. R. Niketić, N. Sakagami-Yoshida, S. Kaizaki, Transition Met. Chem. 2003, 28, 37. 82 S. Grubišić, M. Gruden-Pavlović, S. R. Niketić, S. Kaizaki, N. Sakagami-Yoshida, Inorg. Chem. Commun. 2003, 9, 1180. 83 M. Gruden-Pavlović, S. Grubišić, S. R. Niketić, J. Inorg. Biochem. 2004, 98, 1293. 84 S. Grubišić, S. R. Niketić, D. D. Radanović, U. Rychlewska, B. Warzajtis, Polyhedron, 2005, 24, 1701. 85 S. Grubišić, M. K. Miličić, D. D. Radanović, S. R. Niketić, J. Mol. Struct. 2006, 794, 125. 86 S. Grubišić, D. D. Radanović, U. Rychlewska, B. Warzajtis, N. S. Drašković, M. Đuran, S. R. Niketić, Polyhedron, 2007, 26, 3437. 87 S. Grubišić, M. Gruden-Pavlović, S. R. Niketić, N. Sakagami-Yoshida, S. Kaizaki, J. Coord. Chem. 2007, 60, 851. 88 M. Gruden-Pavlović, S. Grubišić, M. Zlatar, S. R. Niketić, Int. J. Mol. Sci. 2007, 8, 810. 89 S. Grubišić, M. Gruden-Pavlović, D. D. Radanović, M. Perić, S. R. Niketić, J. Mol. Struct. 2009, 919, 54. 90 A. Schaefer, H. Horn, R. Ahlrichs, J. Chem. Phys. 1992, 97, 2571. 91 A. Schaefer, C. Huber, R. Ahlrichs, J. Chem. Phys. 1994, 100, 5829. 92 O. Castillo, A. Luque, P. Roma´n, F. Lioret, M. Julve, Inorg. Chem. 2001, 40, 5526. 93 X. Zhao, D. Liang, S. Liu, C. Sun, R. Cao, C. Gao, Y. Ren, Z. Su, Inorg. Chem. 2008 47, 7133. 94 H. Ōkawa, A. Shigematsu, M. Sadakiyo, T. Miyagawa, K. Yoneda, M. Ohba, H. Kitagawa, J. Am. Chem. Soc. 2009, 131, 13516. 95 I. Muga, J. M. Gutiérrez-Zorrilla, P. Vitoria, A. Luque, M. Insausti, P. Román, F. Lloret, Eur. J. Inorg. Chem. 2000, 30, 2541. 96 S. Decurtins, H. W. Schmalle, H. R. Oswald, A. Linden, J. Ensling, P. Gutlich, A. Hauser, Inorg. Chim. Acta 1994, 65, 216. 97 S. Decurtins, H. W. Schmalle, P. Schneuwly, H. R. Oswald, Inorg. Chem. 1993, 32, 1888. 98 98 M. Perić, M. Zlatar, M. Gruden-Pavlović, S. Grubišić, Monatsh. Chem. 2012, 143, 569. 99 W. Plass, Inorg. Chem. 1997, 36, 10. 100 W. Plass, Inorg. Chem. 1998, 37, 3167. 101 J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, M. R. Pederson, D. J. Sing, C. Fiolhais, Phys. Rev. 1992, 46, 6671. 102 S. Triki, F. Bérézovsky, J. S. Pala, M. T. Garland, Inorg. Chim. Acta 2000, 308, 31. 103 S. Triki, F. Bérézovsky, J. S. Pala, C. J. Gomez-Garcla, E. Coronado, K. Costuas, J. F. Calet, Inorg. Chem. 2001, 40, 5127. 104 X. D. Zhang, Z. Zhao, J. Y. Sun, Y. C. Ma, M. L. Zhu, Acta. Cryst. 2005, 61, 2643. 105 S. Triki, F. Bérézovsky, J. S. Pala, E. Coronado C. J. Gomez-Garcla, J. M. Clemente, A. Riou, P. Molinié, Inorg. Chem. 2000, 39, 3771. 106 D. Armentano, G. De Munno, J. Faus, F. Lloret, M. Julve, Inorg. Chem. 2001, 40, 655. 107C. J. Cramer, D. G. Truhlar, Phys. Chem. Chem. Phys. 2009, 11, 10757. 108 P. J. Hay, J. C. Thibeault, R. Hoffmann, J. Am. Chem. Soc. 1975, 97, 4884. 109 A. D. Becke, J. Chem. Phys. 1986, 84, 4524. 110 J. P. Perdew, Phys. Rev. 1986, 33, 8822. 111 J. P. Perdew, W. Yue, Phys. Rev. B 1986, 33, 8800. 112 O. Oms, J. B. Rota, L. Norel, C. J. Calzado, H. Rousselière, C. Train, V. Robert, Eur. J. Inorg. Chem. 2010, 2010, 5373. 113 J. B. Rota, B. Le Guennic, V. Robert, Inorg. Chem. 2010, 49, 1230. 114 L. Norel, J. B. Rota, L. M. Chamoreau, G. Pilet, V. Robert, C. Train, Angew. Chem. Int. Ed. 2011, 50, 7128. 115 J. B. Rota, C. J. Calzado, C. Train, V. Robert, J. Chem. Phys. 2010, 132, 154702. 116 J. S. Costa, N. A. G. Bandeira, B. Le Guennic, V. Robert, P. Gamez, G. Chastanet, L. O. Frade, L. Gasque, Inorg. Chem. 2011, 50, 5696. 117 B. Le Guennic, N. Ben Amor, D. Maynau, V. Robert, J. Chem. Theory. Comput. 2009, 5 1506. 118 U. Bentrup, K. Harms, W. Massa, J. Pebler, Solid State Sci. 2000, 2, 373. 99 119 J. L. Manson, S. H. Lapidus, P. W. Stephens, P. K. Peterson, K. E. Carreiro, H. I. Southerland, T. Lancaster, S. J. Blundell, A. J. Steele, P. A. Goddard, F. L. Pratt, J. Singleton, Y. Kohama, R. D. McDonald, R. E. Del Sesto, N. A. Smith, J. Bendix, S. A. Zvyagin, J. Kang, C. Lee, M. H. Whangbo, V. S. Zapf, A. Plonczak, Inorg. Chem. 2011, 50, 5990. 120 F. Awwadi, S. F. Haddad, R. D. Willett, B. Twamley, Cryst. Growth. Des. 2010, 10, 158. 121 M. Perić, M. Zlatar, S. Grubišić, M. Gruden-Pavlović, Polyhedron, 2012, 42, 89. 122 J. L. Manson, J. A. Schlueter, K. A. Funk, H. I. Southerland, B. Twamley, T. Lancaster, S. J. Blundell, P. J. Baker, F. L. Pratt, J. Singleton, R. D. McDonald, P. A. Goddard, P. Sengupta, C. D. Batista, L. Ding, M. H. Whangbo, I. Franke, S. Cox, C. Baines, D. Trial, J. Am. Chem. Soc. 2009, 131, 6733. 123 J. L. Manson, M. M. Conner, J. A. Schlueter, A. C. McConnell, H. I. Southerland, I. Malfant, T. Lancaster, S. J. Blundell, M. L. Brooks, F. L. Pratt, J. Singleton, R. D. McDonald, C. Lee, M. H. Whangbo, Chem. Mater. 2008, 20, 7408. 124 J. Tang, J. S. Costa, A. Golobič, B. Kozlevčar, A. Robertazzi, A. V. Vargiu, P. Gamez, J. Reedijk, Inorg. Chem. 2009, 48, 5473. 125 S. G. Kang, K. Nam, K. S. Min, U. Lee, Bull. Korean Chem. Soc. 2011, 32, 1037. 126 D. J. Cárdenas, J. M. Cuerva, M. Alías, E. Bunuel, A. G. Campana, Chem. Eur. J. 2011, 17, 8318. 127 R. T. Butcher, M. M. Turnbull, Inorg. Chem. 2010, 49, 427. 128 A. Shapiro, C. P. Landee, M. M. Turnbull, J. Jornet, M. Deumal, J. J. Novoa, M. A. Robb, W. Lewis, J. Am. Chem. Soc. 2007, 129, 952. 129 R. D. Willett, C. Galeriu, M. M. Turnbull, B. Twamley, Inorg. Chem. 2004, 43, 3804. 130 I. Muga, J. M. Gutiérrez-Zorrilla, P. Vitoria, P. Román, F. Lloret, Polyhdron, 2002, 21, 2631. 131 S. X. Liu, C. Ambrus, S. Dolder, A Neels, S. Decurtins, Inorg. Chem. 2006, 45, 9622. 132 F. Fabrizi de Biani, E. Ruiz, J. Cano, J. J. Novoa, S. Alvarez, Inorg. Chem. 2000, 39, 3221. 100 133 E.J. Baerends, T. Ziegler, J. Autschbach, D. Bashford, A. Bérces, F.M. Bickelhaupt, C. Bo, P.M. Boerrigter, L. Cavallo, D.P. Chong, L. Deng, R.M. Dickson, D.E. Ellis, M. van Faassen, L. Fan, T.H. Fischer, C. Fonseca Guerra, A. Ghysels, A. Giammona, S.J.A. van Gisbergen, A.W. Götz, J.A. Groeneveld, O.V. Gritsenko, M. Grüning, S. Gusarov, F.E. Harris, P. van den Hoek, C.R. Jacob, H. Jacobsen, L. Jensen, J.W. Kaminski, G. van Kessel, F. Kootstra, A. Kovalenko, M.V. Krykunov, E. van Lenthe, D.A. McCormack, A. Michalak, M. Mitoraj, J. Neugebauer, V.P. Nicu, L. Noodleman, V.P. Osinga, S. Patchkovskii, P.H.T. Philipsen, D. Post, C.C. Pye, W. Ravenek, J.I. Rodríguez, P. Ros, P.R.T. Schipper, G. Schreckenbach, J.S. Seldenthuis, M. Seth, J.G. Snijders, M. Solà, M. Swart, D. Swerhone, G. te Velde, P. Vernooijs, L. Versluis, L. Visscher, O. Visser, F. Wang, T.A. Wesolowski, E.M. van Wezenbeek, G. Wiesenekker, S.K. Wolff, T.K. Woo, A.L. Yakovlev, ADF2010.01. 134 S. R. Niketić, Kj. Rasmussen (1977) The Consistent Force Field: A Documentation, Lecture Notes in Chemistry, Springer, Berlin, Germany. 135 F. Neese, Orca, an ab initio, DFT and semiempirical electronic structure package, version 2.8, revision 15. Max-Planck-Institute für Bioanorganische Chemie, Germany, Mulhein, 2009. 136 L. Laaksonen, gOpenMol, Center for Scientific Computing, Espoo, Finland, 2004. 101 12 Prilog Prilog A Perturbaciona teorija u kvantnoj hemiji Prilog B Racah parametri B i C (cm-1), konstante spin-orbitalnog kuplovanja (cm-1), orbitalni redukcioni faktori i LF parametri (1000cm-1 ) dobijeni iz LFDFT-a za seriju ispitivanih model-sistema Prilog C Racah parametri B i C (cm-1), konstante spin-orbitalnog kuplovanja (cm-1), orbitalni redukcioni faktori i LF parametri (1000cm-1 [NiCl ) dobijeni iz LFDFT-a za sintetisani 3(Hdabco)2] + Prilog D Parametri istezanja veza za ispitivane strukture [{Mn(3- Fpcyd)(MeOH)(phen)}2(µ-3-Fpcyd)2], 1 [{Mn(3-Fpcyd)(EtOH)(phen)}2(µ- 3-Fpcyd)2], 2 [{Mn(3- Clpcyd)(MeOH)(phen)}2(µ-3-Clpcyd)2], 3 [{Mn(4- lpcyd)(EtOH)(phen)}2(µ-4-Clpcyd)2 ] 4 Prilog E Parametri savijanja uglova za ispitivane strukture [{Mn(3- Fpcyd)(MeOH)(phen)}2(µ-3-Fpcyd)2], 1 [{Mn(3-Fpcyd)(EtOH)(phen)}2(µ- 102 3-Fpcyd)2], 2 [{Mn(3- Clpcyd)(MeOH)(phen)}2(µ-3-Clpcyd)2], 3 [{Mn(4- lpcyd)(EtOH)(phen)}2(µ-4-Clpcyd)2 Prilog F ] 4 Parametri torzionih uglova za ispitivane strukture [{Mn(3- Fpcyd)(MeOH)(phen)}2(µ-3-Fpcyd)2], 1 [{Mn(3-Fpcyd)(EtOH)(phen)}2(µ- 3-Fpcyd)2], 2 [{Mn(3- Clpcyd)(MeOH)(phen)}2(µ-3-Clpcyd)2], 3 [{Mn(4- lpcyd)(EtOH)(phen)}2(µ-4-Clpcyd)2 ] 4 Prilog G Parametri elektrostatičkih interakcija za ispitivane strukture [{Mn(3- Fpcyd)(MeOH)(phen)}2(µ-3-Fpcyd)2], 1 [{Mn(3-Fpcyd)(EtOH)(phen)}2(µ- 3-Fpcyd)2], 2 [{Mn(3- Clpcyd)(MeOH)(phen)}2(µ-3-Clpcyd)2], 3 [{Mn(4- lpcyd)(EtOH)(phen)}2(µ-4-Clpcyd)2 ] 4 Prilog H Van der Waals parametri za ispitivane strukture [{Mn(3- Fpcyd)(MeOH)(phen)}2(µ-3-Fpcyd)2], 1 [{Mn(3-Fpcyd)(EtOH)(phen)}2(µ- 3-Fpcyd)2], 2 [{Mn(3- Clpcyd)(MeOH)(phen)}2(µ-3-Clpcyd)2], 3 [{Mn(4- lpcyd)(EtOH)(phen)}2(µ-4-Clpcyd)2 ] 4 103 Prilog A Perturbaciona teorija u kvantnoj hemiji U okviru ovog dodatka predstavljene su ukratko osnove perturbacione teorije. Schrödinger- ova jednačina za neperturbovan sistem može se prikazati jednačinom D.1. 00 ˆ φφ EH = D.1 Ukoliko se u izraz unese operator za perturbaciju Schrödinger-ova jednačina se može izraziti na drugačiji način. ( )( ) ( )( ).........ˆˆ 221022102210 φλλφφλλφλλφφλ ++++=+++ EEEVH D. 2 Energija E0 predstavlja energiju nultog reda, odnosno energiju sistema bez perturbacije, dok E1 i E2 predstavljaju korekcije prvog i drugog reda. Funkcija φ0 predstavlja neperturbovanu talasnu funkciju. Da bi se dobio izraz za korekciju prvog reda, neophodno je izvući sve članove sa koeficijentom λ. ( ) ( )011001 ˆˆ φλλφφλφλ EEVH +=+ D.3 Ukoliko se koeficijenti svedu na jedinicu i ako se cela jednačina pomnoži sa talasnom funkcijom φ0 100 ˆ EV =φφ , a zatim integrali po koordinatama prostora dobija se izraz D.4. D.4 104 Shodno tome moguće je dobiti korekciju energije drugog reda, ako se izvuku svi članovi sa λ2 ( ) ( )0221122021222 ˆˆ φλφλφλφλφλ EEEVH ++=+ . D.5 Množenjem celog izraza sa φ0 i integraljenjem, dobija se izraz za korekciju drugog reda. 210 ˆ EV =φφ D.6 Ako se funkcija φ1 napiše kao ciΨi, a zatim ceo izraz D.3 pomnoži sa Ψi 00 ˆ EcVEc iiii =Ψ+ φ uz integraljenje dobija se izraz D.7. D.7 Kombinovanjem izraza D.6 i D.7 korekcija energije drugog reda može se izraziti jednačinom D.8 i ii EE VV E − ΨΨ = 0 00 2 ˆˆ φφ D.8 105 Prilog B Tabela P1-Racah parametri B i C (cm-1), konstante spin-orbitalnog kuplovanja (cm-1), orbitalni redukcioni faktori i LF parametri (1000cm-1 ) dobijeni iz LFDFT-a za seriju ispitivanih model-sistema Struktura B C ξ orf ‹dRxyR׀Â׀dRxyR› ‹dRyzR׀Â׀dRyzR› ‹dRz2R׀Â׀dRz2R› ‹dRxzR׀Â׀dRxzR› ‹dRx2- y2R׀Â׀dRx2- y2R› [CrFR5R]P2- 655,211 3190,558 273 0,8848 -24,4783 -28,8586 -13,5383 -28,8586 - 24,4783 [CrClR5R]P2- 655,211 3190,558 273 0,8640 -24,4783 -28,8586 -13,5383 -28,8586 - 24,4783 [CrBrR5R]P2- 655,211 3190,558 273 0,8308 -24,4783 -28,8586 -13,5383 -28,8586 - 24,4783 [CrIR5R]P2- 396,083 2759,379 273 0,8426 -16,1683 -19,4556 -8,1078 -19,4556 - 16,1683 [CrClR2RFR3R]P2- 585,708 3106,750 273 0,8866 -19,1858 -27,0892 -15,5945 -27,0892 - 19,1858 [VFR5R]P3- 371,109 4200,687 167 0,8739 -52,8143 -55,5356 -44,4666 -55,5357 - 52,8143 [VClR5R]P3- 557,210 2868,871 167 0,9742 -46,5957 -48,3161 -40,2428 -48,3161 - 46,5957 [VBrR5R]P3- 546,460 2803,922 167 0,9352 -43,3308 -44,8798 -37,6797 -44,8798 - 43,3308 106 Tabela P1 nastavak [FeF5] 748,638 3- 3692,281 410 0,9344 -38,8575 -36,0988 -43,6321 -36,0988 - 38,8575 [FeCl5] 614,303 3- 3458,530 410 0,8932 -32,3054 -29,4352 -35,2852 -29,4352 - 32,3054 [FeBr5] 563,714 3- 3402,792 410 0,8660 -30,0278 -27,2759 -32,6733 -27,2759 - 30,0278 [NiF5] 741,618 3- 3648,723 649 0,8590 -84,9269 -82,8656 -91,3452 -82,8656 - 84,9269 [NiCl5] 461,766 3- 3398,325 649 0,7684 -71,0772 -67,8820 -75,7370 -67,8820 - 71,0772 [NiBr5] 344,343 3- 3384,521 649 0,7228 -66,0982 -62,7196 -70,4194 -62,7196 - 66,0982 [NiCl2F3] 622,917 2- 3628,674 649 0,8297 -79,5332 -74,6607 -81,2018 -74,6607 - 79,5332 [NiBr2F3] 584,159 2- 3717,715 649 0,8182 -76,8266 -71,0689 -75,8544 -75,8544 - 76,8266 [NiI2F3] 580,847 2- 3642,660 649 0,8070 -73,5335 -67,7963 -72,7198 -67,7963 - 73,5335 107 Prilog C Tabela P2-Racah parametri B i C (cm-1), konstante spin-orbitalnog kuplovanja (cm-1), orbitalni redukcioni faktori i LF parametri (1000cm-1) dobijeni iz LFDFT-a za sintetisani [NiCl3(Hdabco)2] + Metod CF B C ξ orf ‹dRxyR׀Â׀dRxyR› ‹dRyzR׀Â׀dRyzR› ‹dRz2R׀Â׀dRz2R› ‹dRz2R׀Â׀ dRx2- y2R› ‹dRxzR׀Â׀dRxzR› ‹dRx2- y2R׀Â׀dRx2- y2R› ‹dRz2R׀Â׀ dRx2-y2R› DR3 468,668 2301,230 643 1,0000 0,8704 6,1713 -7,4818 ----- 6,1713 0,8704 ----- LDA CR2R(min) 474,407 2303,931 643 1,0000 6,0267 1,6430 -1,2600 2,8573 6,3718 -5,9913 2,8573 CR2R(TS) 470,607 2270,725 643 1,0000 6,2318 0,1325 -1,2369 4,1208 5,9511 -4,5849 4,1208 DR3 535,472 2841,440 643 1,0000 -0,8520 4,2501 -3,4837 ----- 4,2501 -0,8520 ----- OPBE CR2R(min) 537,494 2830,745 643 1,0000 4,3681 0,2201 -1,6829 1,3665 4,6502 -4,1357 1,3665 CR2R(TS) 529,550 2803,651 643 1,0000 4,5111 -1,0348 -1,4227 2,2852 4,2826 -3,2223 2,2852 DR3 561,277 2192,360 643 1,0000 -9,4791 -5,5394 -13,1511 ----- -5,5394 -9,4791 ----- B3LYP CR2R(min) 568,393 2191,989 643 1,0000 -5,5669 -8,9014 -10,3897 1,3958 -5,3578 -12,7834 1,3958 CR2R(TS) 561,703 2163,741 643 1,0000 -5,4636 -9,7581 -10,3694 2,0376 -5,6711 -12,0136 2,0376 108 Prilog D Tabela P3 – Parametri istezanja veza za ispitivane strukture [{Mn(3- Fpcyd)(MeOH)(phen)}2(µ-3-Fpcyd)2], 1 [{Mn(3-Fpcyd)(EtOH)(phen)}2(µ- 3-Fpcyd)2], 2 [{Mn(3- Clpcyd)(MeOH)(phen)}2(µ-3-Clpcyd)2], 3 [{Mn(4- lpcyd)(EtOH)(phen)}2(µ-4-Clpcyd)2 ] 4 kr (kJ/mol Å2 r) 0 (Å) kr (kJ/mol Å2 r) 0 Mn-N (Å) 1256,03 2,150 Cp-C 921,09 p 1,40 Mn-O 418,68 2,240 Cp 962,95 -H 1,00 Mn-N 628,01 ph 2,300 Cp 628,01 -Cl 1,76 O-C 628,01 sp3 1,390 Cp 628,01 -F 1,34 O-H 628,01 0,870 Nph-C 1046,69 p 1,38 Csp3-C 2653,23 sp3 1,340 Nph -C 1046,69 sp 1,29 Csp3 628,01 -H 0,960 N- C 1046,69 sp 1,17 109 Prilog E Tabela P4 – Parametri savijanja uglova za ispitivane strukture [{Mn(3-Fpcyd)(MeOH)(phen)}2(µ-3-Fpcyd)2], 1 [{Mn(3- Fpcyd)(EtOH)(phen)}2(µ-3-Fpcyd)2], 2 [{Mn(3- Clpcyd)(MeOH)(phen)}2(µ- 3-Clpcyd)2], 3 [{Mn(4-lpcyd)(EtOH)(phen)}2(µ-4-Clpcyd)2 ] 4 kθ (kJ/mol rad2 θ (rad) ) kθ (kJ/mol rad2 θ (rad) ) Mn-Nph -C 125,60 p 2,18 Cp-Cp-N 209,34 ph 2,18 Mn-N-C 1256,03 sp 2,69 Cp-Cp-C 209,34 p 2,18 Csp-Nph -C 3349,40 p 2,09 H-Cp-C 209,34 p 2,18 N-Mn-N 209,34 1,58 Cl-Cp-C 209,34 p 2,18 Nph-Mn-Nph 20,93 2,82 F-Cp-C 209,34 p 2,18 O-Mn-N 209,34 1,49 H-O-C 2052,34 sp3 1,96 O-Mn- N 209,34 ph cis 3,01 H-Csp3 259,24 -O 2,18 O-Mn-N 1046,69 trans 3,04 H-Csp3 192,93 -H 1,86 N-Mn-N 20,93 ph cis 1,84 H-Cp-N 293,07 ph 2,18 N-Mn-N 209,34 trans 2,75 Csp3-Csp3 259,24 -O 1,96 Mn-O-C 753,66 sp3 2,15 Csp3-Csp3 259,24 -H 1,91 Mn-O-H 293,11 1,99 Nph -Csp 2093,38 -N 3,06 Cp-Nph-C 209,34 p 1,83 110 Prilog F Tabela P5 – Parametri torzionih uglova za ispitivane strukture [{Mn(3-Fpcyd)(MeOH)(phen)}2(µ-3-Fpcyd)2], 1 [{Mn(3- Fpcyd)(EtOH)(phen)}2(µ-3-Fpcyd)2], 2 [{Mn(3- Clpcyd)(MeOH)(phen)}2(µ- 3-Clpcyd)2], 3 [{Mn(4-lpcyd)(EtOH)(phen)}2(µ-4-Clpcyd)2 ] 4 kΦ n (kJ/mol) kΦ n (kJ/mol) N-Mn-Nph-C 0,004 sp 12 Cp-Cp-Cp-C 29,307 p -2 O-Mn-N-C 4,187 sp 12 H-Cp-Cp 523,344 -H -2 Mn-N-Csp-N 41,868 ph -2 H-Cp-Cp 314,006 -X -2 Csp-N-Cp 0,029 -N -2 Mn-Nph-Cp 4,187 -H -2 Nph-Mn-O-C 0,029 sp3 3 H-Cp-Cp-C 14,654 p -2 Nph-Mn-Nph-C 20,934 p 3 Nph-Cp-Cp-N 29,307 ph -2 O-Mn-A -C 29,307 p 3 Cp-A -Cp 29,307 -H -2 Cp-Nph-Cp-C 29,307 p -2 H-Csp3 13,942 -O-Mn 3 Mn-Nph-Cp-C 523,344 p -2 H-Csp3-Csp3 0,084 -H 3 Nph-Cp-Cp-C 29,307 p -2 O-Csp3-Csp3 0,167 -H 3 N-Mn-N-C 4,187 sp 3 O-Mn-Nph-C 4,187 sp 3 Nph-Mn-N-C 0,004 sp 3 Mn-Nph-Csp 41,868 -N -2 X=Cl, F 111 Prilog G Tabela P6– Parametri elektrostatičkih interakcija za ispitivane strukture [{Mn(3-Fpcyd)(MeOH)(phen)}2(µ-3-Fpcyd)2], 1 [{Mn(3- Fpcyd)(EtOH)(phen)}2(µ-3-Fpcyd)2], 2 [{Mn(3- Clpcyd)(MeOH)(phen)}2(µ- 3-Clpcyd)2], 3 [{Mn(4-lpcyd)(EtOH)(phen)}2(µ-4-Clpcyd)2 ] 4 esu esu Mn 0,882 N -0,048 Cl -0,130 C 0,401 sp3 F -0,216 Cp, C 0,131 sp O -0,310 H -0,116 112 Prilog H Tabela P7 – Van der Waals parametri za ispitivane strukture [{Mn(3-Fpcyd)(MeOH)(phen)}2(µ-3-Fpcyd)2], 1 [{Mn(3- Fpcyd)(EtOH)(phen)}2(µ-3-Fpcyd)2], 2 [{Mn(3- Clpcyd)(MeOH)(phen)}2(µ- 3-Clpcyd)2], 3 [{Mn(4-lpcyd)(EtOH)(phen)}2(µ-4-Clpcyd)2 ] 4 ε (kJ/mol) r* (Å) ε (kJ/mol) r* Mn-C (Å) 0,54 p 3,24 F-N 0,25 3,42 Mn-Csp 0,54 3 3,10 Cl-H 0,46 3,53 N-C 0,21 p 3,76 F-H 0,21 3,10 N- Csp 0,21 3 3,76 Mn- C 0,59 sp 3,24 N-Mn 0,21 3,72 Mn-O 2,26 3,94 N- C 0,21 sp 3,76 Mn-Mn 0,63 2,60 Cp- C 0,17 p 3,88 Csp3- C 0,17 sp 3,88 Csp3- Csp 0,17 3 3,88 Csp- C 0,17 sp 3,88 Csp3- C 0,17 p 3,88 H- C 0,21 sp 3,34 H-H 0,21 3,00 Cl- C 0,42 sp 3,97 H- Csp 0,21 3 3,34 F- C 0,21 sp 3,54 H- C 0,21 p 3,34 O- C 0,21 sp 3,64 H-Mn 0,54 2,70 Cp- C 0,17 sp 3,88 N-H 0,21 3,32 Cl-N 0,50 3,85 N-N 0,21 3,60 H-O 0,25 3,24 Cl-Cl 1,00 4,00 Cl-O 0,50 3,85 F-F 0,25 3,20 F-O 0,21 3,34 Cl- C 0,42 p 3,97 O- C 0,21 p 3,64 F- C 0,21 p 3,54 O-N 0,21 3,56 Cl- Csp 0,42 3 3,97 O-Mn 0,59 3,08 F- Csp 0,21 3 3,54 O- Csp 1,97 3 3,64 Cl-Mn 0,50 3,28 O-O 0,21 3,48 F-Mn 0,38 2,90 H-H 0,21 3,00 113 Biografija Marko Perić rođen je 25.09.1983. godine u Beogradu, gde je završio osnovnu i srednju školu. Diplomirao je 2007. godine na Hemijskom fakultetu Univerziteta u Beogradu sa prosečnom ocenom 9.5 i ocenom 10 na diplomskom radu. Iste godine upisao je doktorske studije na Hemijskom fakultetu Univerziteta u Beogradu, pri Katedri za Opštu i neorgansku hemiju. Od marta 2007. godine do januara 2011. godine bio je stipendista Ministarstva nauke i tehnološkog razvoja Republike Srbije. Od januara 2011. godine Marko Perić je zaposlen u Centru za hemiju, Institut za hemiju, tehnologiju i metalurgiju (IHTM), Univerzitet u Beogradu. Rezultate svojih istraživanja objavio je u 5 radova1-5 . 1. S. Grubišić, M. Gruden-Pavlović, D. Radanović, M. Perić, S. Niketić, J. Mol. Struc. 2009, 919, 54-58. 2. M. Perić, M. Zlatar, S. Niketić, M. Gruden-Pavlović, S. Grubišić, Monatsh. Chem. 2011, 142, 585-592. 3. M. Perić, M. Zlatar, M. Gruden-Pavlović, S. Grubišić, Monatsh. Chem. 2012, 143, 569- 577. 4. Lj. Andjelković, M. Perić, M. Zlatar, S. Grubišić, M. Gruden-Pavlović, Tetrahedron Lett. 2012, 53. 794-799. 5. M. Perić, M. Zlatar, S. Grubišić, M. Gruden-Pavlović, Polyhedron, 2012, 42, 89-94.