[CP2K-user] Stuck when using BFGS

ma...@gmail.com ma455... at gmail.com
Fri Jun 18 00:44:15 UTC 2021


Hi,

I'm a new user of cp2k v8.2. When I do GEO_OPT using BFGS, it has been 
stuck many times in the step "OPTIMIZATION STEP: 1" as shown below without 
showing any errors. If I don't kill the job manually, the job will just 
keep running without printing any new results. This does not occurs when 
using SCALAPACK instead of ELPA or using LBFGS/CG instead of BFGS. I'm 
wondering is there any way to solve this BFGS issue? Thanks.

---------------------------------------------------------------------------------------------
[image: 微信截图_20210618103738.png]
-----------------------------------------------------------------------------------------------

Regards,
Hongyang
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20210617/97e9c1e2/attachment.htm>
-------------- next part --------------
#Generated by Multiwfn
&GLOBAL
        PROJECT cp2k
        PRINT_LEVEL LOW
        RUN_TYPE  GEO_OPT
&END GLOBAL

&FORCE_EVAL
        METHOD Quickstep
        &SUBSYS
     &CELL
       A     1.1267270000000000E+01    0.0000000000000000E+00    0.0000000000000000E+00
       B    -5.6238762900000001E+00    9.7458406400000008E+00    0.0000000000000000E+00
       C     0.0000000000000000E+00    0.0000000000000000E+00    5.0000000000000007E+01
       PERIODIC  XYZ
       MULTIPLE_UNIT_CELL  1 1 1
     &END CELL
     &COORD
O    9.3431279038092541E+00    1.0950770714940301E+00    4.0039476813519288E+00
O    7.4711918274004088E+00    4.2614833265257683E+00    4.2457746294427521E+00
O    5.6392854967484727E+00    7.5719234950053300E+00    4.0721272687836017E+00
O    1.8562416466350100E+00    1.0316940992690469E+00    4.0858429255945277E+00
O   -3.1124623047520342E-02    4.3807635955833053E+00    3.9726734539672206E+00
O   -1.8214216020822098E+00    7.6454916746580226E+00    4.0000342258769130E+00
O    5.6749681012547573E+00    1.1484144446601381E+00    4.0825295606341969E+00
O    3.8763558026774008E+00    4.3898500581572515E+00    4.0845716161095496E+00
O    1.8319672645438771E+00    7.5028968529451161E+00    4.2423441313837724E+00
Ce   -3.7475582215798431E+00    8.6792523732996898E+00    3.3090769916199019E+00
Ce   -8.2039980593684239E-03    2.1553486186103448E+00    3.3089201853429868E+00
Ce   -1.8617601537769293E+00    5.4330034065263764E+00    3.2933956601785219E+00
Ce    4.7114990974838006E-03    8.6475934278180269E+00    3.2282755367083253E+00
Ce    3.7917950853040385E+00    2.1778630270683346E+00    3.2529835566921768E+00
Ce    1.8881240638893293E+00    5.4151064538052900E+00    3.3097624801348227E+00
Ce    3.7806858521597801E+00    8.6700211793547197E+00    3.3144437018126491E+00
Ce    7.5198918977798730E+00    2.1559224359264073E+00    3.2240963977162189E+00
Ce    5.6378780381887843E+00    5.4389611334034091E+00    3.2373827393742696E+00
O    5.6446842605353229E+00    9.7482011324851037E+00    2.4838178672152775E+00
O    9.3534267264563855E+00    3.2875366712163352E+00    2.4661669745347790E+00
O    7.5291253886657739E+00    6.4666804583784181E+00    2.4894085605989913E+00
O    3.7769691663124485E+00    4.4029953973416823E-02    2.4340206482326048E+00
O    1.9060147134566858E+00    3.2486860728812470E+00    2.4927263003229689E+00
O    3.4782120838758167E-02    6.5116236942470263E+00    2.4186810215629397E+00
O    1.8691781706627804E+00    9.6758159603479363E+00    2.4646838568371314E+00
O    5.6445654114449795E+00    3.2524347317964213E+00    2.4129603604792460E+00
O    3.7167230106473532E+00    6.5464137346705362E+00    2.4803463682864142E+00
O   -3.7483816396908809E+00    8.6607367526560690E+00    8.7779558214418463E-01
O    1.4029657408344588E-03    2.1659308802623354E+00    8.8023430774421330E-01
O   -1.8746269943268190E+00    5.4160434543707181E+00    8.7715395427644827E-01
O    8.4203801565361260E-03    8.6640132755121879E+00    8.7348307036509143E-01
O    3.7586324729197185E+00    2.1691717919024156E+00    8.7308266891602748E-01
O    1.8869325286207601E+00    5.4128540421709221E+00    8.7839012011466366E-01
O    3.7669641788583070E+00    8.6645083836294656E+00    8.7840950052046807E-01
O    7.5155010710345556E+00    2.1570819126447609E+00    8.7525150613329372E-01
O    5.6370429152069494E+00    5.4201443123703203E+00    8.7808195133613054E-01
Ce    1.4895160927608299E-04    9.2726198716391406E-04    7.6806693830516623E-02
Ce   -1.8693891011811101E+00    3.2520429056872668E+00    6.9018655950358887E-02
Ce   -3.7428507260199440E+00    6.5024519524105004E+00    6.1900752042274680E-02
Ce    3.7545074781984753E+00    2.4870650885461441E-03    6.1532462695824609E-02
Ce    1.8815903228373494E+00    3.2528942781697299E+00    7.1912106223871494E-02
Ce    7.9966349569722400E-03    6.4920266665024178E+00    6.2751878251598378E-02
Ce    1.8928133732483137E+00    9.7378357111711455E+00    5.3162820855380863E-02
Ce    5.6365473605648102E+00    3.2482369878671706E+00    4.6646503151263612E-02
Ce    3.7561239594535434E+00    6.5040476132501661E+00    6.3990446590702355E-02
O    9.3913792519538095E+00    1.0827312421606277E+00    4.9236617944757150E+01
O    7.5162225937803804E+00    4.3320296906865554E+00    4.9261612570222589E+01
O    5.6407560853952994E+00    7.5798578177260474E+00    4.9233705865508909E+01
O    1.8780094775763065E+00    1.0845577526591021E+00    4.9259614998027963E+01
O    5.2492428445285134E-03    4.3304533527437457E+00    4.9232562673876267E+01
O   -1.8698464624508639E+00    7.5795360611222247E+00    4.9223120839312401E+01
O    5.6352721802231489E+00    1.0818704188122941E+00    4.9216953448098970E+01
O    3.7587566241189245E+00    4.3320802126149207E+00    4.9221207923096564E+01
O    1.8856955584202886E+00    7.5795500027703540E+00    4.9236514491362570E+01
O    1.4628211883358853E-03    1.1220858524953194E-03    4.7680772153688679E+01
O   -1.8752827821844038E+00    3.2499506735497570E+00    4.7680479160872117E+01
O    7.5182559932923416E+00    6.4955131883518211E+00    4.7680097326199487E+01
O    3.7533473481619772E+00    9.9687544591464065E-04    4.7678314597669370E+01
O    1.8804380960946392E+00    3.2473316628458999E+00    4.7680681591018583E+01
O    7.4599991240358607E-03    6.4965623933478129E+00    4.7677708836642175E+01
O    7.5132863308975919E+00   -1.3773750558984772E-03    4.7677789482772944E+01
O    5.6387569773143484E+00    3.2500848773701589E+00    4.7677383771728771E+01
O    3.7616767713783688E+00    6.4988911637902769E+00    4.7678913375026461E+01
Ce    9.3924265118789361E+00    1.0832281133936656E+00    4.6846977204854483E+01
Ce    7.5175912592794552E+00    4.3322415907863059E+00    4.6845395381734221E+01
Ce    5.6412672248389137E+00    7.5812479865113369E+00    4.6847580533342949E+01
Ce    1.8779714627443176E+00    1.0834571732581431E+00    4.6848224779495126E+01
Ce    4.7639584127831822E-03    4.3307013875322333E+00    4.6845247592980598E+01
Ce   -1.8704032629068565E+00    7.5801627732450561E+00    4.6844612592124264E+01
Ce    5.6342943924855202E+00    1.0815344137013723E+00    4.6842445717604647E+01
Ce    3.7582042928654453E+00    4.3324554074296016E+00    4.6844022215106499E+01
Ce    1.8857628607536494E+00    7.5791278360109624E+00    4.6841073247036626E+01
O   -3.7469734050465102E+00    8.6638634578793461E+00    4.6035241323929931E+01
O    2.8387494090480387E-03    2.1663240239554176E+00    4.6034263817069764E+01
O   -1.8722726990642264E+00    5.4142166630791637E+00    4.6034468177576244E+01
O    9.1807066012225252E-03    8.6627465748787689E+00    4.6034881667230970E+01
O    3.7571842720216493E+00    2.1659598522375929E+00    4.6035544417154291E+01
O    1.8830598992695178E+00    5.4146903883208513E+00    4.6033912846556397E+01
O    3.7647976310375357E+00    8.6635107130300817E+00    4.6034993866788923E+01
O    7.5147520388043283E+00    2.1667237659718803E+00    4.6036136541293537E+01
O    5.6399234824742308E+00    5.4149548306958089E+00    4.6036257848643260E+01
O    1.6933982193357759E+00    4.6269325169307391E+00    6.0061529664787567E+00
H    9.0150887654215339E-01    4.0356847471242565E+00    5.9582586700521842E+00
H    2.0231533669370938E+00    4.5497406427034877E+00    6.9345973321188970E+00
O    3.2678109559070236E+00    2.2967881768945455E+00    5.8442328405360282E+00
H    2.5776958292079164E+00    1.6323795541722526E+00    5.5622318914580164E+00
H    2.8074664464950874E+00    3.1650194057965297E+00    5.7266769733679537E+00
H    4.1959622746477905E+00    8.0854899802893758E+00    6.4276014703239834E+00
O    3.4979602955742126E+00    8.5936343264719159E+00    5.9336227302635747E+00
H    2.7794454180860155E+00    7.9639083603473573E+00    5.6729036374782060E+00
H    6.4455733370629673E+00    4.5551618118633757E+00    5.6683575781085009E+00
O    5.6839728522048567E+00    4.7663547707486975E+00    6.2863350182339550E+00
H    4.8923512679613284E+00    4.4213366073812486E+00    5.8116959217513164E+00
O   -3.5695745074799552E+00    8.3150166109391428E+00    6.0099803041509903E+00
H   -4.4122376376543793E+00    7.9404562938606071E+00    6.3451014476424232E+00
H   -2.9589235791180681E+00    7.5563293580879165E+00    5.8877337177161255E+00
O    2.7823234463467544E-01    7.1633606824924154E+00    6.5667011256975343E+00
H    9.4708380864815955E-01    7.0347995124220333E+00    5.8488901862297267E+00
H   -2.8998686347072401E-02    8.1227936040836060E+00    6.4229651632566265E+00
H   -9.5867127216654435E-01    6.2114911162893574E+00    6.1557193370297991E+00
O   -1.8384581212982885E+00    5.7694826531958956E+00    5.9269057928363074E+00
H   -1.8381118501853677E+00    4.8407271688896660E+00    6.2774289753115493E+00
O   -7.4246434314074072E-02    2.4509840794311959E+00    5.8617412231068942E+00
H    1.3446639357527418E-01    1.9779995280440872E+00    6.7157287076080872E+00
H   -1.0253854933833457E+00    2.7322680232720358E+00    6.0003881821674678E+00
O   -4.6067669049121068E-01    9.5942253852003123E+00    6.2288697167682008E+00
H    5.3837359048715561E+00    4.9534397868788305E-01    5.4892086007138943E+00
H    5.3651888224395314E+00    3.2384009621535942E-01    7.0708162776439289E+00
O    8.7256678026609809E+00    3.2610385790338960E+00    6.4007455017520591E+00
H    8.1566959416053493E+00    3.5337990679521387E+00    5.6219513027639758E+00
H    8.3889192054110957E+00    2.3989705927608251E+00    6.7998313034518096E+00
O    2.8019615164157305E+00    3.3869228562905267E+00    8.2361928256834780E+00
H    3.2381131165190946E+00    2.7625055995902126E+00    7.6081028662536143E+00
H    3.0809903069467022E+00    3.0963476425079031E+00    9.1441137504378407E+00
O   -4.5510047697649919E-01    4.3919613866836782E+00    9.2766007858917785E+00
H   -1.2344398034272066E+00    4.9048715524126916E+00    8.9147767062317307E+00
H   -2.8653107165633063E-01    3.5315189131208724E+00    8.7873718624501880E+00
H    4.5572119414078767E+00    7.3798475052692920E+00    8.2882629041685192E+00
O    5.3011974695921058E+00    7.2857228682577739E+00    7.5983985668240077E+00
H    5.3135292904100631E+00    6.3557577116370698E+00    7.2609594234035555E+00
H    7.3519704381925866E-01    6.9719894527057455E+00    8.2268317294222140E+00
O    9.7051259907557497E-01    6.7264167937546002E+00    9.1627587260749692E+00
H    7.2316603753472730E-01    5.7706870315213976E+00    9.1938191782738201E+00
O    8.1377667750650726E+00    9.5902515971158497E-01    7.3892044357769606E+00
H    8.4571634558449436E+00    1.7683566375435197E-01    6.8812722815855443E+00
H    7.3146978909227895E+00    7.0632104370552329E-01    7.8581439603589507E+00
O    6.2240870337411840E+00    3.5393178429044490E+00    8.7734572602752827E+00
H    6.1808913502803229E+00    4.0118589957083897E+00    7.9080729389276199E+00
H    7.1646277297381413E+00    3.5324604605161936E+00    9.0934030801529122E+00
H    5.6225887588708572E+00    1.9375789895838611E+00    8.6703507535649837E+00
O    5.4737699823663393E+00    9.5007509553964264E-01    8.6445953804568667E+00
H    4.4918740099673888E+00    7.5345608362974470E-01    8.6013711480942128E+00
O    4.0630996145779374E-01    2.0333536723325465E+00    8.4015978848106325E+00
H    3.3821064053370009E-01    1.4102198767361436E+00    9.1893774644599340E+00
H    1.3388703645868090E+00    2.3692127095173521E+00    8.3952658594026257E+00
O   -2.5606040829870418E+00    5.8198826759995317E+00    8.5170196414321815E+00
H   -2.5098117205161858E+00    6.0222253551315612E+00    7.5584536460997187E+00
H   -2.8658987667720974E+00    6.6466787816873865E+00    8.9525986615802609E+00
O    2.9875694475225596E+00    1.2143763682679246E-01    8.3204094103360866E+00
H   -2.7827982139861174E+00    9.4465628956809731E+00    7.4418326051796297E+00
H   -3.0853565372486900E+00    9.2618619837750877E+00    8.9444881128557974E+00
O    9.2257020885944585E-01    3.7547954543599333E+00    1.1525764185714026E+01
H    4.3849282340880080E-01    4.0838912843656372E+00    1.0733055140322529E+01
H    1.8344492659299200E-01    3.4952027804647736E+00    1.2131931562942780E+01
O    3.1819983867537278E+00    2.5965382238731971E+00    1.0718896462514337E+01
H    3.5558771746130406E+00    1.7167900251724371E+00    1.0990299199568005E+01
H    2.2873954073970797E+00    2.7599058825325993E+00    1.1113040646503800E+01
H    3.9158684743407628E+00    7.5699724994600581E+00    1.0336217226068118E+01
O    3.5281577615474800E+00    7.8298685094402876E+00    9.4421387330213218E+00
H    2.6536271637993711E+00    7.3678581232609917E+00    9.3930191181274729E+00
H    5.5677067460062828E+00    4.2562189786256797E+00    1.0197653110455288E+01
O    5.0827912192741396E+00    4.4650402531738260E+00    1.1040473324683996E+01
H    4.2812676170328725E+00    3.8809929588256789E+00    1.0951974661113535E+01
O   -4.1528981667312355E+00    7.8167626956331002E+00    9.6867802646054937E+00
H    6.4541213011161709E+00    7.5014032011338099E+00    9.0252648905436406E+00
H   -4.1076852331399483E+00    7.1599244602795880E+00    1.0430092835997771E+01
O   -8.3132127788350108E-01    7.2786887307869161E+00    1.1279404131969377E+01
H   -2.1581677973832758E-01    7.1467006485180722E+00    1.0509683121530051E+01
H   -9.3701603292453983E-01    8.2657865567802720E+00    1.1339504486731608E+01
H   -2.3436209031980373E+00    6.3690434206087279E+00    1.1420044764256929E+01
O   -3.2405946390474942E+00    5.9465973188694123E+00    1.1525647660528621E+01
H    8.1018355678978118E+00    5.0452051967216196E+00    1.1135805470483314E+01
O    2.2359241213965036E-01    4.9295171424845663E-01    1.0556958994108021E+01
H   -5.0404430463310135E+00    9.3409670580502304E+00    1.0363119435657286E+01
H    7.6990251487007966E-01    8.0668110774873436E-01    1.1316169839012158E+01
O    4.7756313712863809E+00    3.0596536622165188E-01    1.1203210001612984E+01
H    5.3461808749102193E+00    4.8503034960649605E-01    1.0416551713143736E+01
H    5.3439932405514377E+00    5.0654877191831960E-01    1.1991135814135481E+01
O    8.3604286566449488E+00    3.3039835612087689E+00    1.0390676469859077E+01
H   -1.9822221929007733E+00    3.4562596993130903E+00    1.0092839940528208E+01
H    8.3178171988235512E+00    2.3295869459208909E+00    1.0588828896034256E+01
O    2.6405229170231359E+00    3.5279663196527999E+00    1.3772240025714554E+01
H    2.3144300525067578E+00    4.1308218205656244E+00    1.3082561754743049E+01
H    3.6282280254379669E+00    3.6591942890469897E+00    1.3823789187156745E+01
O   -1.0360648724001817E+00    3.4828257764086090E+00    1.3439916730324878E+01
H   -1.9513649622105305E+00    3.8298220978534561E+00    1.3589119239029591E+01
H   -9.7522000247860263E-01    2.6008997368740800E+00    1.3875226266136297E+01
H    3.8337919780373135E+00    6.9397619886165014E+00    1.2409075494009839E+01
O    4.5957733478291489E+00    7.0797038592489354E+00    1.1767898755686200E+01
H    4.8577523897303312E+00    6.1526969959528444E+00    1.1532968379846730E+01
H   -2.4131824953810435E-01    6.5326093178068030E+00    1.2734205774286909E+01
O    9.7970115555477369E-02    6.0725970190180458E+00    1.3551901461761437E+01
H   -2.0857140646785854E-01    5.1367696151678670E+00    1.3521104339632963E+01
O    8.0877278078786787E+00    5.8529498714536488E-01    1.0800133820140049E+01
H    2.9535997485377292E+00    9.7322272394422704E+00    1.0194416412439503E+01
H    8.3494055735440469E+00    2.8683854963770400E-01    1.1688571806821184E+01
O    5.3842067873431958E+00    3.5370389275838705E+00    1.3521479986180386E+01
H    5.4218329583710316E+00    3.8137577151627040E+00    1.2563376006402629E+01
H    6.1580374924013270E+00    4.0203205678375600E+00    1.3911695952103969E+01
H    5.5564748606162331E+00    1.7938913289257612E+00    1.3707773964006384E+01
O    5.7284898629473320E+00    8.1250086526183696E-01    1.3751021254421463E+01
H    4.8446378165108746E+00    3.6097352251707115E-01    1.3850242016549970E+01
O    1.6383430076025198E+00    1.0921707514201520E+00    1.2881978980934790E+01
H    8.4568477397422714E-01    1.0456849216940101E+00    1.3483048795157815E+01
H    2.0870390591317469E+00    1.9329306905921950E+00    1.3165781659301400E+01
O    7.8161924139566290E+00    4.7789478078274046E+00    1.3995827658020385E+01
H    7.8390943699315310E+00    5.4228368783232002E+00    1.3238883435041114E+01
H    8.0336816840320893E+00    5.2597658598712123E+00    1.4828542374342264E+01
O   -2.1502658309682490E+00    9.0915734018355625E+00    1.3823827176869324E+01
H   -2.8989016427042920E+00    9.6291787319219182E+00    1.3448138218059606E+01
H   -2.5304682949510364E+00    8.6137739021701893E+00    1.4578376389047049E+01
O    1.3011924912978936E+00    2.5554341249943731E+00    1.5983318545470235E+01
H    1.6365869861979470E+00    3.1465541670521064E+00    1.5263932665787756E+01
H    1.1676507651025281E+00    3.1342712725545936E+00    1.6786932553164348E+01
O    3.7796618655789778E+00    1.8960375665151044E+00    1.6652511184983123E+01
H    4.2168072131156142E+00    1.0252617616833974E+00    1.6451096161421692E+01
H    2.8577322307300896E+00    1.9519992057097104E+00    1.6282027408337797E+01
H    3.2593591877697894E+00    6.3173943434047226E+00    1.4388834316889238E+01
O    2.7953268561701647E+00    6.7145550691906362E+00    1.3617377480756840E+01
H    1.8551896486993786E+00    6.4087295654003311E+00    1.3629029147247854E+01
H    6.5937881602246762E+00    3.3579231768061133E+00    1.6498365312941235E+01
O    5.7605191124262172E+00    3.7443668659988401E+00    1.6872414765871355E+01
H    5.0455605880719014E+00    3.1019922065819010E+00    1.6615016905462234E+01
O    5.7169092223372608E+00    8.1466278729058690E+00    1.4331615974438996E+01
H    5.4923083496686935E+00    7.8376281710383635E+00    1.3426558844417722E+01
H    5.3057748875270532E+00    7.4605027440984415E+00    1.4921380433853455E+01
O   -2.0287753733319916E-01    6.9932273419342703E+00    1.6021531087016900E+01
H   -1.1181038896630520E-01    6.7210051110795241E+00    1.5063784107964910E+01
H   -2.2693135473718337E-01    7.9918750200288535E+00    1.6021271063624596E+01
H   -1.5549435025885261E+00    5.7820382003345401E+00    1.6346427289484275E+01
O   -2.3788476721759428E+00    5.2738187396355585E+00    1.6559199042602266E+01
H   -2.2334009792435339E+00    4.3243967683197937E+00    1.6357040232668076E+01
O   -4.8633794546670200E-01    1.0046914328319225E+00    1.4577829105970611E+01
H    1.0882276052484997E+01    5.9635207550950723E-03    1.4590308491330459E+01
H    5.9148059262134581E-02    1.3759095802181249E+00    1.5316369164218647E+01
O   -8.7200103779089974E-02    9.6609968979837042E+00    1.6249399338945132E+01
H    5.9427298180340200E+00    2.9925369865024493E-01    1.5432611976951419E+01
H    6.0974452830622949E+00    1.1171983196663403E-01    1.7052499045456319E+01
O    8.1798215178383575E+00    2.6045411720882901E+00    1.6040089915880895E+01
H    8.1775948237327594E+00    2.3739545647015738E+00    1.5097319803076021E+01
H    8.5693916369981977E+00    1.7962144655918799E+00    1.6499709843369398E+01
O    3.1739018976299720E+00    2.1665954207314160E+00    1.9257809876608871E+01
H    3.4191180363207554E+00    2.0742968478991384E+00    1.8289211593589673E+01
H    4.0177392555664593E+00    2.4948667637819839E+00    1.9645634956303457E+01
O    1.2520513161072679E+00    4.0133475826911384E+00    1.8245036108578638E+01
H    3.8921610669019346E-01    4.2638737743867337E+00    1.8686914953184186E+01
H    1.7837678315989516E+00    3.4482991367980613E+00    1.8856553082161383E+01
H    3.6804805254672321E+00    6.2459983382380937E+00    1.6424906020803231E+01
O    4.4479191279438615E+00    6.1424129682308903E+00    1.5785479035534454E+01
H    4.9632097165550837E+00    5.3695496853752056E+00    1.6109951759964350E+01
H    1.4375172570860362E+00    6.6470003644140423E+00    1.6753121505813905E+01
O    2.2382994767181179E+00    6.4271211336829914E+00    1.7284898609862132E+01
H    1.9945792984877897E+00    5.5749314803361472E+00    1.7725295147632977E+01
O    9.0850769472086252E+00    4.9125263822450932E-01    1.7361082769453645E+01
H    3.6780976240965311E+00    9.4435649842976197E+00    1.6841291980425225E+01
H    8.3255516959959746E+00    2.2670741467204358E-01    1.7938068845486075E+01
O    5.9622883633394483E+00    2.8033198238277701E+00    1.9367388906933996E+01
H    5.9376113164680966E+00    3.2613244700749737E+00    1.8471285079045209E+01
H    6.6007439706462607E+00    3.2944605386354642E+00    1.9907754342822663E+01
H    6.4258690126454887E+00    1.0781446797277521E+00    1.9043581464958976E+01
O    6.7235536051148639E+00    2.1924382868455666E-01    1.8640652225688385E+01
H    5.8932503554565063E-01    9.2102353077217920E+00    1.9030572598542474E+01
O   -2.7596845699003771E+00    9.2270108756838933E+00    1.9275844122048920E+01
H   -3.6091767460636355E+00    8.9454125353163612E+00    1.9649028177241501E+01
H    2.9204681941388149E+00    4.6728574197985873E-01    1.9396568408002164E+01
O   -1.1086604446860575E+00    4.9715511815038544E+00    1.9185938325036066E+01
H   -1.6750519410251279E+00    5.0402344041502172E+00    1.8383199825898846E+01
H   -8.9007608922738291E-01    5.9121641575413610E+00    1.9383266568556273E+01
O   -5.3658911185887337E-01    7.7833275961092871E+00    1.9048254217577011E+01
H   -4.6545409927638826E-01    7.5221951645865639E+00    1.8108504588839743E+01
H   -1.4099489713870685E+00    8.2634546669680287E+00    1.9120175087811372E+01
     &END COORD
                #   &VELOCITY #You can set initial atomic velocities in this section
                #   &END VELOCITY
                &KIND O    
                        ELEMENT O 
                        BASIS_SET DZVP-MOLOPT-SR-GTH
                        POTENTIAL GTH-PBE
                &END KIND
                &KIND Ce   
                        ELEMENT Ce
                        BASIS_SET DZVP-MOLOPT-SR-GTH-q12
                        POTENTIAL GTH-PBE-q12
                &END KIND
                &KIND H    
                        ELEMENT H 
                        BASIS_SET DZVP-MOLOPT-SR-GTH
                        POTENTIAL GTH-PBE
                &END KIND
        &END SUBSYS

        &DFT
                BASIS_SET_FILE_NAME  BASIS_MOLOPT
                BASIS_SET_FILE_NAME  BASIS_MOLOPT_UCL
                POTENTIAL_FILE_NAME  POTENTIAL
                WFN_RESTART_FILE_NAME ./cp2k-RESTART.wfn
                CHARGE    0 #Net charge
                MULTIPLICITY    1 #Spin multiplicity
                LSD
                &QS
                        EPS_DEFAULT 1E-10 #This is default. Set all EPS_xxx to values such that the energy will be correct up to this value
                        #EXTRAPOLATION ASPC #Extrapolation for wavefunction during e.g. MD. ASPC is default, PS also be used
                        #EXTRAPOLATION_ORDER 3 #Order for PS or ASPC extrapolation. 3 is default
                &END QS
                &POISSON
                        PERIODIC XYZ #Direction(s) of PBC for calculating electrostatics
                        PSOLVER PERIODIC #The way to solve Poisson equation
                &END POISSON
                &XC
                        &XC_FUNCTIONAL PBE
                        &END XC_FUNCTIONAL
                &END XC
                &MGRID
                        CUTOFF 400
                        REL_CUTOFF 60
                &END MGRID
                &SCF
                        MAX_SCF 20 #Maximum number of steps of inner SCF
                        EPS_SCF 1.0E-05 #Convergence threshold of density matrix of inner SCF
                        SCF_GUESS RESTART #Use wavefunction from WFN_RESTART_FILE_NAME file as initial guess
                        &OT
                                PRECONDITIONER FULL_SINGLE_INVERSE #FULL_SINGLE_INVERSE is also worth to try. FULL_ALL is better but quite expensive for large system
                                MINIMIZER DIIS #CG is worth to consider in difficult cases
                                LINESEARCH 2PNT #1D line search algorithm for CG. 2PNT is default, 3PNT is better but more costly. GOLD is best but very expensive
                        &END OT
                        &OUTER_SCF
                                MAX_SCF 20 #Maximum number of steps of outer SCF
                                EPS_SCF 1.0E-05 #Convergence threshold of outer SCF
                        &END OUTER_SCF
                        &PRINT
                                &RESTART #Use "&RESTART OFF" can prevent generating wfn file
                                        BACKUP_COPIES 0 #Maximum number of backup copies of wfn file
                                &END RESTART
                        &END PRINT
                &END SCF
        &END DFT
&END FORCE_EVAL

&MOTION
        &MD
                ENSEMBLE NVT
                STEPS 200
                TIMESTEP 1.0 #fs. Decrease it properly for high temperature simulation
                TEMPERATURE 298.15 #Initial and maintained temperature (K)
                &THERMOSTAT
                        TYPE NOSE
                &END THERMOSTAT
        &END MD
        &GEO_OPT
                TYPE MINIMIZATION #Search for minimum
                OPTIMIZER BFGS #Can also be CG (more robust for difficult cases) or LBFGS
                &BFGS
                        TRUST_RADIUS 0.2 #Trust radius (maximum stepsize) in Angstrom
                        #     RESTART_HESSIAN T #If read initial Hessian, uncomment this line and specify the file in the next line
                        #     RESTART_FILE_NAME to_be_specified
                &END BFGS
                #&CG
                #        &LINE_SEARCH
                #                TYPE 2PNT
                #        &END LINE_SEARCH
                #&END CG
                MAX_ITER 250 #Maximum number of geometry optimization
                #The following thresholds of geometry convergence are the default ones
                MAX_DR 3E-3 #Maximum geometry change
                RMS_DR 1.5E-3 #RMS geometry change
                MAX_FORCE 4.5E-4 #Maximum force
                RMS_FORCE 3E-4 #RMS force
        &END GEO_OPT
        &PRINT
                &RESTART
                        BACKUP_COPIES 0 #Maximum number of backing up restart file
                &END RESTART
        &END PRINT
&END MOTION
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ????_20210618103738.png
Type: image/png
Size: 8712 bytes
Desc: not available
URL: <https://lists.cp2k.org/archives/cp2k-user/attachments/20210617/97e9c1e2/attachment.png>
-------------- next part --------------
 JOBID| JOBFS                                                  24025363.gadi-pbs
 DBCSR| CPU Multiplication driver                                           XSMM
 DBCSR| Multrec recursion limit                                              512
 DBCSR| Multiplication stack size                                           1000
 DBCSR| Maximum elements for images                                    UNLIMITED
 DBCSR| Multiplicative factor virtual images                                   1
 DBCSR| Use multiplication densification                                       T
 DBCSR| Multiplication size stacks                                             3
 DBCSR| Use memory pool for CPU allocation                                     F
 DBCSR| Number of 3D layers                                               SINGLE
 DBCSR| Use MPI memory allocation                                              F
 DBCSR| Use RMA algorithm                                                      F
 DBCSR| Use Communication thread                                               T
 DBCSR| Communication thread load                                             87
 DBCSR| MPI: My node id                                                        0
 DBCSR| MPI: Number of nodes                                                  48
 DBCSR| OMP: Current number of threads                                         1
 DBCSR| OMP: Max number of threads                                             1
 DBCSR| Split modifier for TAS multiplication algorithm                  1.0E+00


  **** **** ******  **  PROGRAM STARTED AT               2021-06-18 10:31:26.412
 ***** ** ***  *** **   PROGRAM STARTED ON        gadi-cpu-clx-2898.gadi.nci.org
 **    ****   ******    PROGRAM STARTED BY                                hm1876
 ***** **    ** ** **   PROGRAM PROCESS ID                               3076930
  **** **  *******  **  PROGRAM STARTED IN /scratch/ad73/hm1876/cp2k_exercise/6H
                                           2O-NV

 CP2K| version string:                                          CP2K version 8.2
 CP2K| source code revision number:                                  git:310b7ab
 CP2K| cp2kflags: omp libint fftw3 libxc pexsi elpa parallel mpi3 scalapack xsmm
 CP2K|             max_contr=4 plumed2 spglib mkl check_diag libvori libbqb
 CP2K| is freely available from                            https://www.cp2k.org/
 CP2K| Program compiled at                          Thu Jun 3 09:55:52 AEST 2021
 CP2K| Program compiled on                     gadi-cpu-clx-0454.gadi.nci.org.au
 CP2K| Program compiled for                               Linux-x86-64-intel-nci
 CP2K| Data directory path                            /home/561/hm1876/cp2k/data
 CP2K| Input file name                                                  cp2k.inp
 
 GLOBAL| Force Environment number                                              1
 GLOBAL| Basis set file name                                        BASIS_MOLOPT
 GLOBAL| Potential file name                                           POTENTIAL
 GLOBAL| MM Potential file name                                     MM_POTENTIAL
 GLOBAL| Coordinate file name                                      __STD_INPUT__
 GLOBAL| Method name                                                        CP2K
 GLOBAL| Project name                                                       cp2k
 GLOBAL| Run type                                                        GEO_OPT
 GLOBAL| FFT library                                                       FFTW3
 GLOBAL| Diagonalization library                                            ELPA
 GLOBAL| Minimum number of eigenvectors for ELPA usage                        16
 GLOBAL| Orthonormality check for eigenvectors enabled                 5.000E-14
 GLOBAL| Matrix multiplication library                                    SCALAP
 GLOBAL| All-to-all communication in single precision                          F
 GLOBAL| FFTs using library dependent lengths                                  F
 GLOBAL| Grid backend                                                       AUTO
 GLOBAL| Global print level                                                  LOW
 GLOBAL| MPI I/O enabled                                                       T
 GLOBAL| Total number of message passing processes                            48
 GLOBAL| Number of threads for this process                                    1
 GLOBAL| This output is from process                                           0
 GLOBAL| CPU model name             Intel(R) Xeon(R) Platinum 8274 CPU @ 3.20GHz
 GLOBAL| CPUID                                                              1003
 GLOBAL| Compiled for CPUID                                                 1002

 *** HINT in environment.F:872 :: The compiler target flags (x86_avx2)     ***
 *** used to build this binary cannot exploit all extensions of this CPU   ***
 *** model (x86_avx512). Consider compiler target flags as part of FCFLAGS ***
 *** and CFLAGS (ARCH file).                                               ***


 MEMORY| system memory details [Kb]
 MEMORY|                        rank 0           min           max       average
 MEMORY| MemTotal            197732028     197732028     197732028     197732028
 MEMORY| MemFree             173232012     173232012     173232012     173232012
 MEMORY| Buffers                     4             4             4             4
 MEMORY| Cached                3402020       3402020       3402020       3402020
 MEMORY| Slab                  5753828       5753828       5753828       5753828
 MEMORY| SReclaimable          1579588       1579588       1579588       1579588
 MEMORY| MemLikelyFree       178213624     178213624     178213624     178213624

 
 GENERATE|  Preliminary Number of Bonds generated:                             0
 GENERATE|  Achieved consistency in connectivity generation.

 *******************************************************************************
 *******************************************************************************
 **                                                                           **
 **     #####                         ##              ##                      **
 **    ##   ##            ##          ##              ##                      **
 **   ##     ##                       ##            ######                    **
 **   ##     ##  ##   ##  ##   #####  ##  ##   ####   ##    #####    #####    **
 **   ##     ##  ##   ##  ##  ##      ## ##   ##      ##   ##   ##  ##   ##   **
 **   ##  ## ##  ##   ##  ##  ##      ####     ###    ##   ######   ######    **
 **    ##  ###   ##   ##  ##  ##      ## ##      ##   ##   ##       ##        **
 **     #######   #####   ##   #####  ##  ##  ####    ##    #####   ##        **
 **           ##                                                    ##        **
 **                                                                           **
 **                                                ... make the atoms dance   **
 **                                                                           **
 **            Copyright (C) by CP2K developers group (2000-2021)             **
 **                      J. Chem. Phys. 152, 194103 (2020)                    **
 **                                                                           **
 *******************************************************************************


 TOTAL NUMBERS AND MAXIMUM NUMBERS

  Total number of            - Atomic kinds:                                   3
                             - Atoms:                                        261
                             - Shell sets:                                   261
                             - Shells:                                      1200
                             - Primitive Cartesian functions:               1359
                             - Cartesian basis functions:                   3708
                             - Spherical basis functions:                   3216

  Maximum angular momentum of- Orbital basis functions:                        4
                             - Local part of the GTH pseudopotential:          2
                             - Non-local part of the GTH pseudopotential:      3


 SCF PARAMETERS         Density guess:                                   RESTART
                        --------------------------------------------------------
                        max_scf:                                              20
                        max_scf_history:                                       0
                        max_diis:                                              4
                        --------------------------------------------------------
                        eps_scf:                                        1.00E-05
                        eps_scf_history:                                0.00E+00
                        eps_diis:                                       1.00E-01
                        eps_eigval:                                     1.00E-05
                        --------------------------------------------------------
                        level_shift [a.u.]:                                 0.00
                        --------------------------------------------------------
                        Outer loop SCF in use 
                        No variables optimised in outer loop
                        eps_scf                                         1.00E-05
                        max_scf                                               20
                        No outer loop optimization
                        step_size                                       5.00E-01

 BFGS| Use rational function optimization for step estimation:                NO
 BFGS| Use model Hessian for initial guess:                                  YES
 BFGS| Restart Hessian:                                                       NO
 BFGS| Trust radius:                                                       0.378

 *******************************************************************************
 ***                     STARTING GEOMETRY OPTIMIZATION                      ***
 ***                                   BFGS                                  ***
 *******************************************************************************

 Spin 1

 Number of electrons:                                                        564
 Number of occupied orbitals:                                                564
 Number of molecular orbitals:                                               564

 Spin 2

 Number of electrons:                                                        564
 Number of occupied orbitals:                                                564
 Number of molecular orbitals:                                               564

 Number of orbital functions:                                               3216
 Number of independent orbital functions:                                   3216

 Extrapolation method: initial_guess


 SCF WAVEFUNCTION OPTIMIZATION

  ----------------------------------- OT ---------------------------------------
  Minimizer      : DIIS                : direct inversion
                                         in the iterative subspace
                                         using   7 DIIS vectors
                                         safer DIIS on
  Preconditioner : FULL_SINGLE_INVERSE : inversion of 
                                         H + eS - 2*(Sc)(c^T*H*c+const)(Sc)^T
  Precond_solver : DEFAULT
  stepsize       :    0.08000000                  energy_gap     :    0.08000000
  eps_taylor     :   0.10000E-15                  max_taylor     :             4
  ----------------------------------- OT ---------------------------------------

  Step     Update method      Time    Convergence         Total energy    Change
  ------------------------------------------------------------------------------
     1 OT DIIS     0.80E-01    8.1     0.00000529     -2950.2946781928 -2.95E+03

  *** SCF run converged in     1 steps ***


  Electronic density on regular grids:      -1127.9999998265        0.0000001735
  Core density on regular grids:             1127.9999999423       -0.0000000577
  Total charge density on r-space grids:        0.0000001158
  Total charge density g-space grids:           0.0000001158

  Overlap energy of the core charge distribution:               0.00000121103769
  Self energy of the core charge distribution:              -6993.15267366503576
  Core Hamiltonian energy:                                   1889.57626789863343
  Hartree energy:                                            2738.31612832138308
  Exchange-correlation energy:                               -585.03440195885344

  Total energy:                                             -2950.29467819283445

  outer SCF iter =    1 RMS gradient =   0.53E-05 energy =      -2950.2946781928
  outer SCF loop converged in   1 iterations or    1 steps


  Integrated absolute spin density  :                               0.0000000612
  Ideal and single determinant S**2 :                    0.000000       0.000000

 ENERGY| Total FORCE_EVAL ( QS ) energy [a.u.]:            -2950.294685059895528

 --------  Informations at step =     0 ------------
  Optimization Method        =                 BFGS
  Total Energy               =     -2950.2946850599
  Used time                  =               41.096
 ---------------------------------------------------

 --------------------------
 OPTIMIZATION STEP:      1
 --------------------------


More information about the CP2K-user mailing list