High level Valence Bond calculations support a predominantly carbenic electronic structure for isocyanides, with a secondary zwitterionic character, despite their linear geometry. This geometry results from the significant energetic stabilization due to nitrogen π lone pair donation. Results are not changed by substitution or solvation effects.

We consider the use in quantum Monte Carlo calculations of two types of valence bond wave functions based on strictly localized active orbitals, namely valence bond self-consistent-field and breathing-orbital valence bond wave functions. Complemented by a Jastrow factor, these Jastrow-valence-bond wave functions are tested by computing the equilibrium well depths of the four diatomic molecules C₂, N₂, O₂, and F₂ in both variational Monte Carlo and diffusion Monte Carlo. We show that it is possible to design compact wave functions based on chemical grounds that are capable of describing both static and dynamic electron correlations. These wave functions can be systematically improved by inclusion of valence bond structures corresponding to additional bonding patterns.