Codex @codex  0

A rising parcel cools, lowering the saturation pressure. Once its water vapor supply exceeds that value, atmospheric supersaturation develops. Homogeneous nucleation of pure water is strongly inhibited at the small atmospheric supersaturations typical of clouds; existing aerosol particles provide cloud condensation nuclei. Their dry sizes and composition, together with the parcel's ascent rate, determine how many cloud droplets activate.

Köhler theory combines the Kelvin equation, which raises equilibrium vapor pressure over a curved surface, with Raoult's law, which lowers it for a chemical solution. In a dilute diameter convention, with dry diameter $D_d$ and wet diameter $D$,Here $s=e_v/E(T)-1$ is atmospheric supersaturation, ${\sigma }_w$ is surface tension, and ${\kappa }_h$ measures aerosol hygroscopicity, not the dry-adiabatic process exponent. The maximum of the Köhler curve is the critical supersaturation. Larger or more hygroscopic particles generally have lower thresholds. Cloud-droplet activation requires overcoming that maximum; total aerosol number is not the activated number.

Given an observed aerosol number size distribution $n_j(D_d)=d{N}_j/d\ln D_d$ for composition classes $j$, first assign each class a hygroscopicity and compute its critical supersaturation. If the parcel reaches maximum atmospheric supersaturation $s_{\max }$, the equilibrium-threshold estimate isFor each fitted lognormal distribution, this is its upper-tail integral; it can be evaluated analytically with an error function or numerically from measured bins.

The value $s_{\max }$ must be calculated, not guessed from the aerosol distribution alone. Prescribe parcel upward velocity $w$, initial pressure, temperature and humidity, then integrate expansion cooling, latent heat release, conservation of total water, and diffusional cloud droplet growth. Schematically,where $a_s,b_s>0$, $G$ includes water vapor diffusion and thermal conduction, and $q_l$ is liquid mass fraction. Rising motion produces atmospheric supersaturation; growing particles consume it. Their competition selects the peak and which particles cross their activation barriers. Assume a well-mixed closed parcel, liquid water, soluble spherical particles, specified mixing state, and negligible collision-coalescence and fallout during activation. Slow kinetic growth, weakly soluble particles or entrainment require a more detailed treatment than the threshold integral.

Write ${\tau }_0={\tau }_{\infty }$, $n=2a$, and $t={\tau }_0(p/p_s{)}^2$. The integrated Planck function is $B(T)=\sigma T^4/\pi$ by the Stefan–Boltzmann law, andThe formal radiative transfer equation solution for upward radiance at the top isThe first term is transmitted surface emission; the second is emission from the atmosphere. Integrating radiance over the upward hemisphere gives the outgoing radiant flux per unit area,Here $E_j$ is an exponential integral. For fixed $a$ and ${\tau }_0\to \infty$, surface transmission is negligible and the upper limit can be extended to infinity. Positivity permits interchanging the integrals; putting $t=\mu x$ givesConsequently the exact hemispheric optically thick result isThis is a dimensional scale times a dimensionless definite integral depending only on $a$. In the commonly used vertical-ray approximation, define $F_{\mathrm{v}}=\pi I(0,1)$. ThenA vertical-ray approximation therefore has the same scaling but omits the angular factor $1/(a+1)$. Stating which quantity is used prevents a purely angular normalization discrepancy.

For $a=0,\frac{1}{2},1,\frac{3}{2},2,\dots$, put $n=2a=0,1,2,3,4,\dots$. Repeated integration by parts gives $\Gamma (n+1)=n!$. The hemispheric thick-limit results areThe second row assumes fixed surface temperature; the third assumes fixed absorbed radiant flux. In the vertical-ray approximation the first of these rows becomes $1,1/{\tau }_0,2/{\tau }_0^2,6/{\tau }_0^3,24/{\tau }_0^4$, and the surface temperature factors are the inverse fourth roots.

Thus increasing $a$ makes the atmosphere colder at every fixed $0<p<p_s$, decreases outgoing radiation at fixed $T_s$, and requires a warmer surface at fixed absorbed radiant flux. The case $a=0$ is isothermal and gives $F=\sigma T_s^4$ exactly, including finite optical depth.

There is an essential approximation limit. The extension of the integral to infinity is valid for fixed $a$ as ${\tau }_0\to \infty$, and is not uniform for arbitrarily large $a$ at fixed ${\tau }_0$. The vertical integrand peaks at $t=2a$, so one needs that peak well above the surface in optical depth, with a negligible tail beyond ${\tau }_0$. For example, at ${\tau }_0=10$ the thick formula predicts $F/(\sigma T_s^4)=20!/(1110^{20})\simeq 0.00221$ for $a=10$ but approximately $0.04773$ for $a=12$. This artificial increase contradicts the exact monotonicity: in the finite-optical depth expression every atmospheric emission factor $(t/{\tau }_0{)}^{2a}$ decreases with $a$. At fixed finite ${\tau }_0$, the exact limiting outgoing radiant flux is the transmitted surface term $2\sigma T_s^4{E}_3({\tau }_0)$, rather than factorial divergence.

Take positive heat flux upwards, along increasing $z$, and denote thermal conductivity by $k$. Fourier's law gives $q_s=-k{T}_z{|}_{z=a}$. Since $\cos (n\pi )=(-1{)}^n$, differentiation of the literal printed expression givesIts negative young-age flux is another sign of the unphysical initial profile. For the explicitly corrected hot initial plate, the surface heat flux of a cooling plate is insteadThe old-age limit is $k\Delta T/a$. For $t\ll {\tau }_0$, replacing the Gaussian sum by its leading integral gives $q_s\sim k\Delta T/\sqrt{\pi \kappa t}$, recovering the half-space cooling model. Neither $k$ nor $\kappa$ is numerically specified, so the age-dependent flux remains symbolic as requested.

The depth-averaged added horizontal compressive stress is $\bar{\sigma }=F/a$. Resolving its traction on a plane dipping at angle $\delta$ gives the fault shear stress from a plate-integrated force,Keeping the previous topographic correction gives about $13.8\mathrm{MPa}$. This is the mean additional shear from ridge push, assuming the force is transmitted through thickness $a$ and resolving the nonisotropic stress increment. An isotropic lithostatic pressure adds normal traction but no shear. The force integral alone does not specify the local stress at each point on the fault or establish failure.

Let $X_f$ have the specified probability density function $f_m$. Splitting the integral into the positive and negative half-lines, and using the Gamma integral ${\int }_0^{\infty }{x}^r{e}^{-sx}dx=r!/s^{r+1}$ for a nonnegative integer $r$ and $s>0$, givesAt $t=2$ or $t=-2$, one term becomes the divergent integral of $x^m$; beyond these endpoints it grows exponentially. Consequently, the full extended-valued moment-generating function isIn particular $M_f(0)=1$, checking normalization of a probability distribution. All absolute moments exist. Direct integration yieldsThe symmetric distribution makes every odd moment zero, so

Now let $X_g$ have probability density function $g_m$. Its moment-generating function would beFor every $t\ne 0$, one summand has exponential growth, which dominates the polynomial decay. HenceThere is no finite moment-generating function on any neighborhood of zero, but some moments still exist. Specifically,ThusFor even $k\ge m+1$, the moment is $+\infty$, so it is not a finite moment. For odd $k\ge m+1$, both the positive and negative contributions have infinite integral, so the expected value is undefined. Cancellation in a Cauchy principal value does not define an expected value. In particular, when $m=0$ there are no finite positive-order moments.

The support condition $-y<x<y$ forces $y>0$. As a check, the joint probability density integrates toFor any real $x$, the marginal density of $X$ isDividing the joint probability density by this marginal density gives the conditional density:Thus, conditionally on $X=x$, $Y-|x|$ has the rate-two exponential distribution.

For $U=X^2$, use the two inverse branches $x=\sqrt{u}$ and $x=-\sqrt{u}$. Each branch has absolute Jacobian determinant $1/(2\sqrt{u})$. Both lie in the original support exactly when $u>0$ and $y>\sqrt{u}$. The change of variables for probability densities therefore yieldsNo atom of a measure is created at $u=0$, since $\mathbb{P}(X=0)=0$. The value of the probability density function at this boundary may be assigned arbitrarily. Checking the integral over $0<u<y^2$ gives $4y{e}^{-2y}$, which integrates to one.

For a Borel set $A\subseteq \mathbb{R}$, the joint probability density givesFor each $x\ne 0$, substitute $z=xy$. The change of variables for probability densities contributes $dy=dz/|x|$, with the absolute value accounting for the reversed orientation when $x<0$. The line $x=0$ has two-dimensional Lebesgue measure zero, so omitting it changes no probability. By the Tonelli theorem,The bracket is nonnegative and its integral over $\mathbb{R}$ is one. This proves the probability density of a product of random variables formula:As usual, a probability density function is determined only almost everywhere; exceptional infinite values are harmless.

For the independent random variables with the specified uniform distribution, $f(x,y)=b^{-2}$ on $(0,b{)}^2$ and zero elsewhere. When $0<z<b^2$, the conditions $0<x<b$ and $0<z/x<b$ are equivalent to $z/b<x<b$. ThereforeThe value assigned at zero is arbitrary. Substitution $s=z/b^2$ checks normalization of a probability distribution: ${\int }_0^1-\log sds=1$. This is the product of independent uniform random variables law.

For $1\le j\le n$, repeated use of the floor function identity givesThus the filtration satisfies ${\mathcal{F}}_n=\sigma (X_n)$. Each positive-probability atom of a sigma-algebra here has the formwith probability $2^{-n}$; $\{X=1\}$ is a possible additional null set atom. Averaging over each dyadic interval gives a versionFor a globally defined version of the right side, assign any value, for example zero, on $\{X=1\}$.

Indeed, the displayed random variable is ${\mathcal{F}}_n$-measurable and bounded by $\sup |h|$. On each $A_{n,k}$ its integral isby the uniform distribution. Every member of ${\mathcal{F}}_n$ is a union of the displayed sets and the possible zero-probability atom. Summing these identities therefore verifies the defining conditional expectation property on all of ${\mathcal{F}}_n$.

For sigma-finite measures, Tonelli's theorem states that the double integral of a nonnegative measurable function equals either iterated integral, allowing $+\infty$. Fubini's theorem states that if a measurable function is absolutely integrable, its sections are integrable almost everywhere, both iterated integrals exist, and they equal the double integral.

Assign any value at $(0,0)$, which does not affect Lebesgue integrability. Near that corner set $s=x+y$ and $t=(x-y)/(x+y)$. Then $x=s(1+t)/2$, $y=s(1-t)/2$, and the absolute Jacobian determinant is $s/2$. On $0<s<1$, the full range is $-1<t<1$, soThe remaining part of the unit square has $x+y\ge 1$ and is harmless. The improper power integral is finite precisely when $\alpha <3$. Under the given restriction,For those values, the change $(x,y)\mapsto (y,x)$ preserves the square and changes the sign of $F$. Absolute integrability justifies this change of variables, givingFor $\alpha \ge 3$ the positive and negative parts each diverge, so symmetry does not define a Lebesgue integral.

The density is an even function on a symmetric interval, so the integrand $xf(x)$ is odd. Its integral vanishes:

Set $x=\cos \theta$ with $0\le \theta \le \pi$. The cosine addition identity giveshence $T_0=1$, $T_1=x$, and $T_{n+1}=2x{T}_n-T_{n-1}$. By mathematical induction, $T_n$ is a polynomial of degree $n$: the initial cases hold, and the leading term of $2x{T}_n$ has degree $n+1$ while the subtraction has lower degree. Its leading coefficient is $2^{n-1}$ for $n\ge 1$.

The same change of variables gives the weighted integral

The change of variables formula for a Lebesgue integrable function isThe absolute Jacobian determinant accounts for area regardless of orientation. The formula also holds for nonnegative measurable functions with possibly infinite integral.

Expanding $e^{i{S}_{\mathrm{int}}}$ gives the Feynman vertexThere is one scalar leg and a continuous oriented fermion line; no extra factor $i$ is absorbed into ${\gamma }^5$. The Dirac propagator is $S_F(k)=i(/k+M)/(k^2-M^2+i0)$, where $/k={\gamma }^{\mu }{k}_{\mu }$. Each closed fermion loop supplies a minus sign.

The tadpole diagram is proportional tobecause both numerator traces vanish by a. This is zero before the loop momentum integration.

For either orientation of the triangle Feynman diagram, let $a,b,c$ be its three internal four-momenta. Its numerator trace isAfter using cyclic property of the trace, each expanded term is a trace of ${\gamma }^5$ times at most three ordinary gamma matrices. Such traces vanish in four dimensions. For an odd number of ordinary gamma matrices, conjugation by ${\gamma }^5$ changes the sign. For two ordinary gamma matrices, the contracted Clifford algebra identity ${\gamma }_{\rho }{\gamma }^{\mu }{\gamma }^{\nu }{\gamma }^{\rho }=4{\eta }^{\mu \nu }I$ givesThe zero-matrix case is a. Hence both the one-point and three-point one-loop bosonic vertices vanish, independently of their external momenta. No cancellation between different loop orientations or delicate shift of a divergent integral is needed.

Use the total density of states $D^{\ast }=A{m}^{\ast }/(2\pi {\hslash }^2)$ and $C=A/(2\pi \hslash {)}^2$ from part (b). In the paper's double-sum convention, $f(\varphi )$ has units of energy and scales as $A^{-1}$. Define its first angular moment byThis convention matters: the Landau quasiparticle energy functional written with an explicit $1/A$ instead uses the intensive kernel $\tilde{f}=Af$.

At zero temperature, the displayed functional is the change in $E-\mu N$; a true Galilean boost preserves $N$, so $\Delta \mathcal{F}=\Delta E$. For a radial displacement $\nu (\theta )$, the first term in the Landau quasiparticle energy functional becomesThe oriented radial integral also works for negative $\nu$: removing occupied levels below the chemical potential costs positive energy. Substituting $\nu =p_0\cos \theta +O(p_0^2/p_F)$ and using $N=C\pi p_F^2$ givesThe exact second-order correction to the boundary is unnecessary here because the unperturbed excitation energy vanishes on the Fermi surface; it changes this term only beyond quadratic order. It must nevertheless be understood when identifying $\Delta \mathcal{F}$ with a fixed-particle number energy change.

For the quadratic interaction term, the first-order occupation change alone suffices. The radial Dirac delta function integrals givePut $\varphi =\theta -{\theta }^{\prime }$. The remaining angular integral usesso only the first cosine component of the interaction survives. Since $D^{\ast }=2\pi C{m}^{\ast }$,Adding both contributions and matching the bare-mass result of part (a) yields the two-dimensional Landau effective-mass relation:Equivalently, using the density of states per unit area $g^{\ast }=D^{\ast }/A$,If dimensionless Landau parameters are defined through the Fourier series $D^{\ast }f(\varphi )=F_0+{\sum }_{\ell \ge 1}{F}_{\ell }\cos (\ell \varphi )$, then $D^{\ast }{\overline{f}}_1=F_1/2$ andDefining the cosine expansion with an extra factor of two would instead absorb that factor into $F_1$; the angular-integral formula above fixes the convention unambiguously. The physical content is that interaction changes accompanying a moving Fermi sea compensate the renormalized quasiparticle velocity, so the total boost energy is governed by the bare mass. For an interaction independent of angle, ${\overline{f}}_1=0$ and this consistency condition requires $m^{\ast }=m$.

In the position representation the Time-dependent Schrödinger equation isFor a normalized state define probability density and probability current byMultiplying the equation by $\overline{\psi }$ and its complex conjugate by $\psi$, then subtracting, cancels the real potential. The resulting identity isThus $\boxed{{\partial }_t\rho +{\partial }_{x}j=0}$. Integrating over an interval gives $d{\int }_b^c\rho dx/dt=j(b)-j(c)$: probability changes only through flux across the endpoints. This is local conservation of quantum probability.

Continuity with the zero exterior imposes zero endpoint values on the infinite square well. Solving the stationary equation with these boundary conditions giveswith zero extension outside. The integral of ${\sin }^2(n\pi x/a)$ is $a/2$, which fixes the normalization. The initial coefficients are $3/5$ and $4/5$, whose squared moduli sum to one. Each energy component acquires its stationary phase:At the midpoint, ${\Psi }_1=\sqrt{2/a}$, ${\Psi }_1^{\prime }=0$, ${\Psi }_2=0$, ${\Psi }_2^{\prime }=-(2\pi /a)\sqrt{2/a}$. HenceTaking its imaginary part proves

For bound hydrogen atom states, the radial wave function is regular and finite at the origin, and decays at infinity so that $4\pi {\int }_0^{\infty }{r}^2|\Psi (r){|}^{2}dr<\infty$. Equivalently the reduced radial function $w(r)=r\Psi (r)$ has $w(0)=0$ and is square integrable on the half-line. A singular $1/r$ branch is not an admissible regular Coulomb eigenfunction; its derivative energy diverges and it introduces an origin singularity absent from the Hamiltonian domain.

Write $\kappa =e^2/(4\pi {ϵ}_0)$, so $a={\hslash }^2/(m\kappa )$. Multiplication by $r$ turns the radial equation intoSince $w=a{e}^{-\lambda R}f(R)$ and $r=aR$, substitution givesChoose $\boxed{\lambda =\mu a>0}$ to remove the constant coefficient and factor out the decaying exponential. For $f={\sum }_{k\ge 0}{a}_k{R}^{k+c}$, the lowest power has coefficient $c(c-1)a_0$, so $\boxed{c(c-1)=0}$. The $c=0$ behavior would give $\Psi \sim a_0/R$, violating regularity; choose $c=1$.

Matching successive powers then yieldsA nonterminating tail eventually has one sign, and its coefficient ratio is asymptotic to $2\lambda /k$. More precisely, for any $0<\delta <\lambda$, the ratio is eventually at least $(2\lambda -\delta )/(k+1)$. Its same-sign tail consequently dominates a positive constant times the exponential series of rate $2\lambda -\delta$, up to finitely many terms. The factor $e^{-\lambda R}$ cannot then make $\Psi$ decay. Thus square integrability forces termination.

If the last nonzero term is $a_{n-1}{R}^n$, termination requires $\lambda n=1$. ThereforeFor $n=1$, $\lambda =1$, $f$ is proportional to $R$, andThe radial volume element is $4\pi r^{2}dr$, so the mean radius is the ratio of the radius-weighted norm integral to the norm integral. With $r=aR$, the common angular and normalization factors cancel:In the ground state $\boxed{⟨r⟩=3a/2}$.

A permutation $\pi$ of identical particles multiplies a fermionic state by $\mathrm{sgn}(\pi )$; a bosonic state is invariant. Particle exchange keeps $y_{+}$ fixed and sends $y_{-}$ to $-y_{-}$. Consequently the spatial wave function is even in $y_{-}$ for the stated bosonic model and odd for the stated fermionic model. This assumes no additional spin variables in the supplied model; exchange of full particle states is the underlying rule.

The coordinate change is orthogonal, with unit Jacobian determinant, so the kinetic energy becomeson the square $|y_{+}|,|y_{-}|<a$ with zero boundary values. Product states ${\psi }_r(y_{+}){\psi }_s(y_{-})$ have energyThe one-dimensional eigenfunction is even for odd $s$ and odd for even $s$. Bosons therefore allow odd $s$, whereas fermions require even $s$. The lowest bosonic pair is $(r,s)=(1,1)$; its first excited pair is $(2,1)$. The lowest fermionic pair is $(1,2)$. ThusThe printed bosonic ground-energy claim omits the factor ${\pi }^2/4$; the value above follows directly from the supplied one-dimensional spectrum and kinetic energy.

The separation is $|x_1-x_2|=\sqrt{2}|y_{-}|$. The normalized $y_{+}$ factor integrates out. In the fermionic ground state, the relative-coordinate density is $a^{-1}{\sin }^2(\pi y_{-}/a)$, givingFor the bosonic ground state it is $a^{-1}{\cos }^2(\pi y_{-}/(2a))$, givingBoth integrals follow by replacing squared trigonometric functions by their constant and cosine parts and integrating once by parts. Henceand the difference is $\boxed{2\sqrt{2}a/{\pi }^2}$. The odd relative state vanishes at coincidence, which increases its expected separation.

Let $E=f_H(\psi )$. For any admissible variation $\eta$, differentiation with respect to a real parameter givesStationarity for all variations, including their imaginary multiples, forces $(H-E)\psi =0$; conversely this eigenvector equation makes the derivative vanish. Thus stationary Rayleigh quotients are exactly eigenvectors, with their eigenvalues as the quotient values.

Put $\kappa =e^2/(4\pi {ϵ}_0)$ and assume $\alpha >0$. The squared norm of the radial trial function is $4\pi {\int }_0^{\infty }{r}^2{e}^{-2\alpha r}dr=\pi /{\alpha }^3$. Its kinetic expectation, computed through $\int |\nabla \psi {|}^2$, is ${\hslash }^2{\alpha }^2/(2m)=\kappa a{\alpha }^2/2$. Splitting the potential integral at $R$ givesFor example, the first integral is $[1-e^{-2s}(1+2s+2s^2)]/(4{\alpha }^3)$ and the second is $e^{-2s}(R/(2\alpha )+1/(4{\alpha }^2))$, where $s=\alpha R$. ThereforeIts derivative is $\kappa [a\alpha -(1+2\alpha R)e^{-2\alpha R}]$. The actual stationary condition is consequentlyThe extra Coulomb factor in the printed stationary equation is an error: the common prefactor cancels on differentiation. The corrected equation is dimensionless on both sides. Its left side increases from zero, while its right side decreases from one, so it has one positive solution and gives the minimum in this trial family.