The PyProcar Python package plots the band structure and the Fermi surface as a function of site and/or s,p,d,f -projected wavefunctions obtained for each k-point in the Brillouin zone and band in an electronic structure calculation. This can be performed on top of any electronic structure code, as long as the band and projection information is written in the PROCAR format, as done by the VASP and ABINIT codes. PyProcar can be easily modified to read other formats as well. This package is particularly suitable for understanding atomic effects into the band structure, Fermi surface, spin texture, etc. PyProcar can be conveniently used in a command line mode, where each one of the parameters define a plot property. In the case of Fermi-surfaces, the package is able to plot the surface with colors depending on other properties such as the electron velocity or spin projection. The mesh used to calculate the property does not need to be the same as the one used to obtain the Fermi surface. A file with a specific property evaluated for each k-point in a k−mesh and for each band can be used to project other properties such as electron-phonon mean path, Fermi velocity, electron effective mass, etc. Another existing feature refers to the band unfolding of supercell calculations into predefined unit cells.

This article reviews different formulations of the thermodynamic cycles used for the prediction of pK a values, their advantages, and disadvantages with special emphasis on the limitations resulting from the necessity of gas-phase calculations, which allow introducing some difficult cases that motivated alternative strategies. Before introducing the protocols that do not consider gas-phase calculations, the two current opinions available in the literature on the debate about the correct formalism for the calculation of free energies in solution are briefly introduced. Then, the isodesmic proton exchange reaction in solution is reviewed by analyzing its performance on difficult cases for thermodynamic cycles such as carbon acids and amino acids. The pK a values predicted by the isodesmic reaction for common acid species are also reviewed to compare their accuracy results in relation with those of thermodynamic cycles. Linear regressions between experimental pK a values and the calculated free energies obtained with the isodesmic reaction provide expressions for the dependence of the error in the calculated pK a s on the pK a difference between the studied acid and the reference species. Finally, it is shown that linear regressions correct the calculated free energies of the isodesmic reaction, when high constant precision is required in a large pK a range. The deviations from the expected behavior are equivalent to those reported previously for different pK a calculation protocols and are determined by the inaccuracies of continuum solvent models on the interactions with ionic species.

A recent study associate carbon with single photon emitters (SPEs) in hexagonal boron nitride (h-BN). This observation, together with the high mobility of carbon in h-BN suggest the existence of SPEs based on carbon clusters. Here, by means of density-functional theory calculations we studied clusters of substitutional carbon atoms up to tetramers in hexagonal boron nitride. Two different conformations of neutral carbon trimers have zero-point line energies and shifts of the phonon sideband compatible with typical photoluminescence spectra. Moreover, some conformations of two small C clusters next to each other result in photoluminescence spectra similar to those found in experiments. We also showed that vacancies are unable to reproduce the typical features of the phonon sideband observed in most measurements due to the large spectral weight of low-energy breathing modes, ubiquitous in such defects.