nga errors in gradient calculation

[mqiyang]

I have met a problem when trying to calculate the transient spectra of H2, which is basically repeat the result of https://doi.org/10.1021/acs.jctc.5b00473. The input file is:

title "H2"                                                                                                                                                                           
start h2                                                                                                                                                                             
                                                                                                                                                                                     
echo                                                                                                                                                                                 
scratch_dir ./scratch                                                                                                                                                                
permanent_dir ./perm                                                                                                                                                                 
                                                                                                                                                                                                                                                                                                                                                                      
geometry "system" units angstroms noautoz noautosym nocenter                                                                                                                         
  H     0.00000000      0.00000000      0.00000000                                                                                                                                   
  H     0.7299595       0.00000000      0.00000000                                                                                                                                   
end                                                                                                                                                                                  
                                                                                                                                                                                     
                                                                                                                                                                                     
set geometry system                                                                                                                                                                  
                                                                                                                                                                                     
basis spherical print                                                                                                                                                                
  * library 6-31G                                                                                                                                                                    
end                                                                                                                                                                                  
                                                                                                                                                                                     
dft                                                                                                                                                                                  
  xc hfexch                                                                                                                                                                          
end         

tddft                                                                                                                                                                                
  nroots 5                                                                                                                                                                           
  notriplet                                                                                                                                                                          
  civecs                                                                                                                                                                             
  grad                                                                                                                                                                               
    root 1                                                                                                                                                                           
    print debug                                                                                                                                                                      
  end                                                                                                                                                                                
end                                                                                                                                                                                  
task tddft gradient                                                                                                                                                                  
                                                                                                                                                                                     
unset rt_tddft:*                                                                                                                                                                     
rt_tddft                                                                                                                                                                             
  tmax 1000.0                                                                                                                                                                        
  dt 0.1                                                                                                                                                                             
  load density h2.dmat                                                                                                                                                               
  tag "ref"                                                                                                                                                                          
  print dipole                                                                                                                                                                       
 end                                                                                                                                                                                 
task dft rt_tddft                                                                                                                                                                    
                                                                                                                                                                                     
unset rt_tddft:*                                                                                                                                                                     
rt_tddft                                                                                                                                                                             
  tmax 1000.0                                                                                                                                                                        
  dt 0.1                                                                                                                                                                             
  load density h2.dmat                                                                                                                                                               
  tag "kick_x"                                                                                                                                                                       
  print dipole                                                                                                                                                                       
  field "kick"                                                                                                                                                                       
    type delta                                                                                                                                                                       
    polarization x                                                                                                                                                                   
    max 0.0001                                                                                                                                                                       
  end                                                                                                                                                                                
  excite "system" with "kick"                                                                                                                                                        
 end                                                                                                                                                                                 
task dft rt_tddft

在此输入代码...

But the calculation stopped at the TDDFT GRADIENT block, with error of:

[0] Received an Error in Communication: (-999) 0:nga_matmul_patch:  g_b indices out of range :
TDDFT Energy Check(  1) =       0.55611168249130
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI COMMUNICATOR 3 DUP FROM 0
with errorcode -999.

在此输入代码...

Does this problem come from my input file or the code? I'm using nwchem 6.8 and will appreciate any suggestions.

[niri....@gmail.com]

Hi, Thanks for reporting this issue. I was able to reproduce this problem with the 7.0.0 version as well.  We will address it and get back to you.

Best,

-Niri

Niri Govind

PNNL

[mqiyang]

Hi Niri, thank you for your reply. I also met problem when calculating the ESA of H2+. My input file is:

```

title "H2+"
start h2plus

echo
scratch_dir ./scratch
permanent_dir ./perm

geometry "system" units angstroms noautoz noautosym nocenter
H 0.00000000 0.00000000 0.00000000

H 1.04098906 0.00000000 0.00000000
end
charge 1

set geometry system

basis spherical print
* library 6-31G
end

dft
xc hfexch

odft
mult 2
end

tddft
nroots 10

notriplet
civecs
grad
root 1

end
end
task tddft gradient

rt_tddft
tmax 1000.0
dt 0.1

load density h2plus.dmat

tag "ref"
print dipole
end
task dft rt_tddft

```

I got the error message of:

```

[0] Received an Error in Communication: (-973) 0:ngai_get_common:cannot locate region: X vector [1:0 ,1:1 ]:

--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI COMMUNICATOR 3 DUP FROM 0

with errorcode -973.

```

I also tried other 1e- system it gave the same error code. I've tried both nwchem 6.8 and 7 and met the same issue.

Bests,

Mengqi

But the calculation stopped at the TDDFT GRADIENT block, with error of:

[0] Received an Error in Communication: (-999) 0:nga_matmul_patch:  g_b indices out of range :
TDDFT Energy Check(  1) =       0.55611168249130
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI COMMUNICATOR 3 DUP FROM 0
with errorcode -999.

Does this problem come from my input file or the code? I'm using nwchem 6.8 and will appreciate any suggestions.

[Niri Govind]

Hi Mengqi,

We are looking into this bug. Yes, I checked these cases as well yesterday.

To doubly make sure, I have tested the other cases in our QA suite and they run fine and agree with our reference outputs.

This is a special case just with these cases. We will address this issue and check in a fix soon.

Best,

-Niri

--
You received this message because you are subscribed to the Google Groups "NWChem Forum" group.
To unsubscribe from this group and stop receiving emails from it, send an email to nwchem-forum...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/nwchem-forum/ec90e53e-6c39-48f5-83ee-77c208ec8922n%40googlegroups.com.

[Edoardo Aprà]

Commits
https://github.com/nwchemgit/nwchem/commit/7d012039f26ad35673ac4f33b7ebaf74041924f8
https://github.com/nwchemgit/nwchem/commit/674831de9976b0bb97dca0357abdfd4ae283bfc1
https://github.com/nwchemgit/nwchem/commit/bfc02ce293bcf1fdaab37ca5085fcec0e2f71131
address this bug

Commits have been pushed both to the `master` and `hotfix/release-7-0-0` branches