孤立した渦糸周辺での秩序変数

超伝導体中の渦糸の周辺での超伝導秩序変数は、渦糸中心からの距離をrとして、以下のような式で近似されます。（Tinkhamの教科書参照）ここで、νは1に近い定数で、ξはコヒーレンス長です。

これを用いて渦糸周辺での秩序変数の振幅を表す図を作りました。上図は下のスクリプトを実行して作られたPostScriptファイルsingle-vortex.epsからコピーしたものです。

描画には3次元グラフの表面の色づけの項で説明したpm3dを利用しています。さらに、set pm3d hidden3dオプションを用いて、3次元表面にグリッドを描画しています。

# vortexの式
vortex(x,y)=tanh(nu*sqrt(x**2+y**2)/xi)
nu=1.0
xi=0.01


# samplesとisosamplesは必ず奇数にする
set samples 51
set isosamples 51

set xrange [-0.05:0.05]
set yrange [-0.05:0.05]

# 3D表面にgrid（grid不要の場合は hidden3d 1 を消す）
set pm3d depthorder interpolate 1,1 hidden3d 1 

# grid用のlinestyle
set style line 1  linetype 1 linecolor rgb "black" linewidth 0.25 


set view 65,338

unset colorbox
set border 0
unset xtics 
unset ytics 
unset ztics
unset key

# 色使いの指定
set palette defined ( 0 "black", 0.5 "dark-blue" , 1 "cyan")

splot vortex(x,y) with pm3d

pause -1


set out "single-vortex.eps"
set term post color enhanced
replot
set out
set term pop