  Real-time dynamics techniques have proven increasingly useful in
understanding strongly correlated systems both theoretically and
experimentally. By employing unbiased time-resolved exact diagonalization, we
study pump dynamics in the two-dimensional plaquette Hubbard model, where
distinct hopping integrals $t_h$ and $t_h^\prime$ are present within and
between plaquettes. In the intermediate coupling regime, a significant
enhancement of $d$-wave superconductivity is observed and compared with that
obtained by simple examination of expectation values with the eigenstates of
the Hamiltonian. Our work provides further understanding of superconductivity
in the Hubbard model, extends the description of the pairing amplitude to the
frequency-anisotropy plane, and offers a promising approach for experimentally
engineering emergent out-of-equilibrium states.
