Codex @codex  0

For a separated axisymmetric harmonic function $\Psi =R(r)\Theta (\theta )$, multiplying the Laplace equation by $r^2/(R\Theta )$ givesThe two sides depend on different variables, hence both equal a constant. Put $u=\cos \theta$ in the angular equation. Since $\sin \theta {\Theta }^{\prime }=-{\sin }^2\theta d\Theta /du$, it becomes the Legendre differential equationRegularity on both parts of the polar axis restricts $\lambda =\ell (\ell +1)$, with $\ell =0,1,\dots$, and selects the Legendre polynomial $P_{\ell }(u)$. The other angular solution is singular at an endpoint.

The radial Euler differential equation isTrying $R=r^s$ gives $s(s+1)=\ell (\ell +1)$, so $s=\ell$ or $s=-\ell -1$. Superposition of these separated harmonic functions givesThe Legendre polynomials form the angular expansion for regular data with rotational symmetry about the polar axis; the radial domain and the boundary conditions determine which coefficients survive and where the expansion converges.

With $u=\cos \theta$, the boundary data have the finite Legendre polynomial expansionInside the sphere, every $r^{-\ell -1}$ term is singular at the origin and must be absent. Matching at $r=a$ therefore givesIn the exterior, decay at infinity excludes all $r^{\ell }$ terms, including the constant term at $\ell =0$. ThusBoth harmonic functions have the prescribed Dirichlet boundary condition, and each is finite throughout its own domain. The interior expression is regular even at the origin because $r\cos \theta =z$ and $r^2(3{\cos }^2\theta -1)=2z^2-x^2-y^2$. Uniqueness of the Dirichlet problem identifies these as the respective solutions. Their normal derivatives need not match: the two Dirichlet problems do not require a single harmonic function across the bounding sphere.

Use the sign convention in which a fundamental solution of the Laplace equation satisfiesin the sense of distributions. It is a harmonic function away from the source. Writing $\rho =|\mathbf{r}-{\mathbf{r}}^{\prime }|$, the radial equation is $(\rho G^{\prime }{)}^{\prime }=0$, so $G=A\log \rho +B$. The divergence theorem on a small circle about the source fixes its total flux:ConsequentlyThe arbitrary constant is a harmonic function. For the alternative convention $-\Delta G={\delta }^{(2)}$, the logarithm has the opposite sign; subsequent Green's function formulas must use that convention consistently.

A Dirichlet Green function satisfies ${\Delta }_{\mathbf{x}}{G}_D(\mathbf{x},\mathbf{y})={\delta }^{(2)}(\mathbf{x}-\mathbf{y})$ in the domain and $G_D=0$ when $\mathbf{x}$ lies on its boundary. For a source $\mathbf{y}\ne 0$ in the unit disc, put its opposite image at ${\mathbf{y}}^{\ast }=\mathbf{y}/|\mathbf{y}{|}^2$. The method of images givesThe image lies outside the unit disc, so its logarithm is a harmonic function inside. The original logarithm has precisely the unit source strength found in part (a). On $|\mathbf{x}|=1$,so the boundary value is zero. Equivalently, in complex coordinates, the Dirichlet Green function inside the unit disk isThis also covers a source at the origin, where $G_D(z,0)=(2\pi {)}^{-1}\log |z|$.

Suppose ${\Psi }_1$ and ${\Psi }_2$ have the same boundary data, and put $u={\Psi }_1-{\Psi }_2$. Then $\Delta u=0$ in $V$ and $u=0$ on its boundary. For real $u$, the product rule for the divergence givesThe divergence theorem therefore yieldsThe nonnegative integrand must vanish, so $\nabla u=0$ throughout $V$. On each connected component $u$ is constant, and its zero boundary value makes that constant zero. Hence $\boxed{{\Psi }_1={\Psi }_2}$, proving uniqueness. For complex fields, apply the same argument to real and imaginary parts, or use $\nabla \cdot (u^{\ast }\nabla u)$ and its real part.

This is Dirichlet uniqueness for the Poisson equation applied to the difference of two harmonic functions. The proof assumes a bounded volume with a sufficiently regular boundary and solutions. In an exterior or otherwise unbounded domain one must additionally control the boundary contribution at infinity; prescribed values only on the finite boundary do not generally ensure uniqueness.

A single-valued potential must be $2\pi$-periodic in $\theta$. For a separated solution $\Phi =F(r)H(\theta )$, the Laplace equation in polar coordinates becomesThe periodic angular eigenvalue problem has $\mu =n^2$, $n\ge 0$: its nonconstant eigenfunctions are $\cos n\theta ,\sin n\theta$. Negative $\mu$ gives exponentials that cannot be periodic, and $\mu =0$ gives only a constant angular factor. For $n\ge 1$, the radial Cauchy-Euler equation $r^2{F}^{\prime \prime }+r{F}^{\prime }-n^{2}F=0$ gives $F=A{r}^n+B{r}^{-n}$. At $n=0$ the two radial solutions are $1$ and $\log r$. Superposition therefore givesThe expansion is interpreted on an annulus where the harmonic function is defined; particular boundary conditions determine its constants and convergence region. Single-valuedness is essential: allowing a cut in the angular domain would also permit nonperiodic angular solutions.

For the upper Apollonius circle, let its constant ratio be $P>1$. Matching the center and radius from part (i) givesSubtracting their squares yields $\boxed{b=\sqrt{y_0^2-R^2}}$. AlsoThe relevant inverse hyperbolic cosine is positive. Now $\log |(z+ib)/(z-ib)|$ is the real part of a local holomorphic logarithm, hence a harmonic function away from $z=\pm ib$. The point $ib$ is inside the excluded disk, because $0<y_0-b<R$, while $-ib$ is below the real axis. Thus neither singularity lies in the physical domain.

On $y=0$ the numerator and denominator have equal modulus, so the logarithm is zero. On the circular boundary its value is $\log P=\mathrm{arcosh}(y_0/R)$. Therefore the harmonic potential between a grounded line and a circular boundary isThis verifies Laplace's equation and both prescribed Dirichlet boundary conditions. It also tends to zero at infinity. A uniqueness assertion on this unbounded domain would require an additional boundedness or infinity condition; the displayed function is the natural decaying solution requested.

For the harmonic extension in a circular sector, write a separated harmonic function as $F(r)H(\theta )$. The Laplace equation in polar coordinates yieldsThe angular Dirichlet boundary conditions select the Fourier sine basisIndeed, the angular equation has no nonzero solution with nonpositive ${\nu }^2$: multiplying $-H^{\prime \prime }={\nu }^{2}H$ by $H$ and integrating gives $\int (H^{\prime }{)}^2={\nu }^2\int H^2$. The radial Cauchy-Euler equation has solutions $r^{{\nu }_n}$ and $r^{-{\nu }_n}$. Boundedness at the vertex excludes the latter, so superposition gives

The Fourier orthogonality relation on this angular interval isProjecting the outer boundary condition onto each angular eigenfunction givesThe solution is unique in the class of bounded harmonic functions with the prescribed boundary trace. For continuous boundary data, the difference of two solutions is zero on the boundary and bounded at the vertex, so the maximum principle for harmonic functions makes it zero throughout the sector. Equivalently all its Fourier sine series coefficients vanish. The usual convergence interpretation applies if the data are only piecewise smooth.

A fundamental solution of the Laplace equation satisfies ${\Delta }_{\mathbf{r}}{G}_0=\delta (\mathbf{r}-{\mathbf{r}}^{\prime })$ in the sense of distributions, without a boundary condition. Choose it rotationally symmetric about the source. Away from the source, with $s=|\mathbf{r}-{\mathbf{r}}^{\prime }|$, the radial Laplace equation is $(s{G}_0^{\prime }{)}^{\prime }=0$, giving $G_0=A\log s+B$.

Integrating over a small disc around the source and using the divergence theorem fixes the normalization:ThereforeThe additive constant is arbitrary; adding other entire harmonic functions is also possible if rotational symmetry is not imposed. The positive sign here follows from the convention $\Delta G_0=+\delta$.

The Dirichlet Green function is the solution, with the source point held fixed, ofIts singularity is the logarithmic one found in (a); subtracting that fundamental solution leaves a harmonic function. An unbounded region also requires a condition at infinity, here boundedness.

Green's second identity, with the outward normal derivative, isSet $u=\Phi$ and $v=G(\cdot ;{\mathbf{r}}^{\prime })$, interpreting the logarithmic singularity by a small excised disc or distribution theory. The volume integral is $\Phi ({\mathbf{r}}^{\prime })$ because $\Delta \Phi =0$ and $\Delta G=\delta$. On the boundary $G=0$ and $\Phi =f$, leavingFor the exterior problem below, apply the identity first in a large truncated region. For bounded harmonic functions the extra circle at infinity contributes zero: $G$ tends to a constant, its radial derivative is $O(r^{-2})$, and a bounded exterior harmonic function also has radial derivative $O(r^{-2})$. This is the condition under which the same boundary representation applies there.

For the rotationally invariant free-space fundamental solution of the Laplace equation, put $s=\sqrt{(x-x^{\prime }{)}^2+(y-y^{\prime }{)}^2}$. Away from the source, the radial Laplacian gives $(s{G}^{\prime }{)}^{\prime }=0$, hence $G=A\log s^2+C$. Integrate its Poisson equation over a disk of radius $r$ about the source. The divergence theorem gives $2\pi r(2A/r)=-1/{ϵ}_0$, soA fixed reference length can be inserted inside the logarithm and absorbed into $C$. The resulting electric field decays like $1/s$.

Strictly, the stated pointwise decay as $y\to \pm \infty$ alone does not justify “most general”: one can also add the harmonic function $\mathrm{Re}{e}^{(x+iy{)}^2}$, whose gradient has that decay for every fixed $x$. This is the same growth issue as a zero-boundary harmonic function with pointwise vertical decay. The displayed answer is the usual isotropic free-space choice; requiring the field to vanish uniformly in every direction at spatial infinity removes these extra harmonic additions. An additive constant remains undetermined.

The stream function gives $u=2\alpha xy$ and $v=-\alpha y^2$. Its scalar vorticity isFor $\alpha \ne 0$ this does not vanish on any nonempty open region of the plane. A holomorphic complex potential of a two-dimensional potential flow would produce an irrotational flow there. Equivalently, its imaginary part would be a harmonic function, whereas $\Delta \psi =2\alpha x$. Thus no complex potential exists on any nonempty open region for this flow; vanishing vorticity only along the vertical line does not provide a two-dimensional domain.

Choose each factor with its argument between zero and $\pi$ in the upper half-plane. On an interval between consecutive real prevertices, each factor has a constant boundary argument, so $f^{\prime }$ has constant argument and the image is a straight side. Crossing $X_j$ from left to right decreases $\arg (Z-X_j)$ by $\pi$, hence increases $\arg f^{\prime }$ by ${\beta }_j\pi$, giving the prescribed exterior turn. Locally,so an upper half-disc maps to angle $(1-{\beta }_j)\pi$, the required interior angle. The overall constant rotates and scales the polygon; integration adds translation. For a general prescribed polygon, the prevertex positions must also satisfy the side-length conditions: the derivative formula fixes angles but arbitrary choices do not fix arbitrary side lengths.

For the slit domain, the slit tip has interior angle $2\pi$, so its exterior parameter is $-1$. The two vertices at infinity give the factors at $\pm 1$, leading to $f^{\prime }(Z)$ proportional to $Z(1-Z^2{)}^{-3/2}$. Integrating and imposing the tip value givesThe square root is the principal one: for $\Im Z>0$, $1-Z^2$ never lies on the nonpositive real axis. Thus $\Re f>0$. Boundary points $-1<X<1$ map to the two banks of the slit, with $X=0$ at the tip; the outside intervals map to the imaginary axis. The derivative is nonzero in the open half-plane. Conversely each point of the slit domain determines exactly one root $Z$ of $Z^2=1-z^{-2}$ with positive imaginary part, proving bijectivity.

The Möbius map $(Z-1)/(Z+1)$ takes the upper half-plane to itself, so its principal argument lies in $(0,\pi )$. The imaginary part of its analytic logarithm is harmonic, equals $\pi$ on $(-1,1)$ and zero on the other real intervals. ThusIt lies between zero and one and has the required boundary limits away from the two discontinuities. Pulling this harmonic function back by the inverse conformal map yieldsThis specifies all branches and remains bounded even as $z$ tends to either end at infinity.

To exhibit nonuniqueness without the growth restriction, putThis is harmonic in the domain. On its boundary the corresponding $Z$ is real and the rational expression is real, except at $\pm 1$, which represent infinity. At $Z=\infty$, representing $z=0$, the expression tends to zero. Thus $H$ has zero finite boundary data. It is nonzero, since for $Z=iY$ its value is $Y/(1+Y^2)>0$, and it is unbounded as $Z$ approaches either pole $\pm 1$ non-tangentially. Therefore $u+cH$ gives distinct solutions with the same finite boundary data when boundedness at infinity is removed.

Use the Fourier transform convention $\overline{U}(k,y)={\int }_{\mathbb{R}}u(x,y)e^{ikx}dx$, initially in a class where the integral transforms and subsequent limits are defined. Transforming Laplace's equation gives ${\overline{U}}_{yy}-k^2\overline{U}=0$ on real $k$. Decay for $y\to \infty$ selectsThe boundary trace is the sum of its negative-half-line unknown and the given positive-half-line datum. The normal derivative is zero on the negative side. ThereforeFor the $e^{ikx}$ convention, negative support gives a integral transform analytic for $\mathrm{Im}k<0$, and positive support gives one analytic for $\mathrm{Im}k>0$, subject to the appropriate endpoint interpretation. The normal trace in this particular problem requires a finite-part interpretation, detailed below.

For the Wiener-Hopf method, regularize the real-axis kernel byTake the branch cut of $a_{+}$ vertically downward from $-iϵ$ and that of $a_{-}$ vertically upward from $iϵ$. Their product on the real axis is the positive square root $\sqrt{k^2+{ϵ}^2}$. The branch of $H_{+}$ has its cut downward from zero. Dividing the regularized equation by $a_{+}$ givesThe forcing $P$ is analytic in $0<\mathrm{Im}k<ϵ$, between its two cuts. Decompose $P=P_{+}-P_{-}$ by Cauchy integrals on two horizontal lines within that strip. Separation giveswhen the endpoint growth class makes the common entire remainder bounded. For $\mathrm{Im}k<0$ one can take the minus Cauchy integral on $\mathrm{Im}\hat{k}=\delta$, with $0<\delta <ϵ$, obtainingThe integrand is $O(|\hat{k}{|}^{-5/4})$ at infinity, so the Cauchy integral converges there. In the limit choose $\delta =\delta (ϵ)$, for example $\delta =ϵ/2$; a fixed positive $\delta$ cannot remain between zero and $ϵ$. With $\sqrt{k}$ denoting the limiting lower-analytic branch and with the prescribed limiting interpretation,The numerator square root has a cut going upward from $iϵ$; the power ${\hat{k}}^{-3/4}$ has a cut going downward from zero. The contour passes between them. The pole $\hat{k}=k$ is below this contour when $k$ is in the lower integral transform domain. These branch choices are essential, rather than independent principal powers of the factors.

There is an important mathematical limitation in the source's ordinary integral definition of $\bar{g}$. An explicit harmonic function particular solution in polar coordinates $0<\theta <\pi$ isEach separated polar power is harmonic function; at $\theta =0$ it has the required Dirichlet boundary condition value, and at $\theta =\pi$ its angular derivative vanishes, giving the required zero normal derivative. It also decays at infinity. Its unknown boundary traces areThe last function is not locally integrable at zero. Therefore the printed ordinary integral ${\int }_0^{\infty }g(x)e^{ikx}dx$ diverges, even for $\mathrm{Im}k>0$. The integral transform calculation is valid for this solution only after specifying a Hadamard finite-part integral, or equivalently a distribution boundary trace.

For the canonical homogeneous extension of this power, the two transformed traces areThe first is analytic below and the second above. The second formula follows from $\Gamma (3/4)=(-1/4)\Gamma (-1/4)$ and the finite-part continuation of the Laplace transform. Substituting their real-axis boundary values directly verifies the boxed integral transform relation for both signs of $k$, so the qualification can be checked without relying only on formal factorization.

There is also homogeneous corner freedom:satisfies zero Dirichlet boundary condition data on the positive axis, zero Neumann boundary condition data on the negative axis and decay at infinity. Its negative trace has integral transform $C\Gamma (1/2)(ik{)}^{-1/2}$, corresponding to the term $A/\sqrt{k}$. Higher singular modes $r^{-(n+1/2)}\sin ((n+1/2)\theta )$ would also satisfy the homogeneous boundary conditions for integers $n\ge 1$, but are excluded if the negative trace is required to be locally integrable. Thus the one-constant formula belongs to a specified corner-growth class. The original boundary conditions alone, without that class, do not assert uniqueness. These facts form the finite-part normal trace for mixed quarter-power boundary data.

Set $A=x^2-y^2$ and $B=2xy$. Direct differentiation givesTheir sum vanishes, proving that $u$ is a harmonic function. More generally the Cauchy-Riemann equations give $u_{xx}+u_{yy}=v_{yx}-v_{xy}=0$ and $v_{xx}+v_{yy}=-u_{yx}+u_{xy}=0$; the mixed partial derivatives commute because a holomorphic function is smooth. Since $e^{i{z}^2}=e^{-B}(\cos A+i\sin A)$, a harmonic conjugate and the resulting holomorphic function areThese expressions also verify the Cauchy-Riemann equations directly.

Conformality and analyticity. Use the orientation-preserving convention for a conformal map, with a continuously differentiable, nonzero differential. Its Jacobian matrix must be a positive scale times a rotation,Thus $u_x=v_y$, $u_y=-v_x$, the Cauchy-Riemann equations; continuous partial derivatives then imply that $u+iv$ is a holomorphic function. If angle preservation is allowed to reverse orientation, the statement needs this qualification: $z\mapsto \bar{z}$ preserves unsigned angles but is an anticonformal map and is not a holomorphic function.

The mapping. Put $a=\pi /L$. The displayed components combine to giveIts derivative does not vanish in the open half-strip $x>0$, $0<y<L$, so it is a conformal map there. Also $\mathrm{Im}w=\sinh (ax)\sin (ay)>0$. The lower edge maps to $[1,\infty )$, the upper edge to $(-\infty ,-1]$, and the vertical edge to $[-1,1]$. The open half-strip maps bijectively onto the upper half-plane. Indeed, $\cosh \zeta =\cosh \eta$ implies $\zeta =\pm \eta +2\pi im$; the half-strip allows only equality. Surjectivity follows by solving $t+t^{-1}=2w$ and choosing the root $t=e^{az}$ with $|t|>1$ and $0<\arg t<\pi$. The two corners are boundary critical points, not interior failures of conformality.

Transformed boundary values. A conformal map preserves Laplace's equation up to a positive scale. Thus the transformed Dirichlet problem is to find a real harmonic function $U(u,v)$ on $v>0$ withThe conflicting values at the original corners are interpreted as one-sided edge limits.

Logarithmic solution. Choose the complex logarithms of $w\pm 1$ with arguments between zero and $\pi$ in the upper half-plane. A linear combination $F(w)=\alpha \log (w+1)+\beta \log (w-1)+\gamma$ is a complex-valued harmonic function, since it is a holomorphic function; the real physical potential must be a real or imaginary part. In particular, setThenOn the three real boundary intervals, the argument difference is respectively $0,\pi ,0$, which gives the required values. This is the bounded harmonic solution, with $U\to V_0$ at infinity; uniqueness holds in this bounded class by the maximum principle. Without a boundedness or far-field condition, the edge data alone do not give uniqueness: one may add $C\sinh (\pi x/L)\sin (\pi y/L)$.

Back substitution. Sincethe required continuous argument difference is $2\arctan [\sin (ay)/\sinh (ax)]$ in the open half-strip. ConsequentlyAs $y\to 0,L$ with $x>0$, the arctangent tends to zero; as $x\to 0^{+}$ with $0<y<L$, it tends to $\pi /2$. These limits verify every open edge of the Dirichlet problem.

Set $A=x^2-y^2$ and $B=2xy$. Direct differentiation givesTheir sum vanishes, proving that $u$ is a harmonic function. More generally the Cauchy-Riemann equations give $u_{xx}+u_{yy}=v_{yx}-v_{xy}=0$ and $v_{xx}+v_{yy}=-u_{yx}+u_{xy}=0$; the mixed partial derivatives commute because a holomorphic function is smooth. Since $e^{i{z}^2}=e^{-B}(\cos A+i\sin A)$, a harmonic conjugate and the resulting holomorphic function areThese expressions also verify the Cauchy-Riemann equations directly.

Write $Z_t=X_t+i{Y}_t$ and $u=\mathrm{Re}g$. Since $g$ is entire, $u$ is a harmonic function, so $u_{xx}+u_{yy}=0$. The two-dimensional Itô formula givesThe drift vanishes. Stop on successively larger disks to bound the derivatives and define the stochastic integrals without global growth assumptions. ThusThis is the harmonic functions of Brownian motion principle.