Showing 1–17 of 17 results for author: Gavini, V

Quasicrystal bulk and surface energies from density functional theory

Authors: Woohyeon Baek, Sambit Das, Shibo Tan, Vikram Gavini, Wenhao Sun

Abstract: Are quasicrystals stable or metastable? Density functional theory (DFT) is often used to evaluate thermodynamic stability, but quasicrystals are long-range aperiodic and their energies cannot be calculated using conventional ab initio methods. Here, we perform first-principles calculations on quasicrystal nanoparticles of increasing sizes, from which we can directly extrapolate their bulk and surf… ▽ More Are quasicrystals stable or metastable? Density functional theory (DFT) is often used to evaluate thermodynamic stability, but quasicrystals are long-range aperiodic and their energies cannot be calculated using conventional ab initio methods. Here, we perform first-principles calculations on quasicrystal nanoparticles of increasing sizes, from which we can directly extrapolate their bulk and surface energies. Using this technique, we determine with high confidence that the icosahedral quasicrystals ScZn7.33 and YbCd5.7 are ground-state phases--revealing that translational symmetry is not a necessary condition for the T = 0 K stability of inorganic solids. Although we find the ScZn7.33 quasicrystal to be thermodynamically stable, we show on a mixed thermodynamic and kinetic phase diagram that its solidification from the melt is nucleation-limited, which illustrates why even stable materials may be kinetically challenging to grow. Our techniques here broadly open the door to first-principles investigations into the structure-bonding-stability relationships of aperiodic materials. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Tucker tensor approach for accelerating exchange computations in a real-space finite-element discretization of generalized Kohn-Sham density functional theory

Authors: Vishal Subramanian, Sambit Das, Vikram Gavini

Abstract: The evaluation of Fock exchange is often the computationally most expensive part of hybrid functional density functional theory calculations in a systematically improvable, complete basis. In this work, we employ a Tucker tensor based approach that substantially accelerates the evaluation of the action of Fock exchange by transforming 3-dimensional convolutional integrals into a tensor product of… ▽ More The evaluation of Fock exchange is often the computationally most expensive part of hybrid functional density functional theory calculations in a systematically improvable, complete basis. In this work, we employ a Tucker tensor based approach that substantially accelerates the evaluation of the action of Fock exchange by transforming 3-dimensional convolutional integrals into a tensor product of 1-dimensional convolution integrals. Our numerical implementation uses a parallelization strategy that balances the memory and communication bottlenecks, alongside overalapping compute and communication operations to enhance computational efficiency and parallel scalability. The accuracy and computational efficiency is demonstrated on various systems, including Pt clusters of various sizes and a $\text{TiO}_{\text{2}}$ cluster with 3,684 electrons. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: 30 pages, 10 figures

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

Authors: Bikash Kanungo, Jeffrey Hatch, Paul M. Zimmerman, Vikram Gavini

Abstract: 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 t… ▽ More 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)$. △ Less

Submitted 5 November, 2023; v1 submitted 24 May, 2023; originally announced May 2023.

Unconventional Error Cancellation Explains the Success of Hartree-Fock Density Functional Theory for Barrier Heights

Authors: Bikash Kanungo, Aaron D. Kaplan, Chandra Shahi, Vikram Gavini, John P. Perdew

Abstract: Energy barriers, which control the rates of chemical reactions, are seriously underestimated by computationally-efficient semi-local approximations for the exchange-correlation energy. The accuracy of a semi-local density functional approximation is strongly boosted for reaction barrier heights by evaluating that approximation non-self-consistently on Hartree-Fock electron densities, as known for… ▽ More Energy barriers, which control the rates of chemical reactions, are seriously underestimated by computationally-efficient semi-local approximations for the exchange-correlation energy. The accuracy of a semi-local density functional approximation is strongly boosted for reaction barrier heights by evaluating that approximation non-self-consistently on Hartree-Fock electron densities, as known for about 30 years. The conventional explanation is that Hartree-Fock theory yields the more accurate density. This article presents a benchmark Kohn-Sham inversion of accurate coupled-cluster densities for the reaction H$_2$ + F $\rightarrow$ HHF $\rightarrow$ H + HF, and finds a strong, understandable cancellation between positive (excessively over-corrected) density-driven and large negative functional-driven errors (expected from stretched radical bonds in the transition state) within this Hartree-Fock density functional theory. This confirms earlier conclusions [Kaplan et al., J. Chem. Theory Comput. 19, 532--543 (2023)] based on 76 barrier heights and three less reliable, but less expensive, fully-nonlocal density-functional proxies for the exact density. △ Less

Submitted 2 November, 2023; v1 submitted 9 March, 2023; originally announced March 2023.

Comments: Revisions in response to peer comments. Addition of Supplemental Material. Addition of BCCD analysis, further KS inversion analysis

Accelerating self-consistent field iterations in Kohn-Sham density functional theory using a low rank approximation of the dielectric matrix

Authors: Sambit Das, Vikram Gavini

Abstract: We present an efficient preconditioning technique for accelerating the fixed point iteration in real-space Kohn-Sham density functional theory (DFT) calculations. The preconditioner uses a low rank approximation of the dielectric matrix (LRDM) based on Gâteaux derivatives of the residual of fixed point iteration along appropriately chosen direction functions. We develop a computationally efficient… ▽ More We present an efficient preconditioning technique for accelerating the fixed point iteration in real-space Kohn-Sham density functional theory (DFT) calculations. The preconditioner uses a low rank approximation of the dielectric matrix (LRDM) based on Gâteaux derivatives of the residual of fixed point iteration along appropriately chosen direction functions. We develop a computationally efficient method to evaluate these Gâteaux derivatives in conjunction with the Chebyshev filtered subspace iteration procedure, an approach widely used in large-scale Kohn-Sham DFT calculations. Further, we propose a variant of LRDM preconditioner based on adaptive accumulation of low-rank approximations from previous SCF iterations, and also extend the LRDM preconditioner to spin-polarized Kohn-Sham DFT calculations. We demonstrate the robustness and efficiency of the LRDM preconditioner against other widely used preconditioners on a range of benchmark systems with sizes ranging from $\sim$ 100-1100 atoms ($\sim$ 500--20,000 electrons). The benchmark systems include various combinations of metal-insulating-semiconducting heterogeneous material systems, nanoparticles with localized $d$ orbitals near the Fermi energy, nanofilm with metal dopants, and magnetic systems. In all benchmark systems, the LRDM preconditioner converges robustly within 20--30 iterations. In contrast, other widely used preconditioners show slow convergence in many cases, as well as divergence of the fixed point iteration in some cases. Finally, we demonstrate the computational efficiency afforded by the LRDM method, with up to 3.4$\times$ reduction in computational cost for the total ground-state calculation compared to other preconditioners. △ Less

Submitted 4 March, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted in Physical Review B

Roadmap on Electronic Structure Codes in the Exascale Era

Authors: Vikram Gavini, Stefano Baroni, Volker Blum, David R. Bowler, Alexander Buccheri, James R. Chelikowsky, Sambit Das, William Dawson, Pietro Delugas, Mehmet Dogan, Claudia Draxl, Giulia Galli, Luigi Genovese, Paolo Giannozzi, Matteo Giantomassi, Xavier Gonze, Marco Govoni, Andris Gulans, François Gygi, John M. Herbert, Sebastian Kokott, Thomas D. Kühne, Kai-Hsin Liou, Tsuyoshi Miyazaki, Phani Motamarri , et al. (16 additional authors not shown)

Abstract: Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources d… ▽ More Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing. △ Less

Submitted 26 September, 2022; originally announced September 2022.

Comments: Submitted as a roadmap article to Modelling and Simulation in Materials Science and Engineering; Address any correspondence to Vikram Gavini (vikramg@umich.edu) and Danny Perez (danny_perez@lanl.gov)

Accurate approximations of density functional theory for large systems with applications to defects in crystalline solids

Authors: Kaushik Bhattacharya, Vikram Gavini, Michael Ortiz, Mauricio Ponga, Phanish Suryanarayana

Abstract: This chapter presents controlled approximations of Kohn-Sham density functional theory (DFT) that enable very large scale simulations. The work is motivated by the study of defects in crystalline solids, though the ideas can be used in other applications. The key idea is to formulate DFT as a minimization problem over the density operator, and to cast spatial and spectral discretization as systema… ▽ More This chapter presents controlled approximations of Kohn-Sham density functional theory (DFT) that enable very large scale simulations. The work is motivated by the study of defects in crystalline solids, though the ideas can be used in other applications. The key idea is to formulate DFT as a minimization problem over the density operator, and to cast spatial and spectral discretization as systematically convergent approximations. This enables efficient and adaptive algorithms that solve the equations of DFT with no additional modeling, and up to desired accuracy, for very large systems, with linear and sublinear scaling. Various approaches based on such approximations are presented, and their numerical performance demonstrated through selected examples. These examples also provide important insight about the mechanics and physics of defects in crystalline solids. △ Less

Submitted 11 December, 2021; originally announced December 2021.

Li$_x$CoO$_2$ phase stability studied by machine learning-enabled scale bridging between electronic structure, statistical mechanics and phase field theories

Authors: Gregory H. Teichert, Sambit Das, Muratahan Aykol, Chirranjeevi Gopal, Vikram Gavini, Krishna Garikipati

Abstract: Li$_xTM$O$_2$ (TM={Ni, Co, Mn}) are promising cathodes for Li-ion batteries, whose electrochemical cycling performance is strongly governed by crystal structure and phase stability as a function of Li content at the atomistic scale. Here, we use Li$_x$CoO$_2$ (LCO) as a model system to benchmark a scale-bridging framework that combines density functional theory (DFT) calculations at the atomistic… ▽ More Li$_xTM$O$_2$ (TM={Ni, Co, Mn}) are promising cathodes for Li-ion batteries, whose electrochemical cycling performance is strongly governed by crystal structure and phase stability as a function of Li content at the atomistic scale. Here, we use Li$_x$CoO$_2$ (LCO) as a model system to benchmark a scale-bridging framework that combines density functional theory (DFT) calculations at the atomistic scale with phase field modeling at the continuum scale to understand the impact of phase stability on microstructure evolution. This scale bridging is accomplished by incorporating traditional statistical mechanics methods with integrable deep neural networks, which allows formation energies for specific atomic configurations to be coarse-grained and incorporated in a neural network description of the free energy of the material. The resulting realistic free energy functions enable atomistically informed phase-field simulations. These computational results allow us to make connections to experimental work on LCO cathode degradation as a function of temperature, morphology and particle size. △ Less

Submitted 22 April, 2021; v1 submitted 16 April, 2021; originally announced April 2021.

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

Authors: Krishnendu Ghosh, He Ma, Mykyta Onizhuk, Vikram Gavini, Giulia Galli

Abstract: 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 a… ▽ More 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. △ Less

Submitted 30 January, 2021; originally announced February 2021.

All-electron density functional calculations for electron and nuclear spin interactions in molecules and solids

Authors: Krishnendu Ghosh, He Ma, Vikram Gavini, Giulia Galli

Abstract: The interaction between electronic and nuclear spins in the presence of external magnetic fields can be described by a spin Hamiltonian, with parameters obtained from first principles, electronic structure calculations. We describe an approach to compute these parameters, applicable to both molecules and solids, which is based on Density Functional Theory (DFT) and real-space, all-electron calcula… ▽ More The interaction between electronic and nuclear spins in the presence of external magnetic fields can be described by a spin Hamiltonian, with parameters obtained from first principles, electronic structure calculations. We describe an approach to compute these parameters, applicable to both molecules and solids, which is based on Density Functional Theory (DFT) and real-space, all-electron calculations using finite elements (FE). We report results for hyperfine tensors, zero field splitting tensors (spin-spin component) and nuclear quadrupole tensors of a series of molecules and of the nitrogen-vacancy center in diamond. We compare our results with those of calculations using Gaussian orbitals and plane-wave basis sets, and we discuss their numerical accuracy. We show that calculations based on FE can be systematically converged with respect to the basis set, thus allowing one to establish reference values for the spin Hamiltonian parameters, at a given level of DFT. △ Less

Submitted 17 March, 2019; v1 submitted 19 February, 2019; originally announced February 2019.

Journal ref: Phys. Rev. Materials 3, 043801 (2019)

Origins of the transformability of Nickel-Titanium shape memory alloys

Authors: Xian Chen, Colin Ophus, Chengyu Song, Jim Ciston, Sambit Das, Yintao Song, Yuriy Chumlyakov, Andrew Minor, Vikram Gavini, Richard D. James

Abstract: The near equiatomic NiTi alloy is the most successful shape memory alloy by a large margin. It is widely and increasingly used in biomedical devices. Yet, despite having a repeatable superelastic effect and excellent shape-memory, NiTi is very far from satisfying the conditions that characterize the most reversible phase transforming materials. Thus, the scientific reasons underlying its vast succ… ▽ More The near equiatomic NiTi alloy is the most successful shape memory alloy by a large margin. It is widely and increasingly used in biomedical devices. Yet, despite having a repeatable superelastic effect and excellent shape-memory, NiTi is very far from satisfying the conditions that characterize the most reversible phase transforming materials. Thus, the scientific reasons underlying its vast success present an enigma. In this work, we perform rigorous mathematical derivation and accurate DFT calculation of transformation mechanisms to seek previously unrecognized twin-like defects that we term involution domains, and we observe them in real space in NiTi by the aberration-corrected scanning transmission electron microscopy. Involution domains lead to an additional 216 compatible interfaces between phases in NiTi, and we theorize that this feature contributes importantly to its reliability. They are expected to arise in other transformations and to alter the conventional interpretation of the mechanism of the martensitic transformation. △ Less

Submitted 9 October, 2020; v1 submitted 18 January, 2019; originally announced January 2019.

Comments: 16 pages, 4 figures

Journal ref: Phys. Rev. Mater. 2020

Electronic structure study of screw dislocation core energetics in Aluminum and core energetics informed forces in a dislocation aggregate

Authors: Sambit Das, Vikram Gavini

Abstract: We use a real-space formulation of orbital-free DFT to study the core energetics and core structure of an isolated screw dislocation in Aluminum. Using a direct energetics based approach, we estimate the core size of a perfect screw dislocation to be $\approx$ 7 $|{\bf b}|$, which is considerably larger than previous estimates of $1-3~|{\bf b}|$ based on displacement fields. The perfect screw upon… ▽ More We use a real-space formulation of orbital-free DFT to study the core energetics and core structure of an isolated screw dislocation in Aluminum. Using a direct energetics based approach, we estimate the core size of a perfect screw dislocation to be $\approx$ 7 $|{\bf b}|$, which is considerably larger than previous estimates of $1-3~|{\bf b}|$ based on displacement fields. The perfect screw upon structural relaxation dissociates into two Shockley partials with partial separation distances of 8.2 Å~and 6.6 Å~ measured from the screw and edge component differential displacement plots, respectively. Similar to a previous electronic structure study on edge dislocation, we find that the core energy of the relaxed screw dislocation is strongly dependent on macroscopic deformations. Next, we use the edge and screw dislocation core energetics data with physically reasonable assumptions to develop a continuum energetics model for an aggregate of dislocations. Further, we use this model in a discrete dislocation network, and from the variations of the core energy with respect to the nodal positions of the network, we obtain the nodal core force which can directly be incorporated into discrete dislocation dynamics frameworks. We analyze and classify the nodal core force into three different contributions based on their decay behavior. Two of these contributions to the core force, both arising from the core energy dependence on macroscopic deformations, are not accounted for in currently used discrete dislocation dynamics models which assume the core energy to be a constant. Using case studies involving simple dislocation structures, we demonstrate that the contribution to the core force from the core energy dependence on macroscopic deformations can be significant in comparison to the elastic Peach-Koehler force even up to distances of $10-15$ nm between dislocation structures. △ Less

Submitted 16 March, 2017; v1 submitted 30 January, 2017; originally announced January 2017.

Comments: In press- Journal of the Mechanics and Physics of Solids

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

Authors: Balachandran Radhakrishnan, Vikram Gavini

Abstract: 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 clusterin… ▽ More 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. △ Less

Submitted 8 June, 2015; originally announced June 2015.

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

Authors: Sambit Das, Mrinal Iyer, Vikram Gavini

Abstract: 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 electros… ▽ More 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. △ Less

Submitted 12 June, 2015; v1 submitted 23 April, 2015; originally announced April 2015.

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

Authors: Mrinal Iyer, Balachandran Radhakrishnan, Vikram Gavini

Abstract: 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 dislocatio… ▽ More 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. △ Less

Submitted 23 December, 2014; v1 submitted 14 June, 2014; originally announced June 2014.

Comments: 27 pages, 9 figures

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

Authors: Vikram Gavini, Liping Liu

Abstract: 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 th… ▽ More 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. △ Less

Submitted 10 December, 2010; originally announced December 2010.

Comments: 26 pages, 1 figure

Role of the defect-core in energetics of vacancies

Authors: Vikram Gavini

Abstract: 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 dep… ▽ More 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. △ Less

Submitted 4 August, 2009; originally announced August 2009.