def d2_theta(t, alpha):
    """
    theta''(t) = alpha / (alpha + |t|)^2
    Also called phi'' or psi''.
    Nikolova et al 2013, table 1, f3.
    """
    assert alpha > 0
    denom = alpha + abs(t)
    return alpha / (denom*denom)