Ivan a Peter,

 presiel som cez par starsich clankov o tych Si povrchoch (Chadi, Cohen,
 Artacho atd).

 1) rozdiel energii medzi symetricke-asymtricke dimers je velmi maly, takze 
    vypocty v pripade vacsich systemov budu dost narocne 
 
 2) kedze symetricky dimer je efektivne (slaby) antiferromagnet, vlnova funkcia
    bude multi-reference, tj. minimalne dve konfiguracie na dimer

 3) Artacho aj Goddard modelovali povrch cez cluster a vazby okolo saturovali
    vodikmi. Myslel som, ze je to bluf, ze to nemoze fungovat (ze tam
    budu velke zmeny dlzok vazieb, atd). Ale dal som si robotu to otestovat,
    napodiv uz prvy pokus ktory som urobil pre Si9H12 vysiel necakane 
    dobre: dlzky vazieb suhlasia s experimentami, clustre su stabilne a 
    existuju tam aj dve minima jasne korespondujuce sym and asym strukture,
    vid obrazok dole.  

 Vypocitane energie (a.u., BPW91, 631G*):

 sym   -2612.37060292928  
 asym  -2612.37108996626  

    Takze asym je v beznom GGA nizsi v zhode s clankami Chadi, Cohen atd. 
    Otazka teda je ci v skutocnosti je sym nizsi (ako tvrdi Artacho, Goddard 
    a nove experimenty z Tsukuby).

Navrhujem nasledovny postup:

   a) robit radsej s Si clustrami s H terminatormi, nez s periodickymi
      hranicnymi podmienkami. Je na to viacej dovodov. Jeden vychdza
      z toho ze rozdiel bude radu 0.0005 a.u. takze vsetko
      bude musiet byt ozaj "ciste": totalna energie a fluktuacie budu vacsie
      nez pre H-clusters tak to bude challenging. 
      Dalsi: analyza vlnovych funkcii
      musi byt perfektna (na tom budem potom musiet popracovat) inak mozeme
      lahko dostat hluposti. Je toho este viac, ale toto je to klucove, aby
      sa realisticky dali ocakavat vysledky. 

   b) zacat s Si9H12 (alebo podobnym ak mate lepsi navrh) a ten urobit
      podrobne - su tam minimalne dve veci 
      - zistit ci je sym ozaj nizsie nez asym
      - ak je sym nizsie tak potom overit ci je asym vobec lokalne minimum
        (odhadujem ze nebude, ak by bolo, pozriet sa na barieru sym-asym) 

   c) ked budeme mat toto hotove pozriet sa na vacsie clustre s dvomi,
      styrmi, pripadne viacerymi dimermi a tak odhadnut bulk limit,
      konvergencia sa da sledovat aj v GGA takze budeme mat trendy na
      porovnanie

   d) z hladiska buducnosti SiH systemy su perspektivne (silany v chemii,
      clustre v polovodicoch, atd) tazke toto moze byt dobra platforma na
      dalsi projekt

 Takze mi dajte vediet co si o tom myslite.

 H paper citam - komentare budu nasledovat. 

 Teraz komunikacia: Ivan tebe posielam emaily do Tsukuby
                    Peter tebe na ncsa adresu

 Dajte mi vediet ci to funguje

						Lubos

-------------------------------------------------------------------------
Peter,

a. v mojej scratch kniznici na modi4 je si25 directory kde najdes
   gamess behy pre elongated and spherical struktury, gamess pouzivam 
   LEN na vlnove funkcie, nie na geometrie - ako dostat vlnove funckie
   pre tieto clusters je velmi tricky, musim ti o tom napisat viac neskor,
   je to dlhe vysvetlovanie :-(
   (ked urobis bezny run s bazou aka tam je tak to spravidla nekonverguje)

b. geometrie treba robit v ramci DFT

-------------------------------------------------------------------------