 Solution

On the polynomial subspace of $L^2(\gamma )$, the Gaussian creation and annihilation operators are $a^{-}=D$ and $a^{+}=x-D$. Gaussian integration by parts makes them adjoints there and gives $[a^{-},a^{+}]=I$. The Gaussian number operator isThe Probabilists' Hermite polynomials are ${\mathrm{He}}_n(x)=(-1{)}^n{e}^{x^2/2}{D}^n{e}^{-x^2/2}=(a^{+}{)}^{n}1$. Their generating function is $e^{sx-s^2/2}={\sum }_n{\mathrm{He}}_n(x)s^n/n!$, giving $D{\mathrm{He}}_n=n{\mathrm{He}}_{n-1}$. The commutator also yieldsThe orthogonality follows by repeated integration by parts; the leading coefficient is one, so these polynomials span every polynomial. By the allowed density assumption $h_n={\mathrm{He}}_n/\sqrt{n!}$ is an orthonormal basis. Define the closed number operator by $N\sum c_n{h}_n=\sum n{c}_n{h}_n$ on $\sum n^2|c_n{|}^2<\infty$. This diagonal multiplication operator is self-adjoint and positive semidefinite; polynomial truncations show it is the closure of the polynomial operator.

The Ornstein-Uhlenbeck semigroup is $P_t=e^{-tN}$, or $P_{t}f={\sum }_n{e}^{-nt}{c}_n{h}_n$. It also has the Mehler formula for the Ornstein-Uhlenbeck semigroupTo verify the formula, apply its right side to $e^{sx-s^2/2}$: the Gaussian moment-generating function makes the result $e^{s{e}^{-t}x-s^2{e}^{-2t}/2}$, proving the eigenvalue identity on every Hermite polynomial. Positivity and invariance of Gaussian measure give $L^2$ contraction by Jensen inequality, so polynomial density extends the equality to all $L^2(\gamma )$. The coefficient expansion givesFor continuously differentiable $f$ with bounded derivative, differentiate the Mehler expectation by the dominated convergence theorem to get the Ornstein-Uhlenbeck gradient commutation identityThe same identity extends to the Gaussian Sobolev form domain. Mehler's formula also gives pointwise convergence to the Gaussian mean for such $f$, since bounded derivative permits at most linear growth.

The Gaussian Dirichlet energy is the energy form of $N$:For polynomials this follows from $N=(a^{-}{)}^{\ast }{a}^{-}$, and closure extends it to the form domain. For the given continuously differentiable $f$, Gaussian integration by parts yields $⟨f^{\prime },h_{n-1}⟩=\sqrt{n}{c}_n$, so Parseval identity proves the derivative-energy equality directly. Bounded derivative puts $f$ in the Gaussian Sobolev space, hence in that form domain; it need not be in the full operator domain, so $⟨f,Nf⟩$ is not always a legitimate initial definition.

Use the entropy functional ${\mathrm{Ent}}_{\gamma }(h)=\int h\log hd\gamma -(\int hd\gamma )\log \int hd\gamma$. Begin with bounded smooth $h>0$ bounded away from zero, and put $h_t=P_{t}h$. Invariance and Gaussian integration by parts give the Ornstein-Uhlenbeck entropy dissipation identityThe gradient identity and the permitted weighted Cauchy-Schwarz inequality, applied with weight $h$, implyThe entropy tends to zero as $t\to \infty$ by Mehler's formula and bounded convergence. Integrating in time and using invariance yields ${\mathrm{Ent}}_{\gamma }(h)\le \frac{1}{2}\int |h^{\prime }{|}^2/hd\gamma$. Set $h=f^2$ to obtain the sharp Gaussian logarithmic Sobolev inequalityFor the stated possibly unbounded $f$, apply the argument to smooth positive clipped approximations with a common positive lower bound and uniformly bounded derivatives. They can converge pointwise together with their derivatives and have a common linear-growth bound. Since Gaussian measure integrates every polynomial moment, dominated convergence theorem passes both the entropy and the derivative energy to the limit. This also justifies use of the supplied inequality when $2f{f}^{\prime }$ itself is unbounded. The normalization $\Vert f{\Vert }_1=1$ concerns $f$, and does not imply $\int f^{2}d\gamma =1$; the second term in the entropy definition must be retained.