Ab Initio Total Energy Calculations

AB INITIO – literally “from first principles”; a computational approach t… #

Related terms: first‑principles, quantum‑mechanical modeling. In the context of total‑energy calculations, ab initio methods such as Density Functional Theory (DFT) provide the ground‑state energy of a material solely from its atomic composition. Example: Calculating the cohesive energy of silicon using only the atomic number of Si and the crystal lattice. Practical application: Predicting phase stability of alloys without experimental data. Challenge: Computational cost grows rapidly with system size and with the need for high‑accuracy exchange‑correlation functionals.

BRILLOUIN ZONE – the primitive cell in reciprocal space defined by the se… #

Related terms: reciprocal lattice, k‑point sampling. For total‑energy calculations, integration over the Brillouin zone determines electronic properties such as total energy, density of states, and forces. Example: Using a 6×6×6 Monkhorst‑Pack grid to sample the Brillouin zone of a cubic perovskite. Practical application: Accurate band‑structure prediction for semiconductors. Challenge: Converging the energy with respect to k‑point density can be demanding for metallic systems where the Fermi surface is complex.

BROYDEN MIXING – an iterative scheme used to accelerate convergence of th… #

Related terms: SCF convergence, Pulay mixing. In Quantum Espresso and VASP, Broydens method adjusts the mixing coefficients based on previous iterations, reducing the number of SCF steps needed. Example: Setting “mixing_beta = 0.4” and “mixing_ndim = 8” in a VASP INCAR file to invoke Broyden mixing. Practical application: Enabling reliable convergence for large supercells with many atoms. Challenge: Choosing optimal mixing parameters; inappropriate values may cause divergence or oscillatory behavior.

CAR‑PARINELLO MOLECULAR DYNAMICS (CPMD) – a technique that couples the el… #

Related terms: Born‑Oppenheimer MD, ab initio MD. Although not a primary total‑energy method, CPMD generates forces from the same DFT Hamiltonian used in static calculations, providing insight into finite‑temperature phenomena. Example: Simulating proton transfer in water clusters using the CPMD implementation in Quantum Espresso. Practical application: Studying reaction pathways and vibrational spectra at the quantum level. Challenge: The need for a small fictitious electron mass and tiny time steps makes CPMD computationally expensive.

CONVERGENCE TESTS – systematic checks that ensure calculated quantities (… #

Related terms: cutoff convergence, k‑mesh convergence. Performing convergence tests is a prerequisite for reliable total‑energy results; otherwise, reported energies may be artifacts of insufficient numerical precision. Example: Increasing the kinetic‑energy cutoff from 40 Ry to 80 Ry until the total energy changes by less than 1 meV per atom. Practical application: Establishing a reproducible workflow for high‑throughput materials screening. Challenge: Balancing accuracy against computational cost, especially for large supercells or low‑symmetry structures.

DENSITY FUNCTIONAL THEORY (DFT) – the foundational quantum‑mechanical fra… #

Related terms: Hohenberg‑Kohn theorem, Kohn‑Sham equations. DFT is the workhorse for ab initio total‑energy calculations in both Quantum Espresso and VASP. Example: Solving the Kohn‑Sham equations for a NaCl crystal to obtain its equilibrium lattice constant. Practical application: Predicting formation energies, elastic constants, and defect formation energies across a wide range of materials. Challenge: The exact exchange‑correlation functional is unknown; approximate functionals introduce systematic errors such as band‑gap underestimation.

DENSITY OF STATES (DOS) – a spectral function that counts the number of e… #

Related terms: projected DOS, partial DOS. The DOS provides insight into metallic, semiconducting, or insulating behavior and is essential for interpreting photoemission experiments. Example: Plotting the total DOS of graphene and observing the characteristic linear increase near the Dirac point. Practical application: Assessing the suitability of a material for thermoelectric applications based on the shape of the DOS near the Fermi level. Challenge: Accurate DOS requires dense k‑point sampling and careful smearing, especially for materials with sharp features.

EXCHANGE‑CORRELATION FUNCTIONAL (XC) – the term in DFT that accounts for… #

Related terms: LDA, GGA, hybrid functional. The choice of XC functional strongly influences total‑energy predictions, structural parameters, and electronic properties. Example: Using the Perdew‑Burke‑Ernzerhof (PBE) GGA functional to compute the lattice constant of Al, which typically yields a slight overestimation compared with experiment. Practical application: Selecting a hybrid functional such as HSE06 for accurate band‑gap calculations of wide‑gap semiconductors. Challenge: Hybrid functionals increase computational cost dramatically; meta‑GGA and dispersion‑corrected functionals add further complexity.

FERMI LEVEL – the energy at which the probability of electron occupation… #

Related terms: chemical potential, Fermi surface. In total‑energy calculations, the Fermi level determines the occupation of electronic states and therefore influences the total energy and forces. Example: Shifting the Fermi level in a doped silicon calculation to simulate n‑type carriers. Practical application: Aligning band edges of heterojunctions for photovoltaic device modeling. Challenge: Accurate determination of the Fermi level requires fine k‑point meshes and appropriate smearing schemes for metals.

FOURIER TRANSFORM – mathematical operation that converts real‑space quant… #

g., charge density) into reciprocal‑space representations, enabling efficient evaluation of convolution integrals. Related terms: reciprocal space, plane‑wave basis. Both Quantum Espresso and VASP rely on fast Fourier transforms (FFT) to switch between real and reciprocal space during SCF cycles. Example: Using a 64×64×64 FFT grid to represent the electron density of a 2‑atom primitive cell. Practical application: Reducing the computational scaling of the Hartree potential from O(N²) to O(N log N). Challenge: Selecting an FFT grid dense enough to avoid aliasing while keeping memory usage manageable.

GENERALIZED GRADIENT APPROXIMATION (GGA) – a class of XC functionals that… #

Related terms: PBE, PW91, meta‑GGA. GGA functionals are widely used for structural optimization because they often yield better lattice constants and bulk moduli. Example: Optimizing the geometry of MgO with the PBE functional and comparing the resulting lattice parameter with experimental data. Practical application: High‑throughput screening of oxides for catalytic activity where accurate bond lengths are crucial. Challenge: GGA may over‑delocalize electrons, leading to inaccurate magnetic moments in transition‑metal oxides.

GPAW – a real‑space DFT code that can be interfaced with VASP‑type PAW da… #

Related terms: projector‑augmented wave, real‑space grid. Although not the primary focus of the certificate, familiarity with alternative implementations helps students appreciate methodological differences. Example: Reproducing the VASP total energy of bulk Ti using GPAW’s PAW potentials and a 0.18 Å grid spacing. Practical application: Benchmarking the effect of basis‑set choice on formation energies. Challenge: Ensuring consistent pseudopotential or PAW dataset versions across codes to avoid systematic discrepancies.

HARTREE POTENTIAL – the classical electrostatic potential generated by th… #

Related terms: Poisson solver, FFT. In DFT calculations, the Hartree term contributes to the total energy and influences the Kohn‑Sham potential. Example: Using the Poisson solver in Quantum Espresso to obtain the Hartree energy of a water molecule. Practical application: Analyzing charge redistribution in heterostructures by inspecting the Hartree potential profile. Challenge: Accurate evaluation requires a sufficiently fine real‑space grid; coarse grids introduce spurious self‑interaction errors.

HYBRID FUNCTIONAL – an XC functional that mixes a fraction of exact (Hart… #

Related terms: HSE06, PBE0, screened exchange. Hybrid functionals are computationally demanding because they require evaluation of non‑local exchange integrals. Example: Calculating the band gap of GaN with the HSE06 functional, obtaining a value within 0.1 eV of experiment. Practical application: Designing optoelectronic materials where precise band‑edge positions are essential. Challenge: Scaling of exact exchange with system size limits hybrid calculations to a few hundred atoms even on modern supercomputers.

IONIC RELAXATION – the process of minimizing the total energy with respec… #

Related terms: geometry optimization, force convergence. In VASP, ionic steps are driven by the forces obtained from the SCF cycle; in Quantum Espresso, the “vc‑relax” calculation type performs simultaneous cell‑shape and ionic relaxation. Example: Relaxing a 2 × 2 × 2 supercell of Si with a vacancy until forces are below 0.01 eV/Å. Practical application: Determining defect formation energies and migration barriers. Challenge: Converging forces for metallic systems often requires denser k‑point meshes and tighter SCF thresholds.

K‑POINT SAMPLING – discretization of the Brillouin zone into a finite set… #

Related terms: Monkhorst‑Pack grid, Gamma‑centered mesh. Proper k‑point sampling is vital for accurate total energies; insufficient sampling leads to “k‑point noise” and erroneous forces. Example: Employing a 4×4×4 Gamma‑centered grid for a 20‑atom perovskite slab, then testing convergence by increasing to 6×6×6. Practical application: High‑throughput calculations where a standard k‑mesh is applied across many materials. Challenge: Metals require especially dense meshes to resolve the Fermi surface, increasing computational expense.

LDA (LOCAL DENSITY APPROXIMATION) – the simplest XC functional that assum… #

Related terms: Perdew‑Zunger, Ceperley‑Alder. LDA often underestimates lattice constants but can provide accurate bulk moduli for many metals. Example: Computing the equilibrium volume of bulk Al with LDA, obtaining a value ~2 % smaller than experiment. Practical application: Benchmarking more sophisticated functionals; LDA serves as a baseline for error analysis. Challenge: LDA lacks gradient information, leading to poor description of systems with rapidly varying densities such as surfaces and molecules.

MAGNETIC MOMENT – a vector quantity representing the net spin polarizatio… #

Related terms: spin‑polarized DFT, collinear magnetism. Total‑energy calculations can predict magnetic ordering (ferro‑, antiferro‑, ferri‑) by comparing energies of different spin configurations. Example: Calculating the magnetic moment of Fe in bcc Fe using spin‑polarized VASP, obtaining ~2.2 µB per atom. Practical application: Designing magnetic storage materials where high moments and anisotropy are desired. Challenge: Converging magnetic solutions may require careful initialization of the spin density and the use of symmetry‑breaking perturbations.

METAGGA FUNCTIONAL – an advanced class of XC functionals that depend on t… #

Related terms: SCAN, TPSS. Meta‑GGA functionals often improve thermochemical accuracy without the full cost of hybrids. Example: Applying the SCAN functional to compute the formation energy of TiO₂, achieving ≈0.1 eV per formula unit of experimental values. Practical application: High‑throughput screening where a balance between accuracy and speed is needed. Challenge: Some meta‑GGAs are more sensitive to numerical parameters, requiring finer grids and stricter SCF convergence.

MOST PROBABLE CONFIGURATION (MPC) – the atomic arrangement that minimizes… #

Related terms: ground‑state structure, global optimization. In total‑energy studies, the MPC serves as the reference for defect formation energies and phase diagrams. Example: Using a genetic algorithm to locate the lowest‑energy Si‑Ge alloy configuration at 50 % composition. Practical application: Predicting stable crystal structures of novel compounds before synthesis. Challenge: The configurational space grows combinatorially; exhaustive search is impossible, necessitating heuristic methods.

PAW (PROJECTOR‑AUGMENTED WAVE) METHOD – a technique that reconstructs the… #

Related terms: pseudopotential, augmentation charge. Both VASP and Quantum Espresso support PAW datasets, enabling high‑precision total‑energy calculations for transition metals and heavy elements. Example: Employing the PAW‑PBE potential for Fe to compute magnetic moments and compare with all‑electron FLAPW results. Practical application: Modeling materials under extreme pressure where core‑level effects become important. Challenge: Generating reliable PAW datasets for exotic elements; inconsistencies between PAW and norm‑conserving pseudopotentials can lead to systematic errors.

PHONON CALCULATIONS – determination of vibrational frequencies and eigenv… #

Related terms: DFPT, finite‑difference method, phonon dispersion. Within the DFT framework, phonons are obtained either via Density‑Functional Perturbation Theory (DFPT) or by constructing supercells and displacing atoms. Example: Using Quantum Espresso’s PHonon module to compute the phonon band structure of NaCl and identify the acoustic‑optic splitting at the Γ point. Practical application: Evaluating thermodynamic properties (entropy, free energy) and predicting phase transitions. Challenge: Accurate phonon spectra require dense q‑point meshes and well‑converged forces; anharmonic effects are neglected in the harmonic approximation.

PLANE‑WAVE BASIS SET – a set of sinusoidal functions used to expand the e… #

Related terms: energy cutoff, basis‑set convergence. Plane waves are orthogonal and systematically improvable, making them ideal for total‑energy calculations in periodic crystals. Example: Setting “ENCUT = 520 eV” in VASP for a system containing O, Si, and Al to ensure convergence of total energy within 1 meV per atom. Practical application: Uniform treatment of diverse chemical environments within a single computational framework. Challenge: Large cutoffs increase memory and CPU demands; for elements with hard pseudopotentials (e.g., transition metals), cutoffs may exceed 600 eV.

PROJECTED DENSITY OF STATES (PDOS) – the decomposition of the total DOS o… #

Related terms: orbital‑resolved DOS, Mulliken analysis. PDOS is extracted after a converged SCF calculation and often visualized alongside the total DOS. Example: Identifying the dominant Ti‑3d contribution near the Fermi level in TiO₂, confirming its role in photocatalytic activity. Practical application: Designing dopants that introduce states at desired energies for electronic devices. Challenge: The projection depends on the choice of atomic‑like basis; overlapping spheres can lead to double counting.

QUASI‑NEWTON OPTIMIZATION – a family of algorithms (e #

g., BFGS) that approximate the Hessian matrix to efficiently locate the minimum of the total‑energy surface during geometry relaxation. Related terms: conjugate‑gradient, trust‑radius. In VASP, the IBRION = 2 tag activates the BFGS quasi‑Newton method, often converging faster than conjugate‑gradient for well‑behaved systems. Example: Relaxing a 100‑atom slab of Au using the BFGS method, achieving force convergence in 12 ionic steps. Practical application: Rapid structural optimization for large supercells in surface science. Challenge: For highly anharmonic potentials, the approximate Hessian may become inaccurate, leading to oscillatory steps or failure to converge.

RECIPROCAL SPACE – the Fourier‑dual of real space; used to describe perio… #

Related terms: Brillouin zone, G‑vectors. In plane‑wave DFT, many operators (kinetic energy, Hartree potential) are evaluated in reciprocal space because of computational efficiency. Example: Representing the kinetic‑energy operator as |G|²/2 in the plane‑wave expansion. Practical application: Analyzing diffraction patterns to validate calculated lattice parameters. Challenge: Transforming quantities between real and reciprocal space introduces numerical errors if the FFT grid is insufficiently dense.

SCF (SELF‑CONSISTENT FIELD) CYCLE – iterative process that solves the Koh… #

Related terms: charge density mixing, convergence criterion. The SCF cycle is the core of any total‑energy calculation; its stability determines overall computational efficiency. Example: Setting “EDIFF = 1E‑6” in VASP to require total‑energy changes below 1 µeV between SCF steps. Practical application: Ensuring reliable forces for geometry optimization and molecular dynamics. Challenge: Metals, low‑symmetry systems, and poorly chosen initial guesses can cause SCF to stall or converge to metastable solutions.

SMEARING TECHNIQUES – methods that broaden the occupation numbers of elec… #

Related terms: Methfessel‑Paxton, Gaussian, Fermi‑Dirac. Smearing introduces an artificial temperature; the total energy must be extrapolated to zero temperature (e.g., using the “tetrahedron method” for final energies). Example: Applying a 0.02 Ry Methfessel‑Paxton smearing for a Cu bulk calculation, then correcting the energy with the “ISMEAR = -5” tetrahedron option. Practical application: Stabilizing SCF for complex alloys with partially filled d‑bands. Challenge: Excessive smearing can distort the electronic structure; careful selection of smearing width is required.

SPIN‑POLARIZED CALCULATIONS – DFT runs in which spin‑up and spin‑down ele… #

Related terms: collinear magnetism, non‑collinear magnetism. In VASP, setting “ISPIN = 2” activates spin polarization; in Quantum Espresso, the “nspin = 2” flag does the same. Example: Computing the antiferromagnetic ordering of MnO by initializing opposite spins on neighboring Mn atoms. Practical application: Predicting magnetic phase diagrams and Curie temperatures. Challenge: Convergence can be slower for magnetic systems; the initial magnetic moment must be supplied to avoid collapse to a non‑magnetic solution.

SUPERCELL APPROACH – constructing a larger periodic cell that contains de… #

Related terms: defect formation energy, slab model. The supercell must be large enough to minimize spurious interactions between periodic images. Example: Using a 3×3×3 Si supercell to model a vacancy, ensuring that the vacancy‑vacancy distance exceeds 10 Å. Practical application: Calculating migration barriers for diffusion by creating a series of relaxed images (NEB method). Challenge: Larger supercells increase the number of electrons, raising the computational cost and demanding more stringent k‑point reductions.

TIME‑DEPENDENT DFT (TDDFT) – an extension of DFT that treats the electron… #

Related terms: linear response, Casida equation. While not a primary total‑energy method, TDDFT builds on the same ground‑state potentials, allowing seamless transition from static to dynamic simulations. Example: Computing the absorption spectrum of a TiO₂ nanoparticle using the real‑time TDDFT implementation in Quantum Espresso. Practical application: Designing photocatalysts with tailored light‑absorption characteristics. Challenge: The accuracy of exchange‑correlation kernels in TDDFT remains limited for charge‑transfer excitations.

VASP (VIENNA AB‑INITIO SIMULATION PACKAGE) – a commercial plane‑wave DFT… #

Related terms: VASP INCAR, POTCAR, KPOINTS. VASP is widely used in the certificate program for its robustness and extensive documentation. Example: Performing a static self‑consistent calculation of bulk Al with “ENCUT = 400 eV”, “ISMEAR = 1”, and “EDIFF = 1E‑6”. Practical application: High‑throughput materials screening pipelines that automate VASP input generation and post‑processing. Challenge: License restrictions limit distribution of VASP binaries; students must access the software through institutional licenses.

VIBRATIONAL FREE ENERGY – the contribution of lattice vibrations to the H… #

Related terms: quasiharmonic approximation, phonon density of states. Adding vibrational free energy to the static DFT total energy yields temperature‑dependent phase stability predictions. Example: Computing the Gibbs free energy of MgO at 300 K by integrating the phonon DOS and comparing with the free energy of NaCl. Practical application: Constructing temperature‑pressure phase diagrams for geophysical minerals. Challenge: Harmonic approximation neglects anharmonicity; for high‑temperature phases, quasiharmonic or molecular‑dynamics‑based methods are required.

VIBRATIONAL PROPERTIES – properties derived from the phonon spectrum, suc… #

Related terms: phonon dispersion, Grüneisen parameter. Total‑energy calculations supply the force constants needed to evaluate these properties. Example: Using the Phonopy package together with VASP forces to calculate the heat capacity of diamond up to 1000 K. Practical application: Assessing materials for thermal management in electronics. Challenge: Accurate thermal conductivity predictions require inclusion of phonon‑phonon scattering, which goes beyond harmonic total‑energy calculations.

WANNIER FUNCTIONS – localized orbitals constructed from Bloch states that… #

Related terms: maximally localized Wannier functions, Wannier90. In total‑energy workflows, Wannier functions are often generated after an SCF run to study transport properties or to build tight‑binding models. Example: Generating Wannier functions for the Fe 3d bands of bcc Fe and interpolating the Fermi surface for high‑resolution transport calculations. Practical application: Computing anomalous Hall conductivity in magnetic materials. Challenge: Selecting an appropriate initial projection and disentanglement window is non‑trivial, especially for entangled bands.

WAVEFUNCTION CONVERGENCE – the criterion that the change in Kohn‑Sham orb… #

Related terms: energy convergence, charge density convergence. Tight wavefunction convergence is essential for accurate forces and stress tensors. Example: Setting “NELM = 200” and “EDIFF = 1E‑8” in VASP to achieve wavefunction convergence suitable for phonon calculations. Practical application: High‑precision total‑energy differences for reaction barriers. Challenge: Achieving wavefunction convergence can be difficult for metallic systems or when using hybrid functionals.

WANNIER90 – a post‑processing code that constructs maximally localized Wa… #

Related terms: disentanglement, interpolation. The workflow typically involves an SCF run, a non‑self‑consistent band calculation, and then Wannier90 to generate the localized basis. Example: Using Wannier90 to interpolate the band structure of SrTiO₃, reproducing DFT bands with sub‑meV accuracy. Practical application: Enabling efficient Brillouin‑zone integration for large‑scale transport calculations. Challenge: Properly handling entangled bands and selecting frozen windows to avoid spurious dispersion.

ZERO‑POINT ENERGY (ZPE) – the quantum mechanical energy associated with t… #

Related terms: vibrational free energy, harmonic approximation. ZPE corrections are essential when comparing total energies of light‑atom systems (e.g., hydrogen storage materials). Example: Adding the ZPE of H₂ (≈0.27 eV) to the DFT energy of a metal hydride to obtain a more accurate formation enthalpy. Practical application: Predicting reaction enthalpies for catalytic processes involving H₂. Challenge: ZPE depends on the phonon spectrum; any error in force constants propagates directly into the ZPE correction.

ZONE‑CENTER (Γ) POINT – the k‑point at the origin of the Brillouin zone;… #

Related terms: Γ‑point only, supercell calculations. While Γ‑only sampling dramatically reduces computational cost, it may miss important dispersion effects in small cells. Example: Performing a Γ‑only calculation for a 4 × 4 × 4 Si supercell containing a vacancy, then verifying convergence with a sparse k‑mesh. Practical application: Rapid screening of defect configurations before refinement with denser k‑point grids. Challenge: Metallic systems require more than Γ‑point sampling to capture Fermi‑surface features accurately.

ZONE‑BOUNDARY (X, L, K) POINTS – high‑symmetry k‑points located at the ed… #

Related terms: high‑symmetry path, band structure plot. In total‑energy calculations, energies at these points are used to construct dispersion curves that reveal band gaps and effective masses. Example: Calculating the band structure of GaAs along Γ‑X‑L‑Γ and identifying the direct band gap at Γ. Practical application: Evaluating carrier mobility by extracting curvature near the band extrema. Challenge: Accurate band energies require well‑converged SCF runs and sufficiently dense k‑point meshes to avoid artificial band flattening.

ZERO‑TEMPERATURE TOTAL ENERGY – the DFT energy obtained after SCF converg… #

Related terms: ground‑state energy, static calculation. This quantity forms the baseline for constructing free‑energy surfaces and phase diagrams. Example: Reporting the zero‑temperature total energy of NaCl as –7.84 eV per formula unit after full relaxation. Practical application: Comparing relative stabilities of polymorphs on the basis of static DFT energies. Challenge: Neglecting zero‑point vibrational energy can lead to incorrect ordering for light‑atom compounds.

ZERO‑POINT VIBRATION CORRECTION – the adjustment applied to static DFT en… #

Related terms: ZPE, vibrational free energy. The correction is computed from phonon frequencies and is especially important for systems containing hydrogen or other light atoms. Example: Adding a ZPE correction of 0.12 eV to the DFT energy of LiH to improve agreement with experimental formation enthalpy. Practical application: Accurate thermochemical predictions for hydrogen storage alloys. Challenge: Requires a separate phonon calculation; any error in force constants directly affects the correction magnitude.

ZONE‑CENTER SYMMETRY BREAKING – a phenomenon where the Γ‑point calculatio… #

Related terms: symmetry reduction, supercell artifacts. When modeling ferroelectric distortions or charge density waves, a Γ‑only grid may suppress the relevant phonon modes. Example: Observing that a Γ‑only calculation of BaTiO₃ does not show the expected tetragonal distortion, requiring a denser k‑mesh to activate the soft mode. Practical application: Properly capturing structural phase transitions in perovskites. Challenge: Determining the minimal k‑point sampling that preserves the physical instability while keeping computational cost reasonable.

ZONE‑SAMPLING SCHEME – the systematic procedure for selecting k‑points (e #

g., Monkhorst‑Pack, Gamma‑centered) that balances accuracy and efficiency. Related terms: k‑mesh density, symmetry reduction. The choice of scheme influences total‑energy convergence, especially for low‑symmetry or metallic systems. Example: Using a 6×6×6 Monkhorst‑Pack grid for cubic Ti and a 4×4×4 Gamma‑centered grid for a hexagonal Mg slab. Practical application: Standardizing k‑point generation across a materials database to ensure comparable energies. Challenge: Automatic generation tools must respect crystal symmetry to avoid unnecessary k‑point duplication.

ZONE‑CENTER OPTICAL TRANSITIONS – electronic excitations that occur at th… #

Related terms: dipole selection rules, exciton binding energy. In total‑energy calculations, the dipole matrix elements can be extracted from the wavefunctions at Γ to estimate transition probabilities. Example: Computing the oscillator strength of the direct transition in GaN at Γ using the momentum matrix elements from VASP. Practical application: Designing LEDs and laser diodes where Γ‑point transitions dominate emission. Challenge: Accurate transition energies often require beyond‑DFT methods (e.g., GW) to correct the underestimation of band gaps.

ZONE‑BOUNDARY PHONON ANOMALIES – softening of phonon modes at the edge of… #