BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Vikram Gavini

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

ResearcherUnclaimed source record
physics.comp-phcond-mat.mtrl-scimath-phmath.MPmath.NAphysics.chem-phquant-ph

Catalog footprint

What is connected

16works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

16 published item(s)

preprint2026arXiv

Intrinsic ductility enhancement in Mg alloys elucidated via large-scale ab-initio calculations

Magnesium is the lightest structural alloy, yet its practical use is limited by its low ductility. Recent studies suggest ductility enhancement in dilute Mg alloys may stem from favorable solute modification of <c+a> pyramidal I/II screw dislocation core energy difference, activating <c+a> slip via a double cross-slip mechanism. This work conducts large-scale DFT calculations, reaching ~6,000 atoms, of <c+a> dislocation energetics in Mg and Mg-Y/Zn alloys. We find that relative solute strengthening effects on pyramidal I and II screw dislocation glide are crucial for cross-slip enhancement in Mg-Y, in contrast to prior investigations, that find solute-mediated dislocation-core energy modification as the main driver. Our predictions align with single- and poly-crystal experimental results and also capture the transition from pyramidal II to I preferred slip in Mg-Y.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Exact and model exchange-correlation potentials for open-shell systems

The conventional approaches to the inverse density functional theory problem typically assume non-degeneracy of the Kohn-Sham (KS) eigenvalues, greatly hindering their use in open-shell systems. We present a generalization of the inverse density functional theory problem that can seamlessly admit degenerate KS eigenvalues. Additionally, we allow for fractional occupancy of the Kohn-Sham orbitals to also handle non-interacting ensemble-v-representable densities, as opposed to just non-interacting pure-v-representable densities. We present the exact exchange-correlation (XC) potentials for six open-shell systems -- four atoms (Li, C, N, and O) and two molecules (CN and $\text{CH}_2$) -- using accurate ground-state densities from configuration interaction calculations. We compare these exact XC potentials with model XC potentials obtained using non-local (B3LYP, SCAN0) and local/semi-local (SCAN, PBE, PW92) XC functionals. Although the relative errors in the densities obtained from these DFT functionals are of $\mathcal{O}(10^{-3}-10^{-2})$, the relative errors in the model XC potentials remain substantially large -- $\mathcal{O}(10^{-1}-10^0)$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization

We present DFT-FE 1.0, building on DFT-FE 0.6 [Comput. Phys. Commun. 246, 106853 (2020)], to conduct fast and accurate large-scale density functional theory (DFT) calculations (reaching ~ $100,000$ electrons) on both many-core CPU and hybrid CPU-GPU computing architectures. This work involves improvements in the real-space formulation -- via an improved treatment of the electrostatic interactions that substantially enhances the computational efficiency -- as well high-performance computing aspects, including the GPU acceleration of all the key compute kernels in DFT-FE. We demonstrate the accuracy by comparing the ground-state energies, ionic forces and cell stresses on a wide-range of benchmark systems against those obtained from widely used DFT codes. Further, we demonstrate the numerical efficiency of our implementation, which yields $\sim 20 \times$ CPU-GPU speed-up by using GPU acceleration on hybrid CPU-GPU nodes. Notably, owing to the parallel-scaling of the GPU implementation, we obtain wall-times of $80-140$ seconds for full ground-state calculations, with stringent accuracy, on benchmark systems containing ~ $6,000-15,000$ electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Ionic forces and stress tensor in all-electron DFT calculations using enriched finite element basis

The enriched finite element basis -- wherein the finite element basis is enriched with atom-centered numerical functions -- has recently been shown to be a computationally efficient basis for systematically convergent all-electron DFT ground-state calculations. In this work, we present the expressions to compute variationally consistent ionic forces and stress tensor for all-electron DFT calculations in the enriched finite element basis. In particular, we extend the formulation of configurational forces in Motamarri & Gavini (Phys. Rev. B 2018) to the enriched finite element basis and elucidate the additional contributions arising from the enrichment functions. We demonstrate the accuracy of the formulation by comparing the computed forces and stresses for various benchmark systems with those obtained from finite-differencing the ground-state energy. Further, we also benchmark our calculations against Gaussian basis for molecular systems and against the LAPW+lo basis for periodic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Spin-spin interactions in solids from mixed all-electron and pseudopotential calculations $-$ a path to screening materials for spin qubits

Understanding the quantum dynamics of spin defects and their coherence properties requires accurate modeling of spin-spin interaction in solids and molecules, for example by using spin Hamiltonians with parameters obtained from first-principles calculations. We present a real-space approach based on density functional theory for the calculation of spin-Hamiltonian parameters, where only selected atoms are treated at the all-electron level, while the rest of the system is described with the pseudopotential approximation. Our approach permits calculations for systems containing more than 1000 atoms, as demonstrated for defects in diamond and silicon carbide. We show that only a small number of atoms surrounding the defect needs to be treated at the all-electron level, in order to obtain an overall all-electron accuracy for hyperfine and zero-field splitting tensors. We also present results for coherence times, computed with the cluster correlation expansion method, highlighting the importance of accurate spin-Hamiltonian parameters for quantitative predictions of spin dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Tensor-structured algorithm for reduced-order scaling large-scale Kohn-Sham density functional theory calculations

We present a tensor-structured algorithm for efficient large-scale DFT calculations by constructing a Tucker tensor basis that is adapted to the Kohn-Sham Hamiltonian and localized in real-space. The proposed approach uses an additive separable approximation to the Kohn-Sham Hamiltonian and an $L_1$ localization technique to generate the 1-D localized functions that constitute the Tucker tensor basis. Numerical results show that the resulting Tucker tensor basis exhibits exponential convergence in the ground-state energy with increasing Tucker rank. Further, the proposed tensor-structured algorithm demonstrated sub-quadratic scaling with system size for both systems with and without a gap, and involving many thousands of atoms. This reduced-order scaling has also resulted in the proposed approach outperforming plane-wave DFT implementation for systems beyond 2,000 electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Treecode-accelerated Green Iteration for Kohn-Sham Density Functional Theory

We present a real-space computational method called treecode-accelerated Green Iteration (TAGI) for all-electron Kohn-Sham Density Functional Theory. TAGI is based on a reformulation of the Kohn-Sham equations in which the eigenvalue problem in differential form is converted into a fixed-point problem in integral form by convolution with the modified Helmholtz Green's function. In each self-consistent field (SCF) iteration, the fixed-points are computed by Green Iteration, where the discrete convolution sums are efficiently evaluated by a GPU-accelerated barycentric Lagrange treecode. Other techniques used in TAGI include adaptive mesh refinement, Fejér quadrature, singularity subtraction, gradient-free eigenvalue update, and Anderson mixing to accelerate convergence of the SCF and Green Iterations. Ground state energy computations of several atoms (Li, Be, O) and small molecules (H$_2$, CO, C$_6$H$_6$) demonstrate TAGI's ability to efficiently achieve chemical accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

DFT-FE -- A massively parallel adaptive finite-element code for large-scale density functional theory calculations

We present an accurate, efficient and massively parallel finite-element code, DFT-FE, for large-scale ab-initio calculations (reaching $\sim 100,000$ electrons) using Kohn-Sham density functional theory (DFT). DFT-FE is based on a local real-space variational formulation of the Kohn-Sham DFT energy functional that is discretized using a higher-order adaptive spectral finite-element (FE) basis, and treats pseudopotential and all-electron calculations in the same framework, while accommodating non-periodic, semi-periodic and periodic boundary conditions. We discuss the main aspects of the code, which include, the various strategies of adaptive FE basis generation, and the different approaches employed in the numerical implementation of the solution of the discrete Kohn-Sham problem that are focused on significantly reducing the floating point operations, communication costs and latency. We demonstrate the accuracy of DFT-FE by comparing the energies, ionic forces and periodic cell stresses on a wide range of problems with popularly used DFT codes. Further, we demonstrate that DFT-FE significantly outperforms widely used plane-wave codes---both in CPU-times and wall-times, and on both non-periodic and periodic systems---at systems sizes beyond a few thousand electrons, with over $5-10$ fold speedups in systems with more than 10,000 electrons. The benchmark studies also highlight the excellent parallel scalability of DFT-FE, with strong scaling demonstrated on up to 192,000 MPI tasks.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A subquadratic-scaling subspace projection method for large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization

We present a subspace projection technique to conduct large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization. The proposed method treats both metallic and insulating materials in a single framework, and is applicable to both pseudopotential as well as all-electron calculations. The key ideas involved in the method include: (i) employing a higher-order spectral finite-element basis that is amenable to mesh adaption; (ii) using a Chebyshev filter to construct a subspace which is an approximation to the occupied eigenspace in a given self-consistent field iteration; (iii) using a localization procedure to construct a non-orthogonal localized basis spanning the Chebyshev filtered subspace; (iv) using a Fermi-operator expansion in terms of the subspace-projected Hamiltonian represented in the non-orthogonal localized basis to compute relevant quantities like the density matrix, electron density and band energy. We demonstrate the accuracy and efficiency of the approach on benchmark systems involving pseudopotential calculations on metallic aluminum nano-clusters up to 3430 atoms and on insulating alkane chains up to 7052 atoms, as well as all-electron calculations on silicon nano-clusters up to 3920 electrons. The benchmark studies revealed that accuracies commensurate with chemical accuracy can be obtained, and a subquadratic-scaling with system size was observed for the range of materials systems studied. In particular, for the alkane chains---close to linear-scaling is observed, whereas, for aluminum nano-clusters---the scaling is observed to be $\mathcal{O} (N^{1.46})$. For all-electron calculations on silicon nano-clusters, the scaling with the number of electrons is computed to be $\mathcal{O} (N^{1.75})$. Furthermore, significant computational savings have been realized with the proposed approach with respect to reference calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Orbital-free DFT study of the energetics of vacancy clustering and prismatic dislocation loop nucleation in aluminum

In the present work, we conduct large-scale orbital-free DFT calculations to study the energetics of vacancy clustering in aluminum from electronic structure calculations. The simulation domains considered in this study are as large as those containing a million atoms to accurately account for both the electronic structure and long-ranged elastic fields. Our results indicate that vacancy clustering is an energetically favorable mechanisms with positive binding energies for a range of vacancy clusters considered in the present study. In particular, the $19$ vacancy hexagonal cluster lying in $\{111\}$ plane has a very large binding energy with the relaxed atomic structure representative of a prismatic dislocation loop. This suggests that vacancy prismatic loops as small as those formed from 19 vacancies are stable, thus providing insights into the nucleation sizes of these defects in aluminum.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Real-space formulation of orbital-free density functional theory using finite-element discretization: The case for Al, Mg, and Al-Mg intermetallics

We propose a local real-space formulation for orbital-free DFT with density dependent kinetic energy functionals and a unified variational framework for computing the configurational forces associated with geometry optimization of both internal atomic positions as well as the cell geometry. The proposed real-space formulation, which involves a reformulation of the extended interactions in electrostatic and kinetic energy functionals as local variational problems in auxiliary potential fields, also readily extends to all-electron orbital-free DFT calculations that are employed in warm dense matter calculations. We use the local real-space formulation in conjunction with higher-order finite-element discretization to demonstrate the accuracy of orbital-free DFT and the proposed formalism for the Al-Mg materials system, where we obtain good agreement with Kohn-Sham DFT calculations on a wide range of properties and benchmark calculations. Finally, we investigate the cell-size effects in the electronic structure of point defects, in particular a mono-vacancy in Al. We unambiguously demonstrate that the cell-size effects observed from vacancy formation energies computed using periodic boundary conditions underestimate the extent of the electronic structure perturbations created by the defect. On the contrary, the bulk Dirichlet boundary conditions, accessible only through the proposed real-space formulation, which correspond to an isolated defect embedded in the bulk, show cell-size effects in the defect formation energy that are commensurate with the perturbations in the electronic structure. Our studies suggest that even for a simple defect like a vacancy in Al, we require cell-sizes of $\sim 10^3$ atoms for convergence in the electronic structure.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Electronic-structure study of an edge dislocation in Aluminum and the role of macroscopic deformations on its energetics

We employed a real-space formulation of orbital-free density functional theory using finite-element basis to study the defect-core and energetics of an edge dislocation in Aluminum. Our study shows that the core-size of a perfect edge dislocation is around ten times the magnitude of the Burgers vector. This finding is contrary to the widely accepted notion that continuum descriptions of dislocation energetics are accurate beyond 1-3 Burgers vector from the dislocation line. Consistent with prior electronic-structure studies, we find that the perfect edge dislocation dissociates into two Shockley partials with a partial separation distance of 12.8 $Å$. Interestingly, our study revealed a significant influence of macroscopic deformations on the core-energy of Shockley partials. We show that this dependence of the core-energy on macroscopic deformations results in an additional force on dislocations, beyond the Peach-Koehler force, that is proportional to strain gradients. Further, we demonstrate that this force from core-effects can be significant and can play an important role in governing the dislocation behavior in regions of inhomogeneous deformations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Higher-order adaptive finite-element methods for Kohn-Sham density functional theory

We present an efficient computational approach to perform real-space electronic structure calculations using an adaptive higher-order finite-element discretization of Kohn-Sham density-functional theory (DFT). To this end, we develop an a-priori mesh adaption technique to construct a close to optimal finite-element discretization of the problem. We further propose an efficient solution strategy for solving the discrete eigenvalue problem by using spectral finite-elements in conjunction with Gauss-Lobatto quadrature, and a Chebyshev acceleration technique for computing the occupied eigenspace. The proposed approach has been observed to provide a staggering 100-200 fold computational advantage over the solution of a generalized eigenvalue problem. Using the proposed solution procedure, we investigate the computational efficiency afforded by higher-order finite-element discretization of the Kohn-Sham DFT problem. Our studies suggest that staggering computational savings of the order of 1000 fold relative to linear finite-elements can be realized, for both all-electron and local pseudopotential calculations. On all the benchmark systems studied, we observe diminishing returns in computational savings beyond the sixth-order for accuracies commensurate with chemical accuracy. A comparative study of the computational efficiency of the proposed higher-order finite-element discretizations suggests that the performance of finite-element basis is competing with the plane-wave discretization for non-periodic local pseudopotential calculations, and compares to the Gaussian basis for all-electron calculations within an order of magnitude. Further, we demonstrate the capability of the proposed approach to compute the electronic structure of a metallic system containing 1688 atoms using modest computational resources, and good scalability of the present implementation up to 192 processors.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Higher-order adaptive finite-element methods for orbital-free density functional theory

In the present work, we investigate the computational efficiency afforded by higher-order finite-element discretization of the saddle-point formulation of orbital-free density functional theory. We first investigate the robustness of viable solution schemes by analyzing the solvability conditions of the discrete problem. We find that a staggered solution procedure where the potential fields are computed consistently for every trial electron-density is a robust solution procedure for higher-order finite-element discretizations. We next study the numerical convergence rates for various orders of finite-element approximations on benchmark problems. We obtain close to optimal convergence rates in our studies, although orbital-free density-functional theory is nonlinear in nature and some benchmark problems have Coulomb singular potential fields. We finally investigate the computational efficiency of various higher-order finite-element discretizations by measuring the CPU time for the solution of discrete equations on benchmark problems that include large Aluminum clusters. In these studies, we use mesh coarse-graining rates that are derived from error estimates and an a priori knowledge of the asymptotic solution of the far-field electronic fields. Our studies reveal a significant 100-1000 fold computational savings afforded by the use of higher-order finite-element discretization, alongside providing the desired chemical accuracy. We consider this study as a step towards developing a robust and computationally efficient discretization of electronic structure calculations using the finite-element basis.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

An analysis of the field theoretic approach to the quasi-continuum method

Using the orbital-free density functional theory as a model theory, we present an analysis of the field theoretic approach to quasi-continuum method. In particular, by perturbation method and multiple scale analysis, we provide a formal justification for the validity of the coarse-graining of various fields, which is central to the quasi-continuum reduction of field theories. Further, we derive the homogenized equations that govern the behavior of electronic fields in regions of smooth deformations. Using Fourier analysis, we determine the far-field solutions for these fields in the presence of local defects, and subsequently estimate cell-size effects in computed defect energies.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Role of the defect-core in energetics of vacancies

Electronic structure calculations at macroscopic scales are employed to investigate the crucial role of a defect-core in the energetics of vacancies in aluminum. We find that vacancy core-energy is significantly influenced by the state of deformation at the vacancy-core, especially volumetric strains. Insights from the core electronic structure and computed displacement fields show that this dependence on volumetric strains is closely related to the changing nature of the core-structure under volumetric deformations. These results are in sharp contrast to mechanics descriptions based on elastic interactions that often consider defect core-energies as an inconsequential constant. Calculations suggest that the variation in core-energies with changing macroscopic deformations is quantitatively more significant than the corresponding variation in relaxation energies associated with elastic fields. Upon studying the influence of various macroscopic deformations, which include volumetric, uniaxial, biaxial and shear deformations, on the formation energies of vacancies, we show that volumetric deformations play a dominant role in governing the energetics of these defects. Further, by plotting formation energies of vacancies and di-vacancies against the volumetric strain corresponding to any macroscopic deformation, we find that all variations in the formation energies collapse on to a universal curve. This suggests a universal role of volumetric strains in the energetics of vacancies. Implications of these results in the context of dynamic failure in metals due to spalling are analyzed.

0 citations0 reviewsNo code linkedOpen access