Your problem is likely that some observations are exactly at the truncation boundary (i.e. zero).
so your truncated normal model with lb = 0 is not actually a proper solution I think.
I would indeed also suggest not to use the absolute values but the original ones with the signs still existing. If you insist on using the absolute values, you may try out the zero_inflated_lognormal family to handle the zeros in some way.