Nuclear couplings for the Se-Se bonds, 1J(Se, Se), are analyzed on the basis of the molecular orbital (MO) theory. The values are calculated by employing the triple ζ basis sets of the Slater type at the DFT level. 1J(Se, Se) are calculated modeled by MeSeSeMe (1a), which shows the typical torsional angular dependence on ϕ(CMeSeSeCMe). The dependence explains well the observed 1Jobsd (Se, Se) of small values (≤ 64 Hz) for RSeSeR′ (1) (simple derivatives of 1a) and large values (330–380 Hz) observed for 4-substituted naphto[1,8-c, d]-1,2-diselenoles (2) which correspond to symperiplanar diselenides. 1J (Se, Se: 2) becomes larger as the electron density on Se increases. The paramagnetic spin-orbit terms contribute predominantly. The contributions are evaluated separately from each MO (ψi) and each ψi → ψa transition, where ψi and ψa are occupied and unoccupied MO's, respectively. The separate evaluation enables us to recognize and visualize the origin and the mechanism of the couplings.