Tobin’s q and Investment Dynamics
“Tobin’s q is the ratio of the market value of capital to its replacement cost. When q > 1, install more; when q < 1, let capital depreciate.” — James Tobin
Cross-reference: Principles Ch. 12 (investment theory: Jorgenson, Tobin’s q, real options, adjustment costs); Ch. 24.2 (financial accelerator and Hayashi’s theorem) [P:Ch.12, P:Ch.24.2]
13.1 Why Adjustment Costs? Motivating the Model¶
The neoclassical investment theory of Jorgenson (1963) [P:Ch.12.1] determines the optimal capital stock by equating the marginal product of capital to the user cost: . In this frictionless framework, the firm immediately jumps to the optimal capital stock in response to any shock. Observed investment is then simply the instantaneous jump — an impulse function.
This prediction is empirically implausible. Investment data show:
Investment responds gradually to shocks, not in a single jump.
Firms do not continuously adjust capital; there are periods of high investment followed by consolidation.
The correlation between Tobin’s and investment is positive but loose — firms do not instantly equate to 1.
The adjustment cost model explains these patterns by introducing costs of rapid capital accumulation: installing new capital requires not just purchasing it but also disrupting ongoing production, training workers, and reorganizing the production process. These installation costs make it optimal to spread investment over time rather than adjusting instantaneously.
Formally, we assume a strictly convex adjustment cost function that increases in (higher investment is more expensive) and is homogeneous of degree one in — a standard specification that leads to the elegant results of this chapter.
13.2 The Firm’s Optimal Control Problem¶
13.2.1 Setup¶
A representative firm owns capital and chooses investment to maximize the present value of profits. The firm’s optimization problem:
subject to: , given.
where:
is operating profit (maximized over labor), an increasing concave function of .
is gross investment (at purchase cost normalized to 1 per unit).
is the quadratic adjustment cost function (Hayashi, 1982 specification), with the adjustment cost parameter.
is the (constant) discount rate.
is the depreciation rate.
Definition 13.1 (Quadratic Adjustment Costs). The adjustment cost function has the properties:
: zero investment has zero adjustment cost.
: marginal adjustment cost is positive and increasing in the investment rate .
: the marginal adjustment cost is increasing in (strictly convex).
Homogeneous of degree 1 in : .
The homogeneity assumption is crucial for Hayashi’s theorem (Section 13.5).
13.2.2 The Current-Value Hamiltonian¶
Apply the optimal control machinery of Chapter 11. The state variable is ; the control variable is ; the costate variable is the shadow price of an additional unit of capital.
Definition 13.2 (Tobin’s q as Costate Variable). The marginal q is the shadow value of an installed unit of capital — the costate variable of the firm’s optimal control problem:
where is the firm’s value function.
The current-value Hamiltonian:
13.2.3 Pontryagin Conditions¶
Condition 1 (FOC with respect to ):
This immediately gives the investment function:
Investment is zero when , positive when , and negative (disinvestment) when . The parameter controls the sensitivity of the investment rate to : high (large adjustment costs) means investment responds slowly to deviations.
Condition 2 (Costate equation for ):
Wait — let us compute carefully:
Therefore:
Substituting the investment function :
Condition 3 (State equation):
Condition 4 (Transversality):
13.3 The Investment Function and Its Interpretation¶
The investment function is the central result of the model. It has a clean economic interpretation:
When : The shadow price of installed capital exceeds its replacement cost. The firm can increase its value by investing — each dollar of investment generates more than one dollar of firm value. It is optimal to invest until capital accumulation drives down (diminishing returns), reducing back toward 1.
When : The shadow price of installed capital is below its replacement cost. The firm’s installed capital is worth less than it costs to replace. It should allow capital to depreciate without replacement (disinvestment at rate ) until capital falls and rises enough to bring back to 1.
When : The firm is at its optimal capital stock. Net investment is zero; gross investment exactly covers depreciation.
Definition 13.3 (Tobin’s Average q and Marginal q). Average q is the ratio of the firm’s market value to the replacement cost of its capital stock: . Marginal q is the shadow price . These two concepts differ in general, but Hayashi’s theorem (Section 13.5) identifies conditions under which they are equal.
13.4 The Phase Plane¶
The two-equation system forms a 2D autonomous system amenable to phase-plane analysis.
13.4.1 Nullclines¶
nullcline: . This is a horizontal line at . (When , the nullcline is .) Above this line, implies , so is growing. Below, is shrinking.
nullcline: . This is a curve in space. For the steady state :
For : — the standard Jorgenson condition (MPK = user cost) is recovered as a special case when adjustment costs do not interact with depreciation.
13.4.2 Steady State and Saddle-Point Property¶
Theorem 13.1 (Saddle-Point Property of the System). Under standard regularity conditions (), the steady state of the system is a saddle point.
Proof. The Jacobian of with respect to at the steady state:
At the steady state: , so the (1,1) entry is and the (2,2) entry is .
Hmm — let me be careful with signs. , so . Then:
since , , and ... this gives . ✓ The determinant is negative, confirming the saddle-point property.
13.4.3 Impulse Responses¶
The saddle-path structure means that for any shock to (e.g., a TFP improvement that raises at every ), the adjustment follows the saddle path:
Permanent TFP shock: A permanent increase in TFP shifts the nullcline upward (higher at every ), raising the steady-state capital . On impact, jumps up (the shadow value of capital increases) and then gradually declines along the new saddle path as accumulates toward its new steady state. Investment is high throughout the transition and declines as approaches .
Temporary TFP shock: A temporary shock generates a smaller initial jump in (since the higher profits are temporary and the old steady state will eventually be restored) and a smaller investment response. This is a key distinction: the model separates permanent from temporary shocks through their differential impact on the shadow price.
13.5 Hayashi’s Theorem¶
The empirical implementation of Tobin’s theory requires measuring . Marginal (the model’s key variable) is not directly observable. Average (market capitalization divided by replacement cost) is observable. Hayashi (1982) identifies conditions under which the two are equal.
Theorem 13.2 (Hayashi’s Theorem). If: (1) the production function is homogeneous of degree one (constant returns to scale), and (2) the adjustment cost function is homogeneous of degree one in , then marginal equals average :
where is the market value of the firm and is the replacement cost of capital.
Proof. Under CRS in production and HD1 adjustment costs, the firm’s value function is homogeneous of degree one in : . Euler’s theorem for homogeneous functions gives: . Therefore .
Importance: Hayashi’s theorem means we can test the model using stock market data. If the model is correct, the investment rate should be a linear function of Tobin’s average (the equity market capitalization plus debt, divided by replacement cost):
This is an empirically testable linear regression. However, the empirical of this regression is typically quite low (around 0.1–0.2), and explains much less of investment variation than the theory predicts. Proposed explanations: measurement error in (stock prices are noisy); non-convex adjustment costs (lumpy investment); financial frictions (the financial accelerator, [P:Ch.24.2]) that drive a wedge between marginal and the observable .
13.6 Real Options: The Investment Trigger¶
The adjustment cost model treats investment as continuously reversible. In reality, many investment projects are irreversible (once a factory is built, it cannot be un-built without large losses). Real options theory analyzes investment timing when investment is irreversible and payoffs are uncertain.
Definition 13.4 (Real Option). A real option is the right, but not the obligation, to undertake an investment project at a future date. Like a financial call option, a real option has value even when the immediate NPV is negative: waiting preserves the option to invest under more favorable conditions.
13.6.1 Geometric Brownian Motion for Profits¶
Let operating profit follow geometric Brownian motion (GBM):
where is the drift, is the volatility, and is a Wiener process increment.
Itô’s Lemma: For any twice-differentiable function :
13.6.2 The Option Value and the Investment Trigger¶
The firm waits to invest until profits exceed a threshold (the investment trigger). Below , the option value of waiting exceeds the NPV of immediate investment. Above , immediate investment is optimal.
For a perpetual project costing with value (Gordon growth formula):
The option value satisfies the Bellman ODE (derived from Itô’s Lemma and the no-arbitrage condition ):
This is an Euler–Cauchy ODE with power-function solutions . The characteristic equation:
The positive root (required for the option value to be finite as ):
Boundary conditions:
Value matching: (option value equals project value minus cost at the trigger).
Smooth pasting: (optimality condition: option and project value curves are tangent at trigger).
Solving these two conditions for :
Definition 13.5 (Investment Trigger). The optimal investment trigger is:
where is the option value multiplier — always greater than one, reflecting the value of waiting.
Key comparative statics:
Higher (more uncertainty): decreases (the characteristic equation root moves), making larger — more uncertainty raises the trigger and delays investment.
Higher (higher discount rate): rises — higher opportunity cost of waiting raises the trigger.
Higher (faster expected profit growth): falls — faster growth makes the option less valuable relative to the project value.
This is the Bloom (2009) uncertainty channel: higher uncertainty (higher ) raises investment thresholds, causing firms to delay investment during uncertain times. This mechanism explains the sharp investment contraction during the 2008–09 financial crisis and the COVID-19 recession [P:Ch.12.4, P:Ch.40].
13.7 Worked Example: Investment Response to a TFP Shock¶
Cross-reference: Principles Ch. 12.1 (investment empirics, q vs. data) [P:Ch.12.1]
Calibration: , , , with (decreasing returns to K in operating profit — firms are not pure CRS in capital), .
Steady state (baseline):
(investment rate at steady state: ... this suggests I/K = δ at SS, so q*=1+ψδ=1.20 gives I/K = 0.10 = δ. ✓)
(for small corrections)
Permanent 10% TFP shock ():
New MPK condition:
On impact: jumps up. The new nullcline requires a higher , achievable only at lower — but cannot jump. So jumps to the value that puts the economy on the new saddle path at .
APL¶
⎕IO←0
⍝ 1. Parameters
r ← 0.05 ⋄ delta ← 0.10 ⋄ psi ← 2 ⋄ alpha_pi ← 0.7
⍝ 2. Operating profit derivative ( Marginal Product of Capital )
pi_Kp ← {alpha_pi × ⍵ * alpha_pi - 1}
⍝ 3. Steady state
qstar ← 1 + psi × delta
Kstar ← ((r + delta) ÷ alpha_pi) * ÷ alpha_pi - 1
⍝ 4. Integrated ODE System & RK4 Solver
rk4_step ← {
h state ← ⍺ ⍵
sys ← {
K q ← ⍵
dK ← K × (((q - 1) ÷ psi) - delta)
⍝ FIXED: Enforced (A - B) - C to prevent right-to-left sign flipping
⍝ Replaced *2 with explicit multiplication for domain safety
dq ← (((r + delta) × q) - (pi_Kp K)) - ((q - 1) × (q - 1)) ÷ (2 × psi)
⍝ Return as a flat vector
dK , dq
}
k1 ← sys state
k2 ← sys state + (h ÷ 2) × k1
k3 ← sys state + (h ÷ 2) × k2
k4 ← sys state + h × k3
state + (h ÷ 6) × (k1 + (2 × k2) + (2 × k3) + k4)
}
⍝ 5. Setup and Run
K0 ← Kstar × 0.7
q0 ← 1.4
h ← 0.1
T ← 200
⍝ Memory-Safe Accumulator
shoot ← {
path ← ⍵
next ← h rk4_step ⊃⌽path
path , ⊂next
}
path_states ← shoot⍣T ⊢ ⊂K0 q0
⍝ 6. Extract and Format Output
K_path ← {⍵[0]}¨ path_states
q_path ← {⍵[1]}¨ path_states
'--- Tobin q Model Steady State ---'
↑('Capital (K*)' 'Tobin q (q*)') ,¨ (⊂' : ') ,¨ ⍕¨ 4 ⍕¨ Kstar qstar
'--- Convergence Verification (End of Path) ---'
↑('Final K' 'Final q') ,¨ (⊂' : ') ,¨ ⍕¨ 4 ⍕¨ (⊃⌽K_path) (⊃⌽q_path)Python¶
import numpy as np
from scipy.integrate import solve_ivp
from scipy.optimize import brentq
import matplotlib.pyplot as plt
# Parameters
r, delta, psi, alpha_pi = 0.05, 0.10, 2.0, 0.70
# Use np.maximum to prevent K from going negative and causing power errors
def pi_Kp(K, A=1.0):
K_safe = np.maximum(K, 1e-6)
return A * alpha_pi * K_safe**(alpha_pi - 1)
# 1. Steady State
qstar = 1 + psi*delta
Rstar = (r + delta)*qstar - (qstar - 1)**2 / (2*psi)
Kstar = (alpha_pi / Rstar)**(1 / (1 - alpha_pi))
def Kq_system(t, state, A=1.0):
K, q = state
K_safe = np.maximum(K, 1e-6) # Guard against power errors
dK = K_safe*(q-1)/psi - delta*K_safe
dq = (r+delta)*q - pi_Kp(K_safe, A) - (q-1)**2/(2*psi)
return [dK, dq]
# Event to stop integration if K hits 0 or grows too large
def hit_boundary(t, state, A=1.0):
return state[0] - 1e-3
hit_boundary.terminal = True
# 2. Saddle Path via Eigenvectors
# J = [[dK_dot/dK, dK_dot/dq], [dq_dot/dK, dq_dot/dq]]
# dK_dot/dK = (q-1)/psi - delta. At SS, this is 0.
J = np.array([[0, Kstar/psi],
[alpha_pi*(1-alpha_pi)*Kstar**(alpha_pi-2), r]])
eigvals, eigvecs = np.linalg.eig(J)
stable_idx = np.argmin(eigvals)
v = eigvecs[:, stable_idx]
if v[0] < 0: v = -v # Ensure it points right
# 3. Solve Saddle Path (Integrating backward from SS)
eps = 1e-3 # Small perturbation for shooting
back_time = (0, 30) # Reduced time to prevent exploding
sol_r = solve_ivp(lambda t,s: [-x for x in Kq_system(t,s)], back_time,
[Kstar + eps*v[0], qstar + eps*v[1]], events=hit_boundary)
sol_l = solve_ivp(lambda t,s: [-x for x in Kq_system(t,s)], back_time,
[Kstar - eps*v[0], qstar - eps*v[1]], events=hit_boundary)
# 4. TFP Shock (A: 1.0 -> 1.1)
A_new = 1.1
Kstar_new = (alpha_pi*A_new / Rstar)**(1 / (1 - alpha_pi))
# Find the jump: Integrate backward from the NEW steady state to find where it crosses OLD K*
sol_shock = solve_ivp(lambda t,s: [-x for x in Kq_system(t,s, A_new)], (0, 50),
[Kstar_new + eps*v[0], qstar + eps*v[1]],
dense_output=True, events=hit_boundary)
# Find q0 by looking for when K(t) in backward integration reached the original Kstar
try:
t_jump = brentq(lambda t: sol_shock.sol(t)[0] - Kstar, 0, sol_shock.t[-1])
q_jump = sol_shock.sol(t_jump)[1]
except:
q_jump = qstar + 0.05 # Rough fallback
# Forward simulation of the shock
sol_path = solve_ivp(Kq_system, (0, 100), [Kstar, q_jump], args=(A_new,), dense_output=True)
print("--- PERMANENT TFP SHOCK (A: 1.0 -> 1.1) ---")
print(f"New Steady State K*: {Kstar_new:.4f}")
print(f"Initial q jump: {q_jump:.4f}")
print(f"Total q increase: {((q_jump - qstar)/qstar)*100:.2f}%")
print("-" * 30)
# --- Visualization ---
fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(14, 5))
# Phase Diagram
ax1.plot(sol_r.y[0], sol_r.y[1], 'g-', label='Saddle Path')
ax1.plot(sol_l.y[0], sol_l.y[1], 'g-')
ax1.axhline(qstar, color='k', linestyle=':', label='K-nullcline')
ax1.plot(Kstar, qstar, 'ko', label='Old SS')
ax1.plot(Kstar_new, qstar, 'ro', label='New SS')
ax1.set_title("Phase Diagram: Tobin's q")
ax1.set_xlabel("Capital (K)"); ax1.set_ylabel("q")
ax1.legend()
# Time Paths
t_grid = np.linspace(0, 100, 200)
ax2.plot(t_grid, sol_path.sol(t_grid)[0], label='K (Capital)')
ax2_right = ax2.twinx()
ax2_right.plot(t_grid, sol_path.sol(t_grid)[1], 'r', label='q (Tobin\'s q)')
ax2.set_title("Response to TFP Shock")
ax2.set_xlabel("Time")
ax2.legend(loc='upper left'); ax2_right.legend(loc='upper right')
plt.tight_layout()
plt.show()--- PERMANENT TFP SHOCK (A: 1.0 -> 1.1) ---
New Steady State K*: 153.7468
Initial q jump: 1.2500
Total q increase: 4.17%
------------------------------

13.8 Programming Exercises¶
Exercise 13.1 (APL — Saddle Path)¶
Implement the full reverse-shooting algorithm for the system in APL using ⍣ iteration. (a) Find the Jacobian eigenvalues and stable eigenvector at the steady state as a dfn. (b) Perturb from the steady state in both directions and integrate backward using the RK4 dfn from Chapter 3. (c) Trace both branches of the saddle path and export for plotting.
Exercise 13.2 (Python — Adjustment Cost Sensitivity)¶
Compute the investment rate response to a permanent 5% TFP shock for . For each : (a) find the on-impact jump in ; (b) simulate the full convergence path; (c) compute the half-life of the capital stock adjustment. Plot all five convergence paths on a single figure. Verify that higher (larger adjustment costs) leads to slower but more gradual investment.
Exercise 13.3 (Julia — Real Options Trigger)¶
using Roots
function investment_trigger(r, mu, sigma, I_cost)
# Characteristic equation root
beta_plus = (-( mu - sigma^2/2) + sqrt((mu-sigma^2/2)^2 + 2*sigma^2*r)) / sigma^2
# Investment trigger
Pi_star = beta_plus / (beta_plus - 1) * (r - mu) * I_cost
return (beta_plus=beta_plus, Pi_star=Pi_star)
end
r, mu, I_cost = 0.05, 0.02, 1.0
println("Trigger analysis:")
for sigma in [0.10, 0.20, 0.30, 0.40]
res = investment_trigger(r, mu, sigma, I_cost)
println(" σ=$(sigma): β+=$(round(res.beta_plus,digits=3)), Π*=$(round(res.Pi_star,digits=3))")
end
println("\nNPV rule threshold: Π*_NPV = $(r-mu)*I = $((r-mu)*I_cost)")
println("Option value ratio at σ=0.20: $(investment_trigger(r,mu,0.20,I_cost).Pi_star/((r-mu)*I_cost))")Exercise 13.4 (R — Hayashi Theorem Test)¶
Using Compustat data (or simulated data calibrated to U.S. manufacturing), run the Hayashi regression: where is Tobin’s average q. (a) What is the estimated ? What does this imply for ? (b) Add cash flow as an additional regressor. If the Hayashi theorem holds perfectly, cash flow should be insignificant (all relevant information is in ). Is it? (c) Interpret the cash flow coefficient as a measure of financial frictions.
Exercise 13.5 — Irreversibility ()¶
Modify the adjustment cost model to allow for irreversibility: (no disinvestment). Add a KKT multiplier for the constraint with complementary slackness . (a) Show that with irreversibility, can fall below 1 (the firm would like to disinvest but cannot). (b) Characterize the “wait and see” region where and investment is zero. (c) Using the real options trigger formula, explain why irreversibility raises the investment threshold for expansion and creates a symmetric lower threshold for contraction.
Exercise 13.6 — Financial Frictions ()¶
Introduce a financial friction: the firm faces an external finance premium that depends on the gap between and 1. Specifically, the effective discount rate for investment becomes where — external finance becomes more expensive when (financial distress). (a) Modify the costate equation to include the premium. (b) Show that this creates an additional amplification mechanism: a negative TFP shock reduces below 1, raises the external finance premium, further discourages investment, and further reduces — a financial accelerator. (c) Calibrate with and simulate the investment response to a 10% negative TFP shock with and without the financial accelerator.
13.9 Chapter Summary¶
Key results:
The firm’s optimal control problem minimizes investment cost subject to the capital accumulation constraint, with the Hamiltonian .
Pontryagin’s conditions yield the investment function: — investment is positive (negative) when exceeds (falls below) 1, with sensitivity to the gap.
The costate equation combined with the state equation forms a 2D saddle-point system.
The phase plane has a horizontal nullcline at and a curved nullcline; the saddle-point property follows from (as in the RCK model).
Hayashi’s theorem: under CRS production and HD1 adjustment costs, marginal equals average — allowing empirical testing using stock market data.
The real options trigger shows that uncertainty raises the investment threshold above the NPV rule: the ratio is the option value multiplier.
In APL: the system is integrated with the same
rk4_stepdfn as the RCK model; the phase diagram nullclines are computed on a grid usingf¨K_grid(apply function to each element).
Connections forward: The financial accelerator mechanism — investment responds to , which depends on the firm’s balance sheet — is central to Chapter 37’s replication of the Great Recession. Chapter 29 (perturbation methods) shows how the second-order approximation of the value function is needed to capture the precautionary investment response to uncertainty.
Next: Part IV — Discrete-Time Dynamic Models in Macroeconomics