Los circuitos de oscilador de cristal estándar tienen una precisión de +/- 20 ppm a 200 ppm de precisión (ppm = pulsos por millón) según el cristal utilizado. 4 Hz +/- 0.01 Hz es una precisión de 2500 ppm. Entonces ... casi cualquier circuito de oscilador de cristal le dará la precisión que necesita.
El problema que encontrará es que un oscilador de cristal funcionará a una velocidad mucho mayor que 4 Hz. 32.768 KHz es aproximadamente el más bajo que encontrará, con el extremo superior en cientos de MHz (o tal vez un GHz o algo así). Entonces, lo que harías es usar una frecuencia de reloj más alta y dividirla para obtener los 4 Hz que deseas.
Exactamente cómo lo haces todo depende de lo que te resulte cómodo diseñar y construir. Puede usar un oscilador de tipo lata de metal y luego dividirlo con fichas de contador simples. O podría usar un CPLD o FPGA (exceso). O incluso un pequeño microcontrolador para el oscilador de cristal y el divisor.
Yo usaría el enfoque del microcontrolador. Si nunca ha hecho un enfoque de microcontrolador, probablemente sea más difícil hacerlo. Pero es bastante flexible, de baja potencia y tamaño pequeño. Aunque seré el primero en admitir que este enfoque no es para todos.