\section{Introduction}
\label{intro}
\indent The Falicov-Kimball model \cite{originalFK} is perhaps the simplest solid-state model for describing strongly correlated electron systems. The model possesses a Mott-insulator transition, a charge-density-wave ordered phase and can be solved exactly in the limit of infinite dimensions \cite{Brandt}. Initially the model was employed to understand phase transitions in transition-metal oxides~\cite{originalFK}, and has since been used to study a number of strongly correlated systems. Its solution via dynamical mean-field theory has also been reviewed~\cite{jimFK}.

Recently, the Falicov-Kimball model has been used to investigate core-level X-ray photoemission spectroscopy \cite{jimXPS1} \cite{jimXPS2}. It has been known for some time that the localized $f$-electron Green's function in the metallic phase of this model is related to the x-ray edge singularity problem \cite{XRAYEDGE1,XRAYEDGE2} (we only reference papers immediately relevant to this work here; a more complete history of X-ray edge and the Falicov-Kimball model appears elsewhere~\cite{jimXPS2}). In particular, at low temperature, the density of states in the metal displays a power-law divergence with an interaction-dependent power law.

The $f$-electrons, unlike the conduction electrons, possess a non-trivial and highly temperature dependent spectral function. This has been studied previously using the numerical renormalization group \cite{NRG} illustrating the power law divergence of the $f$-electron spectrum at $\omega=0$ in the metallic phase. Here we seek to understand the gapped behavior of the $f$-electron spectral function in the Mott-insulating regime, where NRG approaches are known to have accuracy issues. Particularly, we ask if the $f$-electron gap approaches the (temperature-independent) conduction-electron gap in the limit $T\to 0$ in the insulating phase. This question may seem like it should have an obvious answer, because the $f$-electron dynamics are inherited through their interaction with the conduction electrons. But at nonzero temperature, we clearly see $f$-electron spectral weight within the gap, so the answer to the question is far from obvious.

\section{Formalism}
The spinless Falicov-Kimball model describes two electron species, conduction ($c$) and localized ($f$) electrons, which interact via a local Couloumb repulsion $U$ when occupying the same lattice site $i$. We represent the creation (destruction) of a conduction electron at the site $i$ by the second quantized operator $c^{\dagger}_i$ ($c_i^{\phantom\dagger}$) and of a localized electron by $f^{\dagger}_i$ ($f_i^{\phantom\dagger}$). Assuming a common chemical potential $\mu$ for the two species, the spinless Falicov-Kimball Hamiltonian is given by
\begin{equation}
    \mathcal{H}_{FK} = \sum_{ij}\left (-t_{ij}-\mu\delta_{ij}\right )c^{\dagger}_ic_j^{\phantom\dagger} +E_f \sum_i f^{\dagger}_i f_i^{\phantom\dagger} + U\sum_if^{\dagger}_if_i^{\phantom\dagger}c^{\dagger}_ic_i^{\phantom\dagger}
\end{equation}
where $E_f$ is the $f$-electron site energy, and $t_{ij}$ is the hopping matrix, with $t_{ij}=t=t_{ji}$ for nearest-neighbor sites $i$ and $j$ with the hopping integral $t$ a real constant. We work at half filling ($\mu = U/2$ and $E_f=-U/2$) on the infinite-coordination Bethe lattice in equilibrium.\\

Scaling the hopping energy $t$ with coordination number $Z$ as $t = t^{*}/\sqrt{Z}$ then in the limit $Z \to \infty$ we obtain the noninteracting density of states on the Bethe lattice, given by
\begin{equation}
     \rho_{Z\to\infty}(\epsilon) = \frac{\sqrt{4t^{*2}-\epsilon^2}}{2\pi t^{*2}}
\end{equation}
with $t^{*} = 1$ our energy unit and we restrict to $|\epsilon|\le 2t^*$.

The dynamical mean-field theory (DMFT) approach solves the many-body problem by mapping the lattice problem onto an impurity problem because the self-energy has no momentum dependence. A self-consistent iterative approach is employed to determine the Green's functions, with details given in~\cite{jimFK}. We do not need the full conduction-electron solution here, instead,
we require only the bare time-ordered Green's function for the conduction electrons $G_0(t)$ given by
\begin{equation}
    G_0(t) =-\frac{i}{\pi}\int_{-\infty}^{\infty}d\omega\ \text{Im}\left [\mathcal{G}_0(\omega)\right ] e^{-i\omega t}[f(\omega)-\theta(t)]
\end{equation}
with $f(\omega)$ the Fermi-Dirac distribution,
\begin{equation}
    \mathcal{G}_{0}(\omega) = \frac{1}{\omega+i\delta +\mu- \lambda(\omega+i\delta)}
\end{equation}
the effective medium, and the dynamical mean-field $\lambda(\omega)$ obtained from the DMFT algorithm
\cite{jimFK}. Here $\theta(t)$ is the Heaviside unit-step function.

The conduction electrons behave like noninteracting particles, in that they have a temperature-independent density of states, but also like interacting electrons, since they have a Mott metal-insulator transition at $U=2$.
\begin{figure}[htb]
\centerline{
  \includegraphics[width=1.5in]{kbk_contour.png}}
\caption{The Kadanoff-Baym-Keldysh contour, which starts at time 0, moves along the real time axis to time $t$, back along the real axis to time 0, then proceeds down the imaginary time axis a distance $-\beta=-1/T$.}
\label{fig:kbk}
\end{figure}

We define the greater and lesser $f$-electron Green's functions by
\begin{eqnarray}
    G_f^>(t,t') &=& -i \langle f(t)f^{\dagger}(t') \rangle,\\
G_f^<(t,t') &=&i \langle f^{\dagger}(t') f(t)\rangle,
\end{eqnarray}
 where the angle brackets denote the thermal average, $\langle \ldots\rangle = {\text{Tr}(\exp(-\beta \mathcal{H}_{FK})(...))}/{\text{Tr}\exp(-\beta{\mathcal H}_{FK})} $ The creation and annihilation operators are in the Heisenberg representation. These Green's functions can be determined by selecting $t$, $t'$ on certain branches of the Kadanoff-Baym-Keldysh contour (Fig. \ref{fig:kbk}) and using the contour-ordered Green's function,
 \begin{equation}
     G_f^c(t_c,t_c') = -i \langle {\mathcal T}_{c}\big(f(t_c)f^{\dagger}(t_c')\big)\rangle
 \end{equation}
 where $t_c$, $t_c'$ are two times on the contour (Fig. \ref{fig:kbk}). ${\mathcal T}_c$ is the contour time-ordering operator, arranging the operators such that $t_c'$ lies before $t_c$ on the contour. For example, if we pick $t_c'$ on the upper real-time branch of the contour and $t_c$ on the lower real-time branch, we recover $G_f^>(t,t')$ from the contour-ordered Green's function.
 \\

It has been shown in Ref.~\cite{main1}, that in equilibrium the greater Green's functions for the $f$-electrons take the form of a Toeplitz determinant of a continuous matrix operator over \emph{only the positive time branch of the
contour}
\begin{equation}\label{det}
    G_f^>(t) = -iw_0e^{-i(E_f-\mu)t}\text{Det}_{[0,t]}\big|\delta(t_1-t_2)-UG_0(t_1-t_2)\big|
\end{equation}
where $w_0$ is the average density of sites without an $f$-electron ($w_0 = \frac{1}{2}$ at half filling), and $G_0(t)$ is the bare time-ordered Green's function determined from the dynamical mean-field $\lambda(\omega)$. The symbols $t_1$ and $t_2$ denote the matrix indices of the continuous matrix operator for which we evaluate the determinant; note that both times must fall within the interval $[0,t]$. (There is a similar expression for the lesser Green's function.) To approximate this continuous matrix operator, we discretize it to a conventional matrix and calculate the determinant for three different discretization time steps $\Delta t$. We then perform a second-order Lagrange interpolation to extrapolate to the limit $\Delta t\to 0$. These numerical results are checked for accuracy against known spectral moments of the Green's functions \cite{SpecMom}.

In the limit of large times, an exact analytic formula for the $f$-electron Green's function in equilibrium can be obtained using a factorization technique from complex analysis described by McCoy and Wu \cite{Mccoy} and called the Wiener-Hopf sum approach. It relies on a result for infinite-sized determinants of Toeplitz matrices called Szego's theorem. Further finite-time approximations can be made to improve the short-time agreement of this asymptotically exact result with the determinant calculation given by Eq.~(\ref{det}) \cite{main1}, with slightly different formulas for the case when the interaction energy $U$ lies above or below  $U_{c}$, a critical interaction strength $U_{c} = \sqrt{2}$ on the Bethe lattice where the complex function $C(\omega) = 1-UG_0(\omega)$ goes from no winding around the origin for $U<U_{c}$ to winding once in the clockwise direction (Ind$\,C$=-1) for $U\geq U_{c}$; note that this critical $U_c$ is smaller than the critical $U$ for the Mott transition. In the regime with $U>U_c$, the analytic result for the Toeplitz determinant in Eq.~(\ref{det}) is
\begin{eqnarray}
\label{WH}
    {\rm Det}_{[0,t]} &=& \exp \bigg[\frac{t}{2\pi}\int_{-\infty}^{\infty}\,d\omega\ \ln{\bar{C}(\omega)} + \int_0^{\infty}\,dt'\ t'\bar{g}(t')\bar{g}(-t')\bigg]\nonumber\\
    &\times&\frac{\Delta t}{2\pi}\int_{-\frac{\pi}{\Delta t}}^{\frac{\pi}{\Delta t}} \,d\omega'\ e^{i\omega' t}\frac{\bar{P}(-\omega')}{\bar{Q}(\omega')}
\end{eqnarray}
where $\bar{C}(\omega) = \exp[i\omega\Delta t][1-UG_0(\omega)]$, $\,\bar{g}(t) = \int_{-\infty}^{\infty}\,d\omega$\, \\ $\times\exp[-i\omega t]\ln(\bar{C}(\omega))/2\pi$, and $\bar{P}$   $(\bar{Q})$ are integrals over all positive (negative) time of $\bar{g}(t)$ and satisfy
\begin{equation}
\bar{C}(\omega) = \frac{1}{\bar{P}(\omega)\bar{Q}(-\omega)}.
\end{equation}
This expression represents a significant reduction in computational complexity, and allows us to probe a much wider parameter space when directly calculating the discretized determinant is not possible. Our approach will be to calculate the determinant directly for short times and use the analytic expression for long times, allowing us to obtain the spectral function of the $f$-electrons down to temperatures significantly lower than previous calculations. We need to patch the two solutions together smoothly, as described below. We examine the behavior of the $f$-electrons as they approach their $T=0$ limit in both the metallic phase near the Mott-insulator transition $U=2$ and in the Mott insulating regime.

Finally, we define the local density of states of the $f$-electrons, $A_f(\omega)$. At half-filling, there exists a particle-hole symmetry in our system [$A_f(\omega) = A_f(-\omega)$], and consequently the full $f$-electron density of states can be expressed as a Fourier transform of Im $G^>(t)$ alone \cite{Brandt1992},
\begin{equation}
        A_f(\omega) = -\frac{2}{\pi}\int_0^{\infty} \,dt\ \cos(\omega t)\text{Im }G^>(t).
 \end{equation}

\label{sec:1}

\section{Numerics}

We examine the temperature-dependent dynamics of the $f$-electron spectral function in the metallic regime near the Mott transition as well as in weakly and strongly correlated insulating phases above the Mott transition, which occurs at $U=2$. The spectral function is known to have a power-law like divergence in the metallic phase that disappears as we move into the Mott phase \cite{NRG}. In the Mott insulator, the $f$-electron density of states develops a gap with decreasing temperature \cite{main2}. We explore the character of this gap to see if it is the same as in the temperature-independent conduction electron density of states.

To utilize the Weiner-Hopf technique in Eq.~(\ref{WH}), we perform the discretized matrix calculation of Eq.~(\ref{det}) out to the longest time computationally feasible for three different time steps $\Delta t$ in the ratio $1:2:4$. Next we employ a second-order Lagrange extrapolation to take the limit $\Delta t\to 0$. We use Eq.~(\ref{WH}) to calculate the determinant out to even longer times, and use a weighted blending of the two functions over a time range where the analytic approximation is roughly parallel to the discretized determinant results. This procedure works exceptionally well for high temperature, as shown in Fig.~\ref{fig:blend1}, but as we decrease the temperature, the determinant calculation requires significantly smaller time steps and takes longer to reach the "long-time" regime where our analytic result holds (see Fig. \ref{fig:blend2}). This effect limits our low-temperature calculations, and along with truncating the calculation at a finite time, leads to numerical artifacts near $\omega = 0$ (see Figs. \ref{fig:2DOS}, \ref{fig:25DOS}, \ref{fig:3DOS}).
\begin{figure}[htb]
\centering
  {\includegraphics[width=0.49\textwidth]{u15gf_im_new.png}}

\caption{$G^>(t)$ vs.~$t$ for $U=1.5$. Here we have no zero crossing of the time axis, corresponding to the metallic phase of the model. The inset shows the shorter-time behavior.}
\label{fig:15time}
\end{figure}

\begin{figure}[htb]
\centering{
  \includegraphics[width=0.49\textwidth]{u15dos.png}}

\caption{Spectral function $A_f(\omega)$ vs.~$\omega$. Here we see the well documented \cite{NRG} power law divergence of the orthogonality catastrophe set in as $T\to 0$.}
\label{fig:15w}
\end{figure}

\begin{figure}[htb]
\centering{
  \includegraphics[width=0.49\textwidth]{blend_examp.png}}
\caption{Time-domain plot of Im~$G^>(t)$ for high temperature ($T=1$) in the insulating phase $U\geq 2$. Here the analytic formula works well for all times, and the blending between the two approaches is smooth.}
\label{fig:blend1}
\end{figure}

\begin{figure}[htb]
\centering{
  \includegraphics[width=0.5\textwidth]{bad_blend_examp.png}}
\caption{Time-domain plot of Im~$G^>(t)$ for low temperature $T<0.05$ in the insulating phase $U\geq 2$. Here the direct determinant calculation has not reached its asymptotic limit at the maximum time allowed by our computational resources. Note that both Green's functions oscillate with the same frequency, but these oscillations are damped much more quickly when using Eq.~(\ref{WH}). The appearance of high-frequency oscillations at lower temperatures are characteristic of the $f$-electron's low-temperature dynamics.}
\label{fig:blend2}
\end{figure}

\section{Results}

In the metallic regime, for $U=1.5$ (Fig. \ref{fig:15w}) we see evidence of a power-law divergence as the temperature approaches zero in agreement with \cite{NRG}. In the time domain, we observe a delayed decay towards zero with decreasing temperature. Note there is no zero crossing in Im~$G^>(t)$ (Fig. \ref{fig:15time}), indicating we are still in the metallic phase.

\begin{figure}[htb]
    \centering{
    \includegraphics[width=0.5\textwidth]{u19dos_new.png}}
    \caption{$A_f(\omega)$ vs.~$\omega$ for $U=1.9$. We see power-law-like behavior, with a central kink that sharpens with decreasing temperature.}
    \label{fig:19}
\end{figure}

For $U=1.9$, closer to the Mott transition, we still see the power-law like increase in the spectral function (Fig. \ref{fig:19}), but with the development of a kink in the center of our density of states, representing a precursor to the insulating Mott phase.

\begin{figure}[htb]
\centering{
    \includegraphics[width=0.5\textwidth]{u20dos_new.png}}
    \caption{$A_f(\omega)$ vs $\omega$ at the Mott transition, $U=2$. We see a gradual evolution whereby spectral weight is transferred from higher frequency states to states near $\omega=0$.}
    \label{fig:2DOS}
\end{figure}

At the Mott transition, $U=2$ (Fig. \ref{fig:2DOS}), the $f$-electrons partially fill the insulating gap at high temperature, with a kinked density of states that does not reach zero at the minimum temperature we were able to reach ($T=0.001$). We expect that the density of states should touch zero precisely at $T=0$. This is in contrast with the conduction electrons, whose density of states touches zero at $U=2$ for all temperatures.

\begin{figure}[htb]
\centering{
    \includegraphics[width=0.49\textwidth]{u25dos_new.png}}
    \caption{Temperature evolution of $A_f(\omega)$ vs.~$\omega$ at $U=2.5$ and (inset) the conduction electron density of states over two different ranges of frequency. We see that the gap in the lowest temperature $f$-electron density of states seems to be approaching the gap width of the conduction electron density of states.}
    \label{fig:25DOS}
\end{figure}

For $U=2.5$ (Fig. \ref{fig:2DOS}), we see a similar transfer of weight from high to low frequencies as described in \cite{main2}. Here we see the spectral function becomes gapped around $T=0.1$, but continues to change shape down to $T=0.001$, below which the gap is frozen into place. Notice in the right inset of Fig.~\ref{fig:25DOS} the small region where the density of states becomes slightly negative. This is an artifact of not properly capturing the long-time behavior, likely due to the finite time truncation of some long period oscillations in $G^>(t)$ and possibly the blending of the numerical and asymptotic results.

\begin{figure}[htb]
\centering{
    \includegraphics[width=0.49\textwidth]{u30dos_new.png}}
    \caption{Temperature evolution of $A_f(\omega)$ vs.~$\omega$ at $U=3$; (inset) the conduction-electron density of states over two different frequency ranges. We again see that the lowest temperature $f$-electron density of states seems to be approaching the width of the conduction electron density of states.}
    \label{fig:3DOS}
\end{figure}

Finally we examine the strongly correlated insulator with $U=3$ (Fig. \ref{fig:3DOS}). We see a similar temperature evolution to $U=2$ and $U=2.5$. In this case, the system becomes gapped around $T=0.25$, higher  than $U=2.5$. Below this temperature $A_f(\omega)$ changes more slowly as the low-temperature behavior is frozen in. We observe a general trend of the gap appearing and freezing into place at a higher temperature for more strongly correlated materials. Again we see the irregular behavior near $\omega = 0$ from the finite time truncation. More interesting, however, is the comparison with the conduction-electron density of states. We see from the left inset of Fig.~\ref{fig:3DOS} that at low temperature the width of the gap is nearly indistinguishable between the two different electrons---this strongly suggests that the $T=0$ gap of the $f$-electrons is the same as the conduction-electron gap.

\section{Conclusion}

We extensively studied the properties of the $f$-electron spectra of the Falicov-Kimball model in a number of interesting cases: just below the Mott transition, in the strongly correlated regime, and at temperatures approaching the $T=0$ limit. By using the Weiner-Hopf technique, we were able to examine these more complex situations by obtaining an analytic expression for the long-time behavior which is computationally much more efficient than the direct determinant calculation, because the matrices grow in size with increasing time. Due to the rich temperature-dependent dynamics of the $f$-electrons, we asked if the spectral gap in the $f$-spectrum approaches that of the (temperature-independent) conduction electrons at $T=0$. We found strong evidence that this is the case. This technique could be pushed further by making higher-order finite-time corrections to the asymptotic formulas given in \cite{main1} and by carrying the calculations out to longer times.

\begin{acknowledgements}
We would like to thank A.~M.~Shvaika for valuable conversations and suggestions. This work was supported by
the Department of Energy, Office of Basic Energy Sciences, Division of Materials Sciences and Engineering under Contract No.~DE-SC0019126.
J.~K.~F. was also supported by the McDevitt bequest at Georgetown.
\end{acknowledgements}

\section*{Conflict of interest}
 The authors declare that they have no conflict of interest.

\bibliographystyle{spphys}


\section{introduction}
Let $(A,\mathfrak{m} )$ be a Noetherian local ring of dimension $d$. If $M$ is an $A$-module then let $\ell(M)$ denote its length and $\mu(M)$ the number of its minimal generators. The function $H^{(1)}(A, n) = \ell(A/\mathfrak{m} ^{n+1})$ is called the \emph{Hilbert-Samuel function} of $A$ (with respect to $\mathfrak{m} $). It is well known that there exists a polynomial $P_A(z) \in \mathbb{Q}[z]$ of degree $d$ such that $H^{(1)}(A, n) = P_A(n)$ for all $n \gg 0$.
 The polynomial $P_A(z)$ is called the \emph{Hilbert-Samuel polynomial} of $A$.  We write
$$ P_A(z) = \sum_{i = 0}^{d}(-1)^ie_i(A)\binom{z+ d -i}{d-i}.$$
The integers $e_i(A)$ are called the \emph{$i^{th}$-Hilbert coefficient} of $A$. The zeroth Hilbert coefficient $e_0(A)$ is called the \emph{multiplicity} of $A$.
We set $e_i = e_i(A)$ for all $i$ and $e = e_0$.

The graded ring $G(A) = \bigoplus_{n \geq 0}\mathfrak{m} ^n/\mathfrak{m} ^{n+1}$ is called the associated graded ring of $A$(with respect to $\mathfrak{m} $). It's Hilbert series
$$H_A(z) = \bigoplus_{n \geq 0} \ell(\mathfrak{m} ^n/\mathfrak{m} ^{n+1})z^n  = \frac{h_A(z)}{(1-z)^d} \quad  \text{where} \ h_A(z) \in \mathbb{Z}[z].$$
If $f$ is a polynomial we use $f^{(i)}$ to denote the $i$-th formal derivative of $f$. It is
easy to see that $e_{i} = h_A(z)^{(i)}(1)/i! $ for $0 \leq i \leq d$. It is also
convenient to set
$e_{i} = h_A(z)^{(i)}(1)/i! $ for all $i\geq 0$.  Let $h = \mu(\mathfrak{m} ) - d$ be the embedding co-dimension of $A$. It is easily seen  that $h_A(z) = 1 + hz + \text{higher terms}$. We call $h_A(z)$ the \emph{h-polynomial} of $A$. The function $H(A,n) = \ell(\mathfrak{m} ^n/\mathfrak{m} ^{n+1})$ is called the \emph{Hilbert function}
of $A$ (with respect to $\mathfrak{m} $).

Now assume that $A$ is Cohen-Macaulay. Then Abhyankar proved that $e \geq h + 1$, see \cite{A}. Northcott proved that $e_1 \geq e_0 - 1$, see \cite{N}.
Itoh proved that $e_2 \geq e_1 - e + 1$, see \cite[Theorem 12]{It}.
It was observed by Sally that when the Hilbert-coefficients satisfy border values then it forces $G(A)$ to have high depth and also forces the $h_A(z)$
to have a prescribed form; see \cite{S1}, \cite{S2}, \cite{S3} and \cite{S4}. See the nice survey article \cite{V} for an introduction to this area of research. In particular see \cite{V} for classification of Hilbert functions when $e = h+1, h+2$ and $e_1 = e-1, e, e + 1$. Let $\operatorname{type}(A) = \ell(\operatorname{Ext}^d_A(k, A))$ be the Cohen-Macaulay-type of $A$.
For classification of Hilbert functions when $e = h + 3$ and $\operatorname{type}(A) < h$, see \cite{RV}.

In this paper we describe all Hilbert functions that can possibly occur if $e_1 = e + 2$. It turns out that in the process we also have to describe
Hilbert functions that can occur when $e_2 = e_1 - e + 1 = 3$. It was conjectured by Valla that if $e_2 = e_1 - e + 1$ then $G(A)$ is Cohen-Macaulay \ and
$\deg h_A(z) \leq 2$.  This conjecture is true when $\dim A \leq 1$.
 A counter-example  to this conjecture was found by Wang in the case $e_2 = e_1 - e + 1 = 3$  and $\dim A = 2$ ; see \cite[3.10]{CPR}.
We prove
\begin{theorem}\label{e2}
Let $(A,\mathfrak{m} )$ be a Cohen-Macaulay \ local ring  of dimension $d$. Set $G = G(A)$. If $e_2 = e_1 - e + 1 = 3$ then one of the following cases occur
\begin{enumerate}[\rm (i)]
\item
$e = h + 4$, $G$ is Cohen-Macaulay \  and $h_A(z) = 1 + hz + 3z^2$.
\item
$e = h + 3$, $d \geq 2$, $\operatorname{depth} G = d - 2$, $\operatorname{type}(A) = h$ and

$\displaystyle{h_A(z) = 1 + hz + 3z^3 - z^4.}$
\end{enumerate}
\end{theorem}
 Note that in Theorem \ref{e2} the $h$-polynomial of $A$ completely determines \\ $\operatorname{depth} G(A)$. In section 7 we give examples of Cohen-Macaulay \ local rings
having Hilbert functions as described in \ref{e2}.

As a consequence of Theorem \ref{e2} we can completely classify Hilbert functions that can occur when $e_1 = e + 2$.
\begin{theorem}\label{e1}
Let $(A,\mathfrak{m} )$ be a Cohen-Macaulay \ local ring  of dimension $d$. Set $G = G(A)$. If $e_1 = e + 2$ then one of the following cases occur
\begin{enumerate}[\rm (i)]
\item
$e = h + 2$, $\operatorname{type}(A) = h$, $\operatorname{depth} G = d - 1$ and $h_A(z) = 1 + hz + z^4$.
\item
$e = h + 4$, $G$ is Cohen-Macaulay \  and $h_A(z) = 1 + hz + 3z^2$.
\item
$e = h + 3$, $e_2 = 4$, $\operatorname{depth} G \geq d - 1$ and $h_A(z) = 1 + hz + z^2 + z^3$.
\item
$e = h + 3$, $e_2 = 3$, $d \geq 2$, $\operatorname{depth} G = d - 2$, $\operatorname{type}(A) = h$ and

$\displaystyle{h_A(z) =  1 + hz + 3z^3 - z^4.}$
\end{enumerate}
\end{theorem}
In section 7 we give examples of Cohen-Macaulay \ local rings
having Hilbert functions as described in \ref{e1}.

We now describe in brief the contents of this paper. In section two we describe some preliminary results that we need. In section three we discuss the case of Theorem \ref{e2} when $d = 2$. The most difficult case of Theorem \ref{e2} is the case when $d = 3$. This is done in section four.
In section five we prove Theorem \ref{e2}. We prove Theorem \ref{e1} in section six. Finally in section seven we give examples which illustrate our results.
\section{preliminaries}
 In this section we discuss some preliminaries that we need.
In this paper all rings are Noetherian and all modules are assumed to be finitely generated.
Let $(A,\mathfrak{m} )$ be
 a local ring of dimension $d$ with residue field $k = A/\mathfrak{m} $.

\s If $a$ is a non-zero
 element of $A$ and if $j$ is the largest integer such that $a \in \mathfrak{m} ^{j}$,
then we let $a^*$ denote the image of $a$
 in $\mathfrak{m} ^{j}/\mathfrak{m} ^{j+1}$.

\s An element $x \in \mathfrak{m} $ is said to be  \emph{$A$-superficial} if there exists $c$ and $n_0$ such that
$(\mathfrak{m} ^{n+1} \colon x)\cap \mathfrak{m} ^c  = \mathfrak{m} ^n$ for all $n \geq n_0$. Superficial elements exist when the residue field of $A$ is infinite.
If $\operatorname{depth} A > 0$ and $x$ is $A$-superficial then $x$ is $A$-regular and furthermore $(\mathfrak{m} ^{n+1} \colon x) = \mathfrak{m} ^n$ for all $n \gg 0$.

 \begin{remark}
\label{Pr0}
 If the residue field of $A$ is finite then we resort to the
standard trick to replace $A$ by $A' = A[X]_S$ where $S =  A[X]\setminus \mathfrak{m}  A[X]$. The residue
field of $A'$ is $k(X)$, the field of rational functions over $k$.
Furthermore
\begin{equation*}
H(A',n) = H(A,n) \ \  \forall n \geq 0 \ \  \text{and} \  \   \operatorname{depth} G(A') = \operatorname{depth} G(A).
 \end{equation*}
\end{remark}

\s \label{exist-sup} Assume the residue field of $A$ is infinite. We need to recall the construction of superficial element as it is used in the proof of  Theorem \ref{depth0}.  Set $G = G(A)$ and let
$\mathfrak{M} $ be its maximal homogeneous ideal. Set $V$ to be the $k$-vector space $\mathfrak{m} /\mathfrak{m} ^2$. If $P$ is a prime ideal of $G$ and $P \neq \mathfrak{M} $ then note $P \cap V \neq V$. Then as $k$ is infinite we get
\[
V \neq \bigcup_{P \in \operatorname{Ass}{G}, P\neq \mathfrak{M} } P \cap V.
\]
If $x \in \mathfrak{m} $ is  such that
\[
x^* \in V \setminus \bigcup_{P \in \operatorname{Ass}{G}, P\neq \mathfrak{M} } P \cap V.
\]
then $x$ is $A$-superficial.

\s \label{mod-sup}Let $x \in \mathfrak{m} $ be $A$-superficial and $A$-regular. Set $B = A/(x)$. Also set $b_n = \ell((\mathfrak{m} ^{n+1} \colon x)/\mathfrak{m} ^n)$ and
$b_x(z) = \sum_{n \geq 0}b_n z^n$. Then $b_x(z) \in \mathbb{Z} [z]$. Furthermore we have
\[
h_A(z) = h_B(z) - (1-z)^db_x(z).
\]
So $e_i(A) = e_i(B)$ for $0 \leq i \leq d - 1$. Also $x^*$ is $G(A)$-regular if and only if $b_x(z) = 0$. In this case we have $G(B) = G(A)/(x^*)$.

\s A
sequence $x_1,x_2,\ldots,x_r$ in a local ring $(A,\mathfrak{m} )$ is said to
be an $A$-\emph{superficial sequence}  if $x_1$ is $A$-superficial
 and $x_i$ is  $A/(x_1,\ldots,x_{i-1})$-superficial for
$2\leq i \leq r$.

\s\label{SD} Assume $A$ is Cohen-Macaulay \ and $x_1,\ldots, x_r$ is an $A$-superficial sequence with $r \leq d$. Set $B = A/(x_1,\ldots, x_r)$. Then
\begin{enumerate}
\item
$\operatorname{depth} G(A) \geq r $ if and only if $x_1^*,\ldots, x_r^*$ is $G(A)$-regular. In this case we have $G(B) = G(A)/(x_1^*,\ldots, x_r^*)$.
\item
(\emph{Sally descent:}) $\operatorname{depth} G(A) \geq r + 1$ if and only if $\operatorname{depth} G(B) \geq 1$.
\end{enumerate}

\s \label{SA} Let $(A,\mathfrak{m} )$ be Cohen-Macaulay \ and suppose $x_1,\ldots, x_d$ be a maximal $A$-superficial sequence. Set $J = (x_1,\ldots, x_d)$.
Then $e = h+ 1+ \ell(\mathfrak{m} ^2/J\mathfrak{m} )$.

\s \label{dim1}
 Let  $\dim A = 1$  and $A$ is Cohen-Macaulay. Let $x$ be $A$-superficial.
 Set $\rho_n = \lambda(\mathfrak{m} ^{n+1}/x\mathfrak{m} ^n)$.   We have
   \begin{enumerate} [\rm (1)]
    \item
    $ H(A,n) = e - \rho_n$.
    \item
    If $\deg h_A(z) = s$ then
   $h_A(z) = 1 + \sum_{i = 0}^{s}(\rho_{i-1} - \rho_{i})z^i$
    \item
    $ e_i = \sum_{j\geq i -1}\binom{j}{i-1}\rho_j \geq 0 $ for all $i \geq 1$.
   \end{enumerate}

\s Assume $\operatorname{depth} A > 0$. For $i \geq 1$ set
\[
\widetilde{\mathfrak{m} ^i} = \bigcup_{l \geq 1}(\mathfrak{m} ^{i+l} \colon \mathfrak{m} ^l).
\]
Then $\{ \widetilde{\mathfrak{m} ^n} \}_{n \geq 0}$ is called the Ratliff-Rush filtration of $A$ (with respect to $\mathfrak{m} $).
We have
\begin{enumerate} [\rm (1)]
    \item
$\widetilde{\mathfrak{m} ^i}\widetilde{\mathfrak{m} ^j} \subseteq \widetilde{\mathfrak{m} ^{i+j}}$ for all $ i, j \geq 0$.
\item
$\widetilde{\mathfrak{m} ^i} = \mathfrak{m} ^i$ for all $i \gg 0$.
\item
If $x$ is $A$-superficial then $(\widetilde{\mathfrak{m} ^{i+1}} \colon x) = \widetilde{\mathfrak{m} ^i}$ for all $i \geq 0$.
\item
$\operatorname{depth} G(A) > 0$ if and only if $\widetilde{\mathfrak{m} ^i} = \mathfrak{m} ^i$ for all $i \geq 1$.
\end{enumerate}

\s \label{ratliff-hilb} Let $\operatorname{depth} A > 0$. Let $\widetilde{G}(A) = \bigoplus_{n \geq 0}\widetilde{\mathfrak{m} ^n}/\widetilde{\mathfrak{m} ^{n+1}}$ be the associated graded ring
of the Ratliff-Rush filtration. Let its Hilbert series be
\[
\sum_{n \geq 0}\ell\left( \widetilde{\mathfrak{m} ^n}/\widetilde{\mathfrak{m} ^{n+1}}\right)z^n = \frac{\widetilde{h}_A(z)}{(1-z)^d} \quad  \text{where} \ \widetilde{h}_A(z) \in \mathbb{Z}[z].
\]
Set $r(z) = \sum_{n \geq 0} \ell(\widetilde{\mathfrak{m} ^{n+1}}/\mathfrak{m} ^{n+1})z^n$. Then $r(z) \in \mathbb{Z} [z]$ and
\[
h_A(z) = \widetilde{h}_A(z) + (1-z)^{d+1}r(z).
\]

\s\label{bella} (\cite[Theorem 3]{It}) Let $(A,\mathfrak{m} )$ be a Cohen-Macaulay \ local ring of dimension two.  Suppose  $x,y$ is an $A$-superficial sequence. Set $J = (x,y)$.
Let $\sigma_j = \ell(\widetilde{\mathfrak{m} ^{j+1}}/J\widetilde{\mathfrak{m} ^j})$ for $j \geq 0$. Then
\[
e_1 = \sum_{j \geq 0} \sigma_j  \quad \text{and} \  e_2 = \sum_{j \geq j} j\sigma_j.
\]
Furthermore $\widetilde{h}_A(z) = 1 + \sum_{i \geq 1}(\sigma_{i-1} - \sigma_i)z^i$.

We will need the following result in the proof of Theorem \ref{depth0}.
\begin{proposition}\label{r-3}
Set $G = G(A)$ and $\widetilde{G} = \widetilde{G}(A)$.
 Assume $\widetilde{G}$ is Cohen-Macaulay. Then
\begin{enumerate}[\rm (1)]
\item
$G$ is generalized Cohen-Macaulay.
\item
$\dim G/P = d$ for all minimal primes $P$ of $G$.
\end{enumerate}
\end{proposition}
\begin{proof}
(1) We note that $\widetilde{G}$ is finitely generated as a $G$-module. Furthermore  the natural map $G \rightarrow \widetilde{G}$ induces an exact sequence of $G$-modules
\[
0 \rightarrow U \rightarrow G \rightarrow \widetilde{G} \rightarrow V \rightarrow 0,
\]
where $U, V$ are of finite length. Let $\mathfrak{M} $ be the maximal homogeneous ideal of $G$. If $P$ is a prime ideal in $G$ with $P \neq \mathfrak{M} $ then $G_P \cong \widetilde{G}_P$. It  follows that $G_P$ is Cohen-Macaulay \ for all $P \neq \mathfrak{M} $. So $G$ is generalized Cohen-Macaulay.

(2) Let $P$ be a minimal prime of  $G$. Then $P \neq \mathfrak{M} $. Furthermore $G_P \cong \widetilde{G}_P$. So $P$ is an associate prime of
$\widetilde{G}$. The result now follows from \cite[2.1.2(a)]{BH}.
\end{proof}

\s \label{quot} Assume $(A,\mathfrak{m} )$ is Cohen-Macaulay \ of dimension $d \geq 2$. Let $x_1,\ldots, x_d$ be an $A$-superficial sequence. Set $J = (x_1,\ldots, x_d)$ and $(B,\mathfrak{n} ) = (A/(x_1), \mathfrak{m} /(x_1))$.  If $I$ is an ideal in $A$ then set $\overline{I}$ to be its image in $B$. Notice $\overline{\widetilde{\mathfrak{m} ^i}} \subseteq \widetilde{\mathfrak{n} ^i}$. So we have a natural maps
 $$\eta_i \colon \widetilde{\mathfrak{m} ^{i+1}}/ J\widetilde{\mathfrak{m} ^i}  \rightarrow \widetilde{\mathfrak{n} ^{i+1}}/\overline{J} \widetilde{\mathfrak{n} ^i} \quad \text{for all} \ i \geq 0. $$
We show
\begin{proposition}\label{nn}
(with setup as in \ref{quot})
\begin{enumerate}[\rm (1)]
\item
If $\overline{\widetilde{\mathfrak{m} ^s}} = \widetilde{\mathfrak{n} ^s}$ for some $s$ then $\eta_s$ is injective.
\item
If $\overline{\widetilde{\mathfrak{m} ^j}} = \widetilde{\mathfrak{n} ^j}$ for $j = s, s+1$ then $\eta_s$ is bijective.
\end{enumerate}
\end{proposition}
\begin{proof}
(1) Suppose $\eta_s( p + J \widetilde{\mathfrak{m} ^s}) = 0$. Then $\overline{p} = \sum_{ t = 2}^{d}a_tx_t$ for some $a_t \in \widetilde{\mathfrak{n} ^s}$. As $\overline{\widetilde{\mathfrak{m} ^s}} = \widetilde{\mathfrak{n} ^s}$ there exists $b_t \in \widetilde{\mathfrak{m} ^s}$ with $\overline{b_t} = a_t$ for $t = 2,\ldots, d$. So
\[
p = \sum_{t =2}^{d} b_tx_t +   rx_1 \quad \text{for some} \ r \in A.
\]
We get  $xr \in \widetilde{\mathfrak{m} ^{s+1}}$. So $r \in \widetilde{\mathfrak{m} ^s}$. Thus $p \in J \widetilde{\mathfrak{m} ^s}$. It follows that $\eta_s$ is injective.

(2) As $\overline{\widetilde{\mathfrak{m} ^{s+1}}} = \widetilde{\mathfrak{n} ^{s+1}}$ then clearly $\eta_s$ is surjective. As $\overline{\widetilde{\mathfrak{m} ^s}} = \widetilde{\mathfrak{n} ^s}$ we get by (1) that $\eta_s$ is injective. So $\eta_s$ is bijective.
\end{proof}
The following result is needed in the  proof of Theorem \ref{depth0}.
\begin{lemma}\label{m3J}
Let $(A,\mathfrak{m} )$ be a Cohen-Macaulay \ local ring of dimension three with an infinite residue field. Let $x,y$ be an $A$-superficial sequence. Let $z$ be $A\oplus A/(x,y)$-superficial (note $z$ is $A$-superficial and $x,y,z$ is an $A$-superficial sequence). Set $J = (x,y,z),  (B, \mathfrak{n} ) = (A/(x), \mathfrak{m} /(x)) $ and $\overline{J}
= (y,z)B$. Consider the natural complex
$$ 0 \rightarrow \frac{ (\mathfrak{m} ^3  \colon x)}{\mathfrak{m} ^2} \xrightarrow{\alpha} \frac{ \mathfrak{m} ^3}{J\mathfrak{m} ^2} \xrightarrow{\pi} \frac{\mathfrak{n} ^3}{\overline{J} \mathfrak{n} ^2} \rightarrow 0$$
where $\pi$ is the natural surjection and $\alpha(a + \mathfrak{m} ^2) = xa + J\mathfrak{m} ^2$.  Then the above complex is exact.
\end{lemma}
\begin{proof}
Clearly $\pi$ is surjective. Let $\pi(a + J\mathfrak{m} ^2) = 0$. It follows that $a = y\xi_2 + z \xi_2^\prime + xr$ where $\xi, \xi^\prime \in \mathfrak{m} ^2$ and $r \in A$. Note $r \in (\mathfrak{m} ^3 \colon x)$. Furthermore $\alpha(r + \mathfrak{m} ^2) = a + J\mathfrak{m} ^2$.

It remains to prove $\alpha$ is injective. Suppose $\alpha(t + \mathfrak{m} ^2) = 0$. It follows that
$tx  = ax + by + cz$ where $a,b,c \in \mathfrak{m} ^2$. Set $(C, \mathfrak{q} ) = (A/(z), \mathfrak{m} /(z))$. We note that $J$ is a minimal reduction of $\mathfrak{m} $. So $(x,y)C$ is a
minimal reduction of $\mathfrak{q} $.   In $C$ we get $\overline{t}x = \overline{a}x + \overline{b} y$. As $x,y$ is a regular sequence in $C$ we get
$\overline{t}- \overline{a} = \delta y$ and $\overline{b} = \delta x$ for some $\delta \in A$. As $(x,y)C$ is a minimal reduction of $\mathfrak{q} $ and $\overline{b} \in \mathfrak{q} ^2$, by analyticity (see \cite[Lemma 2]{NR}) we get $\delta \in  \mathfrak{m} $.
Thus $t = a + y\delta + rz$ for some $r \in A$. Note $t \in \widetilde{\mathfrak{m} ^2}$. So $rz \in \widetilde{\mathfrak{m} ^2}$. As $z$ is $A$-superficial we get
$r \in \mathfrak{m} $. It follows that $t \in \mathfrak{m} ^2$. Thus $\alpha$ is injective.
\end{proof}
\section{$e_2 = e_1 - e + 1 = 3$ and $\dim A = 2$}
To prove this result we need the following:
\begin{proposition}\label{dim2}
Let $(A,\mathfrak{m} )$ be a Cohen-Macaulay  \ local ring of dimension two and with an infinite residue field. Let $x,y$ be an $A$-superficial sequence. Set $(B,\mathfrak{n} ) = (A/(x), \mathfrak{m} /(x))$ and $J = (x,y)$.  Assume $e_2 = e_1 - e_0 + 1 = 3$.  Set $G = G(A)$. Then
\begin{enumerate}[\rm(1)]
\item
$\operatorname{depth} G = 0$ or $2$. If $\operatorname{depth} G = 2$ then $h_A(z) = 1 + hz + 3z^2$. In particular $e = h + 4$ if $\operatorname{depth} G = 2$.
\item
If $\operatorname{depth} G = 0$ then
\begin{enumerate}[\rm (a)]
\item
$\widetilde{\mathfrak{m} ^{n+1}}  = J \widetilde{\mathfrak{m} ^n}$ for all $n \geq 2$.
\item
$\widetilde{G}(A)$ is Cohen-Macaulay.

\item
 $\widetilde{\mathfrak{m} ^2} \neq \mathfrak{m} ^2$.

\item
 $\ell(\widetilde{\mathfrak{m} ^2}/\mathfrak{m} ^2) = 1$, $\ell(\mathfrak{m} ^2/J\mathfrak{m} ) = 2$  and $e = h + 3$.
\item
 $\widetilde{h}_A(z) = 1 + (h-1)z + 3z^2$.
\item
$h_B(z) = 1 + hz + z^2 + z^3$.

\item
 $\widetilde{\mathfrak{m} ^j} = \mathfrak{m} ^j$ for all $j \geq 3$. Furthermore $\ell(\mathfrak{m} ^3/J\mathfrak{m} ^2) = 2$ and $\mathfrak{m} ^4 = J\mathfrak{m} ^3$.

\item
$h_A(z) = 1 + hz + 3z^3 - z^4$.

\item
$\operatorname{type}(A) = h$.
\end{enumerate}
\end{enumerate}
\end{proposition}
\begin{proof}
Let $x,y$ be an $A$-superficial sequence  with respect to \ $\mathfrak{m} $.
(1) If $\operatorname{depth} G \geq 1$ then we have $e_2(B) = e_1(B) - e_0(B) + 1$. Using \ref{dim1} it follows that $\mathfrak{n} ^3 = y\mathfrak{n} ^2$ and $\ell(\mathfrak{n} ^2/y\mathfrak{n} )  = 3$.  So by \cite[2.1]{S2} we have that $\operatorname{depth} G_\mathfrak{n} (B) = 1$ and $h_B(z) = 1 + hz + 3z^2$. So by Sally descent we have $\operatorname{depth} G = 2$.
Also $h_A(z) = h_B(z)$.

(2) For $j \geq 0$ set $\sigma_j = \ell(\widetilde{\mathfrak{m} ^{j+1}}/J\widetilde{\mathfrak{m} ^{j}})$. Then by \ref{bella} we have $e_1 = \sum_{j \geq 0} \sigma_j$ and $e_2 = \sum_{j \geq 1}j\sigma_j$. Since $e_2 = e_1 - e_0 + 1$ we get $\sigma_j = 0$ for $j \geq 2$ and
$e_2 = \sigma_1 = \ell(\widetilde{\mathfrak{m} ^{2}}/J\mathfrak{m} ) = 3$.

(a) This follows as $\sigma_j = 0$ for $j \geq 2$.\\
(b) We have $\widetilde{\mathfrak{m} ^2} \cap J = J \mathfrak{m}  $ and for $i \geq 2$ we also have $\widetilde{\mathfrak{m} ^{i+1}} \cap J = J\widetilde{\mathfrak{m} ^i}$ (since $\sigma_j = 0$ for $j \geq 2$). Thus by Valabrega-Valla theorem for filtration's we get that $\widetilde{G}(A)$ is Cohen-Macaulay \ (see \cite[3.5]{HM}).

(c) If $\widetilde{\mathfrak{m} ^2} = \mathfrak{m} ^2$ then as $\sigma_2 = 0$ we get $\widetilde{\mathfrak{m} ^3} = J\mathfrak{m} ^2 $. It follows that $\widetilde{\mathfrak{m} ^3} = \mathfrak{m} ^{3}$. As $\sigma_j = 0$ for $j \geq 2$ inductively one can show that $\widetilde{\mathfrak{m} ^j} = \mathfrak{m} ^j$ for $j \geq 2$. It follows that $\operatorname{depth} G \geq 1$, a contradiction. Therefore $\widetilde{\mathfrak{m} ^2} \neq \mathfrak{m} ^2$.

(d) Note
\[
3 = e_2 = \ell(\widetilde{\mathfrak{m} ^{2}}/J\mathfrak{m} ) = \ell(\widetilde{\mathfrak{m} ^2}/\mathfrak{m} ^2) + \ell(\mathfrak{m} ^2/J\mathfrak{m} ).
\]
It follows that $\ell(\mathfrak{m} ^2/J\mathfrak{m} ) \leq 2$. If $\ell(\mathfrak{m} ^2/J\mathfrak{m} ) \leq 1$ then by  \cite{rv1} or \cite{W} we have $\operatorname{depth} G \geq 1$, a contradiction. Thus $\ell(\mathfrak{m} ^2/J\mathfrak{m} ) = 2$. It follows that $\ell(\widetilde{\mathfrak{m} ^2}/\mathfrak{m} ^2) = 1$. We also have $e = h + 3$ by \ref{SA}.

(e) This follows from \ref{bella}.

(f) As $e_i(B) = e_i$ for $i \leq 1$, \ref{mod-sup} we get $e_1(B) = e_0(B) + 2$.
 Note $e_1(B) = \sum_{j\geq 0}\ell(\mathfrak{n} ^{j+1}/y\mathfrak{n} ^j)$. Also $\ell(\mathfrak{n} ^2/y\mathfrak{n} ) = \ell(\mathfrak{m} ^2/J\mathfrak{m} ) = 2$. It follows that $\ell(\mathfrak{n} ^3/y\mathfrak{n} ^2) = 1$ and $\mathfrak{n} ^4 = y \mathfrak{n} ^3$.
 Thus $e_2(B) =  \sum_{j\geq 1}j\ell(\mathfrak{n} ^{j+1}/y\mathfrak{n} ^j) = 4$. The formula for $h_B(z)$ follows from \ref{dim1}.

(g) By  \ref{mod-sup} we have
\[
e_2 = e_2(B) - \sum_{n \geq 2} \ell((\mathfrak{m} ^{n+1} \colon x)/\mathfrak{m} ^n).
\]
So $ \sum_{n\geq 2} \ell((\mathfrak{m} ^{n+1} \colon x)/\mathfrak{m} ^n) = 1$.  We have an exact sequence, see \cite[p.\ 305]{rv2}
\begin{equation*}
0 \rightarrow \frac{(\mathfrak{m} ^3 \colon x)}{\mathfrak{m} ^2} \rightarrow \frac{\mathfrak{m} ^3}{J\mathfrak{m} ^2} \rightarrow \frac{\mathfrak{n} ^3}{y\mathfrak{n} ^2} \rightarrow 0. \tag{*}
\end{equation*}
If $(\mathfrak{m} ^3 \colon x) = \mathfrak{m} ^2$ then $\ell(\mathfrak{m} ^3/J\mathfrak{m} ^2) = 1$. So by a result due to Huckaba \cite{H} we get $\operatorname{depth} G \geq 1$, a contradiction. So $(\mathfrak{m} ^3 \colon x) \neq \mathfrak{m} ^2$.
It follows that $(\mathfrak{m} ^{n+1} \colon x) = \mathfrak{m} ^n $ for $n \geq 3$.
For all $i \geq 0$ we have an exact sequence, see \cite[2.6]{Pu5},
 \begin{equation*}
  0 \rightarrow \frac{(\mathfrak{m} ^{i+1}\colon x)}{\mathfrak{m} ^{i}} \rightarrow \frac{\widetilde{\mathfrak{m} ^i}}{\mathfrak{m} ^i}  \rightarrow \frac{\widetilde{\mathfrak{m} ^{i+1}}}{\mathfrak{m} ^{i+1}}
 \end{equation*}
 It follows that $\widetilde{\mathfrak{m} ^j} = \mathfrak{m} ^j$ for $j \geq 3$. By (*) we get $\ell(\mathfrak{m} ^3/J\mathfrak{m} ^2) = 2$. Also by (a),  $\widetilde{\mathfrak{m} ^4} = J\widetilde{\mathfrak{m} ^3}$. As $\widetilde{\mathfrak{m} ^j} = \mathfrak{m} ^j$ for $j \geq 3$ we get $\mathfrak{m} ^4 = J\mathfrak{m} ^3$.

(h) This follows from \ref{ratliff-hilb}.

(i) If $\operatorname{type}(A) < h$ then it follows from \cite[2.3]{RV} that $\ell(\mathfrak{m} ^3/J\mathfrak{m} ^2) \leq 1$ and so by a result of Huckaba \cite{H} we get $\operatorname{depth} G \geq 1$ which is a contradiction.
We show that $\operatorname{type}(A) > h$ is also not possible. Set $(C,\mathfrak{q} ) = (A/(x,y) , \mathfrak{m} /(x,y))$. As $\mathfrak{m} ^3 \subseteq J$ we  get $\mathfrak{q} ^3 = 0$.
Also as $\mathfrak{m} ^2 \nsubseteq J$ we get $\mathfrak{q} ^2 \neq 0$. So
$h_C(z) = 1 + hz + 2z^2$.  Note $\operatorname{type}(A) = \operatorname{type}(C)$. As $\mathfrak{q} ^2 \neq 0 $ we get $\operatorname{type}(C) \leq h + 1$. If $\operatorname{type}(C) = h + 1$ then
$\mathfrak{q}  = (a, K)$ where $K \subseteq (0 \colon \mathfrak{q} )$. Then $\mathfrak{q} ^2 = (a^2)$ a contradiction as $\mu(\mathfrak{q} ^2) = 2$.  Thus $\operatorname{type}(C) = h$.
\end{proof}

\section{  $e_2 = e_1 - e + 1 = 3$ and $\dim A = 3$ }
The most difficult case of Theorem \ref{e2} occurs when $\dim A = 3$. We prove it separately.
\begin{theorem}\label{depth0}
Let $(A,\mathfrak{m} )$ be a Cohen-Macaulay \ local ring of dimension $3$ and infinite residue field. Assume $e_2 = e_1 - e + 1 = 3$.  Then one of the following cases occur
\begin{enumerate}[\rm (i)]
\item
$e = h + 4$, $G$ is Cohen-Macaulay \  and $h_A(z) = 1 + hz + 3z^2$.
\item
$e = h + 3$, $d \geq 2$, $\operatorname{depth} G = 1$, $\operatorname{type}(A) = h$ and

$\displaystyle{h_A(z) = 1 + hz + 3z^3 - z^4.}$
\end{enumerate}
\end{theorem}
\begin{proof}
 Let $x$ be any $A$-superficial element. Set $(B,\mathfrak{n} ) = (A/(x), \mathfrak{m} /(x))$.  Then $e_i(A) = e_i(B)$ for $i = 0,1,2$. By \ref{dim2} we get that
$e = h + 4$ or $e = h + 3$. If $e = h + 4$ then $G(B)$ is Cohen-Macaulay. So by Sally descent we get $G(A)$ is Cohen-Macaulay. Furthermore
$h_A(z) = h_B(z) = 1 + hz + 3z^2$. If $e = h + 3$ then $G(B)$ has depth $0$. Note $\operatorname{type}(A) = \operatorname{type}(B) = h$. By \ref{SD} we get $\operatorname{depth} G(A) \leq 1$. If $\operatorname{depth} G(A) = 1$ then $x^*$ is $G(A)$-regular and $G(A)/(x^*) = G(B)$. So $h_A(z) = h_B(z) = 1 + hz + 3z^3 - z^4$ (by \ref{dim2}).

Suppose if possible $\operatorname{depth} G(A) \neq 1$. Then by the above argument we get $\operatorname{depth} G(A) = 0$. We show:

\emph{Claim}(1) : $\ell(\widetilde{\mathfrak{m} ^2}/\mathfrak{m} ^2) = 1$.

We have an exact sequence for all $i \geq 1$ (see \cite[2.9]{Pu5})
\[
 0 \rightarrow \frac{(\mathfrak{m} ^{i+1} \colon x)}{\mathfrak{m} ^i} \rightarrow \frac{\widetilde{\mathfrak{m} ^{i}}}{\mathfrak{m} ^{i}}  \rightarrow \frac{\widetilde{\mathfrak{m} ^{i+1}}}{\mathfrak{m} ^{i+1}}\rightarrow
 \frac{\widetilde{\mathfrak{n} ^{i+1}}}{\mathfrak{n} ^{ i + 1}} \tag{*}
\]
As $\widetilde{\mathfrak{m} } = \mathfrak{m} $ we get an inclusion $\widetilde{\mathfrak{m} ^2}/\mathfrak{m} ^2 \rightarrow \widetilde{\mathfrak{n} ^2}/\mathfrak{n} ^2$. Furthermore as $\widetilde{\mathfrak{n} ^j} = \mathfrak{n} ^j$ for
$j \geq 3$  (see \ref{dim2}.2(g)) we get surjections $\widetilde{\mathfrak{m} ^j}/\mathfrak{m} ^j \rightarrow \widetilde{\mathfrak{m} ^{j+1}}/\mathfrak{m} ^{j+1}$ for $j \geq 2$. Thus $\ell(\widetilde{\mathfrak{m} ^2}/\mathfrak{m} ^2) \leq 1$. Furthermore if $\widetilde{\mathfrak{m} ^2} = \mathfrak{m} ^2$ we get $\widetilde{\mathfrak{m} ^j} = \mathfrak{m} ^j $ for all $j$ and so $\operatorname{depth} G(A) > 0$ which is a contradiction.
Thus $\ell(\widetilde{\mathfrak{m} ^2}/\mathfrak{m} ^2) =1$.

\emph{Claim}(2): $\widetilde{G}(A)$ is Cohen-Macaulay.

In the inclusion $\widetilde{\mathfrak{m} ^2}/\mathfrak{m} ^2 \rightarrow \widetilde{\mathfrak{n} ^2}/\mathfrak{n} ^2$ both modules have same length. So $\overline{\widetilde{\mathfrak{m} ^2}}  = \widetilde{\mathfrak{n} ^2}$. Furthermore as $\widetilde{\mathfrak{n} ^j} = \mathfrak{n} ^j$ for $j \geq 3$ we get $\overline{\widetilde{\mathfrak{m} ^j}}  = \widetilde{\mathfrak{n} ^j}$ for $j \geq 3$.
Thus the Ratilff-Rush filtration on $A$ behaves well with respect to \ superficial element $x$. So
 $\widetilde{G}(A)/(x^*) = \widetilde{G}(B)$. As $\widetilde{G}(B)$ is Cohen-Macaulay \  (see \ref{dim2}.2(b) ) and $x^*$ is $\widetilde{G}(A)$-regular we get $\widetilde{G}(A)$ is Cohen-Macaulay.

\emph{Claim}(3):   $G(A)$ is generalized Cohen-Macaulay.

 By \ref{r-3} we get $G(A)$ is generalized Cohen-Macaulay.

\emph{Claim}(4):  $\widetilde{\mathfrak{m} ^j} = \mathfrak{m} ^j $ for all $j \geq 3$. Furthermore $\ell((\mathfrak{m} ^3 \colon x)/\mathfrak{m} ^2) = 1$.

 As $\widetilde{\mathfrak{n} ^3} = \mathfrak{n} ^3$, by (*) we have an exact sequence
\begin{equation*}
 0 \rightarrow \frac{(\mathfrak{m} ^{3} \colon x)}{\mathfrak{m} ^2} \rightarrow \frac{\widetilde{\mathfrak{m} ^2}}{\mathfrak{m} ^{2}}  \rightarrow \frac{\widetilde{\mathfrak{m} ^{3}}}{\mathfrak{m} ^{3}}\rightarrow 0. \tag{**}
\end{equation*}
Suppose if possible $\widetilde{\mathfrak{m} ^3} \neq \mathfrak{m} ^3$.  Then as $\ell(\widetilde{\mathfrak{m} ^2}/\mathfrak{m} ^2) = 1$ we get
$(\mathfrak{m} ^3 \colon x) = \mathfrak{m} ^2$ and $\ell(\widetilde{\mathfrak{m} ^3}/\mathfrak{m} ^3) = 1$. We note that we chose $x$ to be \emph{any} $A$-superficial element.
As $G(A)$ has depth zero and also as it is generalized Cohen-Macaulay  \ it follows that
\[
\operatorname{Ass}(G(A)) = \{ \mathfrak{M} , P_1, \ldots, P_s \};
\]
where $P_1, \ldots, P_s$ are minimal primes of $G(A)$ and $\mathfrak{M} $ is the maximal homogeneous ideal of $G(A)$. By \ref{r-3} we get $\dim G(A) /P_i =   3$ for all $i$. Let $V$ be the $k$-vector space
$\mathfrak{m} /\mathfrak{m} ^2$.  If $W$ is a $k$ vector  space then we set $\operatorname{vdim} W$
to be its dimension as a $k$-vector space.  If $V \subseteq P_i$ then $P_i = \mathfrak{M} $, a contradiction. In particular $\operatorname{vdim} V \cap P_i \leq \operatorname{vdim} V -1$ for all $i$.

\textit{sub-claim:} $\operatorname{vdim} P_i \cap V \leq \operatorname{vdim} V - 2$ for all $i$.\\
Suppose for some $i$ we have $\operatorname{vdim} P_i \cap V =  \operatorname{vdim} V - 1$. Then there exists $t \in V$ such that $P_i \cap V + k t = V$. So we get
$P_i + G(A)t = \mathfrak{M} $. It follows that $(t)$ generates the maximal homogeneous ideal in $G(A)/P_i$. So $\dim G(A)/P_i \leq 1$ a contradiction.
 Thus we have proved the sub-claim.

As $k$ is infinite there exists $k$-linearly independent vectors $u^*, v^*$ in $V$ such that if $H  = ku^* + kv^*$ then $H \cap P_i = 0$ for all $i$. It follows that $u - \lambda v$ is an $A$-superficial sequence for any unit $\lambda$ (see \ref{exist-sup}). Furthermore $u, v$ are $A$-superficial.

Suppose $\widetilde{\mathfrak{m} ^2} = \mathfrak{m} ^2 + (a)$. Let $x$ be any $A$-superficial element. If $xa \in \mathfrak{m} ^3$ then $a \in (\mathfrak{m} ^3 \colon x) = \mathfrak{m} ^2$ a contradiction. Thus $\overline{xa}$ is a non-zero element in $\widetilde{\mathfrak{m} ^3} /\mathfrak{m} ^3$. In particular $\overline{ua}, \overline{va}$ are non-zero in $\widetilde{\mathfrak{m} ^3} /\mathfrak{m} ^3$. As $\widetilde{\mathfrak{m} ^3} /\mathfrak{m} ^3 \cong k$ we get that there is a unit $\alpha$  such that  $\overline{ua} = \alpha \overline{va}$.
So $(u - \alpha v)a \in \mathfrak{m} ^3$. By construction $u - \alpha v$ is an $A$-superficial element. So $a \in (\mathfrak{m} ^3 \colon  u - \alpha v) = \mathfrak{m} ^2$ a contradiction. Thus $\widetilde{\mathfrak{m} ^3} = \mathfrak{m} ^3$. As we have surjections  $\widetilde{\mathfrak{m} ^j}/\mathfrak{m} ^j \rightarrow \widetilde{\mathfrak{m} ^{j+1}}/\mathfrak{m} ^{j+1}$ for $j \geq 3$ we also get $ \widetilde{\mathfrak{m} ^j} = \mathfrak{m} ^j$ for $j \geq 4$.  By exact sequence (**) we also get $\ell((\mathfrak{m} ^3 \colon x)/\mathfrak{m} ^2 ) = 1$. This proves claim (4)

Let $x,y,z$ be an $A$-superficial sequence. Set $J = (x,y,z)$.

\emph{Claim} (5): $\widetilde{\mathfrak{m} ^{i+1}} = J \widetilde{\mathfrak{m} ^i}$ for $i \geq 2$.  In particular $\mathfrak{m} ^3 = J \widetilde{\mathfrak{m} ^2}$ and $\mathfrak{m} ^4 = J \mathfrak{m} ^3$.

We have the natural map $\eta_i \colon \widetilde{\mathfrak{m} ^{i+1}}/ J \widetilde{\mathfrak{m} ^i} \rightarrow  \widetilde{\mathfrak{n} ^{i+1}}/ (y,z)\widetilde{\mathfrak{n} ^i}$  for all $i \geq 1$. By Claim(1) we have $\overline{\widetilde{\mathfrak{m} ^i}} = \widetilde{\mathfrak{n} ^i}$ for all $i \geq 1$. From \ref{nn} we get that $\eta_i$ is an isomorphism for $i \geq 1$. The result now follows from \ref{dim2}.(2)(a)  and Claim (4).

 \emph{Claim}(6) : $\ell(\mathfrak{m} ^2/J\mathfrak{m} ) = 2$ and  $\ell(\mathfrak{m} ^3/J\mathfrak{m} ^2) = 3$.

As $e = h + 3$ we get by \ref{SA} that $\ell(\mathfrak{m} ^2/J\mathfrak{m} ) = 2$.  By \ref{m3J}, for a suitable choice of $x,y,z$,  we get  a short exact sequence
\[
 0 \rightarrow \frac{ (\mathfrak{m} ^3  \colon x)}{\mathfrak{m} ^2} \xrightarrow{\alpha} \frac{ \mathfrak{m} ^3}{J\mathfrak{m} ^2} \xrightarrow{\pi} \frac{\mathfrak{n} ^3}{\overline{J} \mathfrak{n} ^2} \rightarrow 0.
\]
By \ref{dim2}(2)g and claim (4) we get for a choice of $x,y,z$ we have $\ell(\mathfrak{m} ^3/J\mathfrak{m} ^2) = 3$. By \cite{Pu-r};  $\ell(\mathfrak{m} ^3/J\mathfrak{m} ^2)$ is an invariant of $A$ and does not depend on choice of minimal reduction $J$. Thus  for any $A$-superficial sequence $x,y,z$ we have
$\ell(\mathfrak{m} ^3/J\mathfrak{m} ^2) = 3$.

\emph{Claim}(7):  $S = A/\mathfrak{m}  \oplus \mathfrak{m} /(J + \mathfrak{m} ^2) \oplus \mathfrak{m} ^2/J\mathfrak{m}  \oplus \mathfrak{m} ^3/J\mathfrak{m} ^2$ is a finite dimensional homogeneous $k$-algebra.

In fact $S = G(A)/(x^*, y^*, z^*)$ (use $\mathfrak{m} ^3 \subseteq J\mathfrak{m} $ and $\mathfrak{m} ^4 = J\mathfrak{m} ^3)$.

\emph{Claim}(8):  $S$ violates Macaulay's bound  ( \cite[4.2.10(c) ]{BH})for Hilbert function of $k = A/\mathfrak{m} $-algebras.

Let $f$ be the Hilbert function of $S$. We have by claim (6)
\[
f(2) = 2  = \binom{2}{2} + \binom{1}{1}.
\]
By Macaulay's bound we have
\[
3 = f(3) \leq  f(2)^{<2>} = \binom{3}{3} + \binom{2}{2} = 2,  \quad \text{a contradiction.}
\]
Thus our assumption $\operatorname{depth} G(A) \neq 1$ is incorrect. So $\operatorname{depth} G(A) = 1$.
\end{proof}

\section{Proof of Theorem \ref{e2}}
In this section we give
\begin{proof}[Proof of Theorem \ref{e2}]
We may assume that the residue field of $A$ is infinite.
If $\dim A \leq 1$ then by \cite[6.21]{V} we get $G(A)$ is Cohen-Macaulay  \ and $\deg h_A(z) \leq 2$. It follows that $h_A(z) = 1 + hz + 3z^2$.
When $d = \dim A \geq 2$ we prove the result by induction on $d$.
The case when $\dim A = 2$ is analyzed in \ref{dim2}.
The case when $\dim A = 3$ is analyzed in \ref{depth0}.

 Now assume $d = \dim A \geq 4$ and the result holds for Cohen-Macaulay \ local rings of dimension $d - 1 \geq 3$.

Let $x$ be an $A$-superficial. Set $(B, \mathfrak{n} ) = (A/(x), \mathfrak{m} /(x))$.
We  consider the following two cases:\\
Case (1)  $G(A)$ is Cohen-Macaulay.  So $G(B)$ is Cohen-Macaulay. In this case $h_A(z) = h_B(z) = 1 + hz + 3z^2$.

Case (2) $G(A)$ is not Cohen-Macaulay.  So $G(B)$ is not Cohen-Macaulay \ by Sally descent. In this case we have $e = h + 3$.  By induction hypothesis we get $\operatorname{depth} G(B) = d - 3 \geq 1$. So by Sally descent
we get $\operatorname{depth} G(A)  = d -2$. Thus  $x^*$ is $G(A)$-regular. So
\[
h_A(z) = h_B(z) = 1 + hz + 3z^3 - z^4.
\]
\end{proof}

\section{Proof of Theorem \ref{e1}}
In this section we give
\begin{proof}[Proof of Theorem \ref{e1}]
We may assume $k = A/\mathfrak{m} $ is infinite.
We first consider the case when $d = 0$. Clearly $G(A)$ is Cohen-Macaulay.
Let $h_A(z) = 1 + hz + a_2z^2 + \cdots + a_sz^s$.
We have $e_1 = e +2$. So
\[
h+ \sum_{i = 2}^{s}ia_i = 3+ h + \sum_{i = 2}^{s} a_i
\]
So we get
\[
\sum_{i = 2}^{s}(i - 1)a_i = 3
\]
So the possibilities are
\begin{enumerate}
\item
$a_2 = 3$, $a_j = 0$ for $j \geq 3$. In this case $h_A(z) = 1 + hz + 3z^2$.
\item
$a_2 = a_3 = 1$ and $a_j = 0$ for $j \geq 4$.  In this case $h_A(z) = 1 + hz + z^2 + z^3$.
\end{enumerate}

Next we consider the case $d = 1$.  Let $x$ be $A$-superficial. Set $\rho_n = \lambda(\mathfrak{m} ^{n+1}/x\mathfrak{m} ^n)$.    If $\deg h_A(z) = s$ then
   $h_A(z) = 1 + \sum_{i = 0}^{s}(\rho_{i-1} - \rho_{i})z^i$. Also
 $ e_1 = \sum_{j \geq 0}\rho_j  $.  Note $\rho_0 = e -1$
   In particular
\[
e_1 = e - 1 + \sum_{j \geq 1} \rho_j.
\]
Thus we have
\[
\sum_{j \geq 1} \rho_j = 3
\]
Thus $\rho_1 = 3$ or $2$ or $1$.

If $\rho_1 = 3$ then $\rho_j = 0$ for $j \geq 2$. As $\mathfrak{m} ^3 = x \mathfrak{m} ^2$ we get $G(A)$ is Cohen-Macaulay. Furthermore
$h_A(z) = 1 + hz + 3z^2$. In particular $e = h + 4$.

If $\rho_1 = 2$ then $\rho_2 = 1$ and $\rho_j = 0$ for $j \geq 3$. So we have $h_A(z) = 1 + hz + z^2 + z^3$. In particular we have $e = h +3$.

If $\rho_1 = 1$ then $\rho_j \leq 1$ for $j \geq 2$ (see \cite[3.3]{rv2}).
So we get $\rho_2 = \rho_3 = 1$. So we have $h_A(z) = 1 + hz  + z^4$. In particular $e = h + 2$ and $G(A)$ is not Cohen-Macaulay. It follows from \cite[3.1]{S4} that $\operatorname{type}(A) = h$.

Next we consider the case $d = 2$. Let $x$ be an $A$-superficial. Set $(B, \mathfrak{n} ) = (A/(x), \mathfrak{m} /(x))$.
We have $e_i(A) = e_i(B)$ for $i = 0, 1$.  So we have $e = h + 2$ or $e = h + 3$ or $e = h +4$.

 If $e = h + 4$ then $G(B)$ is Cohen-Macaulay. So by
Sally descent we get $G(A)$ is Cohen-Macaulay.  Furthermore $h_A(z) = h_B(z) = 1 + hz + 3z^2$.

If $e = h + 2$ then by \cite{rv1} or \cite{W} we get $\operatorname{depth} G(A) = 1$. So $x^*$ is $G(A)$-regular. Thus $h_A(z) = h_B(z) = 1 + hz  + z^4$.
We also have
$\operatorname{type}(A) =  \operatorname{type}(B) = h$.

Now consider the case $e = h + 3$.  By Itoh's bound we have $e_2(A) \geq e_1 - e + 1 = 3$.  Set $b_i = \ell((\mathfrak{m} ^{i+1} \colon x)/\mathfrak{m} ^i)$ for $i \geq 0$. Note $e_2(B) = 4$. Then by \ref{mod-sup} we get $e_2(A) = 4 - \sum_{i \geq 0}b_i$. Thus $e_2(A) = 4$ or $3$. \\
Subcase (1) $e_2(A) = 4$. In this case $ (\mathfrak{m} ^{i+1} \colon x) = \mathfrak{m} ^i$ for all $i$. In particular $x^*$ is $G(A)$-regular. So
$\operatorname{depth} G(A) \geq 1$. Also $h_A(z) = h_B(z) = 1 + hz + z^2 + z^3$.\\
Subcase (2) $e_2 = 3$. Note $e_2 = e_1 - e + 1 = 3$. Also $\operatorname{depth} G(A) = 0 $ as $(\mathfrak{m} ^{i+1} \colon x) \neq \mathfrak{m} ^i$ for some $i$. The result now follows from \ref{dim2}.

Now assume $d \geq 3$. Let $ \mathbf{x} = x_1,\ldots, x_{d-2}$ be an $A$-superficial sequence. Set $(B, \mathfrak{n} ) = (A/(\mathbf{x}), \mathfrak{m} /(\mathbf{x}))$. Then
$e_i(A) = e_i(B)$ for $i = 0, 1, 2$. The result now follows from our analysis of $d = 2$ case, Sally descent and Theorem \ref{e2}.
\end{proof}

\section{examples}
In this section we give examples which shows that there exists Cohen-Macaulay \ local rings with Hilbert functions as described in Theorems \ref{e2} and \ref{e1}.
We begin by the following
\s \label{const} \emph{Construction:}  Let $(A,\mathfrak{m} )$ be Cohen-Macaulay. Set $T = A[X_1,\ldots, X_m]$ and $\mathfrak{n}  = (\mathfrak{m} , X_1,\ldots, X_m)$. Set $R = T_\mathfrak{n} $.
Then $G(R) \cong G(A)[X_1^*, \ldots, X_m^*]$.

\begin{example} \label{dance}
 Let $C = k[X,Y,Z]/( X^2, XY, Y^2, (X,Y,Z)^3)$. Then $C$ is Artin local and $h_C(z) = 1 + 3z + 3z^2$. Using \ref{const} we can construct for every $m \geq 1$,  Cohen-Macaulay \ local rings $A_m$ of dimension $m$ with $h = 3$, $e= h + 4 = 7$ and $h_A(z) = h_C(z)$. Also $G(A_m)$ is  Cohen-Macaulay.
\end{example}

\begin{example} \label{e2d0}
The following example of Wang, see \cite[3.10]{CPR},  is of a 2-dimensional Cohen-Macaulay \ local ring with $e_2 = e_1 - e + 1 = 3$ and $\operatorname{depth} G(A)  = 0$.
$$A = k[[ x, y, z, u, v]]/(z^2, zu, zv, uv, yz - u^3,  xz-v^3).$$
We have $h_A(z) = 1 + 3z + 3z^3 - z^4$.
Using \ref{const} we may construct Cohen-Macaulay \ local rings $A_m$ of dimension $m + 2$
with $\operatorname{depth} G(A_m) = m$.
\end{example}

We note that Examples \ref{dance} and \ref{e2d0} yield Cohen-Macaulay \ local rings with Hilbert functions as described in Theorem \ref{e2}.

\begin{example}\label{sally1}
By \cite[4.4]{S4}
\[
R = k[[ t^6, t^7, t^{11}, t^{15}, t^{16}]]
\]
has multiplicity $6$, embedding dimension $5$ and reduction number $4$. So $h_R(z) = 1 + 4z + z^4$. Using \ref{const} we may construct
Cohen-Macaulay \ local rings $A_m$ of dimension $m + 1$ with $\operatorname{depth} G(A_m) = m$ and $h_{A_m}(z) = h_R(z)$.
\end{example}

\begin{example}\label{sally2}
By \cite{S3} we get that if $e = h + 3$ and $A$ is Gorenstein then $G(A)$ is Cohen-Macaulay \ and $h_A(z) = 1 + hz + z^2 + z^3$.
By \cite[p.\ 97]{S3};  $R = k[[X,Y,Z]]/(XZ-YZ,  XZ + Y^3 - Z^2)$  is Gorenstein of dimension one, multiplicity 5 and $h = 2$. So $h_R(z) = 1 + 2z + z^2 + z^3$.  Using \ref{const} we may construct
Cohen-Macaulay \ local rings $A_m$ of dimension $m + 1$ with $G(A_m)$ Cohen-Macaulay \   and $h_{A_m}(z) = h_R(z)$.
\end{example}

\begin{example}\label{tony}
Let $(A,\mathfrak{m} )$ be the two dimensional local ring as in \ref{e2d0}. We assume $k$ is infinite. Let $u$ be an $A$-superficial. Set $B = A/(u)$. Then by
\ref{dim2}(2)(f) we get that  $h_B(z) = 1 + 3z + z^2 + z^3$ and $\operatorname{depth} G(B) = 0$.  Using \ref{const} we may construct
Cohen-Macaulay \ local rings $A_m$ of dimension $m + 1$ with $\operatorname{depth} G(A_m) = m$    and $h_{A_m}(z) = h_B(z)$.
\end{example}

We note that by the examples constructed  yield Cohen-Macaulay \ local rings with Hilbert functions as described in Theorem \ref{e1}. Note we have also constructed Cohen-Macaulay \ local rings $A$ with $h_A(z) = 1 + hz + z^2 + z^3$ such that $G(A)$ is Cohen-Macaulay \ (see \ref{sally2}) and when $\operatorname{depth} G(A) = d - 1$ (see \ref{tony}).