QUANTUM ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). The acronym ESPRESSO stands for opEn Source Package for Research in Electronic Structure, Simulation, and Optimization. It is freely available to researchers around the world under the terms of the GNU General Public License. QUANTUM ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively parallel architectures, and a great effort being devoted to user friendliness. QUANTUM ESPRESSO is evolving towards a distribution of independent and interoperable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes.

No abstract

We compare, through first-principles pseudopotential calculations, the structural, electronic and optical properties of different size silicon nanoclusters embedded in a SiO2 crystalline or amorphous matrix, with that of free-standing, hydrogenated and hydroxided silicon nanoclusters of corresponding size and shape. We find that the largest effect on the opto-electronic behavior is due to the amorphization of the embedded nanocluster. In that, the amorphization reduces the fundamental gap while increasing the absorption strength in the visible range. Increasing the nanocluster size does not change substantially this picture but only leads to the reduction of the absorption threshold, following the quantum confinement rule. Finally, through the calculation of the optical absorption spectra both in a indipendent-particle and many-body approach, we show that the effect of local fields is crucial for describing properly the optical behavior of the crystalline case while it is of minor importance for amorphous systems.

The present investigation reports for the first time a detailed theoretical analysis of the optical absorption spectra of corannulene-based materials using state-of-the-art first-principles many-body GW-BSE theory. The study specifically addresses the nature of optical excitations for predictions regarding suitability for device fabrication. The well-defined structure-correlation relationship in functionalized corannulenes is used in a focused investigation of the predicted optoelectronic properties in both the isolated state and bulk crystals. The findings suggest that the excitonic properties are strongly dependent on the specific substituent group as well as the crystalline arrangement. Arylethynyl-substituted corannulene derivatives are shown to be the most suitable for device purposes.

No abstract

No abstract

One of the most promising techniques used for studying the electronic properties of materials is based on Density Functional Theory (DFT) approach and its extensions. DFT has been widely applied in traditional solid state physics problems where periodicity and symmetry play a crucial role in reducing the computational workload. With growing compute power capability and the development of improved DFT methods, the range of potential applications is now including other scientific areas such as Chemistry and Biology. However, cross disciplinary combinations of traditional Solid-State Physics, Chemistry and Biology drastically improve the system complexity while reducing the degree of periodicity and symmetry. Large simulation cells containing of hundreds or even thousands of atoms are needed to model these kind of physical systems. The treatment of those systems still remains a computational challenge even with modern supercomputers. In this paper we describe our work to improve the scalability of Quantum ESPRESSO [1] for treating very large cells and huge numbers of electrons. To this end we have introduced an extra level of parallelism, over electronic bands, in three kernels for solving computationally expensive problems: the Sternheimer equation solver (Nuclear Magnetic Resonance, package QE-GIPAW), the Fock operator builder (electronic ground-state, package PWscf) and most of the Car-Parrinello routines (Car-Parrinello dynamics, package CP). Final benchmarks show our success in computing the Nuclear Magnetic Response (NMR) chemical shift of a large biological assembly, the electronic structure of defected amorphous silica with hybrid exchange-correlation functionals and the equilibrium atomic structure of height Porphyrins anchored to a Carbon Nanotube, on many thousands of CPU cores.

A first-principles investigation of E centers in vitreous silica (v-SiO 2 ) based on calculations of the electron paramagnetic resonance (EPR) parameters is presented. The EPR parameters are obtained by exploiting the gauge including projector augmented wave method as implemented in the QUANTUM-ESPRESSO package. First, we analyze the EPR parameters of a large number of Si 2 dimers. The g tensor of the Si 2 dimers is shown to possess an average rhombic symmetry and larger g principal values with respect to those observed, e.g., for the E γ center in silica. Furthermore, the g principal values clearly show a linear trend with the Si-Si dimer length. Our results suggest that the Si 2 dimers could correspond to an unidentified paramagnetic center, though occasionally the calculated g principal values of the Si 2 dimer might be compatible with those found experimentally for the E δ center. Next, we generate nondimer configurations by a procedure involving structural relaxations in the subsequent positively charged states. In particular, puckered, unpuckered, doubly puckered, and forward-oriented configurations are generated. The distributions of the calculated EPR parameters of the puckered and unpuckered configurations further support the assignment of the E γ center to an unpaired spin localized at a threefold coordinated silicon dangling bond. Moreover, by analyzing Fermi contacts and g tensors of the puckered and forward-oriented configurations, we suggest the assignment of the E α center to the latter type of configurations. This work also suggests that the differences in the EPR parameters of E α and E γ centers mainly arise from the strained geometry of the silicon dangling bond. In the forward-oriented configurations, one Si-O bond is about 0.2Å longer than the remaining two, whereas in the silicon dangling bond of the puckered and unpuckered configurations, all three bonds have a length of ∼1.6Å each.