I'm currently working through the book Heisenberg's Quantum Mechanics (Razavy, 2010), and am reading the chapter on classical mechanics. I'm interested in part of their derivative of a generalized Lorentz force via a velocity-dependent potential.
I understand the generalized force that they derive from a Lagrangian of the form $L = \frac{1}{2}m|\vec v|^2 - V(\vec r,\vec v,t)$
$$F_i = -\frac{\partial V}{\partial x_i} + \frac{d}{dt}\left(\frac{\partial V}{\partial v_i}\right)$$
Through a series of steps that I still do not quite understand, the author derives the identity for the mixed velocity derivatives of the force:
$$\frac{\partial F_i}{\partial v_j\partial v_k} = 0$$
At this point, "by integrating this equation once" with respect to $v_k$ , they obtain the equation:
$$\frac{\partial F_i}{\partial v_j} = \sum_k \varepsilon_{ijk}B_k(\vec r,t)$$
where $B_k$ is the $k^{th}$ component of a vector function $\vec B$ that does not depend on velocity.
I'm having trouble understanding where this expression for the integral comes into play. The left-hand side clearly comes from the FTC. Were I to perform the integration myself I would do the same and include an arbitrary function
$$\frac{\partial F_i}{\partial v_j}=g(\vec r, v_1,...,v_{k-1}, v_{k+1},..., t)$$
where $g$ is a function that does not depend on $v_k$ explicitly. In this way $\frac{\partial g}{\partial v_k} =0 $ as we need.
I've tried to work out how this function is related to the expression with $B_k$, but I cannot find any source that could point me in the right direction, especially because my best guess for $g$ depends on the other $n-1$ components of the velocity while the author's $\vec B$ vector is a function of position and time only.
Could I have some help understanding what's being done here?
Edit: Additional important context
Additionally, Razavy goes a step further and assumes that the generalized force is independent of acceleration, just like the Lagrangian. Using this assumption, we can take the second condition listed in another related question I asked to form the anti-symmetry relation
$$\frac{\partial F_i}{\partial v_j} =- \frac{\partial F_j}{\partial v_i}$$
And then we can start taking partial derivatives, assuming all these derivatives are continuous. Taking the left side first:
$$ \frac{\partial}{\partial v_k}(LHS)=\frac{\partial^2 F_i}{\partial v_j\partial v_k} = \frac{\partial^2 F_i}{\partial v_k\partial v_j} = \frac{\partial}{\partial v_j}\frac{\partial F_i}{\partial v_k}= \frac{\partial}{\partial v_j}\left(-\frac{\partial F_k}{\partial v_i}\right) = -\frac{\partial^2 F_k}{\partial v_i\partial v_j} $$
So, we can differentiate and swap the top index and a bottom index at the cost of a negative sign. In a similar way, the right hand side can be differentiated
$$\frac{\partial}{\partial v_k}(RHS)=-\frac{\partial^2 F_j}{\partial v_i\partial v_k}=\frac{\partial^2 F_k}{\partial v_i\partial v_j}$$
Thus, We can write: $\frac{\partial}{\partial v_k}(LHS) = -\frac{\partial}{\partial v_k}(RHS)$.
Because $LHS=-RHS$, we have
$$\frac{\partial}{\partial v_k}(LHS) = \frac{\partial^2 F_i}{\partial v_j\partial v_k} = 0$$