Source author record

Stefano de Gironcoli

Stefano de Gironcoli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Compact atomic descriptors enable accurate predictions via linear models

We probe the accuracy of linear ridge regression employing a three-body local density representation derived from the atomic cluster expansion. We benchmark the accuracy of this framework in the prediction of formation energies and atomic forces in molecules and solids. We find that such a simple regression framework performs on par with state-of-the-art machine learning methods which are, in most cases, more complex and more computationally demanding. Subsequently, we look for ways to sparsify the descriptor and further improve the computational efficiency of the method. To this aim, we use both principal component analysis and least absolute shrinkage operator regression for energy fitting on six single-element datasets. Both methods highlight the possibility of constructing a descriptor that is four times smaller than the original with a similar or even improved accuracy. Furthermore, we find that the reduced descriptors share a sizable fraction of their features across the six independent datasets, hinting at the possibility of designing material-agnostic, optimally compressed, and accurate descriptors.

preprint2022arXiv

Data-driven simulation and characterisation of gold nanoparticle melting

The simulation and analysis of the thermal stability of nanoparticles, a stepping stone towards their application in technological devices, require fast and accurate force fields, in conjunction with effective characterisation methods. In this work, we develop efficient, transferable, and interpretable machine learning force fields for gold nanoparticles based on data gathered from Density Functional Theory calculations. We use them to investigate the thermodynamic stability of gold nanoparticles of different sizes (1 to 6 nm), containing up to 6266 atoms, concerning a solid-liquid phase change through molecular dynamics simulations. We predict nanoparticle melting temperatures in good agreement with available experimental data. Furthermore, we characterize the solid-liquid phase change mechanism employing an unsupervised learning scheme to categorize local atomic environments. We thus provide a data-driven definition of liquid atomic arrangements in the inner and surface regions of a nanoparticle and employ it to show that melting initiates at the outer layers.

preprint2022arXiv

Mathematical Analysis and Numerical Approximations of Density Functional Theory Models for Metallic Systems

In this paper, we investigate the energy minimization model of the ensemble Kohn-Sham density functional theory for metallic systems, in which a pseudo-eigenvalue matrix and a general smearing approach are involved. We study the invariance and the existence of the minimizer of the energy functional. We propose an adaptive double step size strategy and the corresponding preconditioned conjugate gradient methods for solving the energy minimization model. Under some mild but reasonable assumptions, we prove the global convergence of our algorithms. Numerical experiments show that our algorithms are efficient, especially for large scale metallic systems. In particular, our algorithms produce convergent numerical approximations for some metallic systems, for which the traditional self-consistent field iterations fail to converge.

preprint2020arXiv

Hubbard-corrected density functional perturbation theory with ultrasoft pseudopotentials

We present in full detail a newly developed formalism enabling density functional perturbation theory (DFPT) calculations from a DFT+$U$ ground state. The implementation includes ultrasoft pseudopotentials and is valid for both insulating and metallic systems. It aims at fully exploiting the versatility of DFPT combined with the low-cost DFT+$U$ functional. This allows to avoid computationally intensive frozen-phonon calculations when DFT+$U$ is used to eliminate the residual electronic self-interaction from approximate functionals and to capture the localization of valence electrons e.g. on $d$ or $f$ states. In this way, the effects of electronic localization (possibly due to correlations) are consistently taken into account in the calculation of specific phonon modes, Born effective charges, dielectric tensors and in quantities requiring well converged sums over many phonon frequencies, as phonon density of states and free energies. The new computational tool is applied to two representative systems, namely CoO, a prototypical transition metal monoxide and LiCoO$_2$, a material employed for the cathode of Li-ion batteries. The results show the effectiveness of our formalism to capture in a quantitatively reliable way the vibrational properties of systems with localized valence electrons.

preprint2020arXiv

The CECAM Electronic Structure Library and the modular software development paradigm

First-principles electronic structure calculations are very widely used thanks to the many successful software packages available. Their traditional coding paradigm is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, from the compiler up, with the exception of linear-algebra and message-passing libraries. This model has been quite successful for decades. The rapid progress in methodology, however, has resulted in an ever increasing complexity of those programs, which implies a growing amount of replication in coding and in the recurrent re-engineering needed to adapt to evolving hardware architecture. The Electronic Structure Library (\esl) was initiated by CECAM (European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure programs and redesign them as free, open-source libraries. They include "heavy-duty" ones with a high degree of parallelisation, and potential for adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by scientists when implementing new ideas. It is a community effort, undertaken by developers of various successful codes, now facing the challenges arising in the new model. This modular paradigm will improve overall coding efficiency and enable specialists (computer scientists or computational scientists) to use their skills more effectively. It will lead to a more sustainable and dynamic evolution of software as well as lower barriers to entry for new developers.

preprint2016arXiv

Crystal Structure Prediction of Molecular Crystals from First Principles: Are we there yet?

Accurate molecular crystal structure prediction is a fundamental goal in academic and industrial condensed matter research and polymorphism is arguably the biggest obstacle on the way. We tackle this challenge in the difficult case of the repeatedly studied, abundantly used aminoacid Glycine that hosts still little-known phase transitions and we illustrate the current state of the field through this example. We demonstrate that the combination of recent progress in structure search algorithms with the latest advances in the description of van der Waals interactions in Density Functional Theory, supported by data-mining analysis, enables a leap in predictive power: we resolve, without prior empirical input, all known phases of glycine, as well as the structure of the previously unresolved $ζ$ phase after a decade of its experimental observation [Boldyreva et al. \textit{Z. Kristallogr.} \textbf{2005,} \textit{220,} 50-57]. The search for the well-established $α$ phase instead reveals the remaining challenges in exploring a polymorphic landscape.

preprint2016arXiv

Isobaric first-principles molecular dynamics of liquid water with nonlocal van der Waals interactions

We investigate the structural properties of liquid water at near ambient conditions using first-principles molecular dynamics simulations based on a semilocal density functional augmented with nonlocal van der Waals interactions. The adopted scheme offers the advantage of simulating liquid water at essentially the same computational cost of standard semilocal functionals. Applied to the water dimer and to ice Ih, we find that the hydrogen-bond energy is only slightly enhanced compared to a standard semilocal functional. We simulate liquid water through molecular dynamics in the NpH statistical ensemble allowing for fluctuations of the system density. The structure of the liquid departs from that found with a semilocal functional leading to more compact structural arrangements. This indicates that the directionality of the hydrogen-bond interaction has a diminished role as compared to the overall attractions, as expected when dispersion interactions are accounted for. This is substantiated through a detailed analysis comprising the study of the partial radial distribution functions, various local order indices, the hydrogen-bond network, and the selfdiffusion coefficient. The explicit treatment of the van der Waals interactions leads to an overall improved description of liquid water.

preprint2016arXiv

Molecular bonding with the RPAx: from weak dispersion forces to strong correlation

In a recent article [Phys. Rev. B 90, 125102 (2014)] we showed that the random phase approximation with exchange (RPAx) gives accurate total energies for a diverse set of systems including the high and low density regime of the homogeneous electron gas, the N$_2$ molecule, and the H$_2$ molecule at dissociation. In this work, we present results for the van der Waals bonded Ar$_2$ and Kr$_2$ dimers and demonstrate that the RPAx gives superior dispersion forces as compared to the RPA. We then show that this improved description is crucial for the bond formation of the Mg$_2$ molecule. In addition, the RPAx performs better for the Be$_2$ dissociation curve at large nuclear separation but, similar to the RPA, fails around equilibrium due to the build up of a large repulsion hump. For the strongly correlated LiH molecule at dissociation we have also calculated the RPAx potential and find that the correlation peak at the bond mid-point is overestimated as compared to the RPA and the exact result. The step feature is missing and hence the delocalisation error is comparable to the RPA. This is further illustrated by a smooth energy versus fractional charge curve and a poor description of the LiH dipole moment at dissociation.

preprint2014arXiv

Correlation energy within exact-exchange ACFD theory: systematic development and simple approximations

We have calculated the correlation energy of the homogeneous electron gas (HEG) and the dissociation energy curves of molecules with covalent bonds from a novel implementation of the adiabatic connection fluctuation dissipation (ACFD) expression including the exact exchange (EXX) kernel. The EXX kernel is defined from first order perturbation theory and used in the Dyson equation of time-dependent density functional theory. Within this approximation (RPAx), the correlation energies of the HEG are significantly improved with respect to the RPA up to densities of the order of $r_s \approx 10$. However, beyond this value, the RPAx response function exhibits an unphysical divergence and the approximation breaks down. Total energies of molecules at equilibrium are also highly accurate but we find a similar instability at stretched geometries. Staying within an exact first order approximation to the response function we use an alternative resummation of the higher order terms. This slight redefinition of RPAx fixes the instability in total energy calculations without compromising the overall accuracy of the approach.

preprint2013arXiv

Hubbard-corrected DFT energy functionals: the LDA+U description of correlated systems

The aim of this review article is to assess the descriptive capabilities of the Hubbard-rooted LDA+U method and to clarify the conditions under which it can be expected to be most predictive. The paper illustrates the theoretical foundation of LDA+U and prototypical applications to the study of correlated materials, discusses the most relevant approximations used in its formulation, and makes a comparison with other approaches also developed for similar purposes. Open "issues" of the method are also discussed, including the calculation of the electronic couplings (the Hubbard U), the precise expression of the corrective functional and the possibility to use LDA+U for other classes of materials. The second part of the article presents recent extensions to the method and illustrates the significant improvements they have obtained in the description of several classes of different systems. The conclusive section finally discusses possible future developments of LDA+U to further enlarge its predictive power and its range of applicability.

preprint2011arXiv

{\it Ab initio} $^{27}Al$ NMR chemical shifts and quadrupolar parameters for $Al_2O_3$ phases and their precursors

The Gauge-Including Projector Augmented Wave (GIPAW) method, within the Density Functional Theory (DFT) Generalized Gradient Approximation (GGA) framework, is applied to compute solid state NMR parameters for $^{27}Al$ in the $α$, $θ$, and $κ$ aluminium oxide phases and their gibbsite and boehmite precursors. The results for well-established crystalline phases compare very well with available experimental data and provide confidence in the accuracy of the method. For $γ$-alumina, four structural models proposed in the literature are discussed in terms of their ability to reproduce the experimental spectra also reported in the literature. Among the considered models, the $Fd\bar{3}m$ structure proposed by Paglia {\it et al.} [Phys. Rev. B {\bf 71}, 224115 (2005)] shows the best agreement. We attempt to link the theoretical NMR parameters to the local geometry. Chemical shifts depend on coordination number but no further correlation is found with geometrical parameters. Instead our calculations reveal that, within a given coordination number, a linear correlation exists between chemical shifts and Born effective charges.

preprint2001arXiv

Role of defects in the electronic properties of amorphous/crystalline Si interface

The mechanism determining the band alignment of the amorphous/crystalline Si heterostructures is addressed with direct atomistic simulations of the interface performed using a hierarchical combination of various computational schemes ranging from classical model-potential molecular dynamics to ab-initio methods. We found that in coordination defect-free samples the band alignment is almost vanishing and independent on interface details. In defect-rich samples, instead, the band alignment is sizeably different with respect to the defect-free case, but, remarkably, almost independent on the concentration of defects. We rationalize these findings within the theory of semiconductor interfaces.