Error in running Vevacious on susyhit output

[mbeekveld]

Dear readers,

I have tried to install Vevacious 1.2.02 on my system, but I am running into problems. Most of them I have solved, but there is one that I cannot seem to understand.

I have tried to run a standard example, using the xml file provided in the attachment. The error that I get is the following:

Setting up the VevaciousRunner object took 0.004461 seconds.
Please enter the option of homotopy method
1. The polyhedral homotopy
2. The classical linear homotopy
The polyhedral homotopy method
The # of paths followed = 9

The # of roots = 5
The # of real roots = 3
The # of paths followed = 9
The # of blow_up = 0 (before endgame)
The running time : 0.00400 s
The solutions are listed in file "data.roots"

Department of Mathematics, Michigan State University
Contact Info:
Tsung-Lin Lee : leetsung@msu.edu
Tien-Yien Li : li@math.msu.edu
Chih-Hsiung Tsai : tsaichih@msu.edu

Finding all the tree-level extrema and parsing the results took 0.020012 seconds.
rm: cannot remove './Vevacious.pyc': No such file or directory
Traceback (most recent call last):
File "./Vevacious.py", line 29, in <module>
EffectivePotential = VPD.LoopAndThermalCorrectedPotential )
File "/scratch/mbeekveld/computation_code_susy/Vevacious/Vevacious/VevaciousParameterDependent.py", line 872, in __init__
temperatureValue = currentTemperature )
File "/scratch/mbeekveld/computation_code_susy/Vevacious/Vevacious/VevaciousParameterDependent.py", line 775, in __init__
self.functionAtOrigin = FunctionToScale( 0.0, 0.0, 0.0 )
File "/scratch/mbeekveld/computation_code_susy/Vevacious/Vevacious/VevaciousParameterDependent.py", line 688, in LoopAndThermalCorrectedPotential
massesSquared = MassesSquared[ 0 ]( a, b )
File "/scratch/mbeekveld/computation_code_susy/Vevacious/Vevacious/VevaciousParameterDependent.py", line 562, in MassesSquared1
massSquaredElementArray.size ), -1 )
TypeError: 'float' object cannot be interpreted as an index
Minimizing the one-loop potential, calculating the tunneling time or its upper bound, and printing the results took 0.554703 seconds.
Warning! "./MyResult.vout" was not produced!
Vevacious finished. Total running time was 0.579334 seconds.


--------------------------------

Then I tried to replace the line

vcs = VPD.Vevacious(
EffectivePotential = VPD.LoopAndThermalCorrectedPotential )
to

vcs = VPD.Vevacious(
EffectivePotential = VPD.TreeLevelPotential )

and put
vcs.shouldDoThermalCalculations = False

in the Vevacious.py file.

And then I get:

Setting up the VevaciousRunner object took 0.00453 seconds.
Please enter the option of homotopy method
1. The polyhedral homotopy
2. The classical linear homotopy
The classical linear homotopy
The # of paths followed = 9
The # of roots = 5
The # of real roots = 1
The # of paths followed = 9
The # of blow_up = 0 (before endgame)
The running time : 0.00000 s
The solutions are listed in file "data.roots"


Department of Mathematics, Michigan State University
Contact Info:
Tsung-Lin Lee : leetsung@msu.edu
Tien-Yien Li : li@math.msu.edu
Chih-Hsiung Tsai : tsaichih@msu.edu


Finding all the tree-level extrema and parsing the results took 0.016434 seconds.
rm: cannot remove './Vevacious.pyc': No such file or directory
trying to minimize { { vd -> ( 21.9363852672 ), vu -> ( 242.707564543 ) }, TreeLevelPotentialValue -> 139560509.683, EffectivePotentialValue -> 139560509.683 } (T = 0.0)
rolled to { { vd -> ( 0.622100029286 ), vu -> ( 4.146051226 ) }, TreeLevelPotentialValue -> 26832.7623409, EffectivePotentialValue -> 26832.7623409 }
trying to minimize { { vd -> ( 0 ), vu -> ( 0 ) }, TreeLevelPotentialValue -> 0.0, EffectivePotentialValue -> 0.0 } (T = 0.0)
rolled to { { vd -> ( 0.0 ), vu -> ( 0.0 ) }, TreeLevelPotentialValue -> 0.0, EffectivePotentialValue -> 0.0 }
DSB vacuum: { { vd -> ( 0.622100029286 ), vu -> ( 4.146051226 ) }, TreeLevelPotentialValue -> 26832.7623409, EffectivePotentialValue -> 26832.7623409 }
Panic vacuum: { { vd -> ( 0.0 ), vu -> ( 0.0 ) }, TreeLevelPotentialValue -> 0.0, EffectivePotentialValue -> 0.0 }
Global minimum found: { { vd -> ( 0.0 ), vu -> ( 0.0 ) }, TreeLevelPotentialValue -> 0.0, EffectivePotentialValue -> 0.0 }

tunnelingCalculator.npoints being given 0
falseVacuumArray = [0.62210003 4.14605123]
trueVacuumArray = [0. 0.]
Traceback (most recent call last):
File "./Vevacious.py", line 117, in <module>
thermalNotQuantum = False )
File "/scratch/mbeekveld/computation_code_susy/Vevacious/Vevacious/VevaciousParameterDependent.py", line 1282, in CalculateAction
quickTunneling = False )
TypeError: fullTunneling() got an unexpected keyword argument 'phi'
Minimizing the one-loop potential, calculating the tunneling time or its upper bound, and printing the results took 0.461104 seconds.
Warning! "./MyResult.vout" was not produced!
Vevacious finished. Total running time was 0.482212 seconds.


Could you help me solve this error? Note that for the example SLHA file modifying the functions in Vevacious.py did solve the problem, but it did not for my slhafile (also in the appendix). This slhafile was generated with via the chain softsusy - feynhiggs - susyhit.

Kind regards,
Melissa

[FStaub]

Hi,

I don't think that the SusyHit output is compatible with your Vevacious file. You are using:

Code: Select all

  <model_file>
  /scratch/mbeekveld/computation_code_susy/Vevacious/Vevacious/models/pure_SLHA/pure_SLHA2_MSSM_JustNormalHiggsVevs.vin
  </model_file>
 

while the SusyHit output is in the SLHA 1 format. There are model files using SLHA1 conventions. However, I must say that I don't know if those were ever used in practice.

Cheers,
Florian

[mbeekveld]

Hi Florian,

I indeed thought this was the issue as well, but I forgot to mention that I tried to modify this to slha1 but that did not solve my problem. Same issue arises.

Cheers,
Melissa

[mbeekveld]

Hi Florian,

I believe it has to do with an upgrade of scipy. I had installed Vevacious before, then everything worked for me just fine (also on susyhit files using the slha1 model file). But then I tried to install Vevacious on another cluster. I had to modify in cosmotransistions the file finiteT.py and add import scipy.misc, otherwise scipy.misc.factorial was not recognized. Since I made this change, the problem appeared. Is this useful?

Cheers,
Melissa

[FStaub]

Hi,

unfortunately, I haven't installed Vevacious myself for some time on a new machine. So, I didn't run in that problem before.
What version of CT do you use?

Cheers,
Florian

[mbeekveld]

I am using CT 2.0.3, downloaded from https://github.com/clwainwright/CosmoTransitions.

Another change I made compared to the previous time I installed Vevacious is that I am now using iminuit instead of pyminuit. Do you know whether that could affected anything?

[smandal123]

Hi Florian,

I was trying to use Vevacious for the model B-L-SM model and I also got the same error " unexpected keyword argument 'phi' " ,

tunnelingCalculator.npoints being given 22
falseVacuumArray = [ 247.654147 2500. ]
trueVacuumArray = [ -7.71967874e-08 3.58859398e+03]
Traceback (most recent call last):
File "./Vevacious.py", line 117, in <module>
thermalNotQuantum = False )
File "/home/smandal/packages/VEVA/Vevacious/bin/VevaciousParameterDependent.py", line 1242, in CalculateAction
quickTunneling = False )
TypeError: fullTunneling() got an unexpected keyword argument 'phi'

Minimizing the one-loop potential, calculating the tunneling time or its upper bound, and printing the results took 3.91793 seconds.

Warning! "/home/smandal/packages/VEVA/Vevacious/MyResult.vout" was not produced!

Vevacious finished. Total running time was 4.06078 seconds.

I dont have any clue what is going wrong here. Looking forward for your help.

Cheers,

Sanjoy

[FStaub]

I am using CT 2.0.3, downloaded from https://github.com/clwainwright/CosmoTransitions.

Another change I made compared to the previous time I installed Vevacious is that I am now using iminuit instead of pyminuit. Do you know whether that could affected anything?

Hi,

Since the error appears already before CT is called, when you use the loop potential, it might be a problem with the minimizer.
So, I worry that there are two sources of problems at the end. Maybe, it's worth to try to contact the CT author regarding the problem with scipy. He has left physics but obviously still works on the code from time to time. Since I can't reproduce the errors on my machine, it's hard to give any good advice.

Cheers
Florian

[FStaub]

Hi Florian,

I was trying to use Vevacious for the model B-L-SM model and I also got the same error " unexpected keyword argument 'phi' " ,

tunnelingCalculator.npoints being given 22
falseVacuumArray = [ 247.654147 2500. ]
trueVacuumArray = [ -7.71967874e-08 3.58859398e+03]
Traceback (most recent call last):
File "./Vevacious.py", line 117, in <module>
thermalNotQuantum = False )
File "/home/smandal/packages/VEVA/Vevacious/bin/VevaciousParameterDependent.py", line 1242, in CalculateAction
quickTunneling = False )
TypeError: fullTunneling() got an unexpected keyword argument 'phi'

Minimizing the one-loop potential, calculating the tunneling time or its upper bound, and printing the results took 3.91793 seconds.

Warning! "/home/smandal/packages/VEVA/Vevacious/MyResult.vout" was not produced!

Vevacious finished. Total running time was 4.06078 seconds.

I dont have any clue what is going wrong here. Looking forward for your help.

Cheers,

Sanjoy

This seems to be another problem. Could you please open a new thread for that. I'll delete this post here the next days.

Florian

[mbeekveld]

Hi Florian.

Thanks, I will contact him. Another question I have is regarding the .vin model files. It extracts HMIX[2] from the slha file provided, but in the example SLHA file CMSSM_CCB.slha.out, there is no HMIX[2] entry. How is this dealt with?

Cheers,
Melissa