[Edited to include the connection with Euler's theorem that
there's no nonconstant $4$-term arithmetic progression of squares; briefly,
since $1$ and $4(N^2+N)+1 = (2N+1)^2$ are always squares, $(1,b,c,(2N+1)^2)$
is such a $4$-term progression, so $b=c=1$ and $N^2+N=0$ for any rational
solution. It seems that the same deals with all $k \geq 3$.]
The Diophantine equations
$$ B^2 = b = \frac43(N^2+N) + 1 $$
$$ C^2 = c = \frac83(N^2+N) + 1 $$
have no simultaneous solution in integers other than the eight trivial
solutions with $N=0,-1$, $B = \pm 1$, $C = \pm 1$. Such a result
can be tricky to prove in general, but here we are somewhat lucky
in that these are the only rational solutions, and moreover
$N=0$ and $N=-1$ are the only rational numbers for which the product $bc$
is a square. This can be proved using Fermat's "descent" method;
the calculations, though elementary, may be somewhat laborious
to carry out by hand, but happily this is no longer necessary
thanks to existing tables and software.
In general, if $P,Q$ are quadratic polynomials such that $PQ$ has
four distinct factors then the simultaneous Diophantine equations
$B^2 = P(N)$, $C^2 = Q(N)$ define an elliptic curve. A fundamental
1929
theorem of Siegel assures that such a curve has only finitely many
integral points. The proof is "ineffective" and does not in general
give an algorithm guaranteed to determine all solutions. Later
theoretical and computational advances do provide such an algorithm,
which is feasible at least for $P,Q$ with small coefficients;
but the resulting proof is very far from elementary, and it can be
hard to predict in any given case how hard it is to give an elementary proof.
In the present case, though,
the elliptic curve already has finitely many rational points,
as does the "isogenous" curve $A^2 = P(N) Q(N)$. We bring this curve
$$
A^2 = \left( \frac43(N^2+N) + 1 \right) \left( \frac83(N^2+N) + 1 \right)
$$
into standard Weierstrass form in the usual way starting from the
rational point $(N,A) = (0,1)$: the Taylor expansion of $A$ about $N=0$
starts $1 + 2N + O(N^2)$, so for $N\neq 0$ we can write
$$
A = 1 + 2N + rN^2
$$
for some rational $r$, and divide the resulting equation by $N^2$ to get
$$
(9r^2-32)N^2 + (36r-64)N + (18r-32) = 0.
$$
This equation is quadratic in $N$, and thus
has rational roots iff its discriminant w.r.t. $N$ is a square.
That discriminant factors as $-72(r-2)r(9r-16)$; taking $r=-2x/9$
and removing a factor $(8/3)^2$ we find that
$$
y^2 = x (x+8) (x+9) = x^3 + 17 x^2 + 72 x
$$
for some rational $x,y$. This elliptic curve turns out to have
conductor $24$, so it already appears in Tingley's "Antwerp Tables"
of modular elliptic curves of conductor at most $200$;
it turns out to have label 24C
here.
We find that it has rank zero, and only four rational points, which
must be the point at infinity and the three "$2$-torsion points" with $y=0$.
Alternatively, we can input [0,17,0,72,0] to Cremona's program mwrank
to find that the curve has rank zero, and then find its torsion points
with gp. Either way, we finish by retracing our steps
to find that each of $r=0,2,16/9$
corresponds to one of the known solutions with $N=0$ or $N=-1$
(in each case a double root of the quadratic in $N$), so we are done.
[Thanks to Will Jagy for
calling my attention to this question.]
Added later: It turns out that this $2$-descent calculation
long predates the Antwerp tables: it is essentially equivalent to
Euler's proof of his theorem that there is no nonconstant
$4$-term arithmetic progression of squares. (See for instance
Keith
Conrad's exposition, which states that Euler proved the result in 1780,
answering a question "first raised by Fermat in 1640".)
Indeed if $b$ and $c$ are squares then
$$
1, \
\frac43(N^2+N)+1, \
\frac83(N^2+N)+1, \
4(N^2+N)+1
$$
is such a progression (the last term being $(2N+1)^2$,
so it is constant by Euler, whence $b=c=1$ and we are done.
jamaicanworm writes in a comment that for general $k \geq 3$
the problem is whether there exists a positive integer $N$
such that $c_i := \frac{4(k-i)}{k}(N^2+N) + 1$ is a square
for each $i=1,2,\ldots,k-1$. These $c_i$ form an arithmetic progression,
and extending it by one term in each direction again yields the squares
$c_k = 1$ and $c_0 = (2N+1)^2$. Hence we have an arithmetic progression
of $k+1$ squares, and again Euler's theorem shows that even if we allow
$N$ to be rational the only examples are the trivial ones with $N^2+N=0$.
Consider the system of $n$ congruences
$x\equiv 0\pmod{(2)(3)}$
$x\equiv -1\pmod{(5)(7)}$
$x\equiv -2 \pmod{(11)(13)}$
and so on, up to
$x\equiv -(n-1)\pmod{(p_{2n-3})(p_{2n-2})}$
where $p_i$ is the $i$-th prime.
By the Chinese Remainder Theorem, this system of congruences has infinitely many positive solutions $x$.
If $x$ is such a solution, then $x$ is divisible by $(2)(3)$, $x+1$ is divisible by $(5)(7)$, and so on.
It follows that none of $x$, $x+1$, and so on up to $x+(n-1)$ is a prime power.
Best Answer
$a,a+1,a+2,a+3,a+4,a+5$, and $a+6$ comprise a complete set of residues modulo $7$,
so we can take $a=0$ without loss of generality.
For $n=1$, the sum is congruent to $0+1+2+3+4+5+6\equiv0\bmod7$.
For $n=2$, the sum is congruent to $0+1+4+2+2+4+1\equiv0\bmod7$.
For $n=3$, the sum is congruent to $0+1+1-1+1-1-1\equiv0\bmod7$.
For $n=4$, the sum is congruent to $0+1+2+4+4+2+1\equiv0\bmod7$.
For $n=5$, the sum is congruent to $0+1+4+5+2+3+6\equiv0\bmod7$.
For $n=6$, the sum is congruent to $0+1+1+1+1+1+1\equiv6\bmod7$.
For $n=6k+m$, the sum is congruent to that for $m$, as a consequence of Fermat's little theorem.
Therefore, the sum is a multiple of $7$ unless $n$ is a multiple of $6$.