Getting error message from NSolve

Question

I'm trying to solve for coefficients $k1$, $k2$, $k3$, $k4$ in the following equation:

$\det\left(\,\begin{bmatrix}k1+k2 & -k2 & 0\\-k2 & k2+k3 & -k3 \\ 0 & -k3 & k3+k4\end{bmatrix} - \begin{bmatrix}\omega^2 & 0 & 0\\0 & \omega^2 & 0 \\ 0 & 0 & \omega^2\end{bmatrix}\,\right)=0$

where $\omega>0$ is such that $\frac{\omega_1}{2\pi}=0.0028$, $\frac{\omega_2}{2\pi}=0.0036$, $\frac{\omega_3}{2\pi}=0.0042$.

I've tried to solve this by substituting $\omega^2$ with Replace, yielding

w == 
  Root[-a b c - a b d - a c d - b c d + 
   (3909 a b + 7818 a c + 11727 b c + 3909 a d + 7818 b d + 
    3909 c d) #1 + (-15280281 a - 30560562 b - 30560562 c - 
    15280281 d) #1^2 + 59730618429 #1^3 &, 1] || 
w == 
  Root[-a b c - a b d - a c d - b c d + 
    (3909 a b + 7818 a c + 11727 b c + 3909 a d + 7818 b d + 
    3909 c d) #1 + (-15280281 a - 30560562 b - 30560562 c - 
    15280281 d) #1^2 + 59730618429 #1^3 &, 2] || 
w == 
  Root[-a b c - a b d - a c d - b c d + 
   (3909 a b + 7818 a c + 11727 b c + 3909 a d + 7818 b d + 
    3909 c d) #1 + (-15280281 a - 30560562 b - 30560562 c - 
    15280281 d) #1^2 + 59730618429 #1^3 &, 3]

upon which I use NSolve:

equations = 
 {Root[
    -a b c - a b d - a c d - b c d + 
      (3909 a b + 7818 a c + 11727 b c + 3909 a d + 7818 b d + 3909 c d) #1 + 
      (-15280281 a - 30560562 b - 30560562 c - 15280281 d) #1^2 + 
      59730618429 #1^3 &, 
    1] == (0.0028*2*Pi)^2,
  Root[
    -a b c - a b d - a c d - b c d + 
      (3909 a b + 7818 a c + 11727 b c + 3909 a d + 7818 b d + 3909 c d) #1 + 
      (-15280281 a - 30560562 b - 30560562 c - 15280281 d) #1^2 + 
      59730618429 #1^3 &, 
    2] == (0.0036*2*Pi)^2,
  Root[
    -a b c - a b d - a c d - b c d + 
      (3909 a b + 7818 a c + 11727 b c + 3909 a d + 7818 b d + 3909 c d) #1 + 
      (-15280281 a - 30560562 b - 30560562 c - 15280281 d) #1^2 + 
     59730618429 #1^3 &, 
    3] == (0.0042*2*Pi)^2,
  a > 0,
  b > 0,
  c > 0,
  d > 0};

coeffs = NSolve[equations, {a, b, c, d}, Reals]

But I obtain the error message "Infinite solution set has dimension at least 1" and obtain {} as the output.

Does this mean that there are no solutions to the problem, or that I am solving this problem incorrectly? Is there a better way to solve this problem?

Answer 1

For the first $\omega$, the following use of FindInstance finds a {k1, k2, k3, k4} tuple that satisfies the equation more or less instantly.

FindInstance[
 Det[{{k1 + k2, -k2, 0}, {-k2, k2 + k3, -k3}, {0, -k3, 
       k3 + k4}} - {{w^2, 0, 0}, {0, w^2, 0}, {0, 0, w^2}}] == 0 /. 
  w -> (2 π 0.0028),
 {k1, k2, k3, k4}]

The solution space for this value of $\omega$ can be found with Reduce, but it doesn't easily simplify to anything obviously and inherently meaningful:

Reduce[Det[{{k1 + k2, -k2, 0}, {-k2, k2 + k3, -k3}, {0, -k3, 
       k3 + k4}} - {{w^2, 0, 0}, {0, w^2, 0}, {0, 0, w^2}}] == 0 /. 
  w -> (2 π 0.0028)]