The difficulty is that wB97M(2) defines the molecular orbitals and Fock matrix with one density functional and SCF energy with another. Reproducing it by exposing all the parameters is a little painful, but possible. It requires running two chained calculations and a bit of post-processing by hand. Here is the input:

```
$molecule
0 1
He
He 1 1.0
$end
$rem
exchange = gen
correlation = rimp2
basis = vdz
aux_basis = rimp2-vdz
lrc_dft = 1
hfk_sr_coef = 15000000
hfk_lr_coef = 100000000
omega = 300
dh = 1
scs = 3
sss_factor = 340960
sos_factor = 340960
nl_correlation = vv10
nl_vv_b = 600
nl_vv_c = 100
nl_vv_scale = 100000
! Uncomment these settings to use rVV10 instead of VV10
!use_rvv10 = 1
!nl_vv_b = 620
!nl_vv_c = 100
!nl_vv_scale = 100000
$end
$xc_functional
k 0.15
x wb97mv 1.0
c wb97mv 1.0
$end
@@@
$molecule
READ
$end
$rem
exchange = gen
basis = vdz
lrc_dft = 1
hfk_sr_coef = 62194000
hfk_lr_coef = 100000000
omega = 300
nl_correlation = vv10
nl_vv_b = 1000
nl_vv_c = 100
nl_vv_scale = 65904
scf_guess = read
scf_max_cycles = 0
$end
$xc_functional
k 0.62194
x wb97m2 1.0
c wb97m2 1.0
$end
```

The first job will produce

```
SCF energy in the final basis set = -5.6280140740
Total energy in the final basis set = -5.6280140740
. . .
Doing a Spin Component Scaled (SCS)-MP2 calculation
Energies are scaled by
Same Spin Scaling factor = 0.340960
Opposite Spin Scaling factor = 0.340960
Components of the RIMP2 correlation energy:
aaaa correlation energy = -0.0000080287 a.u.
abab correlation energy = -0.0088227548 a.u.
bbbb correlation energy = -0.0000080287 a.u.
non-Brillouin singles = -0.0000000000 a.u.
total same-spin energy = -0.0000160573 a.u.
total opposite-spin energy = -0.0088227548 a.u.
Total RIMP2 correlation energy = -0.0088388121 a.u.
RIMP2 total energy = -5.6368528860 a.u.
```

We will only be interested in the total RIMP2 correlation energy: `-0.0088388121`

The second job will produce

```
---------------------------------------
Cycle Energy DIIS Error
---------------------------------------
1 -5.6180455678 2.86E-03 Convergence criterion met
---------------------------------------
SCF time: CPU 0.14 s wall 0.14 s
SCF energy in the final basis set = -5.61804557
Total energy in the final basis set = -5.61804557
```

Here we are interested in the final SCF energy: `-5.6180455678`

.

Total wB97M(2) energy is then `-5.6180455678 + (-0.0088388121) = -5.6268843799`

Compare this with the output produced by a streamlined job

```
$molecule
0 1
He
He 1 1.0
$end
$rem
method = wb97m(2)
basis = vdz
aux_basis_corr = rimp2-vdz
$end
```

```
Nonlocal correlation = 0.0056725632
7 -5.6280140742 5.50E-07 Convergence criterion met
---------------------------------------
wB97M-V orbitals are now ready for the wB97M(2) calculation.
wB97M(2) energy = Exch (0.62194 SRHF + 1.0 LRHF + wB97M(2)_EXCH) +
Corr (wB97M(2)_CORR + 0.34096 MP2 + 0.65904 VV10(b=10))
SCF time: CPU 1.08 s wall 1.09 s
SCF energy in the final basis set = -5.61804601
Total energy in the final basis set = -5.61804601
. . .
Doing a Spin Component Scaled (SCS)-MP2 calculation
Energies are scaled by
Same Spin Scaling factor = 0.340960
Opposite Spin Scaling factor = 0.340960
Components of the RIMP2 correlation energy:
aaaa correlation energy = -0.0000080287 a.u.
abab correlation energy = -0.0088227574 a.u.
bbbb correlation energy = -0.0000080287 a.u.
non-Brillouin singles = -0.0000000000 a.u.
total same-spin energy = -0.0000160573 a.u.
total opposite-spin energy = -0.0088227574 a.u.
Total RIMP2 correlation energy = -0.0088388147 a.u.
RIMP2 total energy = -5.6268848268 a.u.
```

The final result is the same: `-5.6268848268`