4
大学数学基礎解説
文献あり

Reissner-Nordström解

447
14
$$\newcommand{all}[1]{\left\langle#1\right\rangle} \newcommand{blr}[1]{\left[#1\right]} \newcommand{car}[1]{\left\{#1\right\}} \newcommand{di}[0]{\displaystyle} \newcommand{fr}[2]{\frac{#1}{#2}} \newcommand{lr}[1]{\left(#1\right)} \newcommand{ma}[1]{\(\di{#1}\)} $$

【修正履歴】
・16May2023: 公式5の「これらより」の下の式が間違っていたので修正しました
・20Sep2023: $g:=\det(g_{\mu\nu})$とすべきところを誤って$g:=\det(-g_{\mu\nu})$としていたので直しました
・14Oct.2023: 本記事ではMinkowski metricとして(+---)を用いています。一方、「準備」の章で引用しているsubmersionさんの記事では(-+++)を用いています。そのため、この章の作用$S$$R$の係数の符号とエネルギー・運動量テンソル$T_{\mu\nu}$の符号は(-+++)で正しいものになってます。(+---)の場合は この記事 をご参照ください。(submersionさんご指摘ありがとうございます)

はじめに

静的球対称かつ電荷が存在する場合のEinstein方程式の解であるReissner-Nordström解を求めます。

計算することはシンプルです。計量テンソルとして球対称静的な場合を仮定し、エネルギー・運動量テンソルとして電磁場のそれを採用します。これらをEinstein方程式および重力と結合したMaxwell方程式に代入し、物質の分布外部であることを仮定して解を求めます。

本記事はRef.[1]を参考に書いています。そのうち非可換ゲージ場やSkyrme模型を伴うブラックホール解のことを自身の勉強も兼ねて書きたいと思っており、その準備のための記事です。

規約

本記事で採用している規約は以下です。基本的にはRef.[2]の規約と同じです。


  1. 特に断らない限り、繰り返して存在する添字に関して和をとる、いわゆるEinsteinの規約を採用する
  2. $c$:光速、$G$:重力定数、$\epsilon_0$:真空の誘電率、$\mu_0$:真空の透磁率として$c=G=\epsilon_0=\mu_0=1$とする単位系を採用する
  3. $x_0\leftrightarrow t, \ x_1\leftrightarrow r, \ x_2\leftrightarrow \theta, \ x_3\leftrightarrow \phi$とする
  4. 以下「運動方程式」とは、場や計量テンソルによる作用の変分がゼロになる条件式のことを言う(いわゆるEuler-Lagrange eq.)。
  5. $A_{,x}$$A$$x$による微分を表す:$A_{,x}:=\partial_x A$
  6. $A_{;x}$$A$$x$による共変微分($\nabla_x A$)を表す。共変微分はテンソルの種類によってその作用が異なり、以下のように定義される
    • $(0,0)$テンソル: $S_{;\beta}:=S_{,\beta}$
    • $(1,0)$テンソル: $V^\alpha{}_{;\beta}:=V^\alpha{}_{,\beta}+\Gamma^\alpha_{\mu\beta}V^\mu$
    • $(0,1)$テンソル: $P_{\alpha;\beta}:=P_{\alpha,\beta} -\Gamma^\mu{}_{\alpha\beta}P_\mu$
    • $(2,0)$テンソル: $A^{\mu\nu}{}_{;\beta}:=A^{\mu\nu}{}_{,\beta}+A^{\alpha\nu}\Gamma^{\mu}{}_{\alpha\beta}+A^{\mu\alpha}\Gamma^\nu_{\alpha\beta}$
    • $(0,2)$テンソル: $T_{\mu\nu;\beta}:=T_{\mu\nu,\beta}-T_{\alpha\nu}\Gamma^{\alpha}{}_{\mu\beta}-T_{\mu\alpha}\Gamma^{\alpha}{}_{\nu\beta}$
    • $(1,1)$テンソル: $B^\mu{}_{\nu;\beta}:=B^\mu{}_{\nu,\beta}+B^\alpha{}_\nu\Gamma^\mu{}_{\alpha\beta}-B^\mu{}_\alpha\Gamma^\alpha{}_{\nu\beta}$

単位系として1.を用いていることに注意してください。そのため$c,G,\epsilon_0,\mu_0$は式に現れません。

準備

Einstein方程式およびMaxwell方程式を構成するのに必要な計量テンソル、Christoffel記号、Ricciテンソル等に関する基礎的なことに関して述べます。

計量テンソルの仮定

以下の球対称静的な計量テンソルを仮定します:
\begin{align} &ds^2=A(r)dt^2-B(r)dr^2-r^2d\theta^2-r^2\sin^2\theta d\phi^2,\\ \therefore & \ g_{00}=A(r), g_{11}=-B(r), g_{22}=-r^2, g_{33}=-r^2\sin^2\theta \end{align}
これはSchwarzshild解を求める際に使う仮定と同じです。$c=1$としていることに注意してください。

一般相対論・宇宙論では、計量テンソル全体の符号を上記の定義に対して逆にする定義(時間成分が負、空間成分が正)が多いですが、ここではそれとは逆にしました。素粒子物理ではこちらの定義が多いです。

電磁場のエネルギー・運動量テンソル

エネルギー・運動量テンソル(Energy-Momentum Tensor, EMT)の定義として便利なものに、重力を結合させた物質場の作用を計量テンソルで変分して定義するものがあります。MathlogのSubmersionさんの記事 「【相対論】一般相対性理論概説」 (Ref.[3])にあるように、Einstein-Hilbert 作用に物質場の部分${\cal L}_M$を加えた作用
\begin{align} S=\int\left[\frac{1}{2\kappa}R+{\cal L}_M\right]\sqrt{-g}d^4x \end{align}
より運動方程式を求め、これをEinstein eq.と比較することで、EMTが以下のように与えられることがわかります:
\begin{align} T_{\mu\nu} &=\frac{-2}{\sqrt{-g}}\frac{\delta (\sqrt{-g}{\cal L}_M)}{\delta g^{\mu\nu}} \end{align}
ここで$g:=\det(g_{\mu\nu})$です。

電磁場のラグランジアンは
\begin{align} {\cal L}=-\frac{1}{4}F_{\mu\nu}F^{\mu\nu}, \ \ F_{\mu\nu}:=\partial_\mu A_\nu-\partial_\nu A_\mu, \ \ A_\mu\text{はU(1)ゲージ場} \end{align}
であり、上記のEMTを計算すると以下のようになります。

電磁場に対するエネルギー・運動量テンソル

\begin{align} T_{\alpha\beta}=\frac{1}{4}g_{\alpha\beta}F_{\mu\nu}F^{\mu\nu} -g_{\beta\nu}F_{\alpha\mu}F^{\nu\mu} \end{align}


このEMTの表式さえあれば計算には十分なのですが、以下に少し重要な蛇足を付け足しておきます。

計量テンソルの変分による EMTとは別なEMTの定義に、正準エネルギー・運動量テンソル(canonical EMT)があります。これは以下のように与えられます。

正準エネルギー・運動量テンソル

\begin{align} T^\mu{}_\nu= \frac{\delta {\cal L}}{\delta(\partial_\mu A_\rho)}\partial_\nu A_\rho-\delta^\mu{}_\nu {\cal L} \end{align}

canonical EMTは並進変換に対するネーター・カレントです。

さきほど定義したmetricの変分によるEMTは、canonical EMTに比べて良い性質を持ちます:

計量テンソルの変分によるEMTの利点
  1. 計算が簡単
  2. Lorentzの足の入れ替えに対し対称

1.に関しては、canonical EMTは場の変分で定義されているのに対し、計量テンソルの変分で定義されるEMTではそれがないことに由来します。U(1)ゲージ場の場合は計算量はさほど変わりませんが、例えばSkyrme模型と呼ばれる核子の模型におけるEMTの計算の場合、計量を用いて定義したEMTのほうが遥かに計算が簡単です。この例はちょっと特殊かもしれませんが...。

2.は計量テンソルの変分によるEMTが${\cal L}_M$$g^{\mu\nu}$変分で定義されていることからわかります。一方でcanonical EMTを計算すると
\begin{align} T^{\mu\nu}_C=-F^{\mu\rho}\partial^{\nu}A_\rho +\frac{1}{4}F_{\alpha\beta}F^{\alpha\beta}\delta^\mu{}_\nu \end{align}
であり(下についた$C$はcanonicalを表す)、これは$\mu\leftrightarrow \nu$に対して対称ではありません。カレントには全微分項を足す不定性があるので、これに全微分項$\partial_\rho(F^{\mu\rho}A_\nu)$を足すと、計量テンソルの変分によるEMTと同じ結果を得ます。この変形には運動方程式$\partial_\rho F^{\mu\rho}=0$を用いていることに注意してください。更に上記の$T^{\mu\nu}_C$はゲージ不変でもないです。このように、canonical EMTはネーターカレントという"正当性"を持つにも関わらず、あまり嬉しくない性質を持ちます。

ちなみに、canonical EMTに全微分項を足して対称にする系統的な方法が存在し、Belinfante improvementなどと呼ばれます。この方法で求めたEMTはBelinfante tensorとかBelinfante–Rosenfeld stress–energy tensorと呼ばれます(Ref.[4][5])。

各種テンソルの定義

以下時空に関する各種テンソルの定義です。

  • Christoffel記号

    \begin{align} \Gamma^\gamma{}_{\beta\mu}:=\frac{1}{2}g^{\alpha\gamma}\left( g_{\alpha\beta,\mu}+g_{\alpha\mu,\beta}-g_{\beta\mu,\alpha} \right) \end{align}

  • Riemannテンソル

    \begin{align} R^\alpha{}_{\beta\mu\nu}:=\Gamma^\alpha{}_{\beta\nu,\mu}-\Gamma^\alpha{}_{\beta\mu,\nu}+\Gamma^\alpha{}_{\sigma\mu}\Gamma^\sigma{}_{\beta\nu}-\Gamma^\alpha{}_{\sigma\nu}\Gamma^\sigma{}_{\beta\mu} \end{align}

  • Ricciテンソル

    \begin{align} R_{\alpha\beta}:=R^\mu{}_{\alpha\mu\beta}\ \ (=R_{\beta\alpha}) \end{align}

  • Ricciスカラー

    \begin{align} R:=g^{\mu\nu}R_{\mu\nu}=g^{\mu\nu}g^{\alpha\beta}R_{\alpha\mu\beta\nu} \end{align}

Einstein方程式

Einstein eq.は以下。

Einstein eq.

\begin{align} G^{\alpha\beta}&+\Lambda g^{\alpha\beta}=8\pi T^{\alpha\beta}, \\G^{\alpha\beta}&:=R^{\alpha\beta}-\frac{1}{2}g^{\alpha\beta}R\ \ (=G^{\beta\alpha}) \\&(\text{unit}: c=G=1) \end{align}

$\Lambda$は宇宙定数ですが、以下ではこれがゼロの場合を考えます:
\begin{align} R_{\alpha\beta}-\frac{1}{2}g_{\alpha\beta}R&=8\pi T^{\alpha\beta} \end{align}

これを少し書き換えます。上の式の両辺に$g_{\alpha\beta}$をかけると以下を得ます:
\begin{align} R=-8\pi T \end{align}
ここで$T:=T^\alpha{}_\alpha$としました。これを使うとEinstein eq.は以下のように変形できます:
\begin{align} R_{\alpha\beta}=8\pi\left(T_{\alpha\beta}-\frac{1}{2}Tg_{\alpha\beta}\right) \end{align}
電磁場に対するEMTはトレースレス:$T^\mu_\mu=0$であるから、結局

$\Lambda=0$及び$T^\alpha{}_\alpha=0$の場合のEinstein eq.

\begin{align} R_{\alpha\beta}=8\pi T_{\alpha\beta} \end{align}

を得ます。以下ではこの形のEinstein eq.を使います。

重力と結合するMaxwell方程式

曲がった時空における電磁場を記述する方程式は以下です:

重力と結合するMaxwell方程式

\begin{align} \ F^{\mu\nu}{}_{;\nu}=0 \leftrightarrow \partial_\nu F^{\mu\nu}+\Gamma^\mu{}_{\alpha\nu}F^{\alpha\nu}+\Gamma^\nu{}_{\alpha\nu}F^{\mu\alpha}=0 \end{align}

対称性よりfield strengthの成分を制限する

field strengthを対称性により簡単にすることを考えます。ここでは静的な状況を考えているので、磁気モノポールを考えなければ電場のみが存在します。さらに球対称性より電場は$r$のみの関数かつ方向も$r$方向のみです。ゆえにfield strengthは以下の形になります:
\begin{align} F_{\mu\nu}= \begin{pmatrix} 0 & E_r(r) & 0 & 0\\ -E_r(r) & 0 & 0 & 0\\ 0 & 0 & 0 & 0\\ 0 & 0 & 0 & 0 \end{pmatrix} \end{align}

すなわち$F_{\mu\nu}$のゼロでない成分は$F_{01}=-F_{10}=E_r(r)$のみです。

各種テンソルの具体的な計算

Christoffel記号及びRicciテンソルの計算

Christoffel記号およびRicciテンソルを計算するのはまあまあ大変であり、簡便な計算法があるとよいです。Christoffel記号に関しては、Mathlogのyuskaさんの記事 「クリストッフェル記号の簡単な計算方法」 (Ref.[6])にある計算法 −変分によるクリストッフェル記号の計算− は間違いも少なく便利な方法だと思います。

ここではRef.[7]に記載されている、計量テンソルが対角的な場合に使えるChristoffel記号およびRicci テンソルの計算公式を紹介します。以下の公式・定理は、特に断らない限り、任意の次元で成立します。

計量テンソルが対角成分のみを持つ場合のChristoffel記号の性質
  • $\Gamma^\mu{}_{\nu\rho}$は少なくともどれか2つが一致しないとゼロになる。
  • $2\Gamma^\mu_{\sigma\mu}=\partial_\sigma\ln |g_{\mu\mu}|, \ 2\Gamma^\nu_{\mu\mu}=-g^{\nu\nu}\partial_\nu g_{\mu\mu}$
    ($\mu\neq\nu, \mu\neq \sigma$, 繰り返しの添字の和はとらない)
    これらより$$ \begin{aligned}4\Gamma^\nu{}_{\mu\mu}\Gamma^\mu{}_{\nu\nu} &=g^{\nu\nu}(\partial_\nu g_{\mu\mu})g^{\mu\mu}(\partial_\mu g_{\nu\nu})\\ &=(\partial_\mu\ln|g_{\nu\nu}|)(\partial_\nu\ln|g_{\mu\mu}|)\\ &=4\Gamma^\mu{}_{\mu\nu}\Gamma^\nu{}_{\mu\nu}\end{aligned} $$
    が成立する。

この事実を使うことで、以下のRicciテンソルの表式を得ます。

計量テンソルが対角成分のみの場合のRicciテンソルの計算

$\mu\neq\nu$および繰り返しの添字の和はとらないこととする。

  • 対角成分

    \begin{align} 4R_{\mu\mu}&=(\partial_\mu\ln|g_{\mu\mu}|-2\partial_\mu)\partial_\mu\ln \left| \frac{g}{g_{\mu\mu}} \right|\\ &\hspace{3cm}-\sum_{\sigma\neq\mu} \left[ (\partial_\mu\ln|g_{\sigma\sigma}|)^2 + \left( \partial_\sigma\ln\frac{|g|}{g_{\mu\mu}^2} +2\partial_\sigma \right) g^{\sigma\sigma}\partial_\sigma g_{\mu\mu} \right] \end{align}

  • 非対角成分

    \begin{align} 4R_{\mu\nu}&=(\partial_\mu\ln|g_{\nu\nu}|-\partial_\mu)\partial_\nu \ln \left| \frac{g}{g_{\mu\mu}g_{\nu\nu}} \right| +(\mu\leftrightarrow \nu)\\ &\hspace{5cm}-\sum_{\sigma\neq \mu,\nu} \partial_\mu\ln|g_{\sigma\sigma}| \partial_\nu\ln|g_{\sigma\sigma}|\end{align}

$(\mu\leftrightarrow \nu)$は、その手前の項の$\mu,\nu$を入れ替えたものを表します。これを使うと計量テンソルから直接Ricciテンソルを計算できます。

この公式より、計量テンソルが対角的な場合、以下が成立することがわかります:

Ricci tensorが対角的になる条件(など)
  1. Ricciテンソルが対角的になる十分条件は、計量テンソルが1つの座標にしか依存しないこと。例えばRobertson-Walker計量:$ds^2=dt^2-a(t)^2(dx^2+dy^2+dz^2)$の場合、Ricciテンソルは対角的
  2. 計量テンソルが2つの座標$x_\alpha,x_\beta$に依存する場合、Ricciテンソルの非対角成分でゼロでない可能性があるのは$R_{\alpha\beta}$のみ
  3. 2次元ではRicciテンソルは必ず対角的
  4. 公式4の「非対角成分の公式」には2階微分の項があるが、この項は$g_{\mu\nu}$が単一の変数の関数の積からなる場合($g_{\mu\mu}=\prod_{l}f^l_\mu(x_l)$と書ける場合)ゼロになる

今の場合計量テンソルは$r,\theta$に依存するので、2.より非対角で残る可能性があるのは$R_{12} \ (=R_{21})$のみです。計算すればこれがゼロになることがわかるので、Ricciテンソルは対角成分のみが残ります。

Christoffel記号およびRicciテンソルを具体的に計算すると以下を得ます:

Christoffel記号 & Riemannテンソルの具体的な表式
[Christoffel記号]
  • $\displaystyle \Gamma^0_{01}=\Gamma^0_{10}=\frac{A'}{2A}$
  • $\begin{aligned}\Gamma^1_{00}&=\frac{A'}{2B}\\ \displaystyle \Gamma^1_{11}&=\frac{B'}{2B}\\ \displaystyle \Gamma^1_{22}&=-\frac{r}{B}\\ \displaystyle \Gamma^1_{33}&=-\frac{r\sin^2\theta}{B}\end{aligned}$
  • $\begin{aligned} \Gamma^2_{12}&=\Gamma^2_{21}=\frac{1}{r}\\ \displaystyle \Gamma^2_{33}&=-\sin\theta\cos\theta\end{aligned}$
  • $\begin{aligned} \Gamma^3_{13}&=\Gamma^3_{31}=\frac{1}{r}\\ \displaystyle \Gamma^3_{23}&=\Gamma^3_{32}=\cot\theta\end{aligned}$

[Riemannテンソル]
  • $\displaystyle R_{00}=-\frac{A'}{4B} \left( \frac{A'}{A}+\frac{B'}{B} \right) +\frac{A''}{2B}+\frac{A'}{Br}$
  • $\displaystyle R_{11}=\frac{A'}{4A} \left( \frac{A'}{A}+\frac{B'}{B} \right) -\frac{A''}{2A}+\frac{B'}{Br}$
  • $\displaystyle R_{22}=-\frac{r}{2B} \left( \frac{A'}{A}-\frac{B'}{B} \right) -\frac{1}{B}+1$
  • $\displaystyle R_{33} =R_{22}\sin^2\theta$

ここで${}'$$r$による微分を意味する

公式5,6および定理1を使うとこれらを計算することができますが、公式6が有用か否かは場合と慣れと好みによると思います(今の場合Christoffel記号も計算しなければならないし...)。一方公式5と定理1は、計量テンソルが対角的な場合いつでも有用です。

エネルギー・運動量テンソル

電磁場のEMTを具体的に計算すると以下のようになります。

電磁場のエネルギー・運動量テンソル
  • $\displaystyle T_{00}=-\frac{1}{2}AF_{01}F^{01}$
  • $\displaystyle T_{11}=\frac{1}{2}BF_{01}F^{01}$
  • $\displaystyle T_{22}=-\frac{1}{2}r^2F_{01}F^{01}$
  • $\displaystyle T_{33}=T_{22}\sin^2\theta$

方程式を解く

以上の準備により、方程式群の具体的な表式を求めることができます。ここからは方程式群を解くことにより$A(r),B(r)$を決定することを目指します。

これらを決定するために、Einstein eq.と重力の結合したMaxwell eq.を使います。

\begin{align} &\text{Einstein eq.}: \ R_{\alpha\beta}=8\pi T_{\alpha\beta}\\ &\text{Maxwell eq.}: \ F^{\mu\nu}{}_{;\nu}=0 \leftrightarrow \partial_\nu F^{\mu\nu}+\Gamma^\mu{}_{\alpha\nu}F^{\alpha\nu}+\Gamma^\nu{}_{\alpha\nu}F^{\mu\alpha}=0 \end{align}

$AB=\text{const.}$を導く

前章のエネルギー・運動量テンソルをみると、次の式が成立することがわかります:
\begin{align} T_{00}/A+T_{11}/B=0 \end{align}
これにEinstein eq.を使うと
\begin{align} R_{00}/A+R_{11}/B=0\\ \leftrightarrow \frac{1}{rB} \left( \frac{A'}{A}+\frac{B'}{B} \right) =0\\ \leftrightarrow \partial_r \ln(AB)=0\\ \therefore AB=\text{const.} \end{align}
を得ます。この定数を$f$とすれば
\begin{align} AB=f \ (=\text{const.}) \end{align}
となります。

Maxwell方程式から電場を決定する

$F^{\mu\nu}{}_{;\nu}=0$$\mu=0$の場合を考えます。
\begin{align} \partial_\beta F^{0\beta}+\Gamma^0{}_{\mu\beta}F^{\mu\beta} +\Gamma^\beta{}_{\mu\beta}F^{0\mu}=0\\ \leftrightarrow \partial_rF^{01} +\left( \frac{A'}{2A}+\frac{B'}{2B} +\frac{2}{r} \right)F^{01}=0 \end{align}
いま$A'/A+B'/B=0$なので
\begin{align} \partial_rF^{01}+\frac{2}{r}F^{01}=0\\ \therefore F^{01}=\text{const.}/r^2\\ \therefore E_r=\text{const.}/r^2 \end{align}
これはCoulomb力の表式であるから定数は$Q/4\pi$$Q$:$r$内に存在する電荷)であるべきであり
\begin{align} F_{01}=E_r=\frac{Q}{4\pi r^2} \end{align}
を得ます。

Einstein eq.の$2\text{-}2$成分より$A,B$を導く

Einstein方程式の2-2成分
\begin{align} R_{22}=8\pi T_{22} \end{align}
より$A$を決定します。$R_{22}$を具体的に計算すると
\begin{align} R_{22}&=-\frac{4}{2B} \left( \frac{A'}{A}-\frac{B'}{B} \right) -\frac{1}{B}+1\\ &=-\frac{1}{f}\frac{\partial}{\partial r}(rA)+1 \end{align}
ここで$B=f/A$と、これより導かれる$B'=-fA'/A^2$を用いました。よって
\begin{align} \frac{\partial}{\partial r}(rA)=f-\frac{Q^2}{4\pi r^2}\\ \therefore A=f+\frac{C}{r}+\frac{Q^2}{4\pi r^2} \end{align}
を得ます。$C$は積分定数です。

ここで$r\rightarrow \infty$で計量テンソルがMinkowskiに近づくことから$f=1$でなければなりません。また、弱重力極限では$g_{00}\simeq 1-2M/r$となることから($M$$r$内の質量。submersionさんの記事(Ref.[3])の「万有引力の幾何学化」の章の$\Phi$$-M/r$にしたもの)
\begin{align} C=-2M \end{align}
です。

以上から$r_s:=2M,\ r_Q^2:=Q^2/4\pi$として
\begin{align} A&=1-\frac{r_s}{r}+\frac{r_Q^2}{r^2},\\ B&=\frac{1}{A}=\left(1-\frac{r_s}{r}+\frac{r_Q^2}{r^2}\right)^{-1} \end{align}
となります。


最終的に、 Reissner-Nordström解は以下のようになります:

Reissner-Nordström解

\begin{align} g_{\mu\nu}&= \begin{pmatrix} 1-\frac{r_s}{r}+\frac{r_Q^2}{r^2} & 0 & 0 & 0\\ 0 & -\left(1-\frac{r_s}{r}+\frac{r_Q^2}{r^2}\right)^{-1} &0&0\\ 0&0&-r^2&0\\ 0&0&0&-r^2\sin^2\theta \end{pmatrix} \\ {}\\ r_s&:=2M, \ \ r^2_Q:=Q^2/4\pi \end{align}

おしまい。${}_\blacksquare$

参考文献

[2]
Schutz, B. (江里口良二、二間瀬敏史 訳), 相対論入門 下 − 一般相対論 −, 丸善株式会社, 1988
[4]
Belinfante F.J., On the current and the density of the electric charge, the energy, the linear momentum and the angular momentum of arbitrary fields, Physica. 7 (5), 1940, 449
[7]
Win, K. Z., Ricci Tensor of Diagonal Metric, arXiv:gr-qc/9602015, 1996
投稿日:2023514
OptHub AI Competition

この記事を高評価した人

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

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

バッジはありません。

投稿者

bisaitama
bisaitama
137
56038

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中