Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electronic structure theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Electronic structure of nitrogen-doped lutetium hydrides

First-principles density functional theory (DFT) calculations of supercell structures based on N-doped Fm $\bar3$m LuH 3 reveal configurations of Fm $\bar3$m Lu 8 H 23–x N that exhibit novel electronic properties such as flat bands, sharply peaked densities of states (van Hove singularities, vHs), and intersecting Dirac cones near the Fermi energy (E F ). These electronic properties are present when N substitutes H in the octahedral interstices of Fm $\bar3$ m LuH 3 . These structures also exhibit an interconnected metallic hydrogen network, a common feature of high-T c hydride superconductors. Electronic property systematics gives an estimate of T c for one structure that is well above the critical temperatures predicted for structures considered previously. DFT + U has an especially strong effect on one of the structures considered, enhancing the vHs and flat bands near E F . Finally, these results provide a basis for understanding the electronic properties observed for nitrogen-doped lutetium hydride.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange

In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Electronic structure of the kagome compound CaTi 3⁢ Bi 4 using high-field torque magnetometry and density functional theory

Here, we report systematic torque magnetometry measurements to investigate the electronic properties of the newly discovered kagome compound CaTi 3 ⁢Bi 4 . Electrical transport, magnetic susceptibility, and thermal measurements reveal no evidence of a magnetic ground state in this material. Torque data obtained in magnetic fields up to 41.5 T exhibit clear de Haas–van Alphen (dHvA) oscillations, with nine distinct frequencies ranging from 13 to 6164 T. Angular-dependent dHvA measurements show that, with the exception of the lowest frequency (13 T), all observed frequencies nearly follow a 1/cos ⁡𝜃 dependence, where 𝜃 is the angle between the crystallographic 𝑐 axis and the magnetic field direction. This behavior is characteristic of quasi-two-dimensional Fermi-surface sheets with nearly circular cross sections. To further elucidate the electronic structure, we performed density functional theory (DFT) calculations of the band structure and Fermi surface. The calculated bands reveal the presence of multiple Dirac points (DP), flat bands (FB), and van Hove singularities (VHS) near the Fermi level. The resulting Fermi surface consists of several quasi-two-dimensional cylindrical sheets, consistent with the experimentally observed 1/cos⁡ 𝜃 dependence. Notably, the theoretical dHvA frequencies, derived from extremal Fermi-surface cross-sectional areas, agree well with the experimental values and reproduce their angular dependence. Remarkably, all experimentally observed frequencies are captured by the DFT predictions. Pressure-dependent calculations up to 10 GPa show that the electronic features—DP, FB, and VHS—evolve systematically with pressure. In particular, the VHS shifts closer to the Fermi level, demonstrating that pressure acts as an effective tuning parameter in this material. These combined experimental and theoretical results provide a comprehensive understanding of the electronic structure of CaTi 3 ⁢Bi 4 and demonstrate how pressure can be used to tune its key electronic features, offering valuable guidance for exploring related kagome materials.

Shtefiienko, Kyryl [West Texas A & M University, C↗
DOE OSTI2024

A First-Principles Study of the Structural and Thermo-Mechanical Properties of Tungsten-Based Plasma-Facing Materials

Tungsten (W) and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield, and low long-term disposal radioactive footprint. Despite these relevant features, there is a lack of understanding of how the structural and mechanical properties of W-based alloys are affected by the temperature in fusion power plants. In this work, we present a study on the thermo-mechanical properties of five W-based plasma-facing materials. First-principles density functional theory (DFT) calculations are combined with the quasi-harmonic approximation (QHA) theory to investigate the electronic, structural, mechanical, and thermal properties of these W-based alloys as a function of temperature. The coefficient of thermal expansion, temperature-dependent elastic constants, and several elastic parameters, including bulk and Young’s modulus, are calculated. Our work advances the understanding of the structural and thermo-mechanical behavior of W-based materials, thus providing insights into the design and selection of candidate plasma-facing materials in fusion energy devices.

42 ENGINEERING↗
DOE OSTI2023

Influence of cation species on thermal expansion of Y 2 Si 2 O 7 –Gd 2 Si 2 O 7 solid solutions

Mixtures of Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 were synthesized by solid-state reaction at 1600°C and characterized via in situ x-ray diffraction (XRD) to determine their coefficients of thermal expansion (CTE). All solid solutions within the system exhibited the orthorhombic δ-RE 2 Si 2 O 7 (Pna2 1 ) structure. Thermal expansion measurements of Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 correlated well with reported values in literature, and all synthesized solid solutions exhibited CTEs between Y 2 Si 2 O 7 and Gd 2 Si 2 O 7 . Generally, there was a slight decrease in CTE exhibited by the materials with increasing Gd 2 Si 2 O 7 content, with Gd 2 Si 2 O 7 having the lowest CTEs and Y 2 Si 2 O 7 the highest CTEs. Here, the decrease in CTE was attributed to stronger bonds of Gd-O over Y-O, as determined by calculated crystal orbital Hamilton populations using density functional theory. However, such differences were very small and crystal structure was the dominating factor in CTE trends.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Similarities and differences between nickelate and cuprate films grown on a SrTiO 3 substrate

The recent discovery of superconductivity in Sr-doped NdNiO 2 films grown on SrTiO 3 started a novel field within unconventional superconductivity. Therefore, to understand the similarities and differences between nickelate and cuprate layers on the same SrTiO 3 substrate, here based on the density functional theory we have systematically investigated the structural, electronic, and magnetic properties of NdNiO 2 / SrTiO 3 and CaCuO 2 / SrTiO 3 systems. Our results revealed a strong lattice reconstruction in the case of NdNiO 2 / SrTiO 3 , resulting in a polar film, with the surface and interfacial NiO 2 layers presenting opposite displacements. To avoid the “polar catastrophe,” the NiO 2 surface to the vacuum reconstructs as well. However, for CaCuO 2 / SrTiO 3 , the distortions of those same two CuO 2 layers were in the same direction. In addition, we found this distortion to be approximately independent of the studied range of film thickness for the nickelate films. Furthermore, we also observed a two-dimensional electron gas at the interface between NdNiO 2 and SrTiO 3 , caused by the polar discontinuity, in agreement with recent literature. For NdNiO 2 / SrTiO 3 the two-dimensional electron gas extends over several layers, while for CaCuO 2 / SrTiO 3 this electronic rearrangement is very localized at the interface between CaCuO 2 and SrTiO 3 . The electronic reconstruction found at the interface involves a strong occupation of the Ti 3 d x y state. In both cases, there is a significant electronic charge transfer from the surface Ni or Cu layers to the Ti interface layer. The interfacial Ni and Cu layers are hole and electron doped, respectively. By introducing magnetism and electronic correlation, we observed that the d 3 z 2 - r 2 orbital of Ni becomes itinerant while the same orbital for Cu remains doubly occupied, establishing a clear two- vs one-orbital active framework for the description of these systems. Furthermore, we also observed a strong magnetic reconstruction at the NdNiO 2 surface to vacuum layer where magnetism is basically suppressed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Electronic and vibrational properties of bulk Cr 2 Ge 2 Te 6 from first-principles calculations

Here we use density functional theory based ab initio calculations to investigate the structural, vibrational, magnetic, and electronic properties of the layered ferromagnet Cr 2 Ge 2 Te 6 (CGT) that has attracted attention for potential spintronic applications. We optimize the structure for a fixed c/a ratio. Our results are in excellent agreement with experimental data on structure, phonons, and electronic properties. The use of van der Waals interactions and relativistic spin-orbit coupling yields accurate lattice constants and interlayer distances. Simultaneously, the adopted theoretical methods lead to a rigorous description of the vibrational normal modes, as well as the valence bands, that are in excellent agreement with Raman spectroscopy and angle-resolved photoemission spectroscopy spectra, respectively. The magnetic moment is slightly overestimated, and the magnetic anisotropy has the correct sign but is greater in magnitude than that in experiment. We find that inclusion of the on-site Coulomb repulsion on Cr−3d orbitals worsens agreement with experiment, in contrast to previous studies. Our results provide a path toward ab initio analysis of magnetic heterostructures based on CGT.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Analytical derivatives of the individual state energies in ensemble density functional theory. II. Implementation on graphical processing units (GPUs)

Conical intersections control excited state reactivity, and thus, elucidating and predicting their geometric and energetic characteristics are crucial for understanding photochemistry. Locating these intersections requires accurate and efficient electronic structure methods. Unfortunately, the most accurate methods (e.g., multireference perturbation theories such as XMS-CASPT2) are computationally challenging for large molecules. The state-interaction state-averaged restricted ensemble referenced Kohn–Sham (SI-SA-REKS) method is a computationally efficient alternative. The application of SI-SA-REKS to photochemistry was previously hampered by a lack of analytical nuclear gradients and nonadiabatic coupling matrix elements. We have recently derived analytical energy derivatives for the SI-SA-REKS method and implemented the method effectively on graphical processing units. We demonstrate that our implementation gives the correct conical intersection topography and energetics for several examples. Furthermore, our implementation of SI-SA-REKS is computationally efficient, with observed sub-quadratic scaling as a function of molecular size. This demonstrates the promise of SI-SA-REKS for excited state dynamics of large molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Electronic structure of the kagome metal YbTi 3 Bi 4 studied using torque magnetometry

This study investigates the electronic structure of the kagome metal YbTi 3 ⁢Bi 4 using high-field torque magnetometry. The torque signal measured at a maximum field of 41.5 T reveals clear de Haas–van Alphen (dHvA) oscillations with a major frequency peak at F δ ~ 130 T. By rotating the sample at various tilt angles Θ, we observed that F δ exhibits a nearly 1/cos⁡Θ dependence, indicating the presence of a quasi-two-dimensional (2D) Fermi surface (FS) in YbTi 3 ⁢Bi 4 . This argument is further supported by the detection of a forward-leaning, sawtoothlike waveform in the dHvA effect, a hallmark of 2D FS characteristics. Notably, we identified two high-frequency peaks near F χ ~ 1900 T and F λ ~ 5600 T; however, these peaks quickly disappear at Θ greater than 21°. To better understand experimental observations, we computed the electronic band structure and FS using ab initio density-functional theory (DFT). The electronic bands reveal the presence of several Dirac points, flat bands, and van Hove singularities near the Fermi level. Five bands cross the Fermi level and contribute to the FS of this material. The FS comprises cylindrical sheets, with theoretical frequencies from the FS pockets aligning well with the experimental dHvA frequencies. Several FS parameters characterizing F δ were determined by analyzing the temperature and field dependence of the dHvA oscillations using the Lifshitz-Kosevich theory. Further, the detailed electronic properties presented in this work provide critical insights into the electronic structure of YbTi 3 ⁢Bi 4 and other titanium-based kagome compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Mechanistic implications of excited high-spin states, spin–spin coupling, and differential [2Fe–2S] + cluster temperature relaxations in the electron-bifurcating NfnABC from Thermococcus sibiricus

Electron bifurcation (EB) is a mechanism of biological energy transduction in which multiple oxidation–reduction (redox) reactions are thermodynamically coupled within a single enzyme, enabling the enzyme to harness the excess free energy from an exergonic process to drive an endergonic process. Because of this unprecedented chemistry, there is interest to translate EB principles to artificial and bioengineered systems, but a hurdle is that knowledge pertaining to the fundamental design principles of EB enzymes remains scarce. Here, we investigated the fundamental physical and electronic properties of electron transfer sites in a spectroscopically uncharacterized member of the BfuABC family of EB enzymes, the NADH-dependent reduced-ferredoxin:NADP + oxidoreductase from Thermococcus sibiricus (Tsi NfnABC). Cryo-EM structures of Tsi NfnABC previously demonstrated that it contains twelve redox cofactors: two flavins (one FAD and one FMN), eight [4Fe–4S] clusters, and two [2Fe–2S] clusters. The FMN, one [4Fe–4S] cluster, and one [2Fe–2S] cluster comprise the bifurcating active site termed the electron-bifurcating flavobicluster (BF-FBC), which is found in all BfuABC family members. By using electron paramagnetic resonance spectroscopy, we identified spectral signatures originating from interactions between the FMN radical and [4Fe–4S] + cluster in the BF-FBC and observed temperature dependent behavior of the BF-FBC's [2Fe–2S] + cluster indicative of moderately slow spin–lattice relaxation. Additionally, we uncovered numerous spectral features corresponding to half-integer, S > ½ spin states of [4Fe–4S] + clusters, including one attributable to the consequences of lysine-ligation of a [4Fe–4S] cluster unique to NfnABC. We contextualize these findings to electron transfer theory and NfnABC's structure. Our insights further the understanding of how enzymes are designed to exert control over electron transfer to conduct thermodynamically challenging reactions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Ni Nanoparticles on CeO 2 (111): Energetics, Electron Transfer, and Structure by Ni Adsorption Calorimetry, Spectroscopies, and Density Functional Theory

The morphology, interfacial bonding energetics and charge transfer of Ni clusters and nanoparticles on slightly-reduced CeO 2-x (111) surfaces at 100 to 300 K have been studied using single crystal adsorption calorimetry (SCAC), low-energy ion scattering spectroscopy (LEIS), Xray photoelectron spectroscopy (XPS), low energy electron diffraction (LEED) and density functional theory (DFT). The initial heat of adsorption of Ni vapor decreased with the extent of pre-reduction (x) of the CeO 2-x (111), showing that stoichiometric ceria adsorbs Ni more strongly than oxygen vacancies. On CeO 1.95 (111) at 300 K, the heat dropped quickly with coverage in the first 0.1 ML, attributed to nucleation of Ni clusters on stoichiometric steps, followed by the Ni particles spreading onto less favorable terrace sites. At 100 K, the clusters nucleate on terraces due to slower diffusion. Adsorbed Ni monomers are in the +2 oxidation state, and they bind by ~45 kJ/mol more strongly to step sites than terraces. The measured heat of adsorption versus average particle size on terraces is favorably compared to DFT calculations. The Ce 3d XPS lineshape showed an increase in Ce 3+ /Ce 4+ ratio with Ni coverage, providing the number of electrons donated to the ceria per Ni atom. The charge transferred per Ni is initially large but strongly decreases with increasing cluster size for both experiments and DFT, and shows large differences between clusters at steps versus terraces. This charge is localized on the interfacial Ni and Ce atoms in their atomic layers closest to the interface. This knowledge is crucial to understanding the nature of the active sites on the surface of Ni-CeO 2 catalysts for which metal-oxide interactions play a very important role in the activation of O-H and C-H bonds. In conclusion, the changes in these interactions with Ni particle size (metal loading) and the extent of reduction of the ceria help to explain how previously reported catalytic activity and selectivity change with these same structural details.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Structure, Stability, and Electronic Properties of Boron Suboxide: A Density Functional Theory Study

Boron suboxide (B 6 O) is a boron-rich compound derived from the α-rhombohedral boron lattice with extreme hardness and unusual semiconducting properties. In this work, density functional theory (DFT) was used to show that unit cell volume, mechanical strength, band gaps, and thermodynamic stabilities of B 6 O were influenced by the interstitial elements and point defects at the icosahedral sites. While the hexagonal unit cell volume (HUCV) varies with interstitial occupancy, it is the icosahedral defect that weakens the intrinsic bulk modulus of B 6 O. Using the hybrid HSE functional, we confirmed that the perfect B 6 O bulk is a p-type semiconductor with a direct band gap of 2.8 eV. Furthermore, by screening α-boron compounds systematically, we found that a simple octet rule may offer a consistent explanation for the variations in the computed electronic structures. The formation free energies calculated over a wide range of temperatures (0–2500 K) and pressures (0–80 GPa) predict that formations of interstitial defects become favorable only at higher temperatures (ca. 1800 K) in bulk B 6 O lattices. The nudged elastic band (NEB) method was employed to identify the minimum energy pathways for the diffusions of dislocated B and O atoms. The diffusion of icosahedral B atoms has an energy barrier of 0.16 eV. More complex B diffusion paths involving the reorganization of icosahedral boron atoms incur higher barriers (>1 eV). In contrast, the diffusion of interstitial O atoms is facile with a barrier of 0.4 eV. Lastly, successive O insertions into the α-B lattice were performed using DFT to generate a basic understanding of the oxidation process. Furthermore, these calculations provide fundamental atomistic insights into the growth of B 6 O crystals and control of their point defects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Diverse electronic topography in a distorted kagome metal LaTi 3 ⁢Bi 4

Recent reports on a family of kagome metals of the form 𝐿⁢𝑛⁢Ti 3 ⁢Bi 4 (𝐿⁢𝑛 = Lanthanide) have stoked interest due to the combination of highly anisotropic magnetism and a rich electronic structure. The electronic structure near the Fermi level is proposed to exhibit Dirac points and van Hove singularities (VHSs). Here, in this manuscript, we use angle-resolved photoemission spectroscopy measurements in combination with density functional theory calculations to investigate the electronic structure of an interesting kagome metal LaTi 3 ⁢Bi 4 . Our results reveal multiple VHSs with one VHS located in the vicinity of the Fermi level. We clearly observe two flat bands, which originate from the destructive interference of wave functions within the Ti kagome motif. These flat bands and VHSs originate from Ti 𝑑 orbitals and are very responsive to the polarization of the incident beam. We notice a significant anisotropy in the electronic structure, resulting from the breaking of sixfold rotational symmetry in this material. Our findings demonstrate this member of Ti based kagome materials as a promising platform to explore novel emerging phenomena in the wider 𝐿⁢𝑛⁢Ti 3⁢ Bi 4 (𝐿⁢𝑛 = lanthanide) family.

electronic structure↗
DOE OSTI2023

Exceptional magnetic and magnetoelastic behavior of rare-earth non-centrosymmetric Sm 7 Pd 3

Magnetic compounds possessing an intrinsic combination of near-zero magnetization with high magnetic anisotropy are highly desirable for spinronic applications and memory recording. A comprehensive study of Sm 7 Pd 3 binary compound uncovered a unique combination of strong magnetoelastic behavior, very low net magnetization, and exceptionally high magnetic coercivity. The temperature-dependent X-ray synchrotron powder diffraction study indicates the abrupt changes in the compound's lattice parameters at the magnetic ordering temperature of T C =169 K, although the crystal structure remains non-centrosymmetric hexagonal Th 7 Fe 3 -type down to 6 K. Density functional theory calculations confirm high intrinsic magnetocrystalline anisotropy of Sm 7 Pd 3 , which explains the extremely large coercivity of the polycrystalline sample, up to H cr = 130 kOe at 2 K. In conclusion, this discovery brings to life a novel class of highly anisotropic materials that are distinctly different from known spintronic materials, making them interesting future systems for magnetic memory research.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. Here, the nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Full-potential KKR within the removed-sphere method: A practical and accurate solution to the Poisson equation

Here an efficient and accurate generalization of the removed-sphere method (RSM) to solve the Poisson equation for total charge density in a solid with space-filling convex Voronoi polyhedra (VPs) and any symmetry is presented. The generalized RSM avoids the use of multipoles and VP shape functions for cellular integrals, which have associated ill-convergent large, double-internal L sums in spherical-harmonic expansions, so that fast convergence in single-L sums is reached. Our RSM adopts full Ewald formulation to work for all configurations or when symmetry breaking occurs, such as for atomic displacements or elastic constant calculations. The structure-dependent coefficients A L that define RSM can be calculated once for a fixed structure and speed up the whole self-consistent-field procedure. The accuracy and rapid convergence properties are confirmed using two analytic models, including the Coulomb potential and energy. We then implement the full-potential RSM using the Green's function Korringa-Kohn-Rostoker (KKR) method for real applications and compare the results with other first-principle methods and experimental data, showing that they are equally as accurate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗