Your program has terrible numerical scaling. Change units so that all non-zero input data is within a small number of orders of magnitude if 1.
Also, If you have access to Mosek, use it as solver. Otherwise, follow the instructions at CVXQUAD: How to use CVXQUAD’s Pade Approximant instead of CVX’s unreliable Successive Approximation for GP mode, log, exp, entr, rel_entr, kl_div, log_det, det_rootn, exponential cone. CVXQUAD’s Quantum (Matrix) Entropy & Matrix Log related functions