No puedo demostrar la identidad, pero encontré que implica una identidad bastante extraña de suma de caracteres de Dirichlet, así que estoy registrando la implicación aquí por si alguien reconoce esto último.
Me especializaré en el caso donde $2n-1$ es coprimo con $3$.
Prefiero escribir $\epsilon = -\omega^2$, donde $\omega := e^{2\pi i/3} = z^{2n-1}$ es una raíz cúbica de la unidad. Reúne todos los términos en un lado: $$ \sum_{k=1}^{n-1} \frac{z^k + \omega^2 z^{2k} - \omega^2 z^{4k} - \omega z^{5k}}{1-z^{6k}} = 0.$$ A partir de la fácilmente verificable identidad de "serie de Fourier aserrada discreta" $$ \frac{1}{1-z^{6k}} = -\frac{1}{2n-1} \sum_{j=0}^{2n-2} j z^{6jk}$$ para cualquier $k$ no divisible por $2n-1$ (probado multiplicando el denominador, cancelando términos y aplicando la fórmula de la serie geométrica), escribimos la identidad anterior para demostrar como $$ \sum_{j=0}^{2n-2} j \sum_{k=1}^{n-1} (z^{(6j+1)k} + \omega^2 z^{(6j+2)k} - \omega^2 z^{(6j+4)k} - \omega z^{(6j+5)k}) = 0.$$ [Observación: en este punto se podría usar la fórmula de la serie geométrica para sumar en $k$, pero esto nos devuelve a una identidad de complejidad similar a la que estamos tratando de probar originalmente, y no pude avanzar mucho siguiendo ese camino.]
Suponiendo que $2n-1$ es coprimo con $3$, podemos escribir $z = \omega^{2n-1} \zeta$ para alguna raíz primitiva $(2n-1)^{\mathrm{ésima}}$ de la unidad. Luego reducimos a $$ \sum_{j=0}^{2n-2} j\sum_{k=1}^{n-1} (\omega^{(2n-1)k} \zeta^{(6j+1)k} + \omega^{2(2n-1)k+2} \zeta^{(6j+2)k} - \omega^{(2n-1)k+2} \zeta^{(6j+4)k} - \omega^{2(2n-1)k+1} \zeta^{(6j+5)k}) = 0.$$ Desde la teoría de Galois vemos que el coeficiente neto de $\zeta^a$ tendría que ser independiente de $a$ para cada clase de residuo primitiva $a$ módulo $2n-1$. En particular, si $\chi : {\bf Z}/(2n-1){\bf Z} \to S^1$ es un carácter de Dirichlet no principal de período $2n-1$, entonces una identidad necesaria a establecer es $$ \sum_{j=0}^{2n-2} j \sum_{k=1}^{n-1} (\omega^{(2n-1)k} \chi((6j+1)k) + \omega^{2(2n-1)k+2} \chi((6j+2)k) - \omega^{(2n-1)k+2} \chi((6j+4)k) - \omega^{2(2n-1)k+1} \chi((6j+5)k) = 0.$$ El punto de hacer esto es que el lado izquierdo se factoriza como $$ a_1 b_1 + a_2 b_2 - a_3 b_3 - a_4 b_4$$ donde los $a_i, b_i$ son sumas de caracteres $$ a_1 := \sum_{j=0}^{2n-2} j \chi(6j+1)$$ $$ a_2 := \sum_{j=0}^{2n-2} j \chi(6j+2)$$ $$ a_3 := \sum_{j=0}^{2n-2} j \chi(6j+4)$$ $$ a_4 := \sum_{j=0}^{2n-2} j \chi(6j+5)$$ y $$ b_1 := \sum_{k=1}^{n-1} \omega^{(2n-1)k} \chi(k)$$ $$ b_2 := \sum_{k=1}^{n-1} \omega^{2(2n-1)k+2} \chi(k)$$ $$ b_3 := \sum_{k=1}^{n-1} \omega^{(2n-1)k+2} \chi(k)$$ $$ b_4 := \sum_{k=1}^{n-1} \omega^{2(2n-1)k+1} \chi(k)$$ En este punto ocurre un pequeño milagro algebraico (que posiblemente sugiere que se va por el camino correcto): al usar la sustitución $j \mapsto 2n-2-j$ (así como la naturaleza de media cero de $\chi$), se puede verificar que $$ a_4 = -\chi(-1) a_1; \quad a_3 = -\chi(-1) a_2$$ y de manera similar $$ \chi(-1) b_3 = \sum_{k=1-n}^{-1} \omega^{2(2n-1)k+2} \chi(k)$$ $$ \chi(-1) b_4 = \sum_{k=1-n}^{-1} \omega^{(2n-1)k+1} \chi(k) = \sum_{k=n}^{2n-2} \omega^{(2n-1)k} \chi(k)$$ y así llegamos a la curiosa identidad de suma de caracteres $$(\sum_{j=0}^{2n-2} j \chi(6j+1)) (\sum_{k=1}^{2n-2} \omega^{(2n-1)k} \chi(k)) = -(\sum_{j=0}^{2n-2} j \chi(6j+2)) (\sum_{k=1-n}^{n-1} \omega^{2(2n-1)k+2} \chi(k)),$$ que es válida siempre que $2n-1$ sea coprimo con $3$ y $\chi: {\bf Z}/(2n-1){\bf Z} \to S^1$ sea un carácter de Dirichlet no principal. Puedo verificar esta identidad numéricamente, pero no sé por qué es cierta en general; no reconozco ninguna identidad para estas sumas de caracteres, aunque es ligeramente reminiscente de algunas sumas de caracteres que surgen en algunas demostraciones de la reciprocidad cuadrática. [Observación: se podría intentar aplicar la inversión de Fourier / suma de Poisson a las sumas de caracteres anteriores (en el espíritu de la demostración de la desigualdad de Polya-Vinogradov), pero luego comienzan a surgir denominadores similares a los denominadores iniciales $\frac{1}{1-z^{6k}}$ al hacerlo, por lo que se vuelve a una identidad de complejidad comparable a la identidad inicial, aunque al menos sin el extraño rango de suma $\sum_{k=1}^{n-1}$ que no estaba tan naturalmente relacionado con el período $3(2n-1)$ del problema. Esto no necesariamente tiene que ser un movimiento malo, por ejemplo, podría haber alguna simetría o cancelación milagrosa que explotar en la formulación de "denominador" del problema, pero no pude sacar provecho de tal transformación.]