Based on the achievements of the previous three papers of this series, the best working formulas for scattering matrix are obtained for both the Landau-Zener (LZ) and the nonadiabatic tunneling (NT) case: two formulas valid at b 2 ≥0 and b 2 ≤0 in the LZ case, and three formulas valid at b 2 ≤-1, -1≤b 2 ≤1 and b 2 ≥1 in the NT case, where b 2 represents the effective energy. Simple and compact formulas which work far better than the LZ formula are proposed for nonadiabatic transition probability by one passage of crossing point for both the LZ and NT cases. Furthermore, compact expressions are derived, for the first time, for the nonadiabatic tunneling probability at b 2 ≤1, i.e., at energies lower than the bottom of the upper adiabatic potential. All the formulas proposed here can be usefully utilized at any coupling strength, namely the validity range has been very much expanded compared to the previous formulas by employing certain empirical corrections. Besides, these formulas have convenience to enable an extension to general curved potentials.