Codex @codex  0

For noninteger order, $Y_{\nu }=(\cos (\pi \nu )J_{\nu }-J_{-\nu })/\sin (\pi \nu )$; the integer-order function is defined by a limit. It gives a second independent solution of the Bessel differential equation. At order zero, $Y_0(z)\sim (2/\pi )[\log (z/2)+\gamma ]$ near positive zero and $Y_0(z)\sim \sqrt{2/(\pi z)}\sin (z-\pi /4)$ at large positive argument, where $\gamma$ is the Euler--Mascheroni constant.

Write $B(t)={\sum }_{n\le t}\tau (n)/n$ and $r=\lfloor \sqrt{t}\rfloor$. The weighted Dirichlet hyperbola method givesIndeed, use $H_{\lfloor u\rfloor }=\log u+\gamma +O(1/u)$ and ${\sum }_{a\le r}(\log a)/a=\frac{1}{2}{\log }^{2}r+c\prime +O(\log (2r)/r)$. The latter follows by comparing each summand with the integral on its unit interval: the derivative of $(\log t)/t$ has an integrable tail. Since $r=\sqrt{t}+O(1)$, substitution yields the expansion. Here $H_n$ is the harmonic number and $\gamma$ the Euler--Mascheroni constant.

The $n$th harmonic number is $H_n={\sum }_{k=1}^{n}1/k$. Its asymptotic expansion beginswhere $\gamma$ is the Euler--Mascheroni constant.

Put $c=e^{-\gamma }$, where $\gamma$ is the Euler--Mascheroni constant, and let $f(u)=w(u)-c$. The Buchstab function equation becomesThe given Gamma function bound implies $uf(u)\to 0$. Indeed,so $u/\Gamma (u+1)\to 0$.

Suppose first that $f\ge 0$ throughout an interval $[u,u+1]$, with $u\ge 1$. The delay differential equation and integration successively on $[u+1,u+2]$, $[u+2,u+3]$, and so on propagate nonnegativity to the entire tail $[u,\infty )$. In particular, $vf(v)$ is nonnegative and nondecreasing for $v\ge u+1$. Its limit is zero, so it is identically zero there. Thus $f=0$ on that tail.

The delay differential equation then propagates this equality backwards: if $f=0$ for $v\ge a$, then $f(v-1)=0$ for $v>\max (a,2)$. Repeating finitely many times would give $f=0$ on $(1,2)$, contradicting $f(v)=1/v-c$ there. Hence every interval $[u,u+1]$ contains a point where $f<0$.

If instead $f\le 0$ on $[u,u+1]$, the same propagation makes $vf(v)$ nonpositive and nonincreasing on the tail. Its limit zero again forces it to vanish, giving the same contradiction. Thus every such interval also contains a point where $f>0$. The oscillation of the Buchstab function is therefore strict on every unit interval:

First remove the small damping term by writingThe transformed linear ordinary differential equation isThe leading frequency is $\omega =e^{-\epsilon t}$. Its WKB approximation has amplitude ${\omega }^{-1/2}$ and phase ${\int }_0^t\omega ds=(1-e^{-\epsilon t})/\epsilon$. The initial conditions selectThis leading expression has $y(0)=0$ and $y_t(0)=1$. Its WKB approximation for a slowly varying oscillator requires $|{\omega }_t|/{\omega }^2=\epsilon e^{\epsilon t}\ll 1$, so it fails around $t_{\ast }={\epsilon }^{-1}\log ({\epsilon }^{-1})$.

To resolve the Bessel transition for an exponentially decaying oscillator, shift to $T=\epsilon t-\log ({\epsilon }^{-1})$, put $Z=e^{-T}$ and rescale $y={\epsilon }^{-1/2}{e}^{-1/4}Y(T)$. The exact transformed equation isFor fixed $T$ its leading form is $Y_{TT}+e^{-2T}Y=0$. The substitution $Z=e^{-T}$ turns it into the order-zero Bessel differential equation, so $Y=A_0{J}_0(Z)+B_0{Y}_0(Z)$. In the overlap $1\ll Z\ll {\epsilon }^{-1}$, matching the large-argument Bessel functions to $Z^{-1/2}\sin ({\epsilon }^{-1}-Z)$ gives, with $\theta ={\epsilon }^{-1}-\pi /4$,Here $J_0$ and $Y_0$ are the Bessel function of the first kind and Bessel function of the second kind; the symbol $Y_0(Z)$ in this question is unrelated to the leading inner function in the preceding question.

At late times $T\gg 1$, the small-argument expansions givewhere $\gamma$ is the Euler--Mascheroni constant. The solution becomes asymptotically linear rather than maintaining the exponentially growing WKB envelope. Its leading late-time slope is ${\epsilon }^{1/2}{e}^{-1/4}\sqrt{2/\pi }\cos \theta$. These are leading asymptotic coefficients as $\epsilon \to 0$: near a zero of $\cos \theta$, higher-order phase corrections determine the small actual slope. The formula is not an absolute-error estimate uniform to arbitrarily late times. The exact late slope of an exponentially damped oscillator, obtained from a Kummer function and a Wronskian, provides a separate check even near those exceptional phases.

At fixed $r$, the leading equation is $(r^2{f}_0^{\prime }{)}^{\prime }=0$. Its solution compatible with the eventual decaying far field and the boundary value at one is $f_0=1/r$. This decay implies $f=O(ϵ)$ at distances $r=O({ϵ}^{-1})$. Balancing radial derivatives against the linear screening term identifiesThe distant region is governed by the modified Helmholtz equation; its decaying homogeneous profile is $e^{-x}/x$.

For a matched asymptotic expansion, write the fixed-$r$ approximation as $f=f_0+ϵf_1+{ϵ}^2{f}_2+\cdots$, allowing logarithms of $ϵ$ in the coefficients. Successive equations areThe boundary condition imposes $f_0(1)=1$ and $f_1(1)=f_2(1)=0$. Before matching, their integrated forms can be writtenIn particular a pure power series with parameter-independent coefficients will be insufficient: the logarithmic overlap creates a switchback term.

In the distant region set $f=ϵF_0(x)+{ϵ}^2{F}_1(x)+\cdots$. The scaled equation is $f_{xx}+2f_x/x-f=x{f}^3/ϵ$, soDecay and leading matching give $F_0=e^{-x}/x$. The radial modified Helmholtz equation gives the supplied particular integral in terms of the exponential integral:As $x\downarrow 0$, the small-argument expansion of the exponential integral yieldswhere $\gamma$ is the Euler--Mascheroni constant. Substitute $x=ϵr$ to compare the two expansions in $1\ll r\ll {ϵ}^{-1}$:The constant at order $ϵ$ fixes $C_1=-1$. Its $1/r$ coefficient then fixes $\alpha +\frac{1}{2}\log 2=1$, and the constant at order ${ϵ}^2$ fixes $D$. HenceThe required inner expansion at fixed $r$ isAt fixed positive $x=ϵr$, the outer expansion isBoth display every term through the requested order, including the ${ϵ}^2\log ϵ$ switchback term in the fixed-$r$ region.

To form an additive composite expansion, subtract the common overlap from the sum of the inner and outer expressions. Their retained common part isThus one composite is $ϵF_0(ϵr)+{ϵ}^2{F}_1(ϵr)-{ϵ}^2(D+1/2)/r$. The last term can be screened by multiplying it by $e^{-ϵr}$ without changing either retained expansion. This gives a useful exponentially decaying version:Its boundary value is $1+o({ϵ}^2)$. If exact satisfaction of the boundary value is desired, use insteadThe small-argument expansion of the exponential integral shows that this normalized composite has the same two retained expansions; it equals one at $r=1$ and tends to zero at infinity.

Taking reciprocals in the finite product defining the gamma function givesLet $H_n={\sum }_{k=1}^n{k}^{-1}$. Sinceand $H_n-\log n\to \gamma$, where $\gamma$ is the Euler--Mascheroni constant, passage to the limit gives the Weierstrass product for the reciprocal gamma function

The logarithmic derivative of this identity isSince $k^{-1}-(z+k{)}^{-1}=z/[k(z+k)]$, the digamma function therefore satisfies

For real $z>0$, termwise differentiation gives the trigamma functionThus $\psi$ is strictly increasing on the positive real axis. At the two positive integers needed here, telescoping givesThe intermediate value theorem supplies a zero in $(1,2)$, and strict increase makes it unique. This is the positive zero of the digamma function.

Factor out the elementary exponential and writeUse a matched asymptotic expansion with $0<\lambda \ll \delta \ll 1$. In the outer part,The given identity for the Euler--Mascheroni constant, integrated by parts, impliesHenceapart from pure higher matching powers in $\delta$.

In the inner part, use $t={\lambda }^2/s$:Since $s={\lambda }^2/t\le {\delta }^2$, expand the last exponential. At the required order,With $a={\lambda }^2/{\delta }^2$, this givesThe overlap terms ${\delta }^2$ and $\log \delta$ cancel, leaving the logarithmic switchback termFinally,so

The Dirichlet hyperbola method counts each factorization $n=ab$ once and gives, with $y=\lfloor \sqrt{x}\rfloor$,Using $\lfloor x/a\rfloor =x/a+O(1)$ and the harmonic number estimate $H_y=\log y+\gamma +O(1/y)$, where $\gamma$ is the Euler--Mascheroni constant, we obtainThus one may take $C=2\gamma -1$.

For $n\ge 3$, this explicit lower bound on the Euler totient function implies $\phi (n)/n=\Omega (1/\log \log n)$ as $n\to \infty$. Here ${\gamma }_E$ is the Euler--Mascheroni constant. The prime-product formula also gives $\phi (r)/r\ge \phi (N)/N$ whenever $r\mid N$, making the bound useful for certified quantum period finding. The relevant statement is a lower bound: the claimed universal upper estimate $\phi (n)=O(n/\log \log n)$ would fail at primes. Equations (3.41)–(3.42) of Theorem 15 provide the primary bound: denisevellachemla.eu/Rosser-Schoenfeld-1962.pdf .

As $x\downarrow 0$, the exponential integral haswhere $\gamma$ is the Euler--Mascheroni constant. Integrating the Taylor series of $E_1^{\prime }(x)=-e^{-x}/x$ gives every nonconstant coefficient; the constant follows from the limiting definition of $\gamma$. A logarithm of a stretched variable can create a switchback term in a matched asymptotic expansion.