Plane-wave Density Functional Theory Research Articles

Implementation of frozen density embedding in CP2K and OpenMolcas: CASSCF wavefunctions embedded in a Gaussian and plane wave DFT environment.

Most chemical processes happen at a local scale where only a subset of molecular orbitals is directly involved and only a subset of covalent bonds may be rearranged. To model such reactions, Density Functional Theory (DFT) is often inadequate, and the use of computationally more expensive correlated wavefunction (WF) methods is required for accurate results. Mixed-resolution approaches backed by quantum embedding theory have been used extensively to approach this imbalance. Based on the frozen density embedding freeze-and-thaw algorithm, we describe an approach to embed complete active space self-consistent field simulations run in the OpenMolcas code in a DFT environment calculated in CP2K without requiring any external tools. This makes it possible to study a local, active part of a chemical system in a larger and relatively static environment with a computational cost balanced between the accuracy of a WF method and the efficiency of DFT, which we test on environment-subsystem pairs. Finally, we apply the implementation to an oxygen molecule leaving an aluminum (111) surface and a ruthenium(IV) oxide (110) surface.

Generalizing deep learning electronic structure calculation to the plane-wave basis.

Deep neural networks capable of representing the density functional theory (DFT) Hamiltonian as a function of material structure hold great promise for revolutionizing future electronic structure calculations. However, a notable limitation of previous neural networks is their compatibility solely with the atomic-orbital (AO) basis, excluding the widely used plane-wave (PW) basis. Here we overcome this critical limitation by proposing an accurate and efficient real-space reconstruction method for directly computing AO Hamiltonian matrices from PW DFT results. The reconstruction method is orders of magnitude faster than traditional projection-based methods to convert PW results to the AO basis, and the reconstructed Hamiltonian matrices can faithfully reproduce the PW electronic structure, thus bridging the longstanding gap between the AO basis deep learning electronic structure approach and PW DFT. Advantages of the PW methods, such as high accuracy, high flexibility and wide applicability, thus can be all integrated into deep learning electronic structure methods without sacrificing these methods' inherent benefits. This allows for the construction of large-scale and high-fidelity training datasets with the help of PW DFT results towards the development of precise and broadly applicable deep learning electronic structure models.

First-principles study on the mechanism of Hf and Ti doping promoting hydrogen release in ZrCoH3

This paper reveals the mechanism of Hf and Ti doping promoting hydrogen release in ZrCoH3. Different from previous experimental methods, this study employs the plane wave pseudopotential method and density functional theory to predict the relationship between the formation energy and electronic structure of the ZrCoH3-Hf/Ti system, and uses theoretical calculations to guide the modification design of ZrCoH3. Firstly, by calculating the thermodynamic properties of the ZrCoH3-Hf/Ti system and combining the density of states theory, the bonding characteristics and dehydrogenation preference between H atoms and neighboring metal atoms in the ZrCoH3-Hf/Ti system are determined. Additionally, the electronic structure of the ZrCoH3-Hf/Ti system is calculated to elucidate the hybridization of bonding and antibonding orbitals between atoms in the ZrCoH3-Hf/Ti system. Moreover, DCo-d is utilized to measure the strength of antibonding interactions in the ZrCoH3-Hf/Ti system. The main cause of antibonding interactions in the ZrCoH3-Hf/Ti system is attributed to the variation in formation energy, and the change in antibonding interactions is sensitive to the variation in formation energy. Finally, a feasible method to improve the properties of ZrCoH3 is proposed, which involves modifying the electronic structure to increase the partial density of states at the Fermi energy level, making the formation energy of hydrides less negative and forcing the structure to be more unstable. This method has potential application value in the modification of ZrCoH3.

First-principles calculations to investigate phase stability, elastic and thermodynamic properties of TiMoNbX (X=Cr, Ta, Cr and Ta) refractory high entropy alloys

This work uses first-principles calculations to investigate the phase stability, thermophysical and mechanical properties of refractory high entropy alloys (RHEAs) at finite temperatures. On the basis of plane wave quasi-potential and density functional theory, construct the structure model of a solid solution. The TiMoNbX (X = Cr, Ta, Cr and Ta) RHEAs have been determined to preserve a single body-centered cubic solid solution structure by calculations and the equilibrium lattice parameters and elastic modulus are consistent with experimental data obtained by laser cladding, which is combined with TC4 (Ti-6Al-4V) substrate. Using the quasi-harmonic Debye-Grüneisen model, the thermophysical characteristics of three RHEAs are investigated. The Voigt-Reuss-Hill scheme is used for calculating the Young's modulus (E), bulk modulus (B), shear modulus (G), and Poisson's ratio (ν), which indicates that all three RHEAs are ductile materials. Additionally, the modulus and hardness of materials decrease as temperature rises, whereas the properties of TiMoNbX RHEAs are predicted, as the nanoindentation hardness values at room temperature are comparable to, and slightly higher than the calculated values.

Improving the precision of work-function calculations within plane-wave density functional theory

Abstract Work function is a fundamental property of metals and is related to many surface-related phenomena of metals. Theoretically, it can be calculated with a metal slab supercell in density functional theory (DFT) calculations. In this paper, we discuss how the commensurability of atomic structure with the underlying fast Fourier transform (FFT) grid affects the accuracy of work function obtained from plane-wave pseudopotential DFT calculations. We show that the macroscopic average potential, which is an important property in work function calculations under the ‘bulk reference’ method, is more numerically stable when it is calculated with commensurate FFT grids than with incommensurate FFT grids. Due to the stability of the macroscopic average potential, work function calculated with commensurate FFT grids shows better convergence with respect to basis set size, vacuum length and slab thickness of a slab supercell. After we control the FFT grid commensurability issue in our work function calculations, we obtain well-converged work functions for Al, Pd, Au and Pt of (100), (110) and (111) surface orientations. For all the metals considered, the ordering of our calculated work functions of the three surface orientations agrees with experiment. Our findings reveal the importance of the FFT grid commensurability issue, which is usually neglected in practice, in obtaining accurate metal work functions, and are also meaningful to other DFT calculations which can be affected by the FFT grid commensurability issue.

First principles validation of energy barriers in Ni75Al25

Precipitates in nickel-based superalloys form during heat treatment on a time scale inaccessible to direct molecular dynamics simulation, but can be studied using kinetic Monte Carlo (KMC) modelling. This requires reliable values for the barrier energies separating distinct configurations over the trajectory of the system. In this study, we validate vacancy migration barriers found with the Activation-Relaxation Technique nouveau (ARTn) method in partially ordered Ni75Al25 with a monovacancy using published potentials for the atomic interactions against first-principles methods. In a first step, we confirm that the ARTn barrier energies agree with those determined with the nudged elastic band (NEB) method. As the number of atoms used in those calculations is too great for direct ab initio calculations, we cut the cell size to 255 atoms, thus controlling finite size effects. We then use the plane-wave density functional theory code CASTEP and its inbuilt NEB method in the smaller cells. This provides us with a continuous validation chain from first principles to KMC simulations with interatomic potentials (IPs). We evaluate the barrier energies of five further IPs with NEB, demonstrating that none yields values with sufficient reliability for KMC simulations, with some of them failing completely. This is a first step towards quantifying the errors incurred in KMC simulations of precipitate formation and evolution.

Enhanced figure of merit in two-dimensional ZrNiSn nanosheets for thermoelectric applications

A novel two-dimensional (2D) half-Heusler ZrNiSn nanosheet for thermoelectric applications was designed from bulk half-Heusler ZrNiSn through first-principles calculation. Investigation of bulk half-Heusler and 2D nanosheet ZrNiSn was performed with the Quantum Espresso code based on a density functional theory plane wave basis set. Electronic band structure and density of states calculations were used to study the confinement effects. On moving from bulk to 2D a change of structure is observed from face-centered cubic to trigonal due to confinement effects. The semiconducting nature of bulk ZrNiSn is undisturbed while moving to a 2D nanosheet; however, the band gap is widened from 0.46 to 1.3 eV due to the restricted motion of electrons in one direction. Compared with bulk ZrNiSn, 2D nanosheets were found to have a higher Seebeck coefficient a lower thermal conductivity and higher figure of merit, which makes 2D ZrNiSn nanosheets suitable for thermoelectric applications. Atomically thin 2D structures with a flat surface have the potential to form van der Waals heterojunctions, paving the way for device fabrication at the nanoscale level.

Strain-dependent one-dimensional confinement channels in twisted bilayer 1T′−WTe2

The low-symmetry and anistropic lattice of 1T′WTe2 is responsible for the existence of parallel one-dimensional channels in the moiré patterns of twisted bilayers. This gives the opportunity to explore moiré physics of a different nature to that widely observed in twisted bilayers of materials with hexagonal symmetries. Here, we combine plane-wave and linear-scaling density functional theory calculations to describe the electronic properties of twisted bilayer 1T′WTe2. We investigate the strain dependence of both aligned bilayer configurations across a range of stacking geometries, and of the twisted bilayer itself. For a small change in the lattice parameters of the constituent 1T′WTe2 monolayers, we find a very substantial increase in the moiré-induced striped electrostatic potential landscape in the twisted bilayer. We measure a peak-to-trough magnitude >200 meV between parallel linear channels, sufficient to induce Luttinger liquid behavior. Published by the American Physical Society 2024

First principles calculation and corrosion resistance study of Fe60CrxNi40-x medium entropy alloy

In this paper, Fe60CrxNi40-x medium entropy alloys (MEAs) were synthesized through vacuum arc melting. Based on the plane wave pseudo-potential and density functional theory, the solid solution structure is modeled by virtual crystal approximation (VCA). The first principle method was employed to examine how Cr content affects the phase structure, elastic constant, and elastic properties of the Fe60CrxNi40-x. Physical phase analysis of the MEAs by X-ray diffraction reveals that the Fe60CrxNi40-x maintain a single FCC solid solution structure. Meanwhile, the impact of Cr content on MEAs' resistance to corrosion in a Cl- setting was examined and contrasted with 316 L stainless steel. The result indicated that as the Cr content rises, the MEAs exhibit an extended passivation period, a reduced dimensional passivation current, and a more compact and stable passivation film. When considering the aggregate outcomes of pitting potential (Epit), self-corrosion potential (Ecorr), and corrosion current density (icorr) for 316LSS and MEAs, it's evident that Fe60CrxNi40-x in 3.5 wt% NaCl solution exhibits superior corrosion resistance compared to 316 L stainless steel, with Fe3 alloy leading in corrosion resistance among medium entropy alloys.

Enhanced hydrogen adsorption-desorption reversibility found on NiAl alloy: A first-principles study

Improving the hydrogen storage performance of carbon-based materials via transition metals incorporation has been proved experimentally, where control over the hydrogen adsorption kinetics and hydrogen surface diffusion is the key. In this work, we proved theoretically via density functional theory (DFT) how such incorporation promotes better hydrogen storage performance in Ni-based materials supported on graphene, in which static and dynamic behavior of the system were modeled by plane-wave DFT and ab initio molecular dynamics (AIMD), respectively. In addition, the hydrogen storage materials Ni10Al3 and Ni13 supported on graphene, denoted as Ni10Al3/GP and Ni13/GP, are modeled to realize the experimental observation on Ni3Al alloy under high hydrogen pressure. It was found that the Ni10Al3/GP alloy system exhibited better properties in twofold: (1) increased hydrogen adsorption capacity through providing an additional site for H2 molecular adsorption, particularly on the Al atoms, and (2) enhanced hydrogen desorption at ambient conditions observed from reversible H2 adsorption at all coverages with barrierless migration energy. Therefore, the proof from the DFT point-of-view provides valuable insight into exploring other metal alloy systems that can provide additional sites for molecular hydrogen adsorption and improve hydrogen adsorbate mobility.

Effect of Metal Carbides on Hydrogen Embrittlement: A Density Functional Theory Study

This study uses plane wave density functional theory (DFT) to investigate the effect of certain metal carbides (Niobium carbide, Vanadium carbide, Titanium carbide, and Manganese sulfide) on hydrogen embrittlement in pipeline steels. Our results predict that the interaction of hydrogen molecules with these metal carbides occurs in the long range with binding energy varying in the energy window [0.043 eV to 0.70 eV].In addition, our study shows the desorption of H2 molecules from these metal carbides in the chemisorptions. Since atomic state hydrogen interacts with NbC, VC, TiC, and MnS to cause embrittlement, we classified the strength of the hydrogen trapping as TiC + H > VC + H > NbC + H> MnS + H. In addition, our study reveals that the carbon site is a more favorable hydrogen-trapping site than the metal one.

Physical origin of enhanced electrical conduction in aluminum-graphene composites

The electronic and transport properties of aluminum-graphene composite materials were investigated using the ab initio plane wave density functional theory. The interfacial structure is reported for several configurations. In some cases, the face-centered aluminum (111) surface relaxes in a nearly ideal registry with graphene, resulting in a remarkably continuous interface structure. The Kubo–Greenwood formula and space-projected conductivity were employed to study electronic conduction in aluminum single- and double-layer graphene-aluminum composite models. The electronic density of states at the Fermi level is enhanced by the graphene for certain aluminum–graphene interfaces, thus improving electronic conductivity. In double-layer graphene composites, conductivity varies non-monotonically with temperature, showing an increase between 300 and 400 K at short aluminum-graphene distances, unlike the consistent decrease in single-layer composites.

Massively parallel implementation of iterative eigensolvers in large-scale plane-wave density functional theory

The Kohn-sham density functional theory (DFT) is a powerful method to describe the electronic structures of molecules and solids in condensed matter physics, computational chemistry and materials science. However, large and accurate DFT calculations within plane waves process a cubic-scaling computational complexity, which is usually limited by expensive computation and communication costs. The rapid development of high performance computing (HPC) on leadership supercomputers brings new opportunities for developing plane-wave DFT calculations for large-scale systems. Here, we implement parallel iterative eigensolvers in large-scale plane-wave DFT calculations, including Davidson, locally optimal block preconditioned conjugate gradient (LOBPCG), projected preconditioned conjugate gradient (PPCG) and the Chebyshev subspace iteration (CheFSI) algorithms, and analyze the performance of these algorithms in massively parallel plane-wave computing tasks. We adopt a two-level parallelization strategy that combines the message passing interface (MPI) with open multi-processing (OpenMP) parallel programming to handle data exchange and matrix operations in the construction and diagonalization of large-scale Hamiltonian matrix within plane waves. Numerical results illustrate that these iterative eigensolvers can scale up to 42,592 processing cores with high peak performance of 30% on leadship supercomputers to study the electronic structures of bulk silicon systems containing 10,648 atoms. Program summaryProgram Title: Plane wave density functional theory (PWDFT)CPC Library link to program files:https://doi.org/10.17632/c8v2mx5vn4.1Developer's repository link:https://bitbucket.org/berkeleylab/scalesLicensing provisions: BSD 3-clauseProgramming language: C++Nature of problem: PWDFT is used for electronic structure calculations based on Kohn-Sham density functional theory. The key challenge to address is a constrained energy minimization problem, which can also be formulated as a nonlinear eigenvalue problem. MPI/OpenMP-based approaches are employed to provide multi-core acceleration for the study of the chemical and material properties of larger-scale molecules and solids.Solution method: PWDFT implements self-consistent field (SCF) iterations and direct constrained minimization algorithms with various acceleration strategies. It is written in C++ and offers parallel acceleration based on MPI/OpenMP.

Polarization characteristics and structural modifications of Cu nanoparticles under high electric fields

High electric fields affect the diffusion dynamics of atoms on a metal surface, causing biased surface diffusion that possibly leads to the growth of intensively field emitting protrusions and consequent vacuum breakdown (VBD). The scientific understanding of this process, as well as other fundamental VBD initiation mechanisms, is far from complete. Here we investigate the exact atomic behaviour of metal surfaces exposed to extremely high electric fields using density functional theory (DFT). Previous theories describe the field-surface dynamics in terms of the effective dipole moments and polarizability of surface atoms, disregarding higher-order (hyperpolarizability) terms. The validity of this approximation has been evaluated only for electric fields up to 3 GV/m, due to computational limitations of the plane-wave DFT basis used in previous works. In this work, we test the validity of this approximation for a much wider field range, relevant for VBD and field emission (FE), using Cu nanoparticles as our test structures. We find that although such high fields can change the entire structure of Cu nanoparticles, their energetics are described very precisely by the permanent dipole moment and polarizability terms. Thus, we show that neglecting the hyperpolarizability terms is valid even for field values that exceeds the range that is relevant for intense FE and VBD. This work lays a solid foundation for further developing atomic-level simulation models for electric field-induced surface diffusion on metal surfaces and its effects on protrusion growth and VBD initiation.

Development and application of hybrid AIMD/cDFT simulations for atomic-to-mesoscale chemistry.

Many important chemical processes involve reactivity and dynamics in complex solutions. Gaining a fundamental understanding of these reaction mechanisms is a challenging goal that requires advanced computational and experimental approaches. However, important techniques such as molecular simulation have limitations in terms of scales of time, length, and system complexity. Furthermore, among the currently available solvation models, there are very few designed to describe the interaction between the molecular scale and the mesoscale. To help address this challenge, here, we establish a novel hybrid approach that couples first-principles plane-wave density functional theory with classical density functional theory (cDFT). In this approach, a region of interest described by abinitio molecular dynamics (AIMD) interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. cDFT is a robust but efficient mesoscopic approach to accurate thermodynamics of bulk electrolyte solutions over a wide concentration range (up to 2M concentrations). Benchmarking against commonly used continuum models of solvation, such as SMD, as well as experiments, demonstrates that our hybrid AIMD-cDFT method is able to produce reasonable solvation energies for a variety of molecules and ions. With this model, we also examined the solvent effects on a prototype SN2 reaction of the nucleophilic attack of a chloride ion on methyl chloride in the solution. The resulting reaction pathway profile and the solution phase barrier agree well with experiment, showing that our AIMD/cDFT hybrid approach can provide insight into the specific role of the solvent on the reaction coordinate.

Effects of -H and -OH Termination on Adhesion of Si-Si Contacts Examined Using Molecular Dynamics and Density Functional Theory.

The contact between nanoscale single-crystal silicon asperities and substrates terminated with -H and -OH functional groups is simulated using reactive molecular dynamics (MD). Consistent with previous MD simulations for self-mated surfaces with -H terminations only, adhesion is found to be low at full adsorbate coverages, be it self-mated coverages of mixtures of -H and -OH groups, or just -OH groups. As the coverage reduces, adhesion increases markedly, by factors of ∼5 and ∼6 for -H-terminated surfaces and -OH-terminated surfaces, respectively, and is due to the formation of covalent Si-Si bonds; for -OH-terminated surfaces, some interfacial Si-O-Si bonds are also formed. Thus, covalent linkages need to be broken upon separation of the tip and substrate. In contrast, replacing-Hgroups with -OH groups while maintaining complete coverage leads to negligible increases in adhesion. This indicates that increases in adhesion require unsaturated sites. Furthermore, plane-wave density functional theory (DFT) calculations were performed to investigate the energetics of two Si(111) surfaces fully terminated by either -H or -OH groups. Importantly for the adhesion results, both DFT and MD calculations predict the correct trends for the relative bond strengths: Si-O > Si-H > Si-Si. This work supports the contention that prior experimental work observing strong increases in adhesion after sliding Si-Si nanoasperities over each other is due to sliding-induced removal of passivating species on the Si surfaces.

The EFG Rosetta Stone: translating between DFT calculations and solid state NMR experiments.

We present a comprehensive study on the best practices for integrating first principles simulations in experimental quadrupolar solid-state nuclear magnetic resonance (SS-NMR), exploiting the synergies between theory and experiment for achieving the optimal interpretation of both. Most high performance materials (HPMs), such as battery electrodes, exhibit complex SS-NMR spectra due to dynamic effects or amorphous phases. NMR crystallography for such challenging materials requires reliable, accurate, efficient computational methods for calculating NMR observables from first principles for the transfer between theoretical material structure models and the interpretation of their experimental SS-NMR spectra. NMR-active nuclei within HPMs are routinely probed by their chemical shielding anisotropy (CSA). However, several nuclear isotopes of interest, e.g.7Li and 27Al, have a nuclear quadrupole and experience additional interactions with the surrounding electric field gradient (EFG). The quadrupolar interaction is a valuable source of information about atomistic structure, and in particular, local symmetry, complementing the CSA. As such, there is a range of different methods and codes to choose from for calculating EFGs, from all-electron to plane wave methods. We benchmark the accuracy of different simulation strategies for computing the EFG tensor of quadrupolar nuclei with plane wave density functional theory (DFT) and study the impact of the material structure as well as the details of the simulation strategy. Especially for small nuclei with few electrons, such as 7Li, we show that the choice of physical approximations and simulation parameters has a large effect on the transferability of the simulation results. To the best of our knowledge, we present the first comprehensive reference scale and literature survey for 7Li quadrupolar couplings. The results allow us to establish practical guidelines for developing the best simulation strategy for correlating DFT to experimental data extracting the maximum benefit and information from both, thereby advancing further research into HPMs.

FIRST-PRINCIPLES CALCULATION OF STRUCTURAL, ELECTRONIC, AND OPTICAL PROPERTIES OF InP1-XSbX USING WC-mBJ FOR NANOSCALE IR APPLICATIONS

The structural, electronic, and optical characteristics of cubic InP<sub>1-x</sub>Sb<sub>x</sub>(x = 0, 0.25, 0.50, 0.75, 1) ternary alloys were explored using the full-potential linearized augmented plane wave density functional theory approach. The total energy vs. volume optimization, lattice constants, and density of states were investigated for InP<sub>1-x</sub>Sb<sub>x</sub> alloys using exchange correlation function Wu-Cohen generalized gradient approximation (WC-GGA), available with the WIEN2k code. Band structure of the alloys was calculated using TB-mBJ functional to achieve a higher bandgap accuracy. The results of the mBJ experiment are in close agreement to those of the other experimental studies when compared to WC-GGA. Dielectric function and energy loss function were calculated in order to explore optical properties of the alloys. It was noticed that the estimated lattice parameters exhibit reduction when the Sb content is increased. Furthermore, the compositional dependency of the structural, electronic, and optical properties were also reported. For the InP<sub>1-x</sub>Sb<sub>x</sub> alloys, a band gap of less than 1.6 eVwere observed, making it suitable for usage in infrared optoelectronics devices.

CO2 ultrathin film growth on a monolayer of CO2 adsorbed on the NaCl(100) surface: sticking coefficient and IR-optical signatures in the ν3 region.

CO2 ultrathin molecular films were grown onto a preadsorbed monolayer NaCl(100)/p(2 × 1)-CO2 at 40 K. Polarization infrared spectroscopy (PIRS) reveals that so-prepared films have better quality than directly grown films. A sticking probability of 0.74 ± 0.1 was deduced from the integrated IR absorption. The presence of the monolayer doublet in the film spectra suggests a Stranski-Krastanov film growth with locally varying film thicknesses on the surface. In the region of the ν3(12C16O2) band, fine structure was observed between the well-known transverse-optical (TO) and longitudinal optical (LO) bands. Two independent computational models were applied to analyze the nature of the observed fine structure. Both pair potential calculations in combination with a vibrational exciton model as well as plane-wave density functional theory (DFT) in combination with phonon calculations of IR intensities at the Γ-point reveal that a weak mode visible in s-polarization and p-polarization originates from a vibrational film excitation located near the substrate interface. A series of p-polarized weak bands appearing and partly disappearing upon film-growth is assigned to film stacks of unique local thickness.

Enantiopure molecules form apparently racemic monolayers of chiral cyclic pentamers.

Ultra-high vacuum scanning tunneling microscopy (UHV-STM) was used to investigate two related molecules pulse-deposited onto Au(111) surfaces: indoline-2-carboxylic acid and proline (pyrrolidine-2-carboxylic acid). Indoline-2-carboxylic acid and proline form both dimers and C5-symmetric "pinwheel" pentamers. Enantiomerically pure S-(-)-indoline-2-carboxylic acid and S-proline were used, and the pentamer structures observed for both were chiral. However, the presence of apparently equal numbers of 'right-' and 'left-handed' pinwheels is contrary to the general understanding that the chirality of the molecule dictates supramolecular chirality. A variety of computational methods were used to elucidate pentamer geometry for S-proline. Straightforward geometry optimization proved difficult, as the size of the cluster and the number of possible intermolecular interactions produced an interaction potential with multiple local minima. Instead, the Amber force field was used to exhaustively search all of phase space for chemically reasonable pentamer structures, producing a limited number of candidate structures that were then optimized as gas-phase clusters using density functional theory (DFT). The binding energies of the two lowest-energy pentamers on the Au(111) surface were then calculated by plane-wave DFT using the VASP software, and STM images predicted. These calculations indicate that the right- and left-handed pentamers are instead two different polymorphs.