  We present a functional renormalization group analysis of superconductivity
in the ground state of the attractive Hubbard model on a square lattice.
Spontaneous symmetry breaking is treated in a purely fermionic setting via
anomalous propagators and anomalous effective interactions. In addition to the
anomalous interactions arising already in the reduced BCS model, effective
interactions with three incoming legs and one outgoing leg (and vice versa)
occur. We accomplish their integration into the usual diagrammatic formalism by
introducing a Nambu matrix for the effective interactions. From a random-phase
approximation generalized through use of this matrix we conclude that the
impact of the 3+1 effective interactions is limited, especially considering the
effective interactions important for the determination of the order parameter.
The exact hierarchy of flow equations for one-particle irreducible vertex
functions is truncated on the two-particle level, with higher-order self-energy
corrections included in a scheme proposed by Katanin. Using a parametrization
of effective interactions by patches in momentum space, the flow equations can
be integrated numerically to the lowest scales without encountering
divergences. Momentum-shell as well as interaction-flow cutoff functions are
used, including a small external field or a large external field and a
counterterm, respectively. Both approaches produce momentum-resolved order
parameter values directly from the microscopic model. The size of the
superconducting gap is in reasonable agreement with expectations from other
studies.
