영상법

영상법(映像法, method of images)은 라플라스 방정식의 경계값 문제를 좀 더 쉬운 다른 문제로 바꾸어 푸는 방법이다. 영상법의 타당성은 유일성의 정리에 의해 증명된다. 이 때, 바뀐 문제에서는 원래 문제의 전하에 대응하는 가상의 전하를 추가하므로, 이 가상의 전하를 마치 거울에 비치는 영상에 비유한 것이다. 무한한 도체 평면 위에 전하가 있다고 하자. 그렇다면 정전기 유도에 의하여 도체 표면에 전하가 유도된다. 이러한 계에서 공간 모든 위치의 전위를 계산하는 정전기학 문제를 생각해 보자. 도체 표면 위의 전위는 일정하므로, 이를 편의상 0으로 놓자. 이러한 문제는 직접 계산하기 힘들지만, 그림과 같이 평면 반대쪽에 정확하게 같은 위치에 있지만 그 부호가 다른 전하를 가상으로 추가하고, 도체 평면을 없애자. 이 문제는 비교적 쉽게 풀 수 있다. 이 새로운 문제에서는 대칭에 따라 서로 부호가 다른 두 전하 정중앙에서는 전위가 정확히 0이다. 즉, 새로운 문제의 해는 원래 문제의 디리클레 경계 조건 (즉, 평면 위에 전위가 0인 것)을 만족한다. 라플라스 방정식의 경계 조건을 만족시키는 해는 유일하므로, 원래 문제의 해는 새로운 문제의 해와 (도체 위의 부분에서) 같다. 가장 간단한 예는 그림과 같이 하나의 점전하가 있는 경우지만, 임의의 전하 분포에 대해서도 사용할 수 있다. 구면에서도 영상법을 적용할 수 있다. 그림과 같이, 반지름 R {\displaystyle R} 의 도체 구면 안에 점전하 q {\displaystyle q} 가 구 한가운데에서 p {\displaystyle p} 만큼 떨어져 있는 곳에 있다고 하자. 이 경우도 마찬가지로 정전기 유도에 의하여 도체 표면에 디리클레 경계 조건을 적용하여야 한다. 평면의 경우와 마찬가지로, 도체 표면의 전위를 0으로 놓자. 이 문제는 다음과 같이 바꿀 수 있다. 그림과 같이, 구의 중심에서 R 2 / p {\displaystyle R^{2}/p} 만큼 떨어진 곳에 가상의 전하 − q R / p {\displaystyle -qR/p} 를 놓자. 그렇다면 이 새 문제의 해는 원래 문제의 디리클레 경계 조건을 만족한다는 사실을 계산으로 확인할 수 있다. 구면 밖에 전하가 위치해 있는 경우나 점전하 대신 임의의 전하 분포가 있는 경우도 마찬가지로 다룰 수 있다. 선형 유전체는 P = ϵ 0 χ e E {\displaystyle \mathbf {P} =\epsilon _{0}\chi _{e}\mathbf {E} } 를 만족하는 물질이다. 전기 감수율 χ e {\displaystyle \chi _{e}} 를 가지는 유전체 내부에 전하량 q {\displaystyle q} 를 가지는 입자를 위치시킨다고 가정해보자. 입자가 유전체 내부에 배치되는 순간, 구속 전하가 입자가 만드는 전기장에 의해 재배열되고, 구속 전하와 원래 입자에 의한 새로운 전기장이 유도된다. 이때, 유도되는 구속 전하 밀도는 ρ b = − ∇ ⋅ P = − ∇ ⋅ ( ϵ 0 χ e ϵ D ) = − ( χ e 1 + χ e ) ρ f {\displaystyle \rho _{b}=-\nabla \cdot \mathbf {P} =-\nabla \cdot \left(\epsilon _{0}{\frac {\chi _{e}}{\epsilon }}\mathbf {D} \right)=-\left({\frac {\chi _{e}}{1+\chi _{e}}}\right)\rho _{f}} 와 같이 주어진다. 입자에 의한 q f = q {\displaystyle q_{f}=q} 와 유도되는 q b = − χ e 1 + χ e q {\displaystyle q_{b}=-{\frac {\chi _{e}}{1+\chi _{e}}}q} 를 합하면 전체 전하량 q t o t a l = q f + q b = 1 1 + χ e q {\displaystyle q_{total}=q_{f}+q_{b}={\frac {1}{1+\chi _{e}}}q} 을 얻을 수 있다. 또한 상대유전율은 전기 감수율과 ϵ r = 1 + χ e {\displaystyle \epsilon _{r}=1+\chi _{e}} 의 관계를 가지기 때문에 정리하면 q t o t a l = 1 ϵ r q {\displaystyle q_{total}={\frac {1}{\epsilon _{r}}}q} 을 얻을 수 있다. 진공에 놓인 도체 평판과 동일하게, z > 0 {\displaystyle z>0} 에서 유전율 ϵ 1 {\displaystyle \epsilon _{1}} , z < 0 {\displaystyle z<0} 에서 유전율 ϵ 2 {\displaystyle \epsilon _{2}} 을 가지는 유전체가 z = 0 {\displaystyle z=0} 에서 직선 경계를 가진다고 생각하자. 이때 경계면과 d {\displaystyle d} 만큼 떨어진곳에 전하량 q {\displaystyle q} 를 가지를 입자를 위치했을 때 생기는 전위는 영상법을 이용해 계산 가능하다. z > 0 {\displaystyle z>0} 에서의 전위를 z < 0 {\displaystyle z<0} 에 영상 전하 q ′ {\displaystyle q'} 를, z < 0 {\displaystyle z<0} 에서의 전위를 z > 0 {\displaystyle z>0} 에 영상 전하 q ″ {\displaystyle q''} 를 위치시켜 계산하면 V u p = 1 4 π ϵ 1 [ q x 2 + y 2 + ( z − d ) 2 + q ′ x 2 + y 2 + ( z + d ) 2 ] {\displaystyle V_{up}={\frac {1}{4\pi \epsilon _{1}}}\left[{\frac {q}{\sqrt {x^{2}+y^{2}+(z-d)^{2}}}}+{\frac {q'}{\sqrt {x^{2}+y^{2}+(z+d)^{2}}}}\right]} V d o w n = 1 4 π ϵ 2 [ q ″ x 2 + y 2 + ( z − d ) 2 ] {\displaystyle V_{down}={\frac {1}{4\pi \epsilon _{2}}}\left[{\frac {q''}{\sqrt {x^{2}+y^{2}+(z-d)^{2}}}}\right]} 과 같이 식을 얻을 수 있다. 전위는 연속적이고, 경계면에서 자유 전하가 존재하지 않음으로 z = 0 {\displaystyle z=0} 에서 경계조건은 V u p = V d o w n {\displaystyle V_{up}=V_{down}} , D u p ⊥ = D d o w n ⊥ {\displaystyle \mathbf {D} _{up}^{\perp }=\mathbf {D} _{down}^{\perp }} 으로 주어진다. 변위장은 D ⊥ = ϵ E ⊥ = − ϵ ∂ V ∂ n {\displaystyle \mathbf {D} ^{\perp }=\epsilon \mathbf {E} ^{\perp }=-\epsilon {\frac {\partial V}{\partial n}}} 을 만족하므로, 두번째 조건을 ϵ 1 ∂ V u p ∂ z = ϵ 2 ∂ V d o w n ∂ z {\displaystyle \epsilon _{1}{\frac {\partial V_{up}}{\partial z}}=\epsilon _{2}{\frac {\partial V_{down}}{\partial z}}} 와 같이 변위에 대해 바꿀 수 있다. ρ = x 2 + y 2 {\displaystyle \rho ={\sqrt {x^{2}+y^{2}}}} 으로 두면, 경계조건은 q ′ {\displaystyle q'} , q ″ {\displaystyle q''} 에 대해 다음 두 식을 만족한다. 1 4 π ϵ 1 [ q d ρ 2 + d 2 + q ′ d ρ 2 + d 2 ] = 1 4 π ϵ 2 [ q ″ d ρ 2 + d 2 ] {\displaystyle {\frac {1}{4\pi \epsilon _{1}}}\left[{\frac {qd}{\sqrt {\rho ^{2}+d^{2}}}}+{\frac {q'd}{\sqrt {\rho ^{2}+d^{2}}}}\right]={\frac {1}{4\pi \epsilon _{2}}}\left[{\frac {q''d}{\sqrt {\rho ^{2}+d^{2}}}}\right]} ϵ 1 4 π ϵ 1 [ q d ρ 2 + d 2 3 − q ′ d ρ 2 + d 2 3 ] = ϵ 2 4 π ϵ 2 [ q ″ d ρ 2 + d 2 3 ] {\displaystyle {\frac {\epsilon _{1}}{4\pi \epsilon _{1}}}\left[{\frac {qd}{{\sqrt {\rho ^{2}+d^{2}}}^{3}}}-{\frac {q'd}{{\sqrt {\rho ^{2}+d^{2}}}^{3}}}\right]={\frac {\epsilon _{2}}{4\pi \epsilon _{2}}}\left[{\frac {q''d}{{\sqrt {\rho ^{2}+d^{2}}}^{3}}}\right]} 이를 다시 q ′ {\displaystyle q'} , q ″ {\displaystyle q''} 에 대해 정리하면 ϵ 2 q + ϵ 2 q ′ = ϵ 1 q ″ {\displaystyle \epsilon _{2}q+\epsilon _{2}q'=\epsilon _{1}q''} , q − q ′ = q ″ {\displaystyle q-q'=q''} 이다. 따라서, 영상전하는 q ′ = ϵ r 1 − ϵ r 2 ϵ r 1 + ϵ r 2 q {\displaystyle q'={\frac {\epsilon _{r_{1}}-\epsilon _{r_{2}}}{\epsilon _{r_{1}}+\epsilon _{r_{2}}}}q} , q ″ = 2 ϵ r 2 ϵ r 1 + ϵ r 2 q {\displaystyle q''={\frac {2\epsilon _{r_{2}}}{\epsilon _{r_{1}}+\epsilon _{r_{2}}}}q} 과 같이 주어지며, 각 영역에서의 전위는 V u p = q 4 π ϵ 1 [ 1 x 2 + y 2 + ( z − d ) 2 + ϵ r 1 − ϵ r 2 ϵ r 1 + ϵ r 2 1 x 2 + y 2 + ( z + d ) 2 ] {\displaystyle V_{up}={\frac {q}{4\pi \epsilon _{1}}}\left[{\frac {1}{\sqrt {x^{2}+y^{2}+(z-d)^{2}}}}+{\frac {\epsilon _{r_{1}}-\epsilon _{r_{2}}}{\epsilon _{r_{1}}+\epsilon _{r_{2}}}}{\frac {1}{\sqrt {x^{2}+y^{2}+(z+d)^{2}}}}\right]} V d o w n = q 4 π ϵ 0 [ 2 ϵ r 1 + ϵ r 2 1 x 2 + y 2 + ( z − d ) 2 ] {\displaystyle V_{down}={\frac {q}{4\pi \epsilon _{0}}}\left[{\frac {2}{\epsilon _{r_{1}}+\epsilon _{r_{2}}}}{\frac {1}{\sqrt {x^{2}+y^{2}+(z-d)^{2}}}}\right]} 이다. 먼저 입자에 의한 전기장에 의해 유도되는 구속 전하를 계산해 선형 유전체 내에서의 상황을 진공에서의 상황으로 바꾸어 접근할 수 있다. 입자를 유전체 내부에 위치했을 때 전위에 관여하는 전하는 4개이다. 유전체 내부에 위치한 점전하 q {\displaystyle q} 점전하 주위로 유도되는 구속 전하 위쪽 경계면에 유도되는 구속 전하 아래쪽 경계면에 유도되는 구속 전하 1, 2에 의한 총전하량은 q t o t a l = 1 ϵ r 1 q {\displaystyle q_{total}={\frac {1}{\epsilon _{r_{1}}}}q} 으로 주어진다. z > 0 {\displaystyle z>0} 에서의 전위는 3, 4에 의한 영향을 z < 0 {\displaystyle z<0} 영역에 경계면으로부터 d {\displaystyle d} 만큼 위치한 영상 전하 q ′ {\displaystyle q'} 으로 대신해서 다음과 같이 얻을 수 있다. V u p = 1 4 π ϵ 0 [ 1 ϵ r 1 q x 2 + y 2 + ( z − d ) 2 + q ′ x 2 + y 2 + ( z + d ) 2 ] {\displaystyle V_{up}={\frac {1}{4\pi \epsilon _{0}}}\left[{\frac {1}{\epsilon _{r_{1}}}}{\frac {q}{\sqrt {x^{2}+y^{2}+(z-d)^{2}}}}+{\frac {q'}{\sqrt {x^{2}+y^{2}+(z+d)^{2}}}}\right]} 이때, 유전체에 의한 효과를 구속 전하로 치환해 전위를 구할 때 유전체의 효과를 나타내는 ϵ 1 {\displaystyle \epsilon _{1}} 대신 진공에서의 ϵ 0 {\displaystyle \epsilon _{0}} 를 사용한다. 마찬가지로, z < 0 {\displaystyle z<0} 에서의 전위는 다음과 같이 주어진다. V d o w n = 1 4 π ϵ 0 [ 1 ϵ r 1 q x 2 + y 2 + ( z − d ) 2 + q ″ x 2 + y 2 + ( z − d ) 2 ] {\displaystyle V_{down}={\frac {1}{4\pi \epsilon _{0}}}\left[{\frac {1}{\epsilon _{r_{1}}}}{\frac {q}{\sqrt {x^{2}+y^{2}+(z-d)^{2}}}}+{\frac {q''}{\sqrt {x^{2}+y^{2}+(z-d)^{2}}}}\right]} 이때, 대칭적 구조에 의해 영상 전하는 z > 0 {\displaystyle z>0} 과 반대의 지점에 위치하게 된다. ρ = x 2 + y 2 {\displaystyle \rho ={\sqrt {x^{2}+y^{2}}}} 으로 두면, 경계조건은 q ′ {\displaystyle q'} , q ″ {\displaystyle q''} 에 대해 다음 두 식을 만족한다. 1 4 π ϵ 0 [ 1 ϵ r 1 q d ρ 2 + d 2 + q ′ d ρ 2 + d 2 ] = 1 4 π ϵ 0 [ 1 ϵ r 1 q d ρ 2 + d 2 + q ″ d ρ 2 + d 2 ] {\displaystyle {\frac {1}{4\pi \epsilon _{0}}}\left[{\frac {1}{\epsilon _{r_{1}}}}{\frac {qd}{\sqrt {\rho ^{2}+d^{2}}}}+{\frac {q'd}{\sqrt {\rho ^{2}+d^{2}}}}\right]={\frac {1}{4\pi \epsilon _{0}}}\left[{\frac {1}{\epsilon _{r_{1}}}}{\frac {qd}{\sqrt {\rho ^{2}+d^{2}}}}+{\frac {q''d}{\sqrt {\rho ^{2}+d^{2}}}}\right]} ϵ 1 4 π ϵ 0 [ 1 ϵ r 1 q d ρ 2 + d 2 3 − q ′ d ρ 2 + d 2 3 ] = ϵ 2 4 π ϵ 0 [ 1 ϵ r 1 q d ρ 2 + d 2 3 + q ″ d ρ 2 + d 2 3 ] {\displaystyle {\frac {\epsilon _{1}}{4\pi \epsilon _{0}}}\left[{\frac {1}{\epsilon _{r_{1}}}}{\frac {qd}{{\sqrt {\rho ^{2}+d^{2}}}^{3}}}-{\frac {q'd}{{\sqrt {\rho ^{2}+d^{2}}}^{3}}}\right]={\frac {\epsilon _{2}}{4\pi \epsilon _{0}}}\left[{\frac {1}{\epsilon _{r_{1}}}}{\frac {qd}{{\sqrt {\rho ^{2}+d^{2}}}^{3}}}+{\frac {q''d}{{\sqrt {\rho ^{2}+d^{2}}}^{3}}}\right]} 이를 다시 q ′ {\displaystyle q'} , q ″ {\displaystyle q''} 에 대해 정리하면 q ′ = q ″ {\displaystyle q'=q''} , q − ϵ r 1 q ′ = ϵ r 2 ϵ r 1 q + ϵ r 2 q ″ {\displaystyle q-\epsilon _{r_{1}}q'={\frac {\epsilon _{r_{2}}}{\epsilon _{r_{1}}}}q+\epsilon _{r_{2}}q''} 이다. 따라서, 영상 전하는 q ′ = q ″ = ϵ r 1 − ϵ r 2 ϵ r 1 ( ϵ r 1 + ϵ r 2 ) q {\displaystyle q'=q''={\frac {\epsilon _{r_{1}}-\epsilon _{r_{2}}}{\epsilon _{r_{1}}(\epsilon _{r_{1}}+\epsilon _{r_{2}})}}q} 와 같이 주어지며, 각 영역에서의 전위는 V u p = q 4 π ϵ 1 [ 1 x 2 + y 2 + ( z − d ) 2 + ϵ r 1 − ϵ r 2 ϵ r 1 + ϵ r 2 1 x 2 + y 2 + ( z + d ) 2 ] {\displaystyle V_{up}={\frac {q}{4\pi \epsilon _{1}}}\left[{\frac {1}{\sqrt {x^{2}+y^{2}+(z-d)^{2}}}}+{\frac {\epsilon _{r_{1}}-\epsilon _{r_{2}}}{\epsilon _{r_{1}}+\epsilon _{r_{2}}}}{\frac {1}{\sqrt {x^{2}+y^{2}+(z+d)^{2}}}}\right]} V d o w n = q 4 π ϵ 1 [ 2 ϵ r 1 ϵ r 1 + ϵ r 2 1 x 2 + y 2 + ( z − d ) 2 ] {\displaystyle V_{down}={\frac {q}{4\pi \epsilon _{1}}}\left[{\frac {2\epsilon _{r_{1}}}{\epsilon _{r_{1}}+\epsilon _{r_{2}}}}{\frac {1}{\sqrt {x^{2}+y^{2}+(z-d)^{2}}}}\right]} 와 같다. ϵ 1 = ϵ r 1 ϵ 0 {\displaystyle \epsilon _{1}=\epsilon _{r_{1}}\epsilon _{0}} 이므로, 직접 영상법을 이용해 계산한 전위와 구속 전하를 먼저 계산한 이후 영상법을 사용한 전위가 동일한 것을 확인할 수 있다. 영상법을 통해 구한 구조에서 전하량 q {\displaystyle q} 의 입자가 받는 힘은 F = − 1 4 π ϵ 0 q 2 ( 2 d ) 2 z ^ {\displaystyle \mathbf {F} =-{\frac {1}{4\pi \epsilon _{0}}}{\frac {q^{2}}{(2d)^{2}}}\mathbf {\hat {z}} } 으로 주어진다. 2 d {\displaystyle 2d} 의 거리만큼 떨어진 두 실제 입자계의 에너지는 W 0 = ∫ ∞ 2 d F ⋅ d l = 1 4 π ϵ 0 ∫ ∞ 2 d q 2 r 2 d r = − 1 4 π ϵ 0 q 2 2 d {\displaystyle W_{0}=\int _{\infty }^{2d}\mathbf {F} \cdot d\mathbf {l} ={\frac {1}{4\pi \epsilon _{0}}}\int _{\infty }^{2d}{\frac {q^{2}}{r^{2}}}dr=-{\frac {1}{4\pi \epsilon _{0}}}{\frac {q^{2}}{2d}}} 와 같이 주어진다. 영상 전하-전하 쌍의 경우 에너지는 W = ∫ ∞ d F ⋅ d l = 1 4 π ϵ 0 ∫ ∞ d q 2 ( 2 r ) 2 d r = − 1 4 π ϵ 0 q 2 4 d {\displaystyle W=\int _{\infty }^{d}\mathbf {F} \cdot d\mathbf {l} ={\frac {1}{4\pi \epsilon _{0}}}\int _{\infty }^{d}{\frac {q^{2}}{(2r)^{2}}}dr=-{\frac {1}{4\pi \epsilon _{0}}}{\frac {q^{2}}{4d}}} 으로, W = 1 2 W 0 {\displaystyle W={\frac {1}{2}}W_{0}} 를 만족한다. 이는 W = ϵ 0 2 ∫ E 2 d τ {\displaystyle W={\frac {\epsilon _{0}}{2}}\int E^{2}d\tau } 을 생각했을때, 영상법의 경우 도체 내부에서 전기장이 없기 때문에 전하-전하쌍에 비해 에너지가 반으로 준다는 것을 확인할 수 있다. 각 경우에 구한 전위를 바탕으로 공간에 유도되는 전하를 계산했을때, 유도 전하량과 영상법을 사용할 때 임의로 잡은 영상 전하의 전하량이 일치하는 것을 확인할 수 있다. 즉, 영상 전하의 기원은 실제로 유도되는 전하의 분포에서 온다고 볼 수 있다. 영상법을 이용해 구한 전위는 V = 1 4 π ϵ 0 [ q x 2 + y 2 + ( z − d ) 2 − q x 2 + y 2 + ( z + d ) 2 ] {\displaystyle V={\frac {1}{4\pi \epsilon _{0}}}\left[{\frac {q}{\sqrt {x^{2}+y^{2}+(z-d)^{2}}}}-{\frac {q}{\sqrt {x^{2}+y^{2}+(z+d)^{2}}}}\right]} 으로 주어진다. 이때, 실제로 도체 표면에 유도되는 전하량을 구하기 위해 σ = − ϵ 0 ∂ V ∂ n = − ϵ 0 ∂ V ∂ z | z = 0 {\displaystyle \left.\sigma =-\epsilon _{0}{\frac {\partial V}{\partial n}}=-\epsilon _{0}{\frac {\partial V}{\partial z}}\right\vert _{z=0}} 를 사용하면 σ ( ρ ) = − q d 2 π ( ρ 2 + d 2 ) 3 / 2 {\displaystyle \sigma (\rho )={\frac {-qd}{2\pi (\rho ^{2}+d^{2})^{3/2}}}} 를 얻을 수 있다. 따라서 유도되는 총 전하량은 Q = ∫ σ d a = ∫ 0 2 π ∫ 0 ∞ − q d 2 π ( ρ 2 + d 2 ) 3 / 2 ρ d r d ϕ = − q {\displaystyle Q=\int \sigma da=\int _{0}^{2\pi }\int _{0}^{\infty }{\frac {-qd}{2\pi (\rho ^{2}+d^{2})^{3/2}}}\rho drd\phi =-q} 와 같다. 즉, 도체 표면에 실제로 유도되는 전하량은 도체 위에 위치한 전하량과 크기는 같고 부호는 반대이다. 영상법을 이용해 구한 전위는 V u p = q 4 π ϵ 1 [ 1 x 2 + y 2 + ( z − d ) 2 + ϵ r 1 − ϵ r 2 ϵ r 1 + ϵ r 2 1 x 2 + y 2 + ( z + d ) 2 ] {\displaystyle V_{up}={\frac {q}{4\pi \epsilon _{1}}}\left[{\frac {1}{\sqrt {x^{2}+y^{2}+(z-d)^{2}}}}+{\frac {\epsilon _{r_{1}}-\epsilon _{r_{2}}}{\epsilon _{r_{1}}+\epsilon _{r_{2}}}}{\frac {1}{\sqrt {x^{2}+y^{2}+(z+d)^{2}}}}\right]} , V d o w n = q 4 π ϵ 1 [ 2 ϵ r 1 ϵ r 1 + ϵ r 2 1 x 2 + y 2 + ( z − d ) 2 ] {\displaystyle V_{down}={\frac {q}{4\pi \epsilon _{1}}}\left[{\frac {2\epsilon _{r_{1}}}{\epsilon _{r_{1}}+\epsilon _{r_{2}}}}{\frac {1}{\sqrt {x^{2}+y^{2}+(z-d)^{2}}}}\right]} 으로 주어진다. 마찬가지로 각 표면에 유도되는 구속 전하는 σ b = − ϵ 0 [ ∂ V ∂ z | z = 0 + − ∂ V ∂ z | z = 0 − ] = q 4 π d ( ρ 2 + d 2 ) 3 / 2 [ 2 ( ϵ 2 − ϵ 1 ) ϵ 1 ( ϵ 1 + ϵ 2 ) ] {\displaystyle \sigma _{b}=-\epsilon _{0}\left[\left.{\frac {\partial V}{\partial z}}\right\vert _{z=0^{+}}-\left.{\frac {\partial V}{\partial z}}\right\vert _{z=0^{-}}\right]={\frac {q}{4\pi }}{\frac {d}{(\rho ^{2}+d^{2})^{3/2}}}\left[{\frac {2(\epsilon _{2}-\epsilon _{1})}{\epsilon _{1}(\epsilon _{1}+\epsilon _{2})}}\right]} 와 같다. 이때 표면에 유도되는 전하량은 q b s u r f a c e = ϵ 1 − ϵ 2 ϵ 1 ( ϵ 1 + ϵ 2 ) q {\displaystyle q_{b_{surface}}={\frac {\epsilon _{1}-\epsilon _{2}}{\epsilon _{1}(\epsilon _{1}+\epsilon _{2})}}q} 이다. 반사 원리 Tikhonov, A. N.; Samarskii, A. A. (1963). 《Equations of Mathematical Physics》. New York: Dover Publications. ISBN 0-486-66422-8.  더 이상 지원되지 않는 변수를 사용함 (도움말) Taddei, M M; Mendes, T N C; Farina, C. “Subtleties in energy calculations in the image method”. 《European Journal of Physics》 30 (5): 965–972. doi:10.1088/0143-0807/30/5/005.  Purcell, Edward Mills. 《Electricity and Magnetism》 (영어) 2판. McGraw-Hill.  CS1 관리 - 추가 문구 (링크) Griffiths, David J. (1999). 《Introduction to Electrodynamics》 (영어). Addison-Wesley. ISBN 978-0138053260.  Landau, L. D.; E. M. Lifshitz, L. P. Pitaevskii (1960). 《Electrodynamics of Continuous Media》 2판. London: Elsevier. ISBN 978-0-7506-2634-7.  더 이상 지원되지 않는 변수를 사용함 (도움말) CS1 관리 - 추가 문구 (링크) Jackson, J. D. (1998). 《Classical Electrodynamics》 (영어) 3판. New York: John Wiley & Sons. ISBN 978-0-471-30932-1. OCLC 535998. 2013년 8월 21일에 원본 문서에서 보존된 문서. 2012년 8월 25일에 확인함.  CS1 관리 - 추가 문구 (링크) Feynman, Richard; Robert Leighton, Matthew Sands (1989). 《The Feynman Lectures on Physics》. Addison-Wesley. ISBN 0-201-51003-0.  더 이상 지원되지 않는 변수를 사용함 (도움말) Jeans, James (1908). 〈8장〉. 《The Mathematical Theory of Electricity and Magnetism》. Cambridge University Press.  인하대학교 강의 노트 — 영상법