[Octopus-users] about RandomVelocityTemp variable

[Qin Rao]

Hi, 

I got a problem. When I try to set initial velocity using

 RandomVelocityTemp =273.0

then I got the following message:

    ** Address Error **

    End of diagnostics

What is wrong with this variable? Here I attached my 
input file "inp". 

Thanks for help

Jean

-------------- next part --------------
###############################################################################
# 1. Generalities.
###############################################################################
LinBroadening = 0.06
LinMinEnergy =0.0
LinMaxEnergy = 30.0
LinEnergyStep = 0.01

LinEigenvalues = yes 
LinCasida = yes 
LinPetersilka = yes 

verbose = 300
SystemName = "H2"
CalculationMode = 5
Units = "eVA"

###############################################################################
# 2. Species.
###############################################################################
%Species
  "H"  | 1.0079  | 1 | "tm2" | 0 | 0
%

###############################################################################
# 3. Coordinates.
###############################################################################
%Coordinates
 "H" |  0.3714  | 0.0  |  0.0 | no
 "H" | -0.3714  | 0.0  |  0.0 | no
%
AdjustCoordinates = yes
%MainAxis
0 | 0 | 1
%
OutputCoordinates = no

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

RandomVelocityTemp = 273.0

###############################################################################
# 5. Mesh description.
###############################################################################
BoxShape = 1 # a sphere...
Radius  = 8.0
Zlength = 8.0 # not necessary in this case.
Spacing = 0.3
DerivativesSpace = 0 # Real space derivatives (1 for Fourier space).
OrderDerivatives = 4
DoubleFFTParameter = 2.0

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

SpinComponents = 1
NumberkPoints  = 1
ExcessCharge   = 0
ExtraStates    = 2

ElectronicTemperature = 0.0

###############################################################################
# 7. Hamiltonian.
###############################################################################
NonInteractingElectrons = no
ClassicPotential = no
LocalPotentialSpace = 1
NonLocalPotentialSpace = 0
GridRefinement = 0

###############################################################################
# 8. Exchange and Correlation.
###############################################################################

XFunctional = "PBE"
CFunctional = "PBE"

###############################################################################
# 9. SCF
###############################################################################
MaximumIter = 200
ConvAbsDens = 1e-5
ConvRelDens = 0
ConvAbsEnergy = 0
ConvRelEnergy = 0

LCAOStart = yes
SCFinLCAO = no
LCAOMode = 0 (0 for memory intensive, 1 for cpu intensive)

EigenSolver = cg
EigenSolverInitTolerance = 1e-10
EigenSoverFinalTolerance = 1e-14
EigenSolverFinalToleranceIteration = 20
EigenSolverMaxIter = 25
NuberCG = 3

TypeOfMixing = 2
Mixing = 0.3
BroydenNumber = 3

###############################################################################
# 10. Unoccupied states.
###############################################################################
UnoccNumberStates = 5
UnoccConv = 1.e-3
UnoccMaximumIter = 200

###############################################################################
# 11. Time Dependent.
###############################################################################
TDMaximumIter = 180000
TDSaveIter    = 1000
TDEvolutionMethod = 4
TDExponentialMethod = standard
TDTimeStep = 0.003

MoveIons = 3

%TDLasers
1 | 0 | 0 | 2.745  | 1.610 | 2 | 27 | 50.0
%

TDOutputLaser = yes
TDGauge = 1 # length gauge assumed.

#TDOutputAcceleration = yes   # Note: not for Delta field
#TDOutputOccAnalysis = yes    # Note: not for Delta field
TDOutputElEnergy = yes       #Note: Not for Delta field


AbsorbingBoundaries = mask
AbWidth = 1.0

TDDipoleLmax = 1

###############################################################################
# 12. Photo-electron spectrum.
###############################################################################
CalcPES_rc = no
CalcPes_mask = no
%PES_rc_points
0 | 0 | 5
%

###############################################################################
# 13. Geometry optimization.
###############################################################################
GOMethod = 4
GOTolerance = 0.01
GOMaximumIter = 100

###############################################################################
# 14. Function output for visualization.
###############################################################################
OutputKSPotential = no
OutputDensity     = no
OutputWfs         = no
OutputELF         = no
OutputGeometry    = no

OutputPlaneX = no
OutputPlaneY = no
OutputPlaneZ = no
OutputDX = no
OutputEvery = 1000

###############################################################################
# 15. Spectrum calculations.
###############################################################################
SpecTransMode = "sin"
SpecDampMode  = "exp"
SpecDampFactor = 0.1
SpecStartTime = 0.0
SpecEndTime   = -1.0
SpecEnergyStep = 0.01
SpecMinEnergy = 0.0
SpecMaxEnergy = 30.0
HSPolarization = 'x'


###############################################################################
# 16. Varia.
###############################################################################
PoissonSolver = cg
POLStaticField = 0.001