Phosphorus nitride (PN) is the simplest molecule formed solely by phosphorus and nitrogen. It represents an interesting model for materials, where phosphorus is directly attached to nitrogen. Nevertheless, both theoretical and experimental studies often provide an incomplete picture on the structural, electronic, and spectral properties of PN. Theoretical predictions often suffer from insufficient level of theory, incomplete basis set, or from neglecting several effects, for example, zero‐point vibrational correction (ZPVC). Therefore, we performed an extensive benchmark study on structural, electronic, and spectral properties of PN at the Hartree‐Fock, density functional theory (DFT), or even the coupled‐cluster levels. We paid special attention to the basis set effect. We tested three variants of Dunning's aug‐cc‐pVXZ basis sets with the size from double‐ζ to sextuple‐ζ, as well as Jensen's aug‐pc‐n, aug‐pcJ‐n, and aug‐pcSseg‐n basis sets, where n = 1‐4. Obtained energetics, PN distance, dipole moment, vibrational frequencies, and nuclear magnetic resonance (NMR) parameters were extrapolated to the complete basis set limit (CBS) using three‐ or two‐parameter formulas. The 31P NMR shieldings estimated with the aug‐cc‐pVXZ and aug‐cc‐pV(X + d)Z basis sets strongly depend on the basis set size providing scattered convergence patterns toward CBS. The Hartree‐Fock self‐consistent field (HF‐SCF) NMR parameters evinced similar behavior as the coupled‐cluster data. The only smooth convergence was achieved using the aug‐cc‐pCVXZ basis sets that include core‐valence effects. The KT3 functional underestimated the phosphorus CBS shieldings by about 12 ppm compared to coupled cluster with singles and doubles (CCSD) (T). Nevertheless, KT3 unambiguously surpasses the HF‐SCF and CCSD levels that provide 31P shieldings that are lower by about 150 ppm and 24 ppm compared to CCSD(T). The convergence of nitrogen shieldings was regular for all basis set hierarchies and all theoretical methods. Relativistic and vibrational effects on selected properties were also discussed.

We report on nucleus‐independent magnetic shielding (NICS) scans over the centers of six‐ and five‐membered rings in selected metal phthalocyanines (MPc) and fullerene C60 for more accurate characterization of local aromaticity in these compounds. Detailed tests were conducted on model aromatic molecules including benzene, pyrrole, indole, isoindole, and carbazole and subsequently applied to H2Pc, ZnPc, Al(OH)Pc, and CuPc. Similar behavior of three selected magnetic probes, Bq, 3He, and 7Li+, approaching perpendicularly the ring centers, was observed. For better visualization of shielding zone over the centers of aromatic rings, we introduced a simple mathematical procedure: the first and second derivatives of scan curves with respect to magnetic probe position enabled their additional examination. It allowed an easier localization of curve minimum and discrimination between areas in space varying by the magnetic field magnitude and to illustrate local aromaticity of two different kinds of rings in MPc with better resolution. Our results supported earlier reports on very low aromaticity indexes of pyrrole ring incorporated into MPc and significant aromaticity of the central macrocycle. No direct dependence between harmonic oscillator model of aromaticity and NICS was observed. Instead, a correlation between position of scan curve minimum and its magnitude were observed. In addition, the NICS values and 3He chemical shifts in the middle of neutral C60 and C606− anion agreed well with the reported experimental NMR values for He@C60 and He@C606−.

Detailed study of Jensen's polarization-consistent vs. Dunning's correlation-consistent basis set families performance on the extrapolation of raw and counterpoise-corrected interaction energies of water dimer using coupled cluster with single, double, and perturbative correction for connected triple excitations (CCSD(T)) in the complete basis set (CBS) limit are reported. Both 3parameter exponential and 2-parameter inverse-power fits vs. the cardinal number of basis set, as well as the number of basis functions were analyzed and compared with one of the most extensive CCSD(T) results reported recently. The obtained results for both Jensen-and Dunning-type basis sets underestimate raw interaction energy by less than 0.136 kcal/mol with respect to the reference value of − 4.98065 kcal/mol. The use of counterpoise correction further improves (closer to the reference value) interaction energy. Asymptotic convergence of 3-parameter fitted interaction energy with respect to both cardinal number of basis set and the number of basis functions are closer to the reference value at the CBS limit than other fitting approaches considered here. Separate fits of Hartree-Fock and correlation interaction energy with 3-parameter formula additionally improved the results, and the smallest CBS deviation from the reference value is about 0.001 kcal/mol (underestimated) for CCSD(T)/augcc-pVXZ calculations. However, Jensen's basis set underestimates such value to 0.012 kcal/mol. No improvement was observed for using the number of basis functions instead of cardinal number for fitting.

Nuclear shieldings and chemical shifts of 5-fluorocytosine (5FC) were predicted in the gas phase and DMSO solution modeled by polarizable continuum model using B3LYP density functional and revised STO(1M)-3G basis set. For comparison, eight arbitrary selected basis sets including STO-3G and medium-size Pople-type and larger dedicated Jensen-type ones were applied.The former basis sets were significantly smaller, but the calculated structural parameters, harmonic vibrational frequencies, were very accurate and close to those obtained with larger, polarization-consistent ones. The predicted 13 C and 1 H chemical shieldings of 5FC and cytosine, selected as parent molecule, were acceptable (root mean square for 13 C chemical shifts in DMSO of about 5 ppm and less) though less accurate than those calculated with large basis sets, dedicated for prediction of nuclear magnetic resonance parameters.