How to solve dislotwin model errors?

[tutfighting]

I am new to the dislotwin model and currently conducting relevant exploratory calculations. While configuring material properties via the dislotwin_twip-steel.yaml input file, the computation terminates with Error Code 211: material parameter out of bounds: rho.
I have referred to available literature and benchmark datasets to calibrate multiple relevant parameters, yet this out-of-bounds error persists consistently.
Could you kindly assist in troubleshooting the root cause of this issue? My configured material parameters are attached for your reference.

homogenization:
SX:
N_constituents: 1
mechanical: {type: pass}

phase:
Aluminum:
lattice: cF
thermal:
C_p: 490.0
K_11: 54.0
mechanical:
output: [F, P, F_e, F_p, L_p, gamma_sl]
elastic: {C_11: 219000000000.0, C_12: 126000000000.0, C_44: 130000000000.0, type: Hooke}
plastic:
type: dislotwin
output: [rho_mob, rho_dip, gamma_sl, Lambda_sl, tau_pass, gamma_tw, Lambda_tw, f_tr]

Glide

    N_sl: [12]

    b_sl: [2.48e-10] # a/sqrt(2)
    Q_sl: [3.5e-19]
    p_sl: [0.325]
    q_sl: [1.55]
    B: [0.001]
    i_sl: [7.5]
    v_0: [1.4e+3]
    tau_0: [5.5e+8] # adjusted
    D_a: 2.0
    Q_cl: 3.0e-19

    rho_mob_0: [2.8e+12]
    rho_dip_0: [1.0]

    h_sl-sl: [0.122, 0.122, 0.625, 0.07, 0.137, 0.137, 0.122]

Twin

    N_tw: [12]
    b_tw: [1.47e-10] # a_cF/sqrt(6)
    L_tw: 1.91e-7 # 1300 *b_tw
    i_tw: 10.0
    t_tw: [5.0e-8]
    p_tw: [7] # A, adjusted

    h_tw-tw: [1.0, 1.0]
    h_sl-tw: [1.0, 1.0, 1.0]

Twin & Transformation

    T_ref: 293.15
    Gamma_sf: 2.833e-2
    Gamma_sf,T: 1.214e-4
    Gamma_sf,T^2: 1.473e-7

    x_c: 1.0e-9
    V_cs: 1.67e-29

Slip & Twin & Transformation

    D: 1.08e-5