Differences

This shows you the differences between two versions of the page.

Link to this comparison view

documentation:language_reference:objects:spectra:start [2016/09/25 23:02] – created Maurits W. Haverkortdocumentation:language_reference:objects:spectra:start [2016/10/10 09:41] (current) – external edit 127.0.0.1
Line 1: Line 1:
 +====== Spectra ======
  
 +###
 +The Spectra object stores complex spectra as created by the functions CreateSpectra or CreateResonantSpectra. The spectra are stored as complex functions on a discrete mesh. The object furthermore contains several additional information such as the integral.
 +###
 +
 +===== Shared code for examples =====
 +
 +###
 +In order to run the examples for the different functions, properties and methods of the object spectra we define a set of functions specific to spectra. These functions must be loaded before any of the examples provided later on will work. The code creates three spectra stored in the variable //G//
 +  - a spectrum with 9 peaks at energy 1 to 9 in steps of 1 and intensity 0.1<sup>2</sup> to 0.9<sup>2</sup>.
 +  - a spectrum with 4 peaks at energy 2 to 8 in steps of 2 and intensity 0.2<sup>2</sup> to 0.8<sup>2</sup>.
 +  - a spectrum with 5 peaks at energy 1 to 9 in steps of 2 and intensity 0.1<sup>2</sup> to 0.9<sup>2</sup>.
 +
 +The actual code in Quanty is:
 +###
 +
 +<code Quanty definitions.Quanty>
 +NF = 9
 +NB = 0
 +T1 = NewOperator("Create", NF, {0,1,2,3,4,5,6,7,8},{0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9})
 +T2 = NewOperator("Create", NF, {0,1,2,3,4,5,6,7,8},{0.0,0.2,0.0,0.4,0.0,0.6,0.0,0.8,0.0})
 +T3 = NewOperator("Create", NF, {0,1,2,3,4,5,6,7,8},{0.1,0.0,0.3,0.0,0.5,0.0,0.7,0.0,0.9})
 +H = NewOperator("Number", NF, {0,1,2,3,4,5,6,7,8}, {0,1,2,3,4,5,6,7,8}, {1,2,3,4,5,6,7,8,9})
 +psi = NewWavefunction(NF,NB,{{"000000000",1}})
 +G = CreateSpectra(H, {T1,T2,T3}, psi, {{"Emin",-1}, {"Emax",11}, {"Gamma",0.1}})
 +gnuplotHead = [[
 +set autoscale
 +set xtic auto
 +set ytic auto
 +set style line  1 lt 1 lw 1 lc rgb "#000000"
 +set style line  2 lt 1 lw 1 lc rgb "#FF0000"
 +set style line  3 lt 1 lw 1 lc rgb "#00FF00"
 +set style line  4 lt 1 lw 1 lc rgb "#0000FF"
 +set style line  5 lt 1 lw 3 lc rgb "#808080"
 +set style line  6 lt 1 lw 3 lc rgb "#00FFFF"
 +set style line  7 lt 1 lw 3 lc rgb "#FF00FF"
 +set style line  8 lt 1 lw 3 lc rgb "#FFFF00"
 +set xlabel "Energy (arb. units)" font "Times-Roman,12"
 +set ylabel "Intensity (arb. units)" font "Times-Roman,12"
 +set size 1.0, 0.3
 +set terminal postscript portrait enhanced color  "Times-Roman" 8
 +]]
 +</code>
 +
 +###
 +The code above does not create output. We can save the spectra //G// to file and use gnuplot to plot the spectra with the example below.
 +###
 +
 +<code Quanty Example.Quanty>
 +dofile("definitions.Quanty")
 +G.Print({{"file","Spectra.dat"}})
 +gnuplotScript = gnuplotHead .. [[
 +set output "Spectra.ps"
 +plot "Spectra.dat"  using 1:(-$3  ) title 'T1' with lines ls  1,\
 +     "Spectra.dat"  using 1:(-$5+4) title 'T2' with lines ls  2,\
 +     "Spectra.dat"  using 1:(-$7+8) title 'T3' with lines ls  3
 +]]
 +file = io.open("Spectra.gnuplot", "w")
 +file:write(gnuplotScript)
 +file:close()
 +-- call gnuplot to execute the script
 +os.execute("gnuplot Spectra.gnuplot")
 +os.execute("convert -density 1024 Spectra.ps -resize 1024 Spectra.png")
 +</code>
 +
 +###
 +The resulting picture is:
 +###
 +
 +{{:documentation:language_reference:objects:spectra:spectra.png?nolink |}}
 +
 +===== Table of contents =====
 +{{indexmenu>.#2}}

You are here: Quanty - a quantum many body script language » Documentation » Language Reference » Objects » Spectra

Print/export