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.

On the two independent variables, the total derivative is $D_x={\partial }_x+{ϵ}^{-1}{\partial }_X$. Multiplication of $D_x(a{D}_{x}y)=f$ by ${ϵ}^2$ gives exactly the stated multiple-scale equation.

At order one, ${\partial }_X(a{\partial }_X{y}_0)=0$. Integrating twice givesThe periodic reciprocal coefficient has strictly positive mean. Unless $A=0$, the integral has a nonzero linear drift in $X$, violating the periodic cell ansatz and creating secular growth when $X=x/ϵ$. Thus $A=0$ and $y_0=y_0(x)$.

At order $ϵ$,Periodicity requires the derivative to have zero cell average. ConsequentlyAt order ${ϵ}^2$, the equation isThe second flux inside braces is $C(x)$, independent of $X$. Integrating over one period annihilates the first derivative, by periodicity, and yields $C^{\prime }(x)={\int }_0^{1}f(x,X)dX$. HenceThe effective transport coefficient is the harmonic mean of $a$ over the fast period; the source is averaged arithmetically. Smoothness sufficient for these formal differentiations and uniform positivity on the cells are understood.

For fixed nonzero $x$, expand the denominator geometrically:The effective expansion parameter is ${ϵ}^2/x^2$, so the expansion fails when $x=O(ϵ)$. In the inner variable $X=x/ϵ$,for fixed $X$. Expanding the leading outer term toward $x=0$ gives $x^{-2}$ at leading order; expanding the leading inner term for large $X$ and returning to $x$ gives ${ϵ}^{-2}{X}^{-2}=x^{-2}$. Thus the two one-term re-expansions match in $ϵ\ll x\ll 1$. The next common contribution is $-x^{-1}$ when the corresponding next inner order is retained.

On $x\ge \delta$, the relative denominator error is bounded by ${ϵ}^2/{\delta }^2$, henceIntegration by parts yieldsFor the last step split the remaining integral at one: ${\int }_{\delta }^{1}dx/x=-\log \delta$, and replacing $e^{-x}$ by one on that interval has a bounded error.

In the inner interval, $|e^{-x}-1|\le x$ gives an explicit error estimate:The $1/\delta$ terms cancel between regions. Every displayed remainder is $o({ϵ}^{-1})$ whenever $ϵ/\delta \to 0$ and $\delta \to 0$; for example $ϵ({ϵ}^2/{\delta }^3)=(ϵ/\delta {)}^3$. Therefore the leading answer is independent of the arbitrary intermediate cutoff:

Define $y=k^{-1}\dot{x}+x^3/3-x$. Differentiation and the original equation give $\dot{y}=(\lambda -x)/k$, while rearrangement gives $\dot{x}=k(y-x^3/3+x)$. Write $h(x)=x^3/3-x$.

On the fast scale $\tau =kt$,Fast trajectories are horizontal. Their equilibria lie on the cubic $y=h(x)$; the transverse linearization is $-h^{\prime }(x)=1-x^2$, so the outer branches $|x|>1$ attract and the middle branch repels. On a horizontal line, arrows point right when $y>h(x)$ and left when $y<h(x)$.

The source prints the slow time as $T=k/t$. The slow variable needed here is $T=t/k$. With this correction,At leading order $y=h(x)$ andThe reciprocal printed variable would give $-T^2{y}_T=\lambda -x$, rather than the autonomous slow equations used in the stated period formula.

For $-1<\lambda <1$, the left attracting branch runs from $x=-2$ to the fold $x=-1$, where $y=2/3$. The horizontal fast jump at this height lands at $x=2$. The right branch then runs from $2$ to the fold $1$, where $y=-2/3$, followed by a fast jump to $-2$. The equilibrium at $x=\lambda$ is on the repelling middle branch. The two slow segments dominate the period:The jump times are negligible on this $O(k)$ scale. Substituting $u=-x$ in the first integral also gives the compact valueFor $\lambda =1/2$, the negative slow segment takes $k[1-\frac{3}{4}\log (5/3)]$, and the positive one takes $k[2-\frac{3}{4}\log 3]$. Thus $x(t)$ rises slowly from $-2$ to $-1$, jumps to $2$, decreases more slowly to $1$, then jumps back to $-2$; the jumps are nearly vertical in the time trace.

For $\lambda >1$ the fast horizontal field and attracting branches are unchanged. The right slow branch has arrows toward $x=\lambda$ from both directions; the left branch still moves toward $-1$ and jumps right. The equilibrium $(\lambda ,h(\lambda ))$ is stable, so the large relaxation cycle is replaced by approach to this equilibrium. If $\lambda >2$, a landing at $2$ moves right toward $\lambda$, rather than left. The sketches display the representative case $\lambda =3/2$.

Away from the boundary, derivatives of a smoothly varying solution are much smaller than $k^{2}f$ unless the outer solution vanishes. Nonzero boundary data therefore require a rapidly varying layer. If its thickness is $\delta$, the radial second derivative is of order $f/{\delta }^2$; balancing it with $k^{2}f$ gives $\delta =k^{-1}$.

Put $R=k(1-r)$ and write $f(r,\theta )=F(R,\theta )$. The equation becomesFor fixed $R$ and boundary data with angular variation independent of $k$, the leading equation is $F_{RR}=F$. Matching to the zero outer solution removes the growing exponential, and the boundary value at $R=0$ fixes its coefficient:This approximation describes the layer; it is not a uniform expansion in polar coordinates at the origin.

Laplace's method localizes an exponential integral near maxima of its real phase. At an isolated smooth interior maximum $t_0$ with $G^{\prime \prime }(t_0)<0$ and $F(t_0)\ne 0$, Taylor expansion and Gaussian integral giveAn endpoint maximum with nonzero inward phase slope instead gives an inverse-$k$ factor. Degenerate maxima require their first nonzero local phase term; continuous differentiability alone does not guarantee the quadratic formula.

Here $\cos (ix\sin \varphi )=\cosh (x\sin \varphi )$. In its two exponentials, the positive exponential dominates at the interior maximum $\varphi =\pi /2$ of $\sin \varphi$. Since $\sin (\pi /2+s)=1-s^2/2+O(s^4)$,The negative exponential has no comparable positive maximum and is negligible. Consequently, where $kr\to \infty$,In the boundary layer $r=1-R/k$, $r^{-1/2}=1+O(k^{-1})$, so this reduces to $e^{-R}$, agreeing with the formal result when $g=1$.

The Lax-Milgram theorem gives, for a bounded coercive bilinear form on a real Hilbert space and every bounded linear functional, a unique vector representing that functional through the form, with the inverse bound supplied by coercivity. The Fredholm alternative for a compact operator $K$ says that $I+\mu K$ is invertible if its kernel is zero; in general its range is closed and consists exactly of vectors orthogonal to $\ker (I+\mu K{)}^{\ast }$ in the Hilbert-space case.

Uniform ellipticity gives $\sum a_{ij}{\xi }_i{\xi }_j\ge \theta |\xi {|}^2$ for some $\theta >0$. Bounded coefficients make $B_{\lambda }$ bounded on $H^1$, and for $\lambda <-\gamma$,The forcing $v\mapsto \int fv$ is bounded on $H^1$. The Lax-Milgram theorem therefore gives a unique weak solution for every $\lambda <-\gamma$.

When $c=0$, choose the shift ${\lambda }_0=-1$ and let $Kf$ be the unique solution of $B[Kf,v]+\int (Kf)v=\int fv$. Thus $K=(L+I{)}^{-1}$ in the weak sense. The equation $Lu=f$ is equivalent to $u=K(f+u)$, orConversely a solution in $L^2$ of this operator equation lies in $H^1$ because it equals $K(f+u)$, and satisfies the variational equation. The assumed boundedness $K:L^2\to H^1$, followed by the compact Rellich embedding $H^1\hookrightarrow L^2$, proves compactness of $K$ as an $L^2$ operator.

By the supplied kernel assumption, $\ker (I-K)$ consists of constants. Also $K1=1$. Since $K$ is a self-adjoint operator, the Fredholm alternative condition is $Kf\perp 1$, andTherefore a weak solution exists exactly when the mean of $f$ is zero, and solutions differ by a constant. On a disconnected domain the actual kernel would consist of constants on each component and the condition would apply separately on each; the question's given kernel assumption specifies the single-constant case.

For a general coefficient matrix, the stated variational definition encodes zero conormal flux $(aDu)\cdot n=0$. It agrees with ordinary zero normal derivative only under additional conditions on $a$. The argument above uses the weak formulation explicitly defined in the paper, so does not require that stronger boundary identification.

The weak formulation isThe boundary condition is encoded by $u\in H_0^1$. The Sobolev embedding in dimension three gives $H_0^1\hookrightarrow L^6$, hence also $L^4,L^3,L^2$ on the bounded domain. All energy terms and all the displayed pairings are therefore finite.

Pointwise, $t^4/4-t^3/3\ge t^4/8-C$ for a constant $C$, since the difference is a quartic with positive leading coefficient. The Poincaré inequality and Young inequality give $\int fv\ge -\frac{1}{4}\Vert Dv{\Vert }_2^2-C_f$. ConsequentlyThus the infimum is finite and a minimizing sequence is bounded in $H_0^1$. Extract a weakly convergent subsequence there. The Rellich-Kondrachov theorem states that on a bounded $C^1$ domain in dimension three the embedding $H^1\hookrightarrow L^q$ is compact for $q<6$. A subsequence therefore converges strongly in $L^2,L^3,L^4$ to the same limit $u$. The potential terms and forcing converge, while the gradient norm is weakly lower semicontinuous. The weak limit retains zero trace. Hence $I(u)\le \mathrm{liminf}I(v_j)$ and $u$ attains the absolute minimum.

For any $\varphi \in H_0^1$, the function $t\mapsto I(u+t\varphi )$ has a minimum at zero. Differentiating its polynomial terms gives $\int [Du\cdot D\varphi +(u^3-u^2+f)\varphi ]=0$. These differentiations are justified by the Sobolev integrability already noted. This is the weak formulation, proving every energy minimizer is a weak solution.

Choose a smooth compactly supported mollifier of integral one. On every compact subset of $\Omega$, and for sufficiently small $\epsilon$, each translated mollifier is an allowed test, so the assumption gives $w\ast {\rho }_{\epsilon }=0$. The standard mollifier approximation theorem gives $w\ast {\rho }_{\epsilon }\to w$ in local $L^1$. Thus $w=0$ almost everywhere.

For a multi-index $\alpha$, the weak derivative $v=D^{\alpha }u$ is the locally integrable function satisfyingTwo such functions have zero pairing difference against every test, so the first result proves uniqueness. A compactly supported smooth test on an open subset $E$ extends by zero to one on $\Omega$, proving restriction of the derivative to $E$.

Put $w=D^{\beta }u$. For every test $\varphi$ the defining identity applied to $D^{\alpha }\varphi$ givesThus these two distributions are identical; one is represented by a locally integrable function if and only if the other is. Their weak derivatives then agree, proving $D^{\alpha }(D^{\beta }u)=D^{\alpha +\beta }u$ with exactly the asserted existence equivalence.

On $(0,2)$ the displayed piecewise affine function is continuous at one. Integration by parts separately on the two intervals shows its first weak derivative is one on $(0,1)$ and minus one on $(1,2)$; the first boundary terms cancel by continuity. The derivative of that step function is $-2{\delta }_1$. This distribution is not represented by a locally integrable function: tests of height one supported in shrinking neighborhoods of one have pairing $-2$ with it, whereas their integrals against any fixed local $L^1$ function tend to zero. Hence there is a first weak derivative but no locally integrable second weak derivative.

A continuous embedding $X\hookrightarrow Y$ is an injective continuous linear inclusion, so $\Vert x{\Vert }_Y\le C\Vert x{\Vert }_X$ for a fixed constant. On a finite-measure domain, Holder inequality giveswith $1/\infty =0$. For finite $p>q$ apply Holder inequality to $|u{|}^q\cdot 1$ with exponents $p/q$ and $p/(p-q)$; equality of exponents is immediate, and $p=\infty$ follows by integrating the supremum bound.

Weak convergence $u_j\rightharpoonup u$ means $\ell (u_j)\to \ell (u)$ for every bounded linear functional $\ell \in X^{\ast }$. The Riemann-Lebesgue lemma says that the Fourier transform of an integrable function tends to zero at high frequency. Since $L^2(0,1)\subset L^1(0,1)$, it makes $\cos (j\pi x)$ weakly null in $L^2$. Its absolute value is at most one, so $u_j{v}_j\in L^2$. For a test $\varphi \in L^2$,while $\int u_{j}v\varphi \to 0$ by the Riemann-Lebesgue lemma applied to $v\varphi \in L^1$. Hence $u_j{v}_j\rightharpoonup 0$. Taking $w_j=u_j=\cos (j\pi x)$ instead gives $u_j{w}_j=(1+\cos (2j\pi x))/2\rightharpoonup 1/2$, so two weakly null factors need not have a weakly null product.

For ${\bar{z}}_j=j^{1/p}{1}_{(0,1/j)}$ and $p^{\prime }=p/(p-1)$, Holder inequality givesThus ${\bar{z}}_j\rightharpoonup 0$ in $L^p$. In contrast, if $z_j$ had a weak $L^1$ limit $z$, testing on every interval strictly away from zero would give integral zero for $z$, hence $z=0$ almost everywhere by the permitted fact. Testing against the constant function one would give $\int z=1$, a contradiction. Therefore $z_j$ has no weak $L^1$ limit.