[Octopus-users] Segmentation fault with TDOutputOccAnalysis?

[Son, Sang-Kil]

Dear Octopus authors and users,

     I'm a newbie in octopus and I'm trying to test H2+ TD calculation.   
The static calculation of H2+ was okay for both 1D (soft Coulomb  
potential) and 3D (TM potential), but for TD calculation, I've got  
"Segmentation fault" whenever I put TDOutputOccAnalysis=yes option.   
The same error also occurred with TDOutputAcceleration=yes for the 3D  
case.  Is that an installation problem or internal problem, or  
something wrong in my input file?  Just in case, I'll attach inp for 1D  
calculation of H2+.  I installed octopus 1.4 on SGI altix350 and  
Origin2000, and error occurred on both.  Does anybody suffer from the  
same problem?  Thanks in advance.

Son, Sang-Kil
Graduate Student
Department of Chemistry
The University of Kansas


######################################################################## 
#######
# 1. Generalities.
######################################################################## 
#######
verbose = 30
DebugLevel = 15
SystemName = "H2+"
Dimensions = 1
CalculationMode = recipe
Units = "a.u"

######################################################################## 
#######
# 2. Species.
######################################################################## 
#######
%Species
   "H" | 1.0 | 1 | "-1/sqrt(1+x^2)"
%

######################################################################## 
#######
# 3. Coordinates.
######################################################################## 
#######
%Coordinates
   "H"	1.0	0.0	 0.0	yes
   "H"	-1.0	0.0	 0.0	yes
%

######################################################################## 
#######
# 4. Velocities.
######################################################################## 
#######
%Velocities
   "H" | 0.0 | 0.0 | 0.0
   "H" | 0.0 | 0.0 | 0.0
%

######################################################################## 
#######
# 5. Mesh description.
######################################################################## 
#######
BoxShape = sphere
Radius  = 30.0
Spacing = 0.1
DerivativesSpace = real_space
OrderDerivatives = 4
DoubleFFTParameter = 2.0

######################################################################## 
#######
# 6. States.
######################################################################## 
#######

SpinComponents = polarized
ExcessCharge   = 1
ExtraStates    = 0
%Occupations
   1.0
   0.0
%

######################################################################## 
#######
# 7. Hamiltonian.
######################################################################## 
#######
NonInteractingElectrons = yes

######################################################################## 
#######
# 8. Exchange and Correlation.
######################################################################## 
#######
XFunctional = "LDA"
CFunctional = "ZER"

######################################################################## 
#######
# 9. SCF
# 10. Unoccupied states.
######################################################################## 
#######

######################################################################## 
#######
# 11. Time Dependent.
######################################################################## 
#######
TDEvolutionMethod = exp_mid
TDExponentialMethod = split

w0 = 0.08565            # a.u.
T = 5868.7078175680     # a.u. --> 80 o.c.
dt = 0.0733588477196    # a.u.

TDMaximumIter = T / dt
TDTimeStep = dt

shape = 3		# ramp
intensity = 0.0377456783	# 5.0E13 W/cm^2
t0 = T / 2
tau1 = 1467.1769543920	# 20 o.c.
tau0 = T - tau1		# 60 o.c.
%TDLasers
1 | 0 | 0 | intensity | w0 | shape | tau0 | t0 | tau1
%
TDGauge = 1 # length gauge assumed.

AbsorbingBoundaries = mask
ABWidth = 0.4

TDOutputMultipoles = yes
TDOutputLaser = yes
TDOutputAcceleration = yes
#TDOutputOccAnalysis = yes
#TDOccupationalAnalysis = yes

MoveIons = no

######################################################################## 
#######
# 12. Geometry optimization.
# 13. Function output for visualization.
# 14. Spectrum calculations.
# 15. Varia.
######################################################################## 
#######