R crashes when running random effects on longitudinal data

[maggiemmurphy]

I am trying to run ANCOM-BC with the following random effects formula on a longitudinal data set:
foal_ancom_output_refchange = ancombc2(data = foals, assay_name = "counts", tax_level = "Family",
fix_formula = "timepoint_ref", rand_formula = "(1|Animal.ID)",
p_adj_method = "holm", prv_cut = 0.10, lib_cut = 1000, s0_perc = 0.05,
group = "timepoint_ref", struc_zero = TRUE, neg_lb = TRUE,
alpha = 0.05, n_cl = 15, verbose = TRUE,
global = TRUE, pairwise = FALSE,
dunnet = TRUE, trend = FALSE,
iter_control = list(tol = 1e-5, max_iter = 200,verbose = T),
em_control = list(tol = 1e-5, max_iter = 200),
lme_control = lme4::lmerControl(optimizer = "nlminbwrap"),
mdfdr_control = list(fwer_ctrl_method = "holm", B = 100),
trend_control = FALSE)

It usually crashes or takes forever and then eventually comes out with convergence errors or crashes the RStudio interface I use in our HPRC. I believe I have the most up to date packages (ANCOM 2.4.0) and am using R4.3.1.
I have tried running the data without the random effects and it finishes just fine. However, including the random effects, even with a smaller subset of samples, crashing continues and errors keep popping up.