diff --git a/QA/tests/rt_tddft_restart/na1.nw b/QA/tests/rt_tddft_restart/na1.nw new file mode 100644 index 0000000000..98374cd834 --- /dev/null +++ b/QA/tests/rt_tddft_restart/na1.nw @@ -0,0 +1,73 @@ +## +## Real-time TDDFT restart example. +## +## File 1 of 2. +## +## Na atom under dummy continuous wave excitation. +## +## Example consists of two input decks: +## +## 1: Runs simulation to completion for a short tmax +## and saves restart files along the way. +## +## 2: Runs another simulation with longer tmax starting +## from previously saved restart file. +## +## How to generate dipole moment: +## (a) na1.nw => na1.nwo => z1.dat +## (b) na2.nw => na2.nwo => z2.dat +## (c) cat z1.dat z2.dat > zboth.dat +## +## Practically you would likely keep the old .nwo file +## and re-run the same input deck but with "load vectors" added. +## +## Note: Do NOT use the top-level "restart" option, instead +## "start" a new run every time and let the RT-TDDFT read in +## the self-contained run.rt_restart file +## +title "RT-TDDFT restarts" + +echo +scratch_dir ./scratch +permanent_dir ./perm + +start run +echo + +charge 0 + +geometry "system" units angstroms noautosym nocenter noautoz + Na 0.0 0.0 0.0 +end + +set geometry "system" + +basis + Na library 6-31G +end + +dft + odft + mult 2 + xc b3lyp +end +task dft energy + +rt_tddft + dt 0.2 + tmax 20.0 + + nrestarts 10 + + print dipole field energy charge + + field "driver" + type cw + max 0.0001 + frequency 0.086 + polarization z + end + + excite "system" with "driver" +end +task dft rt_tddft diff --git a/QA/tests/rt_tddft_restart/na1.nwo.gz b/QA/tests/rt_tddft_restart/na1.nwo.gz new file mode 100644 index 0000000000..2281513068 Binary files /dev/null and b/QA/tests/rt_tddft_restart/na1.nwo.gz differ diff --git a/QA/tests/rt_tddft_restart/na2.nw b/QA/tests/rt_tddft_restart/na2.nw new file mode 100644 index 0000000000..505c3bf720 --- /dev/null +++ b/QA/tests/rt_tddft_restart/na2.nw @@ -0,0 +1,75 @@ +## +## Real-time TDDFT restart example. +## +## File 2 of 2. +## +## Na atom under dummy continuous wave excitation. +## +## Example consists of two input decks: +## +## 1: Runs simulation to completion for a short tmax +## and saves restart files along the way. +## +## 2: Runs another simulation with longer tmax starting +## from previously saved restart file. +## +## How to generate dipole moment: +## (a) na1.nw => na1.nwo => z1.dat +## (b) na2.nw => na2.nwo => z2.dat +## (c) cat z1.dat z2.dat > zboth.dat +## +## Practically you would likely keep the old .nwo file +## and re-run the same input deck but with "load vectors" added. +## +## Note: Do NOT use the top-level "restart" option, instead +## "start" a new run every time and let the RT-TDDFT read in +## the self-contained run.rt_restart file +## +title "RT-TDDFT restarts" + +echo +scratch_dir ./scratch +permanent_dir ./perm + +start run +echo + +charge 0 + +geometry "system" units angstroms noautosym nocenter noautoz + Na 0.0 0.0 0.0 +end + +set geometry "system" + +basis + Na library 6-31G +end + +dft + odft + mult 2 + xc b3lyp +end +task dft energy + +rt_tddft + dt 0.2 + tmax 60.0 + + load restart + + nrestarts 10 + + print dipole field energy charge + + field "driver" + type cw + max 0.0001 + frequency 0.086 + polarization z + end + + excite "system" with "driver" +end +task dft rt_tddft diff --git a/QA/tests/rt_tddft_restart/na2.nwo.gz b/QA/tests/rt_tddft_restart/na2.nwo.gz new file mode 100644 index 0000000000..e69f821583 Binary files /dev/null and b/QA/tests/rt_tddft_restart/na2.nwo.gz differ diff --git a/QA/tests/rt_tddft_restart/z1.dat b/QA/tests/rt_tddft_restart/z1.dat new file mode 100644 index 0000000000..5eb4096565 --- /dev/null +++ b/QA/tests/rt_tddft_restart/z1.dat @@ -0,0 +1,126 @@ +# ====================================== +# NWChem Real-time TDDFT output parser +# ====================================== +# +# Runtime options +# --------------- +# compare = False +# target = dipole +# tolerance = 1e-05 +# delim = +# polarization = z +# header = True +# zero = False +# tag = +# geom = system +# spin = total +# +# File 1 +# ------ +# Filename: "na1.nwo" +# Number of data points: 99 +# Time range: [0.0: 19.6] +# +#----------------------------------------- +# Time [au] dipole [au] +#----------------------------------------- +# + 0.0000000000e+00 3.4962030177e-16 + 2.0000000000e-01 1.9662207310e-08 + 4.0000000000e-01 1.1774927641e-07 + 6.0000000000e-01 3.7185657923e-07 + 8.0000000000e-01 8.5809644936e-07 + 1.0000000000e+00 1.6506436874e-06 + 1.2000000000e+00 2.8214469837e-06 + 1.4000000000e+00 4.4401435391e-06 + 1.6000000000e+00 6.5742435734e-06 + 1.8000000000e+00 9.2895189251e-06 + 2.0000000000e+00 1.2650417247e-05 + 2.2000000000e+00 1.6720395179e-05 + 2.4000000000e+00 2.1562191829e-05 + 2.6000000000e+00 2.7238050485e-05 + 2.8000000000e+00 3.3809803444e-05 + 3.0000000000e+00 4.1338779719e-05 + 3.2000000000e+00 4.9885646807e-05 + 3.4000000000e+00 5.9510336642e-05 + 3.6000000000e+00 7.0272082321e-05 + 3.8000000000e+00 8.2229514767e-05 + 4.0000000000e+00 9.5440841742e-05 + 4.2000000000e+00 1.0996418336e-04 + 4.4000000000e+00 1.2585802719e-04 + 4.6000000000e+00 1.4318163601e-04 + 4.8000000000e+00 1.6199528596e-04 + 5.0000000000e+00 1.8236034071e-04 + 5.2000000000e+00 2.0433917809e-04 + 5.4000000000e+00 2.2799489515e-04 + 5.6000000000e+00 2.5339074072e-04 + 5.8000000000e+00 2.8058937937e-04 + 6.0000000000e+00 3.0965216402e-04 + 6.2000000000e+00 3.4063849811e-04 + 6.4000000000e+00 3.7360528205e-04 + 6.6000000000e+00 4.0860649833e-04 + 6.8000000000e+00 4.4569305367e-04 + 7.0000000000e+00 4.8491290263e-04 + 7.2000000000e+00 5.2631132615e-04 + 7.4000000000e+00 5.6993124096e-04 + 7.6000000000e+00 6.1581352480e-04 + 7.8000000000e+00 6.6399735666e-04 + 8.0000000000e+00 7.1452047486e-04 + 8.2000000000e+00 7.6741924740e-04 + 8.4000000000e+00 8.2272858750e-04 + 8.6000000000e+00 8.8048184236e-04 + 8.8000000000e+00 9.4071071518e-04 + 9.0000000000e+00 1.0034451945e-03 + 9.2000000000e+00 1.0687135172e-03 + 9.4000000000e+00 1.1365422718e-03 + 9.6000000000e+00 1.2069566880e-03 + 9.8000000000e+00 1.2799810110e-03 + 1.0000000000e+01 1.3556388404e-03 + 1.0200000000e+01 1.4339534106e-03 + 1.0400000000e+01 1.5149478183e-03 + 1.0600000000e+01 1.5986451163e-03 + 1.0800000000e+01 1.6850681641e-03 + 1.1000000000e+01 1.7742392483e-03 + 1.1200000000e+01 1.8661795957e-03 + 1.1400000000e+01 1.9609088633e-03 + 1.1600000000e+01 2.0584446068e-03 + 1.1800000000e+01 2.1588017695e-03 + 1.2000000000e+01 2.2619923211e-03 + 1.2200000000e+01 2.3680251370e-03 + 1.2400000000e+01 2.4769060690e-03 + 1.2600000000e+01 2.5886381118e-03 + 1.2800000000e+01 2.7032216430e-03 + 1.3000000000e+01 2.8206547464e-03 + 1.3200000000e+01 2.9409335378e-03 + 1.3400000000e+01 3.0640523645e-03 + 1.3600000000e+01 3.1900038426e-03 + 1.3800000000e+01 3.3187788063e-03 + 1.4000000000e+01 3.4503662248e-03 + 1.4200000000e+01 3.5847530658e-03 + 1.4400000000e+01 3.7219241176e-03 + 1.4600000000e+01 3.8618618752e-03 + 1.4800000000e+01 4.0045465803e-03 + 1.5000000000e+01 4.1499563909e-03 + 1.5200000000e+01 4.2980675955e-03 + 1.5400000000e+01 4.4488548548e-03 + 1.5600000000e+01 4.6022914884e-03 + 1.5800000000e+01 4.7583497458e-03 + 1.6000000000e+01 4.9170009441e-03 + 1.6200000000e+01 5.0782154275e-03 + 1.6400000000e+01 5.2419624084e-03 + 1.6600000000e+01 5.4082097481e-03 + 1.6800000000e+01 5.5769236646e-03 + 1.7000000000e+01 5.7480683822e-03 + 1.7200000000e+01 5.9216058240e-03 + 1.7400000000e+01 6.0974954529e-03 + 1.7600000000e+01 6.2756942636e-03 + 1.7800000000e+01 6.4561568642e-03 + 1.8000000000e+01 6.6388356380e-03 + 1.8200000000e+01 6.8236810052e-03 + 1.8400000000e+01 7.0106417368e-03 + 1.8600000000e+01 7.1996652019e-03 + 1.8800000000e+01 7.3906974855e-03 + 1.9000000000e+01 7.5836834087e-03 + 1.9200000000e+01 7.7785664862e-03 + 1.9400000000e+01 7.9752887957e-03 + 1.9600000000e+01 8.1737907549e-03 diff --git a/QA/tests/rt_tddft_restart/z2.dat b/QA/tests/rt_tddft_restart/z2.dat new file mode 100644 index 0000000000..a05ba77c94 --- /dev/null +++ b/QA/tests/rt_tddft_restart/z2.dat @@ -0,0 +1,228 @@ +# ====================================== +# NWChem Real-time TDDFT output parser +# ====================================== +# +# Runtime options +# --------------- +# compare = False +# target = dipole +# tolerance = 1e-05 +# delim = +# polarization = z +# header = True +# zero = False +# tag = +# geom = system +# spin = total +# +# File 1 +# ------ +# Filename: "na2.nwo" +# Number of data points: 201 +# Time range: [19.6: 59.6] +# +#----------------------------------------- +# Time [au] dipole [au] +#----------------------------------------- +# + 1.9600000000e+01 8.1737907551e-03 + 1.9800000000e+01 8.3740108865e-03 + 2.0000000000e+01 8.5758856783e-03 + 2.0200000000e+01 8.7793495377e-03 + 2.0400000000e+01 8.9843348023e-03 + 2.0600000000e+01 9.1907718166e-03 + 2.0800000000e+01 9.3985890930e-03 + 2.1000000000e+01 9.6077135336e-03 + 2.1200000000e+01 9.8180706142e-03 + 2.1400000000e+01 1.0029584473e-02 + 2.1600000000e+01 1.0242177929e-02 + 2.1800000000e+01 1.0455772473e-02 + 2.2000000000e+01 1.0670288192e-02 + 2.2200000000e+01 1.0885643632e-02 + 2.2400000000e+01 1.1101755666e-02 + 2.2600000000e+01 1.1318539453e-02 + 2.2800000000e+01 1.1535908518e-02 + 2.3000000000e+01 1.1753774888e-02 + 2.3200000000e+01 1.1972049306e-02 + 2.3400000000e+01 1.2190641522e-02 + 2.3600000000e+01 1.2409460651e-02 + 2.3800000000e+01 1.2628415478e-02 + 2.4000000000e+01 1.2847414655e-02 + 2.4200000000e+01 1.3066366790e-02 + 2.4400000000e+01 1.3285180462e-02 + 2.4600000000e+01 1.3503764114e-02 + 2.4800000000e+01 1.3722025830e-02 + 2.5000000000e+01 1.3939873043e-02 + 2.5200000000e+01 1.4157212285e-02 + 2.5400000000e+01 1.4373949000e-02 + 2.5600000000e+01 1.4589987399e-02 + 2.5800000000e+01 1.4805230362e-02 + 2.6000000000e+01 1.5019579446e-02 + 2.6200000000e+01 1.5232934985e-02 + 2.6400000000e+01 1.5445196182e-02 + 2.6600000000e+01 1.5656261166e-02 + 2.6800000000e+01 1.5866027004e-02 + 2.7000000000e+01 1.6074389723e-02 + 2.7200000000e+01 1.6281244292e-02 + 2.7400000000e+01 1.6486484544e-02 + 2.7600000000e+01 1.6690003116e-02 + 2.7800000000e+01 1.6891691486e-02 + 2.8000000000e+01 1.7091440120e-02 + 2.8200000000e+01 1.7289138706e-02 + 2.8400000000e+01 1.7484676463e-02 + 2.8600000000e+01 1.7677942549e-02 + 2.8800000000e+01 1.7868826568e-02 + 2.9000000000e+01 1.8057219026e-02 + 2.9200000000e+01 1.8243011692e-02 + 2.9400000000e+01 1.8426097831e-02 + 2.9600000000e+01 1.8606372350e-02 + 2.9800000000e+01 1.8783731784e-02 + 3.0000000000e+01 1.8958074102e-02 + 3.0200000000e+01 1.9129298404e-02 + 3.0400000000e+01 1.9297304602e-02 + 3.0600000000e+01 1.9461993136e-02 + 3.0800000000e+01 1.9623264701e-02 + 3.1000000000e+01 1.9781020019e-02 + 3.1200000000e+01 1.9935159723e-02 + 3.1400000000e+01 2.0085584352e-02 + 3.1600000000e+01 2.0232194385e-02 + 3.1800000000e+01 2.0374890234e-02 + 3.2000000000e+01 2.0513572216e-02 + 3.2200000000e+01 2.0648140532e-02 + 3.2400000000e+01 2.0778495211e-02 + 3.2600000000e+01 2.0904535965e-02 + 3.2800000000e+01 2.1026162040e-02 + 3.3000000000e+01 2.1143272147e-02 + 3.3200000000e+01 2.1255764534e-02 + 3.3400000000e+01 2.1363537142e-02 + 3.3600000000e+01 2.1466487890e-02 + 3.3800000000e+01 2.1564515104e-02 + 3.4000000000e+01 2.1657518118e-02 + 3.4200000000e+01 2.1745397911e-02 + 3.4400000000e+01 2.1828057668e-02 + 3.4600000000e+01 2.1905403254e-02 + 3.4800000000e+01 2.1977343593e-02 + 3.5000000000e+01 2.2043790875e-02 + 3.5200000000e+01 2.2104660526e-02 + 3.5400000000e+01 2.2159870974e-02 + 3.5600000000e+01 2.2209343334e-02 + 3.5800000000e+01 2.2253001070e-02 + 3.6000000000e+01 2.2290769627e-02 + 3.6200000000e+01 2.2322576068e-02 + 3.6400000000e+01 2.2348348825e-02 + 3.6600000000e+01 2.2368017607e-02 + 3.6800000000e+01 2.2381513387e-02 + 3.7000000000e+01 2.2388768390e-02 + 3.7200000000e+01 2.2389716074e-02 + 3.7400000000e+01 2.2384291136e-02 + 3.7600000000e+01 2.2372429465e-02 + 3.7800000000e+01 2.2354067994e-02 + 3.8000000000e+01 2.2329144436e-02 + 3.8200000000e+01 2.2297597056e-02 + 3.8400000000e+01 2.2259364532e-02 + 3.8600000000e+01 2.2214385891e-02 + 3.8800000000e+01 2.2162600539e-02 + 3.9000000000e+01 2.2103948484e-02 + 3.9200000000e+01 2.2038370803e-02 + 3.9400000000e+01 2.1965810256e-02 + 3.9600000000e+01 2.1886211939e-02 + 3.9800000000e+01 2.1799523918e-02 + 4.0000000000e+01 2.1705697840e-02 + 4.0200000000e+01 2.1604689437e-02 + 4.0400000000e+01 2.1496458784e-02 + 4.0600000000e+01 2.1380970287e-02 + 4.0800000000e+01 2.1258192503e-02 + 4.1000000000e+01 2.1128097861e-02 + 4.1200000000e+01 2.0990662281e-02 + 4.1400000000e+01 2.0845864719e-02 + 4.1600000000e+01 2.0693686764e-02 + 4.1800000000e+01 2.0534112408e-02 + 4.2000000000e+01 2.0367127951e-02 + 4.2200000000e+01 2.0192721965e-02 + 4.2400000000e+01 2.0010885298e-02 + 4.2600000000e+01 1.9821611142e-02 + 4.2800000000e+01 1.9624895125e-02 + 4.3000000000e+01 1.9420735275e-02 + 4.3200000000e+01 1.9209131828e-02 + 4.3400000000e+01 1.8990086961e-02 + 4.3600000000e+01 1.8763604529e-02 + 4.3800000000e+01 1.8529689796e-02 + 4.4000000000e+01 1.8288349174e-02 + 4.4200000000e+01 1.8039590075e-02 + 4.4400000000e+01 1.7783421010e-02 + 4.4600000000e+01 1.7519851924e-02 + 4.4800000000e+01 1.7248894649e-02 + 4.5000000000e+01 1.6970563451e-02 + 4.5200000000e+01 1.6684875673e-02 + 4.5400000000e+01 1.6391852394e-02 + 4.5600000000e+01 1.6091518956e-02 + 4.5800000000e+01 1.5783905228e-02 + 4.6000000000e+01 1.5469045677e-02 + 4.6200000000e+01 1.5146979275e-02 + 4.6400000000e+01 1.4817749255e-02 + 4.6600000000e+01 1.4481402678e-02 + 4.6800000000e+01 1.4137989954e-02 + 4.7000000000e+01 1.3787564457e-02 + 4.7200000000e+01 1.3430182300e-02 + 4.7400000000e+01 1.3065902199e-02 + 4.7600000000e+01 1.2694785414e-02 + 4.7800000000e+01 1.2316895839e-02 + 4.8000000000e+01 1.1932300193e-02 + 4.8200000000e+01 1.1541068184e-02 + 4.8400000000e+01 1.1143272527e-02 + 4.8600000000e+01 1.0738988844e-02 + 4.8800000000e+01 1.0328295491e-02 + 4.9000000000e+01 9.9112733066e-03 + 4.9200000000e+01 9.4880052376e-03 + 4.9400000000e+01 9.0585759344e-03 + 4.9600000000e+01 8.6230714809e-03 + 4.9800000000e+01 8.1815793194e-03 + 5.0000000000e+01 7.7341883247e-03 + 5.0200000000e+01 7.2809890038e-03 + 5.0400000000e+01 6.8220738718e-03 + 5.0600000000e+01 6.3575379876e-03 + 5.0800000000e+01 5.8874795052e-03 + 5.1000000000e+01 5.4120000901e-03 + 5.1200000000e+01 4.9312051754e-03 + 5.1400000000e+01 4.4452040910e-03 + 5.1600000000e+01 3.9541100260e-03 + 5.1800000000e+01 3.4580397654e-03 + 5.2000000000e+01 2.9571132595e-03 + 5.2200000000e+01 2.4514531857e-03 + 5.2400000000e+01 1.9411846032e-03 + 5.2600000000e+01 1.4264346880e-03 + 5.2800000000e+01 9.0733255025e-04 + 5.3000000000e+01 3.8400921180e-04 + 5.3200000000e+01 -1.4340221883e-04 + 5.3400000000e+01 -6.7476627879e-04 + 5.3600000000e+01 -1.2099449297e-03 + 5.3800000000e+01 -1.7487974217e-03 + 5.4000000000e+01 -2.2911802181e-03 + 5.4200000000e+01 -2.8369470211e-03 + 5.4400000000e+01 -3.3859489801e-03 + 5.4600000000e+01 -3.9380350649e-03 + 5.4800000000e+01 -4.4930524662e-03 + 5.5000000000e+01 -5.0508469213e-03 + 5.5200000000e+01 -5.6112629654e-03 + 5.5400000000e+01 -6.1741441111e-03 + 5.5600000000e+01 -6.7393328747e-03 + 5.5800000000e+01 -7.3066705961e-03 + 5.6000000000e+01 -7.8759971224e-03 + 5.6200000000e+01 -8.4471504982e-03 + 5.6400000000e+01 -9.0199667152e-03 + 5.6600000000e+01 -9.5942795142e-03 + 5.6800000000e+01 -1.0169920281e-02 + 5.7000000000e+01 -1.0746718133e-02 + 5.7200000000e+01 -1.1324500217e-02 + 5.7400000000e+01 -1.1903092093e-02 + 5.7600000000e+01 -1.2482318115e-02 + 5.7800000000e+01 -1.3062001771e-02 + 5.8000000000e+01 -1.3641965994e-02 + 5.8200000000e+01 -1.4222033334e-02 + 5.8400000000e+01 -1.4802025925e-02 + 5.8600000000e+01 -1.5381765267e-02 + 5.8800000000e+01 -1.5961071945e-02 + 5.9000000000e+01 -1.6539765343e-02 + 5.9200000000e+01 -1.7117663356e-02 + 5.9400000000e+01 -1.7694582137e-02 + 5.9600000000e+01 -1.8270336007e-02