In this paper, we focus on solving the vector scheduling problem with submodular penalties on parallel machines. We are given
n
jobs and
m
parallel machines, where each job is associated with a
d
-dimensional vector. Each job can either be rejected, incurring a rejection penalty, or accepted and processed on one of the
m
parallel machines. The objective is to minimize the sum of the maximum load overall dimensions of the total vector for all accepted jobs, along with the total penalty for rejected jobs. The penalty is determined by a submodular function. Our main work is to design a
2
−
1
/
m
min
r
,
d
-approximation algorithm to solve this problem. Here,
r
denotes the maximum ratio of the maximum load to the minimum load on the
d
-dimensional vectors among all jobs.