# Stationarization#

The previous chapters derive all the equations necessary to solve for the steady-state and nonsteady-state equilibria of this model. However, because labor productivity is growing at rate $$g_y$$ as can be seen in the firms’ production function (40) and the population is growing at rate $$\tilde{g}_{n,t}$$ as defined in (5), the model is not stationary. Different endogenous variables of the model are growing at different rates. We have already specified three potential budget closure rules (91), (92), and (94) using some combination of government spending $$G_t$$ and transfers $$TR_t$$ that stationarize the debt-to-GDP ratio.

Table 3 lists the definitions of stationary versions of all the endogenous variables. Variables with a $$\:\,\hat{}\,\:$$’’ signify stationary variables. The first column of variables are growing at the productivity growth rate $$g_y$$. These variables are most closely associated with individual variables. The second column of variables are growing at the population growth rate $$\tilde{g}_{n,t}$$. These variables are most closely associated with population values. The third column of variables are growing at both the productivity growth rate $$g_y$$ and the population growth rate $$\tilde{g}_{n,t}$$. These variables are most closely associated with aggregate variables. The last column shows that the interest rates $$r_t$$, $$r_{p,t}$$ and $$r_{gov,t}$$, and household labor supply $$n_{j,s,t}$$ are already stationary.

Table 3 Stationary variable definitions. Note: The interest rate $$r_t$$ in firm first order condition is already stationary because $$Y_{m,t}$$ and $$K_{m,t}$$ grow at the same rate and $$p_{m,t}$$ is stationary. Household labor supply $$n_{j,s,t}\in[0,\tilde{l}]$$ is stationary.#

Sources of growth

$$e^{g_y t}$$

$$\tilde{N}_t$$

$$e^{g_y t}\tilde{N}_t$$

Not growing

$$\hat{b}_{j,s,t}\equiv \frac{b_{j,s,t}}{e^{g_y t}}$$

$$\hat{\omega}_{s,t}\equiv\frac{\omega_{s,t}}{\tilde{N}_t}$$

$$\hat{Y}_{m,t}\equiv\frac{Y_{m,t}}{e^{g_y t}\tilde{N}_t}$$

$$n_{j,s,t}$$

$$\hat{bq}_{j,s,t}\equiv \frac{bq_{j,s,t}}{e^{g_y t}}$$

$$\hat{L}_{m,t}\equiv\frac{L_{m,t}}{\tilde{N}_t}$$

$$\hat{K}_{m,t}\equiv\frac{K_{m,t}}{e^{g_y t}\tilde{N}_t}$$

$$r_t$$

$$\hat{c}_{j,s,t}\equiv \frac{c_{j,s,t}}{e^{g_y t}}$$

$$\hat{BQ}_{j,t}\equiv\frac{BQ_{j,t}}{e^{g_y t}\tilde{N}_t}$$

$$r_{p,t}$$

$$\hat{c}_{i,j,s,t}\equiv \frac{c_{i,j,s,t}}{e^{g_y t}}$$

$$\hat{C}_{i,t}\equiv\frac{C_{i,t}}{e^{g_y t}\tilde{N}_t}$$

$$r_{gov,t}$$

$$\hat{tr}_{j,s,t}\equiv \frac{tr_{j,s,t}}{e^{g_y t}}$$

$$\hat{K}_{g,m,t}\equiv\frac{K_{g,m,t}}{e^{g_y t}\tilde{N}_t}$$

$$r_{K,t}$$

$$\hat{ubi}_{j,s,t}\equiv\frac{ubi_{j,s,t}}{e^{g_y t}}$$

$$\hat{TR}_t\equiv\frac{TR_t}{e^{g_y t}\tilde{N}_t}$$

$$p_{i,t} \equiv \frac{\tilde{p}_{i,t}}{\tilde{p}_{M,t}}$$

$$\hat{T}_{j,s,t}\equiv \frac{T_{j,s,t}}{e^{g_y t}}$$

$$\hat{UBI}_t\equiv\frac{UBI_t}{e^{g_y t}\tilde{N}_t}$$

$$p_t \equiv \frac{\tilde{p}_t}{\tilde{p}_{M,t}}$$

$$\hat{w}_t\equiv \frac{w_t}{e^{g_y t}}$$

$$\hat{D}_t\equiv\frac{D_t}{e^{g_y t}\tilde{N}_t}$$

$$p_{m,t} \equiv \frac{\tilde{p}_{m,t}}{\tilde{p}_{M,t}}$$

The usual definition of equilibrium would be allocations and prices such that households optimize (19), (30), (31), and (32), firms optimize (43) and (44), and markets clear (105), (109), (110), (115), (116), and (118). In this chapter, we show how to stationarize each of these characterizing equations so that we can use our fixed point methods described in Sections Steady-state solution method and Stationary non-steady-state solution method of Chapter Equilibrium to solve for the equilibria in the steady-state and transition path equilibrium definitions.

## Stationarized Household Equations#

The stationary versions of the household industry-specific goods preferences and demand equations are obtained by dividing both sides of the equations by the productivity growth rate $$e^{g_y t}$$,

(119)#$\hat{c}_{j,s,t} \equiv \prod_{i=1}^I \left(\hat{c}_{i,j,s,t} - \hat{c}_{min,i,t}\right)^{\alpha_i} \quad\forall j,s,t \quad\text{with}\quad \sum_{i=1}^I\alpha_i=1$
(120)#$\hat{c}_{i,j,s,t} = \alpha_i\left(\frac{p_{i,t}}{p_t}\right)^{-1}\hat{c}_{j,s,t} + \hat{c}_{min,i,t} \quad\forall i,j,s,t$
(121)#$\begin{split} \hat{c}_{min,i,t} \equiv \begin{cases} \frac{c_{min,i}}{e^{g_y t}} \quad\text{for}\quad t < T \\ \frac{c_{min,i}}{e^{g_y T}} \quad\text{for}\quad t \geq T \end{cases} \quad\forall i\end{split}$

where (119) is the stationarized Stone-Geary consumption aggregator for composite consumption and (120) is the stationarized household demand for the composite consumption good. The composite price aggregation equation (20) is already stationary.

Note that the only way to stationarize the consumption aggregator (119) and consumption demand (120) is to divide $$c_{min,i}$$ by the growth rate $$e^{g_y t}$$. However, $$c_{min,i}$$ is already stationary. It is constant for each $$m$$. Therefore, the version of $$\hat{c}_{min,i,t}$$ divided by $$e^{g_y t}$$ would be changing over time (nonstationary) for $$g_y\neq 0$$. For this reason, we define $$\hat{c}_{min,i,t}$$ in (121) as being constant after the steady-state period $$T$$ at whatever value it reaches at that period. In most cases with $$g_y>0$$, that value will be close to zero. But we use $$\bar{c}_{min,i} = c_{min,i}/e^{g_y T}$$ from (121) as the steady-state value of $$c_{min,i}$$.

The stationary version of the household budget constraint (22) is found by dividing both sides of the equation by $$e^{g_y t}$$. For the savings term $$b_{j,s+1,t+1}$$, we must also multiply by $$e^{g_y(t+1)}$$, which leaves an $$e^{g_y} = \frac{e^{g_y(t+1)}}{e^{g_y t}}$$ in front of the stationarized variable.

(122)#$\begin{split} p_t\hat{c}_{j,s,t} + &\sum_{i=1}^I (1 + \tau^{c}_{i,t})p_{i,t}\hat{c}_{min,i} + e^{g_y}\hat{b}_{j,s+1,t+1} = \\ &(1 + r_{p,t})\hat{b}_{j,s,t} + \hat{w}_t e_{j,s} n_{j,s,t} + \\ &\quad\quad\zeta_{j,s}\frac{\hat{BQ}_t}{\lambda_j\hat{\omega}_{s,t}} + \eta_{j,s,t}\frac{\hat{TR}_{t}}{\lambda_j\hat{\omega}_{s,t}} + \hat{ubi}_{j,s,t} - \hat{T}_{j,s,t} \\ &\quad\forall j,t\quad\text{and}\quad s\geq E+1 \quad\text{where}\quad \hat{b}_{j,E+1,t}=0\quad\forall j,t\end{split}$

Because total bequests $$BQ_t$$ and total government transfers $$TR_t$$ grow at both the labor productivity growth rate and the population growth rate, we have to multiply and divide each of those terms by the economically relevant population $$\tilde{N}_t$$. This stationarizes total bequests $$\hat{BQ}_t$$, total transfers $$\hat{TR}_t$$, and the respective population level in the denominator $$\hat{\omega}_{s,t}$$.

We stationarize the Euler equations for labor supply (30) by dividing both sides by $$e^{g_y(1-\sigma)}$$. On the left-hand-side, $$e^{g_y}$$ stationarizes the wage $$\hat{w}_t$$ and $$e^{-\sigma g_y}$$ goes inside the parentheses and stationarizes consumption $$\hat{c}_{j,s,t}$$. On the right-and-side, the $$e^{g_y(1-\sigma)}$$ terms cancel out.

(123)#$\begin{split} &\frac{\hat{w}_t e_{j,s}}{p_t}\bigl(1 - \tau^{mtrx}_{s,t}\bigr)(\hat{c}_{j,s,t})^{-\sigma} = \chi^n_{s}\biggl(\frac{b}{\tilde{l}}\biggr)\biggl(\frac{n_{j,s,t}}{\tilde{l}}\biggr)^{\upsilon-1}\Biggl[1 - \biggl(\frac{n_{j,s,t}}{\tilde{l}}\biggr)^\upsilon\Biggr]^{\frac{1-\upsilon}{\upsilon}} \\ &\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\forall j,t, \quad\text{and}\quad E+1\leq s\leq E+S \\\end{split}$

We stationarize the Euler equations for savings (31) and (32) by dividing both sides of the respective equations by $$e^{-\sigma g_y t}$$. On the right-hand-side of the equation, we then need to multiply and divide both terms by $$e^{-\sigma g_y(t+1)}$$, which leaves a multiplicative coefficient $$e^{-\sigma g_y}$$.

(124)#$\begin{split} \frac{(\hat{c}_{j,s,t})^{-\sigma}}{p_t} &= e^{-\sigma g_y}\Biggl[\chi^b_j\rho_s(\hat{b}_{j,s+1,t+1})^{-\sigma} + \\ &\qquad\qquad\quad \beta_j\bigl(1 - \rho_s\bigr)\left(\frac{1 + r_{p,t+1}\bigl[1 - \tau^{mtry}_{s+1,t+1}\bigr] - \hat{MTR}^w_{j,s+1,t+1}}{p_{t+1}}\right)(\hat{c}_{j,s+1,t+1})^{-\sigma}\Biggr] \\ &\qquad\qquad\qquad\qquad\qquad\qquad\qquad\forall j,t, \quad\text{and}\quad E+1\leq s\leq E+S-1 \\\end{split}$
(125)#$\frac{(\hat{c}_{j,E+S,t})^{-\sigma}}{p_t} = e^{-\sigma g_y}\chi^b_j(\hat{b}_{j,E+S+1,t+1})^{-\sigma} \quad\forall j,t \quad\text{and}\quad s = E+S$

Where $$\hat{MTR}^w_{j,s,t} = \left( \frac{h^{w}p_{w}\hat{b}_{j,s,t}}{(\hat{b}_{j,s,t}h^{w}+m^{w})}\left[2 - \frac{h^{w}p_{w}\hat{b}_{j,s,t}}{(\hat{b}_{j,s,t}h^{w}+m^{w})}\right]\right)$$

## Stationarized Firms Equations#

The nonstationary production function (40) for each industry can be stationarized by dividing both sides by $$e^{g_y t}\tilde{N}$$. This stationarizes output $$\hat{Y}_{m,t}$$ on the left-hand-side. Because the general CES production function is homogeneous of degree 1, $$F(xK,xK_g,xL) = xF(K,K_g,L)$$, the right-hand-side of the production function is also stationarized by dividing by $$e^{g_y t}\tilde{N}_t$$.

(126)#$\begin{split} \begin{split} \hat{Y}_{m,t} &= F(\hat{K}_{m,t}, \hat{K}_{g,m,t}, \hat{L}_{m,t}) \\ &\equiv Z_{m,t}\biggl[(\gamma_m)^\frac{1}{\varepsilon_m}(\hat{K}_{m,t})^\frac{\varepsilon_m-1}{\varepsilon_m} + (\gamma_{g,m})^\frac{1}{\varepsilon_m}(\hat{K}_{g,m,t})^\frac{\varepsilon_m-1}{\varepsilon_m} + ... \\ &\qquad\qquad\qquad (1-\gamma_m-\gamma_{g,m})^\frac{1}{\varepsilon_m}(\hat{L}_{m,t})^\frac{\varepsilon_m-1}{\varepsilon_m}\biggr]^\frac{\varepsilon_m}{\varepsilon_m-1} \quad\forall m,t \end{split}\end{split}$

Notice that the growth term multiplied by the labor input drops out in this stationarized version of the production function. We stationarize the nonstationary profit function (42) in the same way, by dividing both sides by $$e^{g_y t}\tilde{N}_t$$.

(127)#$\begin{split} \hat{PR}_{m,t} &= (1 - \tau^{corp}_{m,t})\Bigl[F(\hat{K}_{m,t},\hat{K}_{g,m,t},\hat{L}_{m,t}) - \hat{w}_t \hat{L}_{m,t}\Bigr] - ... \\ &\qquad\qquad\quad \bigl(r_t + \delta_{M,t}\bigr)\hat{K}_{m,t} + \tau^{corp}_{m,t}\delta^\tau_{m,t}\hat{K}_{m,t} + \tau^{inv}_{m,t}\delta_{M,t}\hat{K}_{m,t} \quad\forall m,t\end{split}$

The firms’ first order equation for labor demand (43) is stationarized by dividing both sides by $$e^{g_y t}$$. This stationarizes the wage $$\hat{w}_t$$ on the left-hand-side and cancels out the $$e^{g_y t}$$ term in front of the right-hand-side. To complete the stationarization, we multiply and divide the $$\frac{Y_{m,t}}{e^{g_y t}L_{m,t}}$$ term on the right-hand-side by $$\tilde{N}_t$$.

(128)#$\hat{w}_t = p_{m,t}(Z_{m,t})^\frac{\varepsilon_m-1}{\varepsilon_m}\left[(1-\gamma_m-\gamma_{g,m})\frac{\hat{Y}_{m,t}}{\hat{L}_{m,t}}\right]^\frac{1}{\varepsilon_m} \quad\forall m,t$

It can be seen from the firms’ first order equation for capital demand (44) that the interest rate is already stationary. If we multiply and divide the $$\frac{Y_{m,t}}{K_{m,t}}$$ term on the right-hand-side by $$e^{g_y t}\tilde{N}_t$$, those two aggregate variables become stationary. In other words, $$Y_{m,t}$$ and $$K_{m,t}$$ grow at the same rate and $$\frac{Y_{m,t}}{K_{m,t}} = \frac{\hat{Y}_{m,t}}{\hat{K}_{m,t}}$$.

(129)#$r_t = (1 - \tau^{corp}_{m,t})p_{m,t}(Z_{m,t})^\frac{\varepsilon_m-1}{\varepsilon_m}\left[\gamma_m\frac{\hat{Y}_{m,t}}{\hat{K}_{m,t}}\right]^\frac{1}{\varepsilon_m} - \delta_{M,t} + \tau^{corp}_{m,t}\delta^\tau_{m,t} + \tau^{inv}_{m,t}\delta_{M,t} \quad\forall m,t$

A stationary version of the firms’ gross revenue attributed to each factor of production (45) is found by dividing both sides of the equation by $$e^{g_y t}\tilde{N}_t$$.

(130)#$\hat{Y}_{m,t} = MPK_{m,t}\hat{K}_{m,t} + MPK_{g,m,t}\hat{K}_{g,m,t} + \hat{MPL}_{m,t}\hat{L}_{m,t} \quad\forall m,t$

Note that this implies that both the marginal product of private capital $$MPK_{m,t}$$ and the marginal product of public capital $$MPK_{g,m,t}$$ are already stationary, as seen in (46) and (48). However, we see in (47) that the marginal product of labor is growing at rate $$e^{g_y t}$$ because of its relationship to the wage $$w_t$$. The division of both sides of (45) by $$e^{g_y t}\tilde{N}_t$$ gives us a stationarized marginal product of labor $$\hat{MPL}_{m,t}$$ and a stationarized labor demand $$\hat{L}_{m,t}$$.

Using the derivation of firm profits when firms are optimizing in (49) and the expressions for optimized stationary revenue (130), we can show the stationary equation for firm profits when firms are optimizing. As before, stationary profits are positive when stationary public capital is positive $$\hat{K}_{g,m,t}>0$$.

(131)#$\hat{PR}_{m,t} = (1 - \tau^{corp}_{m,t})p_{m,t}MPK_{g,m,t}\hat{K}_{g,m,t} \quad\forall m,t$

Using the derivation from (50) and (51) in Chapter Firms, we can stationarize the terms in the right-hand-side of the expression for $$r_{K,t}$$ by multiplying and dividing the quotient in the last term by $$e^{g_y t}\tilde{N}_t$$. This implies that the interest rate paid out by the financial intermediary on private capital $$r_{K,t}$$ is stationary, whether the variables on the right-hand-side are non-stationary in (51) or stationarized as in (132).

(132)#$r_{K,t} = r_t + \frac{\sum_{m=1}^M(1 - \tau^{corp}_{m,t})p_{m,t}MPK_{g,m,t}\hat{K}_{g,m,t}}{\sum_{m=1}^M\hat{K}_{m,t}} \quad\forall t$

## Stationarized Government Equations#

Each of the tax rate functions $$\tau^{etr}_{s,t}$$, $$\tau^{mtrx}_{s,t}$$, and $$\tau^{mtry}_{s,t}$$ is stationary. The total tax liability function $$T_{j,s,t}$$ is growing at the rate of labor productivity growth $$g_y$$ This can be see by looking at the decomposition of the total tax liability function into the effective tax rate times total income (23). The effective tax rate function is stationary, and household income is growing at rate $$g_y$$. So household total tax liability is stationarized by dividing both sides of the equation by $$e^{g_y t}$$.

(133)#$\begin{split} \hat{T}_{js,t} &= \tau^{etr}_{s,t}(\hat{x}_{j,s,t}, \hat{y}_{j,s,t})\left(\hat{x}_{j,s,t} + \hat{y}_{j,s,t}\right) \qquad\qquad\qquad\qquad\:\:\:\forall j,t \quad\text{and}\quad s\geq E+1 \\ &= \tau^{etr}_{s,t}(\hat{w}_t e_{j,s}n_{j,s,t}, r_{p,t}\hat{b}_{j,s,t})\left(\hat{w}_t e_{j,s}n_{j,s,t} + r_{p,t}\hat{b}_{j,s,t}\right) \quad\forall j,t \quad\text{and}\quad s\geq E+1\end{split}$

We can stationarize the simple expressions for total government spending on household transfers $$TR_t$$ in (75) and on public goods $$G_t$$ in (84) by dividing both sides by $$e^{g_y t}\tilde{N}_t$$,

(134)#$p_t \hat{Y}_t \equiv \sum_{m=1}^M p_{m,t} \hat{Y}_{m,t} \quad\forall t$
(135)#$\hat{TR}_t = g_{tr,t}\:\alpha_{tr}\: p_t \hat{Y}_t \quad\forall t$
(136)#$\hat{G}_t = g_{g,t}\:\alpha_{g}\: p_t \hat{Y}_t \quad\forall t$

where the time varying multipliers $$g_{g,t}$$ and $$g_{tr,t}$$, respectively, are defined in (145) and (146) below. These multipliers $$g_{g,t}$$ and $$g_{tr,t}$$ do not have a $$\:\,\hat{}\,\:$$’’ on them because their specifications (91) and (92) that are functions of nonstationary variables are equivalent to (145) and (146) specified in stationary variables.

We can stationarize the expression for total government revenue $$Rev_t$$ in (82) by dividing both sides of the equation by $$e^{g_y t}\tilde{N}_t$$.

(137)#$\begin{split} \hat{Rev}_t &= \underbrace{\sum_{m=1}^M\Bigl[\tau^{corp}_{m,t}\bigl(p_{m,t}\hat{Y}_{m,t} - \hat{w}_t\hat{L}_t\bigr) - \tau^{corp}_{m,t}\delta^\tau_{m,t}\hat{K}_{m,t} - \tau^{inv}_{m,t}\hat{I}_{m,t}\Bigr]}_{\text{corporate tax revenue}} \\ &\qquad + \underbrace{\sum_{s=E+1}^{E+S}\sum_{j=1}^J\lambda_j\hat{\omega}_{s,t}\tau^{etr}_{s,t}\left(\hat{x}_{j,s,t},\hat{y}_{j,s,t}\right)\bigl(\hat{x}_{j,s,t} + \hat{y}_{j,s,t}\bigr)}_{\text{household tax revenue}} \\ &\quad + \underbrace{\sum_{s=E+1}^{E+S}\sum_{j=1}^J\sum_{i=1}^I\lambda_j\omega_{s,t}\tau^{c}_{i,t}p{i,t}\hat{c}_{i,j,s,t}}_{\text{consumption tax revenue}} \\ &\quad + \underbrace{\sum_{s=E+1}^{E+S}\sum_{j=1}^J\lambda_j\omega_{s,t}\tau^{w}_{t}\hat{b}_{j,s,t}}_{\text{wealth tax revenue}} \quad\forall t\end{split}$

Every term in the government budget constraint (83) is growing at both the productivity growth rate and the population growth rate, so we stationarize it by dividing both sides by $$e^{g_y t}\tilde{N}_t$$. We also have to multiply and divide the next period debt term $$D_{t+1}$$ by $$e^{g_y(t+1)}\tilde{N}_{t+1}$$, leaving the term $$e^{g_y}(1 + \tilde{g}_{n,t+1})$$.

(138)#$e^{g_y}\left(1 + \tilde{g}_{n,t+1}\right)\hat{D}_{t+1} + \hat{Rev}_t = (1 + r_{gov,t})\hat{D}_t + \hat{G}_t + \hat{I}_{g,t} + \hat{Pensions}_t + \hat{TR}_t + \hat{UBI}_t \quad\forall t$

The stationarized versions of the rule for total government infrastructure investment spending $$I_{g,t}$$ in (85) and the rule for government investment spending in each industry in (85) are found by dividing both sides of the respective equations by $$e^{g_y t}\tilde{N}_t$$.

(139)#$\hat{I}_{g,t} = \alpha_{I,t}\: p_t\hat{Y}_t \quad\forall t$
(140)#$\hat{I}_{g,m,t} = \alpha_{I,m,t}\: \hat{I}_{g,t} \quad\forall m,t$

The stationarized version of the law of motion for the public capital stock in each industry $$K_{g,m,t}$$ in (87) is found by dividing both sides of the equation by $$e^{g_y t}\tilde{N}_t$$ then multiply and divide the $$K_{g,m,t+1}$$ term on the left-hand-side by $$e^{g_y(t+1)}\tilde{N}_{t+1}$$, leaving the term $$e^{g_y}(1 + \tilde{g}_{n,t+1})$$ in the denominator of the right-hand-side.

(141)#$\hat{K}_{g,m,t+1} = \frac{(1 - \delta_g)\hat{K}_{g,m,t} + \hat{I}_{g,m,t}}{e^{g_y}(1 + \tilde{g}_{n,t+1})} \quad\forall m,t$

Stationary aggregate universal basic income expenditure is found in one of two ways depending on how the individual UBI payments $$ubi_{j,s,t}$$ are modeled. In Section Universal basic income of Chapter Government, we discuss how UBI payments to households $$ubi_{j,s,t}$$ can be growth adjusted so that they grow over time at the rate of productivity growth or non-growth adjusted such that they are constant overtime. In the first case, when UBI benefits are growth adjusted and growing over time, the stationary aggregate government UBI payout $$\hat{UBI}_t$$ is found by dividing (88) by $$e^{g_y t}\tilde{N}_t$$. In the second case, when UBI benefits are constant over time and not growing with productivity, the stationary aggregate government UBI payout $$\hat{UBI}_t$$ is found by dividing (88) by only $$\tilde{N}_t$$.

(142)#$\begin{split} \hat{UBI}_t = \begin{cases} \sum_{s=E+1}^{E+S}\sum_{j=1}^J \lambda_j\hat{\omega}_{s,t} \hat{ubi}_{j,s,t} \quad\forall t \quad\text{if}\quad ubi_{j,s,t} \:\:\text{is growth adjusted} \\ \sum_{s=E+1}^{E+S}\sum_{j=1}^J \lambda_j\hat{\omega}_{s,t} ubi_{j,s,t} \quad\forall t \quad\text{if}\quad ubi_{j,s,t} \:\:\text{is not growth adjusted} \end{cases}\end{split}$

The expression for the interest rate on government debt $$r_{gov,t}$$ in (89) is already stationary because every term on the right-hand-side is already stationary. The net return on capital, $$r_{K,t}$$ is also stationary as shown in (132). The expression for the return to household savings $$r_{p,t}$$ in (104) is equivalent to its stationary representation because the same macroeconomic variables occur linearly in both the numerator and denominator.

(143)#$r_{p,t} = \frac{r_{gov,t}\hat{D}_{t} + r_{K,t}\hat{K}_{t}}{\hat{D}_{t} + \hat{K}_{t}} \quad\forall t \quad\text{where}\quad \hat{K}_t \equiv \sum_{m=1}^M \hat{K}_{m,t}$

The long-run debt-to-GDP ratio condition is also the same in both the nonstationary version in (90) as well as the stationary version below because the endogenous side is a ratio of macroeconomic variables that are growing at the same rate, with the exception of already stationary $$p_t$$.

(144)#$\frac{\hat{D}_t}{p_t\hat{Y}_t} = \alpha_D \quad\text{for}\quad t\geq T$

The three potential budget closure rules (91), (92), and (94) are the last government equations to stationarize. In each of the cases, we simply divide both sides by $$e^{g_y t}\tilde{N}_t$$.

(145)#$\begin{split}\begin{split} &\hat{G}_t = g_{g,t}\:\alpha_{g}\: p_t\hat{Y}_t \\ &\text{where}\quad g_{g,t} = \begin{cases} 1 \qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\:\:\text{if}\quad t < T_{G1} \\ \frac{e^{g_y}\left(1 + \tilde{g}_{n,t+1}\right)\left[\rho_{d}\alpha_{D}p_t\hat{Y}_{t} + (1-\rho_{d})\hat{D}_{t}\right] - (1+r_{gov,t})\hat{D}_{t} - \hat{TR}_{t} - \hat{I}_{g,t} - \hat{UBI}_t + \hat{Rev}_{t}}{\alpha_g p_t\hat{Y}_t} \quad\text{if}\quad T_{G1}\leq t<T_{G2} \\ \frac{e^{g_y}\left(1 + \tilde{g}_{n,t+1}\right)\alpha_{D}p_t\hat{Y}_{t} - (1+r_{gov,t})\hat{D}_{t} - \hat{TR}_{t} - \hat{I}_{g,t} - \hat{UBI}_t + \hat{Rev}_{t}}{\alpha_g p_t\hat{Y}_t} \qquad\qquad\quad\,\text{if}\quad t \geq T_{G2} \end{cases} \\ &\text{and}\quad g_{tr,t} = 1 \quad\forall t \end{split}\end{split}$

or

(146)#$\begin{split}\begin{split} &\hat{TR}_t = g_{tr,t}\:\alpha_{tr}\: p_t\hat{Y}_t \\ &\text{where}\quad g_{tr,t} = \begin{cases} 1 \qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\text{if}\quad t < T_{G1} \\ \frac{e^{g_y}\left(1 + \tilde{g}_{n,t+1}\right)\left[\rho_{d}\alpha_{D}p_t\hat{Y}_{t} + (1-\rho_{d})\hat{D}_{t}\right] - (1+r_{gov,t})\hat{D}_{t} - \hat{G}_{t} - \hat{I}_{g,t} - \hat{UBI}_t + \hat{Rev}_{t}}{\alpha_{tr} p_t\hat{Y}_t} \quad\text{if}\quad T_{G1}\leq t<T_{G2} \\ \frac{e^{g_y}\left(1 + \tilde{g}_{n,t+1}\right)\alpha_{D}p_t\hat{Y}_{t} - (1+r_{gov,t})\hat{D}_{t} - \hat{G}_{t} - \hat{I}_{g,t} - \hat{UBI}_t + \hat{Rev}_{t}}{\alpha_{tr} p_t\hat{Y}_t} \qquad\qquad\quad\,\text{if}\quad t \geq T_{G2} \end{cases} \\ &\text{and}\quad g_{g,t} = 1 \quad\forall t \end{split}\end{split}$

or

(147)#$\begin{split}\begin{split} &\hat{G}_t + \hat{TR}_t = g_{trg,t}\left(\alpha_g + \alpha_{tr}\right)p_t\hat{Y}_t \quad\Rightarrow\quad \hat{G}_t = g_{trg,t}\:\alpha_g\:p_t\hat{Y}_t \quad\text{and}\quad \hat{TR}_t = g_{trg,t}\:\alpha_{tr}\:p_t\hat{Y}_t \\ &\text{where}\quad g_{trg,t} = \begin{cases} 1 \qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\quad\text{if}\quad t < T_{G1} \\ \frac{e^{g_y}\left(1 + \tilde{g}_{n,t+1}\right)\left[\rho_{d}\alpha_{D}p_t\hat{Y}_{t} + (1-\rho_{d})\hat{D}_{t}\right] - (1+r_{gov,t})\hat{D}_{t} - \hat{I}_{g,t} - \hat{UBI}_t + \hat{Rev}_{t}}{\left(\alpha_g + \alpha_{tr}\right)p_t\hat{Y}_t} \quad\text{if}\quad T_{G1}\leq t<T_{G2} \\ \frac{e^{g_y}\left(1 + \tilde{g}_{n,t+1}\right)\alpha_{D}p_t\hat{Y}_{t} - (1+r_{gov,t})\hat{D}_{t} - \hat{I}_{g,t} - \hat{UBI}_t + \hat{Rev}_{t}}{\left(\alpha_g + \alpha_{tr}\right)p_t\hat{Y}_t} \qquad\qquad\quad\,\text{if}\quad t \geq T_{G2} \end{cases} \end{split}\end{split}$

### Stationarized Pension System Equations#

#### Stiationarized Notional Defined Contributions Equations#

The stationarized NDC pension amount is given by:

(148)#$\begin{split}\hat{\theta}_{j,u,t+u-s}= \begin{cases} 0, & \text{if}\ u < R \\ \biggl[\sum_{s=E}^{R-1}\tau^{p}_{t}\frac{\hat{w}_{t}}{e^{g_{y}(u-s)}}e_{j,s}n_{j,s,t}(1 + g_{NDC,t})^{R-s-1}\biggr]\delta_{R, t}, & \text{otherwise} \end{cases}\end{split}$

The stationarized derivative of the pension amount it slightly simpler since it involved only current period wages. We give the derivation first.

The FOC for the choice of labor supply is given by:

$\begin{split}\begin{split} \biggl(\frac{1}{1 + \tau^{c}_{j,s,t}}\biggr)&\biggl(w_{t}e_{j,s} - \frac{\partial T_{j,s,t}}{\partial n_{j,s,t}}\biggr)c^{-\sigma}_{j,s,t} + \sum_{u=R}^{E+S}\beta^{u-s}\prod_{v=s}^{u}(1-\rho_{v})\frac{\partial \theta_{j,u,t+u-s}}{\partial n_{j,s,t}}c^{-\sigma}_{j,u,t+u-s}\biggl(\frac{1}{1+\tau^{c}_{j,u,t+u-s}}\biggr) \\ & = MDU_l(n_{j,s,t})e^{g_{y}t(1-\sigma)} \end{split}\end{split}$

where we now pull the growth factor out of the marginal disutility of labor term to aid in the exposition of the stationarization. To stationarize this equation, we divide both sides through by $$e^{g_{y}t(1-\sigma)}$$.

$\begin{split}\begin{split} \biggl(\frac{1}{1 + \tau^{c}_{j,s,t}}\biggr)&\biggl(\frac{w_{t}}{e^{g_{y}t}}e_{j,s} - \frac{\partial T_{j,s,t}}{\partial n_{j,s,t}}\biggr)\frac{c^{-\sigma}_{j,s,t}}{e^{g_{y}t(-\sigma)}} + \sum_{u=R}^{E+S}\beta^{u-s}\prod_{v=s}^{u}(1-\rho_{v})\frac{\partial \theta_{j,u,t+u-s}}{\partial n_{j,s,t}e^{g_{y}t}}\frac{c^{-\sigma}_{j,u,t+u-s}}{e^{g_{y}t(-\sigma)}}\biggl(\frac{1}{1+\tau^{c}_{j,u,t+u-s}}\biggr) \\ & = MDU_l(n_{j,s,t})\frac{e^{g_{y}t(1-\sigma)}}{e^{g_{y}t(1-\sigma)}} \end{split}\end{split}$

Which we can write as:

$\begin{split}\begin{split} \biggl(\frac{1}{1 + \tau^{c}_{j,s,t}}\biggr)&\biggl(\hat{w}_{t}e_{j,s} - \frac{\partial T_{j,s,t}}{\partial n_{j,s,t}}\biggr)\hat{c}^{-\sigma}_{j,s,t} + \sum_{u=R}^{E+S}\beta^{u-s}\prod_{v=s}^{u}(1-\rho_{v})\frac{\partial \hat{\theta}_{j,u,t+u-s}}{\partial n_{j,s,t}}\hat{c}^{-\sigma}_{j,u,t+u-s}e^{g_{y}(u-s)(-\sigma)}\biggl(\frac{1}{1+\tau^{c}_{j,u,t+u-s}}\biggr) \\ & = MDU_l(n_{j,s,t}) \end{split}\end{split}$

where $$\frac{\partial \hat{\theta}_{j,u,t+u-s}}{\partial n_{j,s,t}}$$ is given by:

(149)#$\begin{split}\frac{\partial \theta_{j,u,t+u-s}}{\partial n_{j,s,t}} = \begin{cases} \tau^{p}_{t}\hat{w}_{t}e_{j,s}(1+g_{NDC,t})^{u - s}\delta_{R,t}, & \text{if}\ s<R-1 \\ 0, & \text{if}\ s \geq R \\ \end{cases}\end{split}$

#### Stationarized Defined Benefits Equations#

Stationarized pension amount:

(150)#$\hat{\theta}_{j,u,t+u-s} = \biggl[\frac{\sum_{s=R-ny}^{R-1}\frac{\hat{w}_{t}}{e^{g_{y}(u-s)}}e_{j,s}n_{j,s,t}}{ny}\biggr]\times Cy \times \alpha_{DB}, \ \ \forall u \geq R$

Stationarized pension amount derivative:

(151)#$\begin{split} \frac{\partial \hat{\theta}_{j,u,t+u-s}}{\partial n_{j,s,t}} = \begin{cases} 0 , & \text{if}\ s < R - Cy \\ \hat{w}_{t}e_{j,s}\alpha_{DB}\times \frac{Cy}{ny}, & \text{if}\ R - Cy <= s < R \\ 0, & \text{if}\ s \geq R \\ \end{cases}\end{split}$

#### Stationarized Points System Equations#

Stationarized pension amount:

(152)#$\hat{\theta}_{j,u,t+u-s} =\sum_{s=E}^{R-1}\frac{\hat{w}_{t}}{e^{g_{y}(u-s)}}e_{j,s}n_{j,s,t}v_{t}, \ \ \forall u \geq R$

Stationarized pension amount derivative:

(153)#$\begin{split}\frac{\partial \hat{\theta}_{j,u,t+u-s}}{\partial n_{j,s,t}} = \begin{cases} \hat{w}_{t}e_{j,s} v_{t}, & \text{if}\ s < R \\ 0, & \text{if}\ s \geq R \\ \end{cases}\end{split}$

## Stationarized Market Clearing Equations#

The labor market clearing equation (105) is stationarized by dividing both sides by $$\tilde{N}_t$$.

(154)#$\sum_{m=1}^M \hat{L}_{m,t} = \sum_{s=E+1}^{E+S}\sum_{j=1}^{J} \hat{\omega}_{s,t}\lambda_j e_{j,s}n_{j,s,t} \quad \forall t$

Total savings by domestic households $$B_t$$ from (106) is stationarized by dividing both sides by $$e^{g_y t}\tilde{N}_t$$. The $$\omega_{s,t-1}$$ terms on the right-hand_side require multiplying and dividing by $$\tilde{N}_{t-1}$$, which leads to the division of $$1 + \tilde{g}_{n,t}$$.

(155)#$\hat{B}_t \equiv \frac{1}{1 + \tilde{g}_{n,t}}\sum_{s=E+2}^{E+S+1}\sum_{j=1}^{J}\Bigl(\hat{\omega}_{s-1,t-1}\lambda_j b_{j,s,t} + i_s\hat{\omega}_{s,t-1}\lambda_j \hat{b}_{j,s,t}\Bigr) \quad \forall t$

And the total domestic savings constraint (107) is stationarized by dividing both sides by $$e^{g_y t}\tilde{N}_t$$.

(156)#$\hat{K}^d_t + \hat{D}^d_t = \hat{B}_t \quad \forall t$

The stationarized law of motion for foreign holdings of government debt (108) and the government debt market clearing condition (109), respectively, are solved for by dividing both sides by $$e^{g_y t}\tilde{N}_t$$.

(157)#$e^{g_y}\bigl[1 + \tilde{g}_{n,t+1}\bigr]\hat{D}^{f}_{t+1} = \hat{D}^{f}_{t} + \zeta_{D}\Bigl(e^{g_y}\bigl[1 + \tilde{g}_{n,t+1}\bigr]\hat{D}_{t+1} - \hat{D}_{t}\Bigr) \quad\forall t$
(158)#$\hat{D}_t = \hat{D}^d_t + \hat{D}^f_t \quad\forall t$

The private capital market clearing equation (110) is stationarized by dividing both sides by $$e^{g_y t}\tilde{N}_t$$, as is the expression for excess demand at the world interest rate (111) and the exogenous expression for foreign private capital flows (112).

(159)#$\hat{K}_t = \hat{K}^d_t + \hat{K}^f_t \quad\forall t \quad\text{where}\quad \hat{K_t} \equiv \sum_{m=1}^M \hat{K}_{m,t}$
(160)#$\hat{ED}^{K,r^*}_t \equiv \hat{K}^{r^*}_t - \hat{K}^d_t \quad\forall t \quad\text{where}\quad \hat{K}^{r^*}_t \equiv \sum_{m=1}^M \hat{K}^{r^*}_{m,t}$
(161)#$\hat{K}^{f}_t = \zeta_{K}\hat{ED}^{K,r^*}_t \quad\forall t$

We stationarize the goods market clearing equations for the first $$M-1$$ industries (115) and for the $$M$$th industry (116) by dividing both sides by $$e^{g_y t}\tilde{N}_t$$. On the right-hand-side, we must multiply and divide the $$K^d_{t+1}$$ term and the $$D^f_{t+1}$$ term, respectively, by $$e^{g_y(t+1)}\tilde{N}_{t+1}$$ leaving the coefficient $$e^{g_y}(1+\tilde{g}_{n,t+1})$$.

(162)#$\hat{Y}_{m,t} = \hat{C}_{m,t} \quad\forall t \quad\text{and}\quad m=1,2,...M-1$
(163)#$\begin{split} \hat{Y}_{M,t} &= \hat{C}_{M,t} + \hat{I}_{M,t} + \hat{I}_{g,t} + \hat{G}_t + r_{p,t} \hat{K}^f_t + r_{p,t}\hat{D}^f_t ... \\ &\quad - \Bigl(e^{g_y}\bigl[1 + \tilde{g}_{n,t+1}\bigr]\hat{K}^f_{t+1} - \hat{K}^f_t\Bigr) - \Bigl(e^{g_y}\bigl[1 + \tilde{g}_{n,t+1}\bigr]\hat{D}^f_{t+1} - \hat{D}^f_t\Bigr) \quad\forall t\end{split}$

where

(164)#$\hat{C}_{m,t} \equiv \sum_{i=1}^{I}\sum_{s=E+1}^{E+S}\sum_{j=1}^{J}\hat{\omega}_{s,t}\lambda_j \pi_{i,m} \hat{c}_{i,j,s,t} \quad\forall m,t$

and

(165)#$\begin{split} \hat{I}_{M,t} &\equiv e^{g_y}\bigl(1 + \tilde{g}_{n,t+1}\bigr)\sum_{m=1}^M \hat{K}_{m,t+1} - (1 - \delta_{M,t})\sum_{m=1}^M \hat{K}_{m,t} \quad\forall t \\ &= e^{g_y}\bigl(1 + \tilde{g}_{n,t+1}\bigr)\hat{K}_{t+1} - (1 - \delta_{M,t})\hat{K}_t \\ &= e^{g_y}\bigl(1 + \tilde{g}_{n,t+1}\bigr)(\hat{K}^d_{t+1} + \hat{K}^f_{t+1}) - (1 - \delta_{M,t})(\hat{K}^d_t + \hat{K}^f_t)\end{split}$

We stationarize the law of motion for total bequests $$BQ_t$$ in (118) by dividing both sides by $$e^{g_y t}\tilde{N}_t$$. Because the population levels in the summation are from period $$t-1$$, we must multiply and divide the summed term by $$\tilde{N}_{t-1}$$ leaving the term in the denominator of $$1+\tilde{g}_{n,t}$$.

(166)#$\hat{BQ}_{t} = \left(\frac{1+r_{p,t}}{1 + \tilde{g}_{n,t}}\right)\left(\sum_{s=E+2}^{E+S+1}\sum_{j=1}^J\rho_{s-1}\lambda_j\hat{\omega}_{s-1,t-1}\hat{b}_{j,s,t}\right) \quad\forall t$