Search NASA⌕ Search

SEARCH · Search NASA

Results for “quantum espresso”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Advanced capabilities for materials modelling with Quantum ESPRESSO

Quantum ESPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the-art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudopotential and projector-augmented-wave approaches. Quantum ESPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement their ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Improved treatment of exact exchange in Quantum ESPRESSO

Here, we present an algorithm and implementation for the parallel computation of exact exchange in Quantum ESPRESSO (QE) that exhibits greatly improved strong scaling. QE is an open-source software package for electronic structure calculations using plane wave density functional theory, and supports the use of local, semi-local, and hybrid DFT functionals. Wider application of hybrid functionals is desirable for the improved simulation of electronic band energy alignments and thermodynamic properties, but the computational complexity of evaluating the exact exchange potential limits the practical application of hybrid functionals to large systems and requires efficient implementations. We demonstrate that existing implementations of hybrid DFT that utilize a single data structure for both the local and exact exchange regions of the code are significantly limited in the degree of parallelization achievable. We present a band-pair parallelization approach, in which the calculation of exact exchange is parallelized and evaluated independently from the parallelization of the remainder of the calculation, with the wavefunction data being efficiently transformed on-the-fly into a form that is optimal for each part of the calculation. For a 64 water molecule supercell, our new algorithm reduces the overall time to solution by nearly an order of magnitude.

97 MATHEMATICS AND COMPUTING↗

Real-time time-dependent density functional theory in Quantum Espresso and related codes

A set of subroutines that interact with an existing computer program, Quantum Espresso (QE), that takes a set of electronic orbitals from QE and propagates them forward in time in response to a number of perturbations. These perturbations include the motion of ions and a wide spectrum of applied electromagnetic fields (e.g., static magnetic fields, slowly-varying electric fields, x-rays, etc.). From the time evolution of the electronic orbitals comes time-varying values of physical observables like the electronic charge and current densities and atomic forces. These are further post-processed to yield properties of interest like stopping powers, conductivities, and the dynamic structure factor. The electronic orbitals are described in a plane wave basis and the equations of motion are carried out using a unitary or approximately unitary time propagation algorithm. Related work will include subroutines for interfacing this capability, as well as others that already exist within QE, with a fork of QE that is more extensible and modular.SAND2019-12611 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

Baczewski, Andrew↗

A Performance Study of Quantum ESPRESSO’s PWscf Code on Multi-core and GPU Systems

We describe the porting of PWscf (Plane-Wave Self Consistent Field), a key component of the Quantum ESPRESSO open-source suite of codes for materials modeling, to GPU systems using CUDA Fortran. Kernel loop directives (CUF kernels) have been extensively used in order to have a single source code for both CPU and GPU implementations. The results of the GPU version have been carefully validated and the performance of the code on several GPU systems (both x86 and POWER8 based) has been compared with traditional Intel multi-core (CPU only) systems. This current GPU version can reduce the time-to-solution by an average factor of 2–3 running two different input cases widely used as benchmarks on small and large high performance computing systems.

CUDA Fortran↗

Nexus: a modular workflow management system for quantum simulation codes

The management of simulation workflows is a significant task for the individual computational researcher. Automation of the required tasks involved in simulation work can decrease the overall time to solution and reduce sources of human error. A new simulation workflow management system, Nexus, is presented to address these issues. Nexus is capable of automated job management on workstations and resources at several major supercomputing centers. Its modular design allows many quantum simulation codes to be supported within the same framework. Current support includes quantum Monte Carlo calculations with QMCPACK, density functional theory calculations with Quantum Espresso or VASP, and quantum chemical calculations with GAMESS. Users can compose workflows through a transparent, text-based interface, resembling the input file of a typical simulation code. A usage example is provided to illustrate the process.

97 MATHEMATICS AND COMPUTING↗

eQE: An open‐source density functional embedding theory code for the condensed phase

Abstract In this work, we present the main features and algorithmic details of a novel implementation of the frozen density embedding (FDE) formulation of subsystem density functional theory (DFT) that is specifically designed to enable ab initio molecular dynamics (AIMD) simulations of large‐scale condensed‐phase systems containing 1000s of atoms. This code (available at http://eqe.rutgers.edu ) has been given the moniker of embedded Quantum ESPRESSO (eQE) as it is a generalization of the open‐source Quantum ESPRESSO (QE) suite of programs. The strengths of eQE reside in a hierarchical parallelization scheme that allows for an efficient and fully self‐consistent treatment of the electronic structure (via the addition of an additional DIIS extrapolation layer) while simultaneously exploiting the inherent symmetries and periodicities in the system (via sampling of subsystem‐specific first Brillouin zones and utilization of subsystem‐specific basis sets). While bulk liquids and molecular crystals are two classes of systems that exemplify the utility of the FDE approach (as these systems can be partitioned into weakly interacting subunits), we show that eQE has significantly extended this regime of applicability by outperforming standard semilocal Kohn–Sham DFT (KS‐DFT) for large‐scale heterogeneous catalysts with quite different layer‐specific electronic structure and intrinsic periodicities. eQE features very favorable strong parallel scaling for a model system of bulk liquid water composed of 256 water molecules, which allows for a significant decrease in the overall time to solution when compared to KS‐DFT. We show that eQE achieves speedups greater than one order of magnitude ( ) when performing AIMD simulations of such large‐scale condensed‐phase systems as: (1) molecular liquids via bulk liquid water represented by 1024 independent water molecules (3072 atoms with a 25.3× speedup over KS‐DFT), (2) polypeptide/biomolecule solvation via ( gly ) 6 solvated in (H 2 O) 395 (1230 atoms with a 38.6× speedup over KS‐DFT), and (3) molecular crystals via a 3 × 3 × 3 periodic supercell of pentacene (1940 atoms with a 12.0× speedup over KS‐DFT). These results represent a significant improvement over the current state‐of‐the‐art and now enable subsystem DFT‐based AIMD simulations of realistically sized condensed‐phase systems of interest throughout chemistry, physics, and materials science.

Genova, Alessandro↗

RIKEN TRIP Magnets Database

This dataset contains ab-initio calculation results for the temperature-dependent anomalous Hall conductivity, the anomalous Nernst effect, and the Seebeck coefficient. All calculations are based on ab-inito Quantum Espresso (PWSCF v.6.3) + Wannier90 (v.3.0.0). The dependence on carrier doping is also calculated. For all calculations a ferromagnetic order has been assumed, which might not correspond to the true ground state of the system. Tabulated values for the magnetic moments and essential input files for Quantum Espresso are available for download as attachments. This project has been supported by the RIKEN Transformative Research Innovation Platform (TRIP), Use Case: Many-body Electron Systems.

36 MATERIALS SCIENCE↗

Evolutionary optimization of PAW data-sets for accurate high pressure simulations

We examine the challenge of performing accurate electronic structure calculations at high pressures by comparing the results of all-electron full potential linearized augmented-plane-wave calculations, as implemented in the WIEN2k code, with those of the projector augmented wave (PAW) method, as implemented in Quantum ESPRESSO or Abinit code. In particular, we focus on developing an automated and consistent way of generating transferable PAW data-sets that can closely produce the all electron equation of state defined from zero to arbitrary high pressures. The technique we propose is an evolutionary search procedure that exploits the ATOMPAW code to generate atomic data-sets and the Quantum ESPRESSO software suite for total energy calculations. We demonstrate different aspects of its workability by optimizing PAW basis functions of some elements relatively abundant in planetary interiors. In addition, we introduce a new measure of atomic data-set goodness by considering their performance uniformity over an extended pressure range. - Highlights: • An attempt made to achieve the AE-FLAPW accuracy, at the favorable computational efficiency of the PAW approach for an outspread pressure range. • Provided an automatic and consistent route for developing efficient and transferable PAW data-sets. • Explicated an improved measure for the goodness of data-sets by including uniformity of the data-set performance over an extended pressure range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Eigensolver performance comparison on Cray XC systems

Summary Hermitian (symmetric) eigenvalue solvers are the core constituents of electronic structure, quantum‐chemistry, and other HPC applications such as Quantum ESPRESSO, VASP, CP2K, and NWChem to name a few. Our understanding of the performance of symmetric eigenvalue algorithms on various hardware is clearly important to the quantum chemistry or condensed matter physics community but in fact goes beyond that community. For instance, big data analytics is increasingly utilizing eigenvalues solvers, in the study of randomized singular value decomposition (SVD) or principal component analysis (PCA). Noise, vibration, and harshness (NVH) is another field where fast and efficient eigenvalue solvers are required. Most eigenvalue solver packages feature numerous different parameters which can be tuned for performance, eg, the number of nodes, number of total ranks, the decomposition of the matrix, etc. In this paper, we investigate the performance of different packages as well as the influence of these knobs on the solver performance.

Cook, Brandon↗

Many-body electronic structure of LaScO 3 by real-space quantum Monte Carlo

We present real space quantum Monte Carlo (QMC) calculations of the scandate LaSc0 3 that proved to be challenging for traditional electronic structure approaches due to strong correlation effects resulting in inaccurate band gaps from DFT and GW methods when compared with existing experimental data. Besides calculating an accurate QMC band gap corrected for supercell size biases and in agreement with numerous experiments, we also predict the cohesive energy of the crystal using the standard fixed-node QMC without any empirical or non-variational parameters. We show that promotion (optical) gap and fundamental gap agree with each other illustrating a clear absence of significant excitonic effects in the ideal crystal. We obtained these results in perfect consistency in two independent tracks that employ different basis sets (plane wave vs. localized gaussians), different codes for generating orbitals (QUANTUM ESPRESSO vs. CRYSTAL), different QMC codes (QMCPACK vs. QWALK) and different high-accuracy pseudopotentials (ccECPs vs. Troullier-Martins) presenting the maturity and consistency of QMC methodology and tools for studies of strongly correlated problems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

IrRep: Symmetry eigenvalues and irreducible representations of ab initio band structures

Here, we present IrRep – a Python code that calculates the symmetry eigenvalues of electronic Bloch states in crystalline solids and the irreducible representations under which they transform. As input it receives bandstructures computed with state-of-the-art Density Functional Theory codes such as VASP, Quantum Espresso, or Abinit, as well as any other code that has an interface to Wannier90. Our code is applicable to materials in any of the 230 space groups and double groups preserving time-reversal symmetry with or without spin-orbit coupling included, for primitive or conventional unit cells. This makes IrRep a powerful tool to systematically analyze the connectivity and topological classification of bands, as well as to detect insulators with non-trivial topology, following the Topological Quantum Chemistry formalism: IrRep can generate the input files needed to calculate the (physical) elementary band representations and the symmetry-based indicators using the [CheckTopologicalMat: https://www.cryst.ehu.es/cgi-bin/cryst/programs/magnetictopo.pl] routine of the Bilbao Crystallographic Server. It is also particularly suitable for interfaces with other plane-waves based codes, due to its flexible structure.

97 MATHEMATICS AND COMPUTING↗

Temperature dependent structural behavior of Mg{sub 70}Zn{sub 30}: A molecular dynamics study

Mg{sub 70}Zn{sub 30} is simplest of metallic glasses having good mechanical properties, biocompatibility with human body and is biodegradable, so it is widely studied material in its class. We present in this paper first principle calculations showing the change in the structure of the system with temperature. The computations are performed using PWSCF tool of Quantum ESPRESSO package, which uses density functional theory for calculations of electronic ground state properties and dynamical structures. With increasing temperature the number of peaks in RDFs reduce in number leaving only a single peak,the system thus makes a transition from non-crystalline solid to liquid state within the range of 300–1000K. First coordination number was found to decrease with rising temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗

Ab Initio Calculations of New Full Heusler Alloys Rh{sub 2}ZrX (X = Al, Ga, In, Si, Ge, Sn)

In recent years, Heusler alloys have attracted extensive attention because of their usefulness in spintronic devices. In this paper, the structural, electronic, and magnetic properties of full Heusler alloys Rh{sub 2}ZrX (X = Al, Ga, In, Si, Ge, Sn) have been studied using Quantum Espresso software package based on density functional theory. Calculations have been done by GGA and GGA+U methods. The investigations on the structural properties of the alloys showed that the Rh{sub 2}ZrSi alloy had the most stable structure (E{sub coh} = − 1.807 eV) among the six compounds reviewed. The results of GGA calculation on electronic properties of Rh{sub 2}ZrX alloys show that while the Rh{sub 2}ZrX (X = Al, Ga, In) alloys had metallic nature due to overlapping between conduction and valance band at Fermi level for both spin directions, the Rh{sub 2}ZrX (X = Si, Ge, Sn) alloys had half-metal characteristic. Also, Rh{sub 2}ZrSi and Rh{sub 2}ZrGe alloys had 100% spin polarization. By studying the magnetic properties, it was found that the Rh{sub 2}ZrAl and Rh{sub 2}ZrSi alloys were antiferromagnetic, while the Rh{sub 2}ZrGa, Rh{sub 2}ZrIn, Rh{sub 2}ZrGe, and Rh{sub 2}ZrSn alloys were ferromagnetic materials. Also, the total spin magnetic moment of alloys was found in accordance with the Slater-Pauling rule. The results obtained using GGA+U calculation are slightly different from the results of GGA calculation.

36 MATERIALS SCIENCE↗

Electronic and Magnetic Properties of Mn-doped and (Mn,C)-codoped w-AlN with the Presence of N Vacancy

In this paper, we discover the structural, electronic, and magnetic properties of Mn-doped (6.25 and 12.5%) and Mn-C-codoped wurtzite AlN with the presence of nitrogen vacancy. These properties have been investigated using density functional theory implemented in Quantum ESPRESSO package. AlN doped and codoped were optimized by using Broyden-Fletcher-Goldfarb-Shanno (BFGS) approximation. The electronic and magnetic properties are determined by the total and partial density of states. The presence of impurities and nitrogen vacancy induces magnetic properties, indicating that these materials are good candidates for electromagnetic and spintronics applications. Our results are in good agreement with experimental and other theoretical studies. We can conclude from this study that this material can be used as diluted magnetic semiconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Exploring the role of electronic structure on photo-catalytic behavior of carbon-nitride polymorphs

A fully self-consistent density-functional theory (DFT) with improved functionals is used herein to provide a comprehensive account of structural, electronic, and optical properties of C 3 N 4 polymorphs. Using our recently developed van Leeuwen-Baerends (vLB) corrected local-density approximation (LDA), we implemented LDA + vLB within full-potential N th -order muffin-tin orbital (FP-NMTO) method and show that it improves structural properties and band gaps compared to semi-local functionals (LDA/GGA). We demonstrate that the LDA + vLB predicts band-structure and work-function for well-studied 2D-graphene and bulk-Si in very good agreement with experiments, and more exact hybrid functional (HSE) calculations as implemented in the Quantum-Espresso (QE) package. The structural and electronic-structure (band gap) properties of C 3 N 4 polymorphs calculated using FP-NMTO-LDA + vLB is compared with more sophisticated hybrid-functional calculations. We also perform detailed investigation of photocatalytic behavior using QE-HSE method of C 3 N 4 polymorphs through work-function, band (valence and conduction) position with respect to water reduction and oxidation potential. Our results show $γ$-C 3 N 4 as the best candidate for photocatalysis among all the C 3 N 4 polymorphs but it is dynamically unstable at ‘zero’ pressure. We show that $γ$-C 3 N 4 can be stabilized under hydrostatic-pressure, which improves its photocatalytic behavior relative to water reduction and oxidation potentials.

36 MATERIALS SCIENCE↗