if N < N_lim3:
                    # Ogston-Zimm
                    L = contour_length(Nbp, b)  # (m)
                    Rg = radius_gyration(L, lp)  # (m)
                    D0 = free_solution(kB, temperature, eta, Rg)  # (m^2/s)