Through the interaction of Higgs with heavy quarks in standard model, we have systematically studied and predicted the production of excited doubly heavy baryons based on non-relativistic QCD theory. The decay widths, differential distributions, and major theoretical uncertainties of the excited doubly heavy baryons via the process $$H \rightarrow \langle QQ'\rangle [n] \rightarrow \Xi _{QQ'}+ \bar{Q'} + \bar{Q}$$
H
→
⟨
Q
Q
′
⟩
[
n
]
→
Ξ
Q
Q
′
+
Q
′
¯
+
Q
¯
are discussed in detail. The spin and color quantum number of the intermediate P-wave diquark state $$\langle QQ'\rangle [n]$$
⟨
Q
Q
′
⟩
[
n
]
can be $$\langle cc\rangle [^{1}P_{1}]_{\varvec{\bar{3}}}$$
⟨
c
c
⟩
[
1
P
1
]
3
¯
, $$\langle cc\rangle [^{3}P_{J}]_{\varvec{6}}$$
⟨
c
c
⟩
[
3
P
J
]
6
, $$\langle bc\rangle [^{1}P_{1}]_{\varvec{\bar{3}}/ \textbf{6}}$$
⟨
b
c
⟩
[
1
P
1
]
3
¯
/
6
, $$\langle bc\rangle [^{3}P_{J}]_{\varvec{\bar{3}}/ \textbf{6}}$$
⟨
b
c
⟩
[
3
P
J
]
3
¯
/
6
, $$\langle bb \rangle [^{1}P_{1}]_{\varvec{\bar{3}}}$$
⟨
b
b
⟩
[
1
P
1
]
3
¯
and $$\langle bb\rangle [^{3}P_{J}]_{\textbf{6}}$$
⟨
b
b
⟩
[
3
P
J
]
6
, with $$J=0, 1, 2$$
J
=
0
,
1
,
2
. The contributions from all summed P-wave states can be about $$3.05\%$$
3.05
%
, $$3.23\%$$
3.23
%
and $$2.19\%$$
2.19
%
of the S-wave states for the production of $$\Xi _{cc}$$
Ξ
cc
, $$\Xi _{bc}$$
Ξ
bc
and $$\Xi _{bb}$$
Ξ
bb
, accordingly. Therefore, there will be about 0.41$$\times 10^4$$
×
10
4
events of $$\Xi _{cc}$$
Ξ
cc
, 6.35$$\times 10^4$$
×
10
4
events of $$\Xi _{bc}$$
Ξ
bc
and 0.28$$\times 10^4$$
×
10
4
events of $$\Xi _{bb}$$
Ξ
bb
produced per year at the HL-LHC, and a smaller number of events would be produced at the CEPC or ILC but with a cleaner background to be measured by the experiments.