Source author record

Jose M. Soler

Jose M. Soler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.mtrl-sci physics.comp-ph physics.chem-ph cond-mat.other

Catalog footprint

What is connected

11works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2020arXiv

SIESTA: recent developments and applications

A review of the present status, recent enhancements, and applicability of the SIESTA program is presented. Since its debut in the mid-nineties, SIESTA's flexibility, efficiency and free distribution has given advanced materials simulation capabilities to many groups worldwide. The core methodological scheme of SIESTA combines finite-support pseudo-atomic orbitals as basis sets, norm-conserving pseudopotentials, and a real-space grid for the representation of charge density and potentials and the computation of their associated matrix elements. Here we describe the more recent implementations on top of that core scheme, which include: full spin-orbit interaction, non-repeated and multiple-contact ballistic electron transport, DFT+U and hybrid functionals, time-dependent DFT, novel reduced-scaling solvers, density-functional perturbation theory, efficient Van der Waals non-local density functionals, and enhanced molecular-dynamics options. In addition, a substantial effort has been made in enhancing interoperability and interfacing with other codes and utilities, such as Wannier90 and the second-principles modelling it can be used for, an AiiDA plugin for workflow automatization, interface to Lua for steering SIESTA runs, and various postprocessing utilities. SIESTA has also been engaged in the Electronic Structure Library effort from its inception, which has allowed the sharing of various low level libraries, as well as data standards and support for them, in particular the PSML definition and library for transferable pseudopotentials, and the interface to the ELSI library of solvers. Code sharing is made easier by the new open-source licensing model of the program. This review also presents examples of application of the capabilities of the code, as well as a view of on-going and future developments.

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.

preprint2019arXiv

Band Unfolding Made Simple

We present a simple view on band unfolding of the energy bands obtained from supercell calculations. It relies on the relationship between the local density of states in reciprocal space (qLDOS) and the fully unfolded band structure. This provides an intuitive and valid approach not only for periodic, but also for systems with no translational symmetry. By refolding into the primitive Brillouin zone of the pristine crystal we recover the conventional unfolded bands. We implement our algorithm in the SIESTA package and apply it to defects on Si and graphene.

preprint2016arXiv

Optimization of an exchange-correlation density functional for water

We describe a method, that we call data projection onto parameter space (DPPS), to optimize an energy functional of the electron density, so that it reproduces a dataset of experimental magnitudes. Our scheme, based on Bayes theorem, constrains the optimized functional not to depart unphysically from existing ab initio functionals. The resulting functional maximizes the probability of being the \correct" parametrization of a given functional form, in the sense of Bayes theory. The application of DPPS to water sheds new light on why density functional theory has performed rather poorly for liquid water, on what improvements are needed, and on the intrinsic limitations of the generalized gradient approximation to electron exchange and correlation. Finally, we present tests of our water-optimized functional, that we call vdW-DF-w, showing that it performs very well for a variety of condensed water systems.

preprint2011arXiv

Ab initio energetics and kinetics study of H_2 and CH_4 in the SI Clathrate Hydrate

We present ab initio results at the density functional theory level for the energetics and kinetics of H_2 and CH_4 in the SI clathrate hydrate. Our results complement a recent article by some of the authors [G. Roman-Perez et al., Phys. Rev. Lett. 105, 145901 (2010)] in that we show additional results of the energy landscape of H_2 and CH_4 in the various cages of the host material, as well as further results for energy barriers for all possible diffusion paths of H_2 and CH_4 through the water framework. We also report structural data of the low-pressure phase SI and the higher-pressure phases SII and SH.

preprint2010arXiv

Density, structure and dynamics of water: the effect of Van der Waals interactions

It is known that ab initio molecular dynamics (AIMD) simulations of liquid water, based on the generalized gradient approximation (GGA) to density functional theory (DFT), yield structural and diffusive properties in reasonable agreement with experiment only if artificially high temperatures are used in the simulations. The equilibrium density, at normal conditions, of DFT water has been recently shown by Schmidt et al. [J. Phys. chem. B, 113, 11959 (2009)] to be underestimated by different GGA functionals for exchange and correlation, and corrected by the addition of interatomic pair potentials to describe van derWaals (vdW) interactions. In this contribution we present a DFTAIMD study of liquid water using several GGA functionals as well as the van der Waals density functional (vdW-DF) of Dion et al. [Phys. Rev. Lett. 92, 246401(2004)]. As expected, we find that the density of water is grossly underestimated by GGA functionals. When a vdW-DF is used, the density improves drastically and the experimental diffusivity is reproduced without the need of thermal corrections. We analyze the origin of the density differences between all the functionals. We show that the vdW-DF increases the population of non-H-bonded interstitial sites, at distances between the first and second coordination shells. However, it excessively weakens the H-bond network, collapsing the second coordination shell. This structural problem is partially associated to the choice of GGA exchange in the vdW-DF. We show that a different choice for the exchange functional is enough to achieve an overall improvement both in structure and diffusivity.

preprint2010arXiv

Stability, Adsorption and Diffusion of CH4, CO2 and H2 in Clathrate Hydrates

We present a study of the adsorption and diffusion of CH4, CO2 and H2 molecules in clathrate hydrates using ab initio van der Waals density functional formalism [Dion et al. Phys. Rev. Lett. 92, 246401 (2004)]. We find that the adsorption energy is dominated by van der Waals interactions and that, without them, gas hydrates would not be stable. We calculate the maximum adsorption capacity as well as the maximum hydrocarbon size that can be adsorbed.The relaxation of the host lattice is essential for a good description of the diffusion activation energies, which are estimated to be of the order of 0.2, 0.4, and 1.0 eV for H2, CO2, and CH4, respectively.

preprint2008arXiv

Efficient implementation of a van der Waals density functional: Application to double-wall carbon nanotubes

We present an efficient implementation of the van der Waals density functional of Dion et al [Phys. Rev. Lett. 92, 246401 (2004)], which expresses the nonlocal correlation energy as a double spacial integral. We factorize the integration kernel and use fast Fourier transforms to evaluate the selfconsistent potential, total energy, and atomic forces, in N log(N) operations. The resulting overhead in total computational cost, over semilocal functionals, is very moderate for medium and large systems. We apply the method to calculate the binding energies and the barriers for relative translation and rotation in double-wall carbon nanotubes.

preprint2008arXiv

Optimal Fourier filtering of a function that is strictly confined within a sphere

We present an alternative method to filter a distribution, that is strictly confined within a sphere of given radius $r_c$, so that its Fourier transform is optimally confined within another sphere of radius $k_c$. In electronic structure methods, it can be used to generate optimized pseudopotentials, pseudocore charge distributions, and pseudo atomic orbital basis sets.

preprint2001arXiv

Efficient index handling of multidimensional periodic boundary conditions

An efficient method is described to handle mesh indexes in multidimensional problems like numerical integration of partial differential equations, lattice model simulations, and determination of atomic neighbor lists. By creating an extended mesh, beyond the periodic unit cell, the stride in memory between equivalent pairs of mesh points is independent of their position within the cell. This allows to contract the mesh indexes of all dimensions into a single index, avoiding modulo and other implicit index operations.

preprint1999arXiv

Linear-scaling ab-initio calculations for large and complex systems

A brief review of the SIESTA project is presented in the context of linear-scaling density-functional methods for electronic-structure calculations and molecular-dynamics simulations of systems with a large number of atoms. Applications of the method to different systems are reviewed, including carbon nanotubes, gold nanostructures, adsorbates on silicon surfaces, and nucleic acids. Also, progress in atomic-orbital bases adapted to linear-scaling methodology is presented.