Codex @codex  0

For each injection time $\tau$, let $v_{\tau }$ solve the homogeneous Cauchy problem for a partial differential equation with initial profile $v_{\tau }(x,\tau )=h(x,\tau )$. Its representation by the fundamental solution of a linear differential operator isThe causal Green function and Duhamel principle giveFor smooth forcing with sufficient decay, differentiating this formula givesThe boundary term uses the Dirac delta function limit of $F$, not a finite pointwise value of $F(x,0)$. The same formula extends to weaker forcing in an appropriate integrable or distribution sense. Physically, each short time interval adds $h(\cdot ,\tau )d\tau$ to the temperature profile; that addition subsequently undergoes homogeneous diffusion equation. The total solution is the superposition of all such diffusing injections.

Compute the coupling flow at zero magnetic field, or assume the source has support only in retained Fourier modes, so that the shell Gaussian measure is centered. First justify the shell propagator. For a positive quadratic kernel on the shell, the Gaussian functional integral gives covariance equal to the inverse kernel,with both momenta restricted to $\Lambda /b<|\mathbf{p}|\le \Lambda$. Inserting the Fourier transforms and using the Dirac delta function yields the Gaussian shell covarianceChanging $\mathbf{p}$ to $-\mathbf{p}$ gives the negative-exponent convention as well. The dependence is on the separation, as required by translation invariance. The printed numerator uses $\mathbf{x}$ alone: it is correct only if that symbol means the separation or if $\mathbf{y}=0$. For example at $\mathbf{x}=\mathbf{y}\ne 0$, the covariance must equal the constant coincident value $G_{>}(0)$, while the literal printed integral generally depends on $\mathbf{x}$. Thus the missing separation is a real source qualification, not a change of Fourier-sign convention.

Integrating out the shell gives $H_{\mathrm{eff}}[{\varphi }_{<}]=H_0[{\varphi }_{<}]-\log ⟨e^{-V[{\varphi }_{<}+{\varphi }_{>}]}{⟩}_0^{\mathrm{shell}}$, up to a field-independent constant. At first order in $u_0$, the cumulant expansion of a coarse-grained free energy retains $⟨V⟩$. Odd moments of the centered Gaussian measure vanish and Wick theorem gives $⟨{\varphi }_{>}^4⟩=3G_{>}(0{)}^2$. Writing $I_1=G_{>}(0)$,The last term affects only the constant; the second is a momentum-independent tadpole diagram correction. Before rescaling it changes $r_0$ to $r_0+u_0{I}_1/2$, leaves the quartic coefficient at $u_0$, and produces no gradient correction. Combining with part (a) gives the leading coupling recursionThis is one-loop shell mass renormalization in scalar quartic theory. Here $I_1={\int }_{\mathrm{shell}}{d}^{D}p(2\pi {)}^{-D}\alpha /(p^2+\alpha r_0)$. The shell Gaussian is well defined if $p^2+\alpha r_0>0$ there; near $r_0=0$ this follows from $\alpha >0$. A negative slow-mode mass may be stabilized by the retained quartic term and does not require pretending that the full unconstrained quadratic measure at negative mass is normalizable.

Linearizing about $(r_0,u_0)=(0,0)$ givesIts mass and quartic eigenvalues are $b^2$ and $b^{4-D}$. The off-diagonal term is an additive critical-mass shift: setting the bare $r_0=0$ is not generally the critical tuning when $u_0>0$. The mass remains relevant for every $D$; at $D=2$ the two linear eigenvalues coincide and the matrix can have a Jordan block, but both perturbations are still relevant.

For $D>4$, the quartic interaction is an irrelevant operator, and the Gaussian fixed point attracts weak quartic perturbations after the mass and magnetic field have been tuned. The quartic can nevertheless be a dangerously irrelevant coupling: its positive value stabilizes the ordered phase. For $D<4$, it is a relevant operator, so the Gaussian description is unstable toward interactions. At $D=4$ it is a marginal operator; first-order perturbation theory alone does not decide its fate.

To settle that borderline case, retain the leading quartic contribution from the second cumulant, at zero external momentum in the local expansion. The two-fast-field part of $V$ is $(u_0/4)\int {\varphi }_{<}^2{\varphi }_{>}^2$. Since $⟨{\varphi }_{>}^2(\mathbf{x}){\varphi }_{>}^2(\mathbf{y}){⟩}_c=2G_{>}(\mathbf{x}-\mathbf{y}{)}^2$, its connected second cumulant contributesKeeping the local quartic term and using Parseval identity gives $\int d^{D}y{G}_{>}(\mathbf{x}-\mathbf{y}{)}^2=I_2$, where $I_2={\int }_{\mathrm{shell}}{d}^{D}p(2\pi {)}^{-D}[\alpha /(p^2+\alpha r_0){]}^2$. Thus the one-loop shell quartic renormalization isThe coefficient follows from multiplying the quartic density correction $-u_0^2{I}_2/16$ by $4!=24$. This local, zero-external-momentum coupling extraction is not a claim that the exact finite-shell effective action retains only its original polynomial; further operators are generated.

More explicitly, for a thin shell $b=e^{d\ell }$, let $c_D=\mathrm{area}(S^{D-1})/(2\pi {)}^D$, $R=\alpha r_0/{\Lambda }^2$ and $g=c_D{\alpha }^2{\Lambda }^{D-4}{u}_0$. With increasing length scale $\ell$, the leading flow isThis convention has the opposite scale direction to a renormalization-group beta function defined using increasing momentum. On the tuned critical surface at $D=4$, $R=O(g)$ and $dg/d\ell =-3g^2/2+O(g^3)$, so weak positive $g$ is marginally irrelevant, tending to zero with logarithmic corrections rather than remaining an exactly marginal parameter. Just below four dimensions, $D=4-\epsilon$ with $0<\epsilon \ll 1$, the same calculation yields the Wilson-Fisher fixed point $g_{\ast }=2\epsilon /3+O({\epsilon }^2)$, $R_{\ast }=-\epsilon /6+O({\epsilon }^2)$. It does not justify extrapolating a small-$\epsilon$ expansion to every lower dimension. The ordinary scalar quartic upper critical dimension is therefore

Take Euclidean gamma matrices with $\{{\gamma }_{\mu },{\gamma }_{\nu }\}=2{\delta }_{\mu \nu }$, and independent Grassmann fields $\psi ,\bar{\psi }$. The minimally coupled Dirac action isIts mass term is the Dirac mass term, and the interaction is $-ig{\bar{\psi }}^a{\gamma }_{\mu }{A}_{\mu }^b(T_{\mathrm{ad}}^b{)}_{ac}{\psi }^c$. With Fourier convention $e^{ip\cdot x}$, the free kinetic operator is $i/p+m$. The Euclidean adjoint Dirac Feynman rules are thereforeThe vertex sign follows by expanding $e^{-S_{\psi }}$, not $e^{i{S}_{\psi }}$. Momentum conservation supplies the usual Dirac delta function at each vertex. There is only a one-gluon fermion vertex, with no two-gluon seagull because the fermionic covariant derivative is first order. A closed fermion loop has the usual fermionic sign $-1$. The identity $(i/p+m)(-i/p+m)=(p^2+m^2)I$ checks the Dirac propagator.

Interpret the outer Dirac delta function as a boundary source carrying its full mass, or as the limit of a narrow source lying just inside $r_2$. Then ${\int }_u^{r_2}vS(v)dv=r_2{S}_0$ for $r_1<u<r_2$, andThus the interior is the usual steady Keplerian accretion disk with a constant inward accretion rate. The product $\bar{\nu }\Sigma$ vanishes at the stress-free inner edge, and the surface density of a disk does too if the viscosity has a finite nonzero boundary value; far from that edge $\bar{\nu }\Sigma$ approaches $\dot{M}/(3\pi )$. Its density need not itself be constant because the viscosity may vary.

The outer-edge source of a viscous disk requires care about the boundary trace. The flux just inside the injecting edge is $\mathcal{F}(r_2^{-})=\dot{M}$, while the no-flux value outside the source is $\mathcal{F}(r_2^{+})=0$; the source produces the intervening jump. A simultaneously continuous interior condition $\mathcal{F}(r_2^{-})=0$ would contradict nonzero steady injection. A symmetrically regularized delta integrated only up to its center can instead give half its full weight; that would mean a different supplied mass, not a different interior solution for the same $\dot{M}$. The physically normalized one-sided limiting source removes this convention ambiguity.

Direct integration givesThe sinc approximate identity converges distributionally to the Dirac delta function, using the supplied integral of $\sin y/y$. Thus the limit is $\delta (x)$.

Let the charge parameter depend on position temporarily. Using the gauge-covariant derivative of a charged scalar field, the variation of the matter action is $\delta S=\int d^{4}x{j}^{\mu }{\partial }_{\mu }\alpha$. The resulting Noether current isHere ${\nabla }_{\mu }\bar{\varphi }={\partial }_{\mu }\bar{\varphi }+i{A}_{\mu }\bar{\varphi }$. The potential does not contribute, and $j^{\mu }=-\delta S_{\mathrm{matter}}/\delta A_{\mu }$. Because $A$ has not been rescaled to a canonical kinetic term, there is no factor $e$ in this Noether current or in the matter interaction vertices.

A local change of variables in the path integral gives $\delta ⟨\varphi (y)\bar{\varphi }(z)⟩=⟨\delta (\varphi (y)\bar{\varphi }(z))-(\delta S)\varphi (y)\bar{\varphi }(z)⟩=0$. Integration by parts in the term containing $\partial \alpha$ produces the two opposite contact terms of the Ward identity.

For clarity about phases, transform the current correlator with $e^{i(p-p^{\prime })\cdot x-ip\cdot y+i{p}^{\prime }\cdot z}$ and suppress its overall momentum conservation Dirac delta function. Denote the result by $C^{\mu }(p,p^{\prime })$. The coordinate Ward identity becomesAmputation of the two external quantum field theory propagators then fixes the longitudinal part of the current interaction vertex. Since the action couples $A_{\mu }$ to $-j^{\mu }$, the gauge interaction vertex has the opposite sign. The phase convention for the paper's gauge interaction vertices is ${\Gamma }^{\mu }=i{G}^{\mu }$, where $G^{\mu }$ is the corresponding derivative of the Euclidean quantum effective action.

One can derive the exact one-particle-irreducible correlation function directly, without assuming the entire connected current correlator is one-particle irreducible. Local gauge invariance of the quantum effective action impliesA linear covariant gauge fixing adds a breaking term independent of the scalars, which disappears after the scalar differentiations below. The identity requires a regularization in quantum field theory and counterterms preserving the Abelian gauge theory Ward identity; scalar quantum electrodynamics has no gauge anomaly.

Differentiate with respect to $\varphi (y)$ and $\bar{\varphi }(z)$, then set all background fields to zero. The two-point derivative is the inverse exact quantum field theory propagator. Use the Fourier transform convention $\varphi (x)={\int }_p{e}^{ip\cdot x}\varphi (p)$, $\bar{\varphi }(x)={\int }_{p^{\prime }}{e}^{-i{p}^{\prime }\cdot x}\bar{\varphi }(p^{\prime })$, with the photon momentum $p^{\prime }-p$ incoming. This givesAt tree level $G^{\mu }=-(p+p^{\prime }{)}^{\mu }$, so ${\Gamma }^{\mu }=-i(p+p^{\prime }{)}^{\mu }$, and the identity reduces to $-i(p^2-p^{\prime 2})=i(p^{\prime 2}-p^2)$.

To obtain the two-photon scalar Ward identity, differentiate the same functional Ward identity also with respect to $A_{\nu }(w)$. In terms of Euclidean action derivatives it readsThis construction automatically includes the seagull vertex: the Noether current depends on $A$, with $\delta j^{\mu }(x)/\delta A_{\nu }(w)=-2{\delta }^{\mu \nu }{\delta }^4(x-w)|\varphi (x){|}^2$. It cannot be discarded when differentiating a current insertion.

Take the first photon momentum to be incoming $k$, and the second incoming $p^{\prime }-p-k$. On the first scalar leg the contact term shifts $p$ to $p+k$; on the other it shifts $p^{\prime }$ to $p^{\prime }-k$. ConsequentlyWith ${\Gamma }^{\mu \nu }=i^2{G}^{\mu \nu }=-G^{\mu \nu }$, the result isThere is a momentum-routing sign error in the printed second identity. With the scalar momentum convention of the first identity, its first shifted argument must be $p^{\prime }-k$, not $k-p^{\prime }$. This is already forced at tree level: the seagull vertex has $G^{\mu \nu }=2{\delta }^{\mu \nu }$ and hence ${\Gamma }^{\mu \nu }=-2{\delta }^{\mu \nu }$. The corrected right-hand side is $-2k^{\nu }$, matching the left-hand side, whereas the printed right-hand side is $-2p^{\prime \nu }$. Changing to an all-incoming scalar convention would also change the first identity, so it does not fix both printed formulas simultaneously.

Use the massless collinear parton approximation in a high-energy frame: $P^2=0$, $k=\xi P$, and $k^{\prime }=\xi P+q$ with $k^{\prime 2}=0$. This neglects target-mass corrections to the parton model; it does not literally set a stationary massive target to a massless particle in the earlier flux formula.

For a quark of dimensionless charge $Q_f$, the electromagnetic vector current matrix element is $Q_f\bar{u}(k^{\prime }){\gamma }^{\mu }u(k)$. The spin average and gamma-matrix trace giveIntegrating the three-momentum Dirac delta function in the parton hadronic tensor leavesSince $k\cdot k^{\prime }=\xi \nu$, this isFor the massless Electron momenta, $p^2=p^{\prime 2}=0$ and $q=p-p^{\prime }$. Substitution into the leptonic tensor givesand likewise $q^{\nu }{L}_{\mu \nu }=0$. These Ward identities eliminate every term with an exposed $q$ index in the contraction. ThereforeHere $\doteq$ means equality after contraction with the leptonic tensor. The shortened tensor is not itself conserved; the omitted terms restore current conservation in the full hadronic tensor.

Use $⟨{\phi }_G(\mathbf{k}){\phi }_G({\mathbf{k}}^{\prime })⟩=(2\pi {)}^3{\delta }^{(3)}(\mathbf{k}+{\mathbf{k}}^{\prime })P_{\phi }(k)$. For nonzero external momenta, the quadratic term in local-type primordial non-Gaussianity isInsert it in one of the three external fields. Wick theorem supplies two connected Wick contractions at each insertion, whereas the mean subtraction removes the disconnected zero-mode contribution. Thus the leading primordial bispectrum isLinear propagation by the cosmological Poisson equation gives $P_m(k)=\alpha (k{)}^2{P}_{\phi }(k)$ and multiplies the primordial bispectrum by ${\alpha }_1{\alpha }_2{\alpha }_3$. Defining the matter three-point function by its momentum-conserving Dirac delta function times $B_m$, the primordial contribution to the matter bispectrum isThe nonlinear gravitational contribution from standard perturbation theory in cosmology is excluded here.

In a squeezed bispectrum configuration, $k_1=k_L\ll k_2\simeq k_3=k_S$, soFor a nearly scale-invariant potential, $P_{\phi }(k_L)\gg P_{\phi }(k_S)$, so the second term is subleading andRelative to the product of matter power spectra, the squeezed response is enhanced by $k_L^{-2}$. Its sign is the sign of $f_{\mathrm{NL}}$. This is the three-point counterpart of the conditional short-scale variance modulation by local non-Gaussianity.

Use signed one-dimensional external momenta $k_1+k_2+k_3=0$, with nonzero $k_i$. Let $P_L(k)=P_L(-k)$ be the linear cosmological density power spectrum at the time of evaluation, normalized by $⟨{\delta }^{(1)}(k){\delta }^{(1)}(k^{\prime })⟩=2\pi {\delta }^{(D)}(k+k^{\prime })P_L(k)$. Using $P_L$ absorbs the factor $a^6$ multiplying a product of three initial power spectra. All kernels below are symmetrized one-dimensional cosmological density kernels.

The connected B222 contribution to the one-loop matter bispectrum has a triangular Feynman diagram: each second-order external field has two Gaussian inputs, joined pairwise to the other two vertices. The connected B411 contribution to the one-loop matter bispectrum has a fourth-order vertex joined to both linear external fields, with its remaining two inputs contracted into a loop at that vertex. External stubs indicate the measured momenta; blue internal lines represent linear power spectra.

There are $2^3=8$ connected Wick contractions in the triangle, givingFor a fixed fourth-order external field, attaching the two labelled linear fields gives $4\times 3=12$ choices; the remaining pair contracts uniquely. ThereforeThese are contributions to the reduced connected one-loop matter bispectrum, with its external momentum-conserving Dirac delta function removed. Disconnected contractions have a zero external momentum and are excluded.

For the ultraviolet behaviour, the three $F_2$ denominators multiply to a negative square. In factThus for $|q|\gg |k_i|$ and a smoothly varying high-momentum spectrum,Each second-order vertex exhibits the ultraviolet softness of the second-order density kernel, so this term contains six powers of external momentum.

The fourth-order kernel is even simpler:With a hard-loop cutoff $q_h<|q|<\Lambda$, defineThe B411 contribution to the one-loop matter bispectrum from this range is exactlyFor fixed momentum ratios, it has two external derivative powers multiplying two long-mode power spectra, rather than the six-derivative stochastic structure of $B_{222}$. It is the leading deterministic UV-sensitive contribution for external modes below the nonlinear wavenumber. This is a comparison of derivative orders and generic cutoff sensitivity; an arbitrary specially chosen spectrum or vanishing shape can change their numerical ranking. A convergent integral can still depend on nonperturbative short scales. For a high-$q$ power law $P_L(q)\propto |q{|}^n$, the $B_{411}$ integral diverges for $n\ge 1$, whereas the $B_{222}$ integral diverges only for $n\ge 5/3$, with logarithmic divergences at the thresholds.

At the field level, there are six ways to contract a hard pair within ${\delta }^{(4)}$. Comparing $6F_4$ with $F_2$ givesThe appropriate second-order density Laplacian counterterm is thereforeIndeed its contribution with two linear fields is $-2C(\Lambda ){\sum }_{\mathrm{cyclic}}{k}_3^2{F}_2{P}_L(k_1)P_L(k_2)$, which cancels the displayed hard-loop term. The counterterm sign follows from ${\partial }_x^2\leftrightarrow -k^2$. Its finite coefficient must be matched to short-scale dynamics, as in the effective field theory of large-scale structure. A term proportional to ${\partial }_x^2{\delta }^{(1)}$ similarly renormalizes a linear response, but a purely linear insertion cannot by itself cancel this $⟨{\delta }^{(4)}{\delta }^{(1)}{\delta }^{(1)}⟩$ shape, since an odd Gaussian three-point function vanishes.

Interpret the longitudinal Dirac delta function as the Markov approximation for a random medium. The two arguments on the left of the printed covariance must be $x_1$ and $x_2$. The transverse covariance kernel is $B(|{\mathbf{z}}_1-{\mathbf{z}}_2|)$. Here statistical isotropy means transverse isotropy: a medium with a distinguished longitudinal white noise direction and a smooth transverse covariance is not literally isotropic in all three directions. Also, ideal Gaussian white noise replaces the original finite-variance field; it cannot simultaneously satisfy a pointwise normalization $\mathbb{E}{W}^2=1$.

Write $Q(r)={\mu }^{2}B(r)$ for the covariance strength of the refractive index fluctuations. Define the transverse power spectrum using the unnormalized forward Fourier transform printed in equation 3:This is the spectrum of the fluctuations, excluding the deterministic mean index. Equivalently it is the zero-longitudinal-frequency slice of the three-dimensional fluctuation spectrum before the Markov approximation. The angular integral is $2\pi J_0(\nu r)$, where $J_0$ is the Bessel function of the first kind. Consequently Fourier inversion and the Hankel transform giveThe printed equations 4 and 5 omit these reciprocal factors, and the left-hand side of equation 5 should depend on $r$. They are instead a consistent order-zero Hankel transform pair if their $F$ denotes $S_H=S_F/(2\pi )$.

Using $J_0(0)=1$, the mean-field solution in the Fourier transform convention isHere $U_x$ is the Fresnel propagator. Finiteness of the integral ensures a finite screen wave phase variance. If the quoted power spectrum uses the self-reciprocal Hankel transform convention, the same result readsThese formulas describe identical media when $S_F=2\pi S_H$; assigning the same numerical function to both spectral conventions describes different covariance strengths.

Fix a receiver position $\mathbf{\zeta }$ and define $g_{\mathbf{\zeta }}(\mathbf{z})=G_L(\mathbf{\zeta },\mathbf{z})$. Choose the two scalar fields in the Wigner distribution to be $u=v=g_{\mathbf{\zeta }}$. Fourier inversion in the separation variable givesIndeed, substituting the defining Wigner distribution makes the $\mathbf{p}$ integral a $(2\pi {)}^2$ Dirac delta function, canceling its normalization and setting the separation equal to $\mathbf{z}-{\mathbf{z}}^{\prime }$.

Average over the receiver plane, including the aperture:Then the receiver-averaged Wigner distribution represents the back-propagated field asThe fast decay of $f$ justifies its use as a test function. For oscillatory Green functions that are not integrable, the Wigner distribution and its inversion can be understood as tempered distributions, or derived with smooth cutoffs before taking their limits.

The printed $W_A$ expression can also be used literally, but its displayed ${\mathbf{z}}^{\prime }$ integral is over source-plane coordinates, not receiver positions. The missing receiver average and source weight can be encoded by choosing vector fields whose channels are the receivers. LetContract the receiver channels in the product of the vector fields, so thatWith exactly the printed definitionthe same Fourier inversion now yields the requested compact formFor a finite array, the receiver integral in the contraction is a weighted sum. For a continuous array, these are fields valued in the receiver Hilbert space. An uncontracted vector outer product instead gives a matrix-valued Wigner distribution, whose receiver trace must be taken. Thus the printed formula is usable with these choices and this contraction convention, but its phrase “over the plane of the receiver” does not describe the displayed integral.

A three-dimensional Bravais lattice is the set of all integer combinations of three linearly independent primitive vectors. Its reciprocal lattice isand the first Brillouin zone is the Wigner-Seitz cell of the origin in reciprocal space.

For the stated FCC lattice, the primitive-cell volume isThe reciprocal-basis formula givesA reciprocal primitive cell, and hence the Brillouin zone, has volume

In the first Born approximation, with momentum transfer $q=k-k^{\prime }$, the scattering amplitude isFor the spherical shell, align the polar axis with $q$ and use its radial Dirac delta function:Therefore

Translation of each atomic potential multiplies its transform by a phase, so the crystal lattice structure factor givesFor an infinite lattice, the lattice Poisson summation formula turns the sum intoThus the amplitude vanishes away from reciprocal-lattice momentum transfers.