3

複素平面上の格子点に零点をもつ関数

195
0

三角関数の無限積を拡張した以下のような関数について考えました。
f(z)=z(n,m)(0,0)(1zn+im)
追記:間違いを訂正しました

背景

複素関数は2次元の静電場を表すことを知り、いろいろな電荷分布における電場を考えました。位置aにある電荷Qaの作る電場は
w=(Qa/2π)ln(za)
という関数で表され、複数の電荷が作る電場は

w=(Qa/2π)ln(za)+(Qb/2π)ln(zb)+(Qc/2π)ln(zc)+ ..=(1/2π)ln((za)Qa(zb)Qb(zc)Qc ..)
と表せます。ここでz=x+iyは空間座標を表し、w=ϕ+iν は電位と等電位線にそった位置を表します。これから電気力線や等電位線を求めるためには、これの逆関数を求め、z=f1(w)の形にする必要があります。

直線に並んだ電荷の列が作る電場

対数の中をg(z)とすると2πw=lng(z) という形になります。電荷を1とし実軸状に間隔1で並べたものを考えると
g(z)=z(z1)(z+1)(z2)(z+2)(z3)(z+3) ..
になります。ここでlng(z)+Cのように定数項を足すことは、電位の基準をずらすことと同じです。電位の基準はどこにとってもいいので、g(z)に好きな定数を掛けることができます。なので
g(z)=πz(1z)(1+z)(1z/2)(1+z/2)(1z/3)(1+z/3) ..=sin(πz)
を考えます。
2πw=lnsin(πz)
より
z=arcsin(e2πw)π
これを描画すると以下のようになります。
並んだ電荷が作る電場 並んだ電荷が作る電場
複素関数が多価のときもプログラムの関数は一つの主値を返すよう作られているので、別の分岐の関数を手動で求めて足さないといけません。arcsin(x)は実軸にπずらすだけなので簡単です。

正と負が交互に並んだ電荷の列が作る電場

g(z)=zz24z21z216z29z236z225 ..
定数倍はノーカンなので
g(z)=z1z2/41z21z2/161z2/91z2/361z2/25 ..=sin(πz/2)cos(πz/2)=tan(πz/2)
逆関数は
z=2arctan(e2πw)π
正負が交互に並んだ場合 正負が交互に並んだ場合

本題

ここまでくると、格子状にならんだ電荷が作る電場を知りたくなります。ここで冒頭の式がでてくるわけです。
f(z)=z(n,m)(0,0)(1zn+im)
まずnについての掛け算を先にします。下図の青の方向です。
格子 格子
まず以下の式
sinπ(ab)sinπa=(ab)(1(ab))(1+(ab))(1(ab)/2)(1(ab)/2) ..a(1a)(1+a)(1a/2)(1+a/2) ..(1)=(1ba)(1ba1)(1ba+1)(1ba2)(1ba+2) ..a=im,b=zを代入して
sinπ(imz)sinπim=(1zim)(1zim1)(1zim+1)(1zim2)(1zim+2) ..
を得ます。m=0 のとき
sinπz=πz(1z)(1+z)(1z/2)(1+z/2)(1z/3)(1+z/3) ..
であることをふまえmについての積をとると
f(x)=sinπzπm0sinπ(imz)sinπim
siniz=isinhzと加法定理を使うと
sinπ(imz)sinπim=cosπzicothπm sinπz
mにマイナスをかけると複素共役になるので
f(x)=sinπzπm=1(cos2πz+coth2πm sin2πz)=sinπzπm=1(1+sin2πzsinh2πm)
となります。今度はmについての掛け算を先にしてみます。(1)a=in,b=ziを代入し
sinπ(iniz)sinπni=(1zn)(1zni)(1zn+i)(1zn2i)(1zn+2i) ..
これから同様の計算をし
f(z)=sinhπzπn=1(cosh2πzcoth2πn sinh2πz)=sinhπzπn=1(1sinh2πzsinh2πn)

おかしい

この2つの式はかける順番を変えただけで同じはずですが、これらの関数には縦横の"向き"が違っています。そもそも縦横対称の格子なのに向きが生じる時点でおかしいです。実際
f1(z)=sinπzπm=1(1+sin2πzsinh2πm)
f2(z)=sinhπzπn=1(1sinh2πzsinh2πn)
とするとf1(z)Im(z)で発散して、Re(z) では発散しません。またf1(iz)=if2(z)となり向きが90度回転します。ある程度で打ち切った部分積をMathematicaの描画関数で描かせると以下のようになります。
!FORMULA[44][1264277066][0] sinπzπm=1(1+sin2πzsinh2πm)
"正解"を知るため9×9の格子についての積を描画すると以下のようになります。
!FORMULA[46][-1847066442][0]の格子による!FORMULA[47][1126158767][0] 9×9の格子によるf(z)

似たようなケースは例えば以下の級数でも起こります。
11/2+1/31/4+1/5 ..=ln2
これは以下のように足し方を変えると違う値になります。
(1+1/3)1/2+(1/5+1/7)1/4+(1/9+1/11)1/6+ ..=(3/2)ln2
今回もx方向とy方向で均等に掛けないために違いが生まれたようです。

均等な掛け方を考える

方法1

なるべく均等にかけるため下図のように4つをまとめて掛けたらどうなるでしょうか
4つづつまとめる 4つづつまとめる
(1±xn±im)
の形の項をすべてかけると
(1x2(n+im)2)(1x2(nim)2)
という形になります。ここで
cos2πbcos2πa1cos2πa=sin((a+b)π)sin(aπ)sin((ab)π)sin(aπ)=((1+ba)(1+ba1)(1+ba+1)(1+ba2)(1+ba+2) ..)((1ba)(1ba1)(1ba+1)(1ba2)(1ba+2) ..)=(1b2a2)(1b2(a1)2)(1b2(a+1)2)(1b2(a2)2)(1b2(a+2)2) ..
これにa=im,b=zをいれると
cos2πzcos2πim1cos2πim=(1+z2m2)(1z2(1im)2)(1z2(1+im)2)(1z2(2im)2)(1z2(2+im)2) ..
また
cos2πzcos2πim1cos2πim=cos2πzcosh2πm1cosh2πm=11cosπz1coshπm=1+sin2(πz)sinh2(πm)
m=0のときはsinπz/πをかければよく
f(x)=sinπzπm=1(1+sin2(πz)sinh2(πm))
前と同じ形になりました。はじめ計算ミスをしていて違うことをしているつもりでしたが、結局同じことをしていました。失礼しました。

方法2

均等に二重カウントする 均等に二重カウントする
上図の同じ色のラインをまとめてかけていくと、すべての点で二重カウントされることがわかります。この計算は結果的にさきほどのf1(z),f2(z)の積になるので
f(x)2=sinπzπsinhπzπm=1(1+sin2πzsinh2πm)(1sinh2πzsinh2πm)
よって
f(x)=sinπzπsinhπzπm=1(1+sin2πzsinh2πm)(1sinh2πzsinh2πm)
これを描画すると以下のようになりました。
!FORMULA[62][2000857207][0] sinπzπsinhπzπm=1(1+sin2πzsinh2πm)(1sinh2πzsinh2πm)
ルートを取ったので正負が逆になっているところがありますが、正しく分岐をとれば9×9の格子のものと十分に一致することがわかりました。しかしこの関数は周期的な模様はあるのですが、x方向、y方向ともに完全な周期性はありません。

まとめ

f(z)=z(n,m)(0,0)(1zn+im)
はかける順番により以下のような違う関数に近づく

  • はじめにnについてかける場合
    f1(z)=sinπzπm=1(1+sin2πzsinh2πm)
  • はじめにmについてかける場合
    f2(z)=sinhπzπn=1(1sinh2πzsinh2πn)
  • 均等にかける場合
    f3(z)=f1(x)f2(x)
    これらどの関数もxy両方向の完全な周期性はない

物理として解く

周期的な電位分布ϕ(x,y)は以下のようにフーリエ級数展開できます。
ϕ(x,y)=n,m=cn,me2πi(nx+my)
ここで
cn,m=1/21/2dx1/21/2dy ϕ(x,y)e2πi(nx+my)
ここでx:(1/2,1/2)y:(1/2,1/2)の領域を基本領域としました。これを用いてポアソンの方程式2ϕ=ρ(x)を解きます。ところで複素関数で電場を表す場合、ϕ+iν=f(x+iy)となりますが、
コーシー・リーマンの方程式ϕx=νy,ϕy,νxより2ϕ=0が成り立ちます。これは正則な点では成立しますが、対数分岐点で右辺が有限になります。実はガウスの定理とコーシーの積分定理は関係していて領域Sの境界をSとし
S2ϕ dS=Sf(z)dz
です。なので2ϕが有限な点に1/zの積分にあたる対数分岐点があるわけです。
同じことを言ってるのがポアソンの方程式で、電荷がある点だけ2ϕが有限になります。これを使うと
2ϕ(x,y)=ρ(x,y)=4π2n,m(n2+m2)cn,me2πi(nx+my)
となります。後に考えたことですが、fが全面で正則だと左辺が0になり、これを満たすのはc0,0以外0となるときなので、楕円関数のリウヴィルの第一定理と同じことだと思います。
この式にe2πi(Nx+My)をかけx:(1/2,1/2)y:(1/2,1/2)の領域で積分します。
1/21/2dx1/21/2dy ρ(x,y)e2iπ(Nx+My)=4π2(N2+M2)cN,M
ρ(x,y)=4π2δ(x,y)とします。すると
cn,m=1(n2+m2)
c0,0=となりますが、定数項なので一次元のときと同じように、電位の基準を調節すれば取り除けます。よって
ϕ(x,y)=(n,m)(0,0)e2πi(nx+my)n2+m2=n=1m=12n2+m2(cos(2π(nx+my))cos(2π(nxmy)))+n=12n2cos(2nπx)+n=12n2cos(2nπy)
部分和をプロットしてみると以下のようになります。
!FORMULA[96][-916878351][0] ϕ(x,y)
これに対応する複素関数はあるのかということなのですが、e2πi(nx+my)nx+myの部分がdxdyの正方形を一般には歪めるので正則な操作になりません。ただ和をとるときに(n,m)の組をうまくとればいけそうな気もします。

楕円関数

いろいろ関連分野について調べたところ、楕円関数がこれに関係あることがわかりました。複素平面で格子の周期性をもつ関数を楕円関数とよぶそうです。(これのことだったのか!!!)
今回得た2つの関数はテータ関数θ1 とヴァイエルシュトラスのσ函数で表せることがわかりました。
まずテータ関数θ1(v,τ)
θ1(v,τ)=1in=(1)nq(n+1/2)2zn+1/2  (q=eπiτ,z=e2πiv)
とします。すると
sinπzπk=1(1+(sinπzsinhπk)2)=eπ/3θ1(z,i)
となるようです。そして
g2=(Γ(1/4)44π)2,g3=0
としたときのヴァイエルシュトラスのσ函数で
sinπzπsinhπzπm=1(1+sin2πzsinh2πm)(1sinh2πzsinh2πm)=σ(z)
となります。そして
eπ/3θ1(z,i)=σ(z)eπz2/2
だそうです。ここらへんは後に勉強しようと思います。

数学に"訂正"された!

格子電荷の電場に対応する正解の複素関数はヴァイエルシュトラスのσ函数のようです。これの逆関数は求まらなそうなので、Mathematicaの等高線を書かせる関数でlnσ(z)の実数部分の等高線を描画します。
!FORMULA[112][-1955731343][0]) Re(lnσ(z))
すきまに負電荷が挟まってる!!!
現実の結晶とかはこうなっているので、数学のほうが現実に詳しいってことでしょうか?驚きです。

ありがとうございました。

投稿日:202425
更新日:202425
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

17世紀の数学を学び始めました。 https://www.17centurymaths.com/ このサイト素晴らしい。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. 背景
  2. 本題
  3. 物理として解く
  4. 楕円関数
  5. 数学に"訂正"された!