First-principles calculations of electron interactions in materials have seen rapid progress in recent years, with electron-phonon (e-ph) interactions being a prime example. However, these techniques use large matrices encoding the interactions on dense momentum grids, which reduces computational efficiency and obscures interpretability. For e-ph interactions, existing interpolation techniques leverage locality in real space, but the high dimensionality of the data remains a bottleneck to balance cost and accuracy. Here we show an efficient way to compress e-ph interactions based on singular value decomposition (SVD), a widely used matrix and image compression technique. Leveraging (un)constrained SVD methods, we accurately predict material properties related to e-ph interactions—including charge mobility, spin relaxation times, band renormalization, and superconducting critical temperature—while using only a small fraction (1%–2%) of the interaction data. These findings unveil the hidden low-dimensional nature of e-ph interactions. Furthermore, they accelerate state-of-the-art first-principles e-ph calculations by about 2 orders of magnitude without sacrificing accuracy. Our Pareto-optimal parametrization of e-ph interactions can be readily generalized to electron-electron and electron-defect interactions, as well as to other couplings, advancing quantitative studies of condensed matter.
more »
« less
Ab initio electron-defect interactions using Wannier functions
Abstract Computing electron–defect (e–d) interactions from first principles has remained impractical due to computational cost. Here we develop an interpolation scheme based on maximally localized Wannier functions (WFs) to efficiently computee–d interaction matrix elements. The interpolated matrix elements can accurately reproduce those computed directly without interpolation and the approach can significantly speed up calculations ofe–d relaxation times and defect-limited charge transport. We show example calculations of neutral vacancy defects in silicon and copper, for which we compute thee–d relaxation times on fine uniform and random Brillouin zone grids (and for copper, directly on the Fermi surface), as well as the defect-limited resistivity at low temperature. Our interpolation approach opens doors for atomistic calculations of charge carrier dynamics in the presence of defects.
more »
« less
- PAR ID:
- 10154073
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- npj Computational Materials
- Volume:
- 6
- Issue:
- 1
- ISSN:
- 2057-3960
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The interactions between solute atoms and crystalline defects such as vacancies, dislocations, and grain boundaries are essential in determining alloy properties. Here we present a general linear correlation between two descriptors of local electronic structures and the solute-defect interaction energies in binary alloys of body-centered-cubic (bcc) refractory metals (such as W and Ta) with transition-metal substitutional solutes. One electronic descriptor is the bimodality of thed-orbital local density of states for a matrix atom at the substitutional site, and the other is related to the hybridization strength between the valancesp-andd-bands for the same matrix atom. For a particular pair of solute-matrix elements, this linear correlation is valid independent of types of defects and the locations of substitutional sites. These results provide the possibility to apply local electronic descriptors for quantitative and efficient predictions on the solute-defect interactions and defect properties in alloys.more » « less
-
A<sc>bstract</sc> We study line defects in the fermionic CFTs in the Gross-Neveu-Yukawa universality class in dimensions 2< d <4. These CFTs may be described as the IR fixed points of the Gross-Neveu-Yukawa (GNY) model ind= 4 −ϵ, or as the UV fixed points of the Gross-Neveu (GN) model, which can be studied using the largeNexpansion in 2< d <4. These models admit natural line defects obtained by integrating over a line either the scalar field in the GNY description, or the fermion bilinear operator in the GN description. We compute the beta function for the defect RG flow using both the epsilon expansion and the largeNapproach, and find IR stable fixed points for the defect coupling, thus providing evidence for a non-trivial IR DCFT. We also compute some of the DCFT observables at the fixed point, and check that theg-function associated with the circular defect is consistent with theg-theorem for the defect RG flow.more » « less
-
A<sc>bstract</sc> We study a surface defect in the free and criticalO(N) vector models, defined by adding a quadratic perturbation localized on a two-dimensional subspace of thed-dimensional CFT. We compute the beta function for the corresponding defect renormalization group (RG) flow, and provide evidence that at long distances the system flows to a nontrivial defect conformal field theory (DCFT). We use epsilon and largeNexpansions to compute several physical quantities in the DCFT, finding agreement across different expansion methods. We also compute the defect free energy, and check consistency with the so-calledb-theorem for RG flows on surface defects.more » « less
-
Nonadiabatic molecular dynamics (NA-MD) is a powerful tool to model far-from-equilibrium processes, such as photochemical reactions and charge transport. NA-MD application to condensed phase has drawn tremendous attention recently for development of next-generation energy and optoelectronic materials. Studies of condensed matter allow one to employ efficient computational tools, such as density functional theory (DFT) and classical path approximation (CPA). Still, system size and simulation timescale are strongly limited by costly ab initio calculations of electronic energies, forces, and NA couplings. We resolve the limitations by developing a fully machine learning (ML) approach in which all the above properties are obtained using neural networks based on local descriptors. The ML models correlate the target properties for NA-MD, implemented with DFT and CPA, directly to the system structure. Trained on small systems, the neural networks are applied to large systems and long timescales, extending NA-MD capabilities by orders of magnitude. We demonstrate the approach with dependence of charge trapping and recombination on defect concentration in MoS2. Defects provide the main mechanism of charge losses, resulting in performance degradation. Charge trapping slows with decreasing defect concentration; however, recombination exhibits complex dependence, conditional on whether it occurs between free or trapped charges, and relative concentrations of carriers and defects. Delocalized shallow traps can become localized with increasing temperature, changing trapping and recombination behavior. Completely based on ML, the approach bridges the gap between theoretical models and realistic experimental conditions and enables NA-MD on thousand-atom systems and many nanoseconds.more » « less
An official website of the United States government
