 Solution

A normalization is important here. I use the covolume-one theta function of a fractional ideal, for which the requested limit is $1$, and also give the formula for the unnormalized Gaussian theta sum.

Let $n=[K:\mathbb{Q}]$, let $d_v=[K_v:\mathbb{R}]\in \{1,2\}$, and choose one field embedding ${\sigma }_v$ for each Archimedean place. The Minkowski embedding of a number field identifieswith a real inner-product space havingIts Euclidean Lebesgue measure is $d\mu ={\prod }_{v\text{ real}}d{x}_v{\prod }_{v\text{ complex}}2d\mathrm{Re}{x}_{v}d\mathrm{Im}{x}_v$. With this convention, the covolume of a fractional ideal lattice $j(\mathfrak{b})$ iswhere $D_K$ is the field discriminant and $N(\mathfrak{b})$ is the positive absolute norm of a fractional ideal. Thus the Euclidean lattice ${\Lambda }_{\mathfrak{b}}=C_{\mathfrak{b}}^{-1/n}j(\mathfrak{b})$ has covolume $1$. DefineThis Gaussian theta sum converges absolutely whenever every $y_v>0$.

The trace dual of a fractional ideal isHere ${\mathfrak{D}}_{K/\mathbb{Q}}^{-1}$ is the inverse different. Since $N({\mathfrak{D}}_{K/\mathbb{Q}})=|D_K|$, we have $C_{{\mathfrak{b}}^{\vee }}=C_{\mathfrak{b}}^{-1}$. There is a subtle distinction between the trace pairing and the positive inner product: the dual lattice of $j(\mathfrak{b})$ is $\overline{j({\mathfrak{b}}^{\vee })}$, where the bar conjugates the complex coordinates and fixes the real ones. Indeed,Consequently ${\Lambda }_{\mathfrak{b}}^{\ast }=\overline{{\Lambda }_{{\mathfrak{b}}^{\vee }}}$. Coordinatewise complex conjugation preserves the weighted squared lengths in the Gaussian theta sum.

Here are the precise analytic formulas used in the proof. For a Schwartz function on $K_{\infty }$, take the Fourier transform to beFor a full Euclidean lattice $\Lambda$ of covolume $C$, the Poisson summation formula for a Euclidean lattice isFor $f_y(x)=\exp (-\pi {\sum }_v{d}_v{y}_v|x_v{|}^2)$, the Gaussian Fourier transform, applied in orthonormal real coordinates, givesMore generally, $\widehat{\exp (-\pi x^{T}Ax)}(z)=(\det A{)}^{-1/2}\exp (-\pi z^T{A}^{-1}z)$ for a real positive-definite matrix $A$ that is symmetric. Each complex coordinate contributes two real coordinates, which explains its exponent $d_v=2$.

Apply the Poisson summation formula for a Euclidean lattice to ${\Lambda }_{\mathfrak{b}}$, whose covolume is $1$, and use the preceding description of its dual lattice. The anisotropic theta functional equation isFor comparison, the unnormalized Gaussian theta sumhas the functional equationand its corresponding limit is $C_{\mathfrak{b}}^{-1}$ rather than $1$. Explicitly, our normalization is $\Theta (y,\mathfrak{b})=\theta (C_{\mathfrak{b}}^{-2/n}y,\mathfrak{b})$.

For the small-parameter asymptotic of a lattice theta sum, $y\to 0$ means that every coordinate tends to zero. In $\Theta (y^{-1},{\mathfrak{b}}^{\vee })$, the zero lattice vector contributes $1$. Every nonzero vector contributes a term tending to zero. Once every $y_v\le 1$, all these terms are bounded by the summable Gaussian theta sum ${\sum }_{\lambda \in {\Lambda }_{{\mathfrak{b}}^{\vee }}}{e}^{-\pi |\lambda {|}^2}$. The dominated convergence theorem therefore gives $\Theta (y^{-1},{\mathfrak{b}}^{\vee })\to 1$. Using the anisotropic theta functional equation,The condition that every coordinate tends to zero matters; $\Vert y\Vert \to 0$ alone would not justify this argument.