Search NASA⌕ Search

SEARCH · Search NASA

Results for “time dependent 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Photophysical Properties and Electronic Structure of Hydroporphyrin Dyads Exhibiting Strong Through-Space and Through-Bond Electronic Interactions

Electronic interactions between tetrapyrroles are utilized in natural photosynthetic systems to tune the light-harvesting and energy/charge-transfer processes in these assemblies. Such interactions also can be employed to tailor the electronic properties of tetrapyrrolic dyads and larger arrays for use in materials science and biomedical research. Here we have utilized static and time-resolved optical spectroscopy to characterize the optical absorption and emission properties of a set of chlorin and bacteriochlorin dyads with varied degrees of through-bond (TB) and through-space (TS) interactions between the constituent macrocycles. The dyads consist of two chlorins or two bacteriochlorins joined by a linker that utilizes a triple-double-triple bond (enediyne) motif in which the double bond portion is a ester-substituted ethylene or o-phenylene unit. The photophysical studies are coupled with density functional theory (DFT) calculations to probe the ground-state molecular orbital (MO) characteristics of the dyads and time-dependent calculations (TDDFT) to elucidate excited-state properties. The latter include electronic characteristics of the singlet excited state manifold and the absorption transitions to these states from the electronic ground state. Comparison of the MO and calculated spectral properties of each dyad with the linker present versus disrupted (by eliminating the double-bond portion) gives insight into the relative contributions of TB versus TS interactions to the electronic properties of the dyads. The results show that the TB and TS contributions are additive (constructively interfere), which is not always the case for molecular dyads. Most of the dyads have shorter lifetimes of the lowest singlet excited state compared to the parent monomer, which derives from increased S 1 → S 0 internal conversion. The enhancement is greater for the dyads in benzonitrile than in toluene. Furthermore, the studies provide insights into the nature of the electronic interactions between the constituents in the tetrapyrrole arrays and how these interactions dictate the spectral properties and excited-state decay characteristics.

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

An Impedance Humidity Sensor Based on CVD Grown WSe 2 2D Films

2D materials-based planar devices have been sought after for gas and humidity sensor applications owing to their high sensitivity due to their ultrathin nature and increased surface area. Herein, 2D WSe 2 films-based planner devices were fabricated to evaluate their performance for potential application in humidity sensors. The structure and morphology of the WSe 2 sensing films were studied by Raman spectroscopy, photoluminescence, optical microscopy, and scanning electron microscopy (SEM) techniques. The relative humidity-dependent sensing performances of the 2D WSe 2 films were evaluated using current–voltage (I–V), relative impedance, and complex total impedance (electrochemical impedance spectroscopy, EIS) techniques. This sensor showed a very stable and repeatable performance investigated over the period of 18 months. The response and recovery times of the WSe 2 -based impedance sensors were found to be 40 and 30 s, respectively. The WSe 2 sensor exhibits a 2 ± 1% minimum hysteresis at a lower humidity level (50% RH) and around 8 ± 2% maximum hysteresis at relative humidity level (70% RH). This value of hysteresis is in the acceptable range during the initial investigation of any material to explore its potential for application in humidity sensors. The density functional theory (DFT) analysis was performed with the interaction from water molecules of the WSe 2 sample for a better understanding of the sensing mechanism of devices.

2D- WSe2↗
DOE OSTI2021

Solid-Density Ion Temperature from Redshifted and Double-Peaked Stark Line Shapes

Heβ spectral line shapes are important for diagnosing temperature and density in many dense plasmas. This work presents Heβ line shapes measured with high spectral resolution from solid-density plasmas with minimized gradients. The line shapes show hallmark features of Stark broadening, including quantifiable redshifts and double-peaked structure with a significant dip between the peaks; these features are compared to models through a Markov chain Monte Carlo framework. Line shape theory using the dipole approximation can fit the width and peak separation of measured line shapes, but it cannot resolve an ambiguity between electron density n e and ion temperature T i , since both parameters influence the strength of quasistatic ion microfields. Here a line shape model employing a full Coulomb interaction for the electron broadening computes self-consistent line widths and redshifts through the monopole term; redshifts have different dependence on plasma parameters and thus resolve the n e – T i ambiguity. The measured line shapes indicate densities that are 80–100% of solid, identifying a regime of highly ionized but well-tamped plasma. This analysis also provides the first strong evidence that dense ions and electrons are not in thermal equilibrium, despite equilibration times much shorter than the duration of x-ray emission; cooler ions may arise from nonclassical thermalization rates or anomalous energy transport. Furthermore, the experimental platform and diagnostic technique constitute a promising new approach for studying ion-electron equilibration in dense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Dynamical structure factors of warm dense matter from time-dependent orbital-free and mixed-stochastic-deterministic density functional theory

Abstract We present the first calculations of the inelastic part of the dynamical structure factor (DSF) for warm dense matter (WDM) using time-dependent orbital-free density functional theory (TD-OF-DFT) and mixed-stochastic-deterministic (mixed) Kohn Sham TD-DFT (KS TD-DFT). WDM is an intermediate phase of matter found in planetary cores and laser-driven experiments, where the accurate calculation of the DSF is critical for interpreting x-ray Thomson scattering measurements. Traditional TD-DFT methods, while highly accurate, are computationally expensive, motivating the exploration of TD-OF-DFT and mixed TD-KS-DFT as more efficient alternatives. We applied these methods to experimentally measured WDM systems, including solid-density aluminum and beryllium, compressed beryllium, and carbon–hydrogen mixtures. Our results show that TD-OF-DFT requires a dynamical kinetic energy potential in order to qualitatively capture the plasmon response. Additionally, it struggles with capturing bound electron contributions. In contrast, mixed TD-KS-DFT offers greater accuracy in distinguishing bound and free electron effects, aligning well with experimental data, though at a higher computational cost. This study highlights the trade-offs between computational efficiency and accuracy, demonstrating that TD-OF-DFT remains a valuable tool for rapid scans of parameter space, while mixed TD-KS-DFT should be preferred for high-fidelity simulations. Our findings provide insight into the future development of DFT methods for WDM and suggest potential improvements for TD-OF-DFT.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular Simulation of Lithium Carbonate Reactive Vapor–Liquid Equilibria Using a Deep Potential Model

We developed a first-principles machine learning model for the reactive vapor–liquid phase behavior of molten Li 2 CO 3 . The model was trained on ab initio electronic density functional theory data using the Deep Potential (DP) methodology, and its accuracy was evaluated by comparing model predictions of density and viscosity to experimental measurements. Direct coexistence simulations with the DP model over time scales of tens of nanoseconds were used to observe equilibrium dissociation of Li 2 CO 3 into CO 2 residing primarily in the vapor phase and Li 2 O which remains dissolved in the liquid. The simulations covered a range of temperatures, overall system sizes, and vapor-to-liquid volume ratios. Results were analyzed in terms of the observed chemical composition of the liquid and vapor phases, product structure, and CO 2 partial pressures. In addition, we calculated equilibrium constants for the dissociation reaction by assuming ideal-solution behavior for the liquid. As expected on the basis of thermodynamic arguments and prior experiments for this system, the observed partial pressure of CO 2 in the gas phase depends on both the temperature and the ratio of vapor to liquid volumes, while the calculated equilibrium constants only depend on temperature. DP model predictions for the equilibrium constant of the reaction are generally consistent with the available experimental measurements. Furthermore, the present study establishes the validity of the DP methodology for the description of reactive, multiphase equilibria from first principles, with possible applications to many other systems of scientific and technological interest even in the absence of relevant experimental measurements.

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

Ring currents modulate optoelectronic properties of aromatic chromophores at 25 T

The properties of organic molecules can be influenced by magnetic fields, and these magnetic field effects are diverse. They range from inducing nuclear Zeeman splitting for structural determination in NMR spectroscopy to polaron Zeeman splitting organic spintronics and organic magnetoresistance. A pervasive magnetic field effect on an aromatic molecule is the aromatic ring current, which can be thought of as an induction of a circular current of π-electrons upon the application of a magnetic field perpendicular to the π-system of the molecule. While in NMR spectroscopy the effects of ring currents on the chemical shifts of nearby protons are relatively well understood, and even predictable, the consequences of these modified electronic states on the spectroscopy of molecules has remained unknown. In this work, we find that photophysical properties of model phthalocyanine compounds and their aggregates display clear magnetic field dependences up to 25 T, with the aggregates showing more drastic magnetic field sensitivities depending on the intermolecular interactions with the amplification of ring currents in stacked aggregates. These observations are consistent with ring currents measured in NMR spectroscopy and simulated in time-dependent density functional theory calculations of magnetic field-dependent phthalocyanine monomer and dimer absorption spectra. Here, we propose that ring currents in organic semiconductors, which commonly comprise aromatic moieties, may present new opportunities for the understanding and exploitation of combined optical, electronic, and magnetic properties.

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

Quantum dot polarisation converter in an optomechanical cavity

We propose a scheme of a quantum photon polarisation converter, which is based on controlled electron – photon – phonon transitions in a hybrid semiconductor nanostructure. This structure consists of a GaAs/InAs quantum dot (QD) that has a parallelepiped shape and contains a single electron, and an optomechanical microcavity (MC) based on a photonic crystal (PC) that supports two orthogonally polarised photonic modes and one mechanical (phonon) mode. Within the framework of the microscopic theory, the QD and MC performance characteristics are found. Populations of states of the system as functions of time and its parameters are calculated. The principal possibility of photon polarisation conversion using transitions in a five-level resonance scheme for coherent (single-photon) and steady-state (subphoton) regimes is shown. The MC optical and mechanical spectra are simulated, and the PC structure parameters are selected to ensure the efficient operation of the converter. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Revisiting superconductivity in the extended one-band Hubbard model: Pairing via spin and charge fluctuations

The leading superconducting instabilities of the two-dimensional extended repulsive one-band Hubbard model within spin-fluctuation pairing theory depend sensitively on electron density, band, and interaction parameters. We map out the phase diagrams within a random-phase-approximation spin- and charge-fluctuation approach, and find that while B 1g (d x 2 –y 2 ) and B 2g (d xy ) pairing dominates in the absence of repulsive longer-range Coulomb interactions V NN , the latter induces pairing in other symmetry channels, including, e.g., A 2g (g-wave), nodal A 1g (extended s-wave), or nodal E u (p-wave) spin-triplet superconductivity. At the lowest temperatures, transition boundaries in the phase diagrams between symmetry-distinct spin-singlet orders generate complex time-reversal symmetry broken superpositions. By contrast, we find that boundaries between singlet and triplet regions are characterized by first-order transitions. Finally, motivated by recent photoemission experiments, we have determined the influence of an additional explicitly attractive nearest-neighbor interaction, V NN < 0, on the superconducting gap structure. Furthermore, depending on the electronic filling, such an attraction boosts E u (p-wave) spin-triplet or B 1g (d x 2 –y 2 ) spin-singlet ordering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Operando probing dynamic migration of copper carbonyl during electrocatalytic CO2 reduction

Single crystals and shape-controlled nanocrystals are well known to exhibit facet-dependent catalytic properties. However, few studies have investigated how those nanocrystals evolve and (de)activate during reactions, calling for the development of nanoscale time-resolved operando methods. In this context, we have designed Cu nanocubes as a model system to elucidate the underlying driving force of dynamic nanocatalyst reconstruction during the CO2 reduction reaction (CO2RR). Operando electrochemical liquid-cell scanning transmission electron microscopy (EC-STEM) and synchrotron-based X-ray spectroscopy reveal the size- and potential-dependent complete transformation from (100)-oriented Cu@Cu2O nanocubes to polycrystalline metallic Cu nanograins under CO2RR conditions. In addition, machine learning-assisted operando four-dimensional STEM reveals that large Cu nanograins derived from nanocubes form mainly crystalline domains, while their smaller counterparts are more amorphous due to faster evolution kinetics. In situ Raman spectroscopy and density functional theory calculations suggest that CO drives the ejection of single Cu atoms, resulting in few-nanometre Cu clusters and the surface migration of highly mobile copper carbonyl (Cu–CO) species. Combined, these multimodal operando methods and theoretical approaches pave the way for understanding the complex structural evolution of energy-related nanocatalysts under electrochemical conditions.

Yang, Yao↗
DOE OSTI2021

Application of transport-based metric for continuous interpolation between cryo-EM density maps

Cryogenic electron microscopy (cryo-EM) has become widely used for the past few years in structural biology, to collect single images of macromolecules "frozen in time". As this technique facilitates the identification of multiple conformational states adopted by the same molecule, a direct product of it is a set of 3D volumes, also called EM maps. To gain more insights on the possible mechanisms that govern transitions between different states, and hence the mode of action of a molecule, we recently introduced a bioinformatic tool that interpolates and generates morphing trajectories joining two given EM maps. This tool is based on recent advances made in optimal transport, that allow efficient evaluation of Wasserstein barycenters of 3D shapes. As the overall performance of the method depends on various key parameters, including the sensitivity of the regularization parameter, we performed various numerical experiments to demonstrate how MorphOT can be applied in different contexts and settings. Finally, we discuss current limitations and further potential connections between other optimal transport theories and the conformational heterogeneity problem inherent with cryo-EM data.

3D shapes↗
DOE OSTI2020

Microscopic Theory of Plasmons in Substrate-Supported Borophene

Here, we compute the dielectric properties of freestanding and metal-supported borophene from first-principles time-dependent density functional theory. We find that both the low- and high-energy plasmons of borophene are fully quenched by the presence of a metallic substrate at borophene–metal distances smaller than ≃9 Å. Based on these findings, we derive an electrodynamic model of the interacting, momentum-dependent polarizability for a two-dimensional metal on a model metallic substrate, which quantitatively captures the evolution of the dielectric properties of borophene as a function of metal–borophene distance. Applying this model to a series of metallic substrates, we show that maximizing the plasmon energy detuning between borophene and substrate is the key material descriptor for plasmonic performance.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

High-Temperature Thermodynamics of Uranium from Ab Initio Modeling

We present high-temperature thermodynamic properties for uranium in its γ phase (γ-U) from first-principles, relativistic, and anharmonic theory. The results are compared to CALPHAD modeling. The ab initio electronic structure is obtained from density-functional theory (DFT) that includes spin–orbit coupling and an added self-consistent orbital-polarization (OP) mechanism for more accurate treatment of magnetism. The first-principles method is coupled to a lattice dynamics scheme that is used to model anharmonic lattice vibrations, namely, Self-Consistent Ab Initio Lattice Dynamics (SCAILD). The methodology can be summarized in the acronym DFT + OP + SCAILD. Upon thermal expansion, γ-U develops non-negligible magnetic moments that are included for the first time in thermodynamic theory. The all-electron DFT approach is shown to model γ-U better than the commonly used pseudopotential method. In addition to CALPHAD, DFT + OP + SCAILD thermodynamic properties are compared with other ab initio and semiempirical modeling and experiments. Our first-principles approach produces Gibbs free energy that is essentially identical to CALPHAD. The DFT + OP + SCAILD heat capacity is close to CALPHAD and most experimental data and is predicted to have a significant thermal dependence due to the electronic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1991

Particle motion in the tail current sheet

Theory of particle motion in current sheets is reviewed. For small, approximately constant normal magnetic field, Bz, particles oscillate about the current sheet and 'live' within the sheet for one-half gyroperiod based on Bz. This lifetime replaces the mean collision time in the Lorentzian conductivity and thus gives rise to the concept of an inertial (or gyro-) conductivity. A substorm model by Coroniti utilizes this conductivity to allow reconnection to proceed without anomalous processes, due to wave-particle interactions. Chaotic particle orbits may at times be important to the dynamics, depending on parameters such as particle energy, current sheet thickness, and field line curvature. A current sheet model with neutral line predicts a ridge structure and asymmetries in the distribution function. Ion distributions near the plasma sheet boundary layer, during the CDAW 6 interval, are consistent with the model predictions. In recent studies by Mitchell et al. and Williams et al., the major current carriers during the growth phase of a substorm were found to be adiabatic electrons not more than 1 keV, but just before a current disruption event, the tail current was mainly carried by energetic ions undergoing current sheet oscillation.

Speiser, T. W.↗
DOE OSTI2023

Non-Linear Optics at Twist Interfaces in h-BN/SiC Heterostructures

Understanding the emergent electronic structure in twisted atomically thin layers has led to the exciting field of twistronics. However, practical applications of such systems are challenging since the specific angular correlations between the layers must be precisely controlled and the layers have to be single crystalline with uniform atomic ordering. Here, in this work, an alternative, simple, and scalable approach is suggested, where nanocrystalline two-dimensional (2D) film on 3D substrates yields twisted-interface-dependent properties. Ultrawide-bandgap hexagonal boron nitride (h-BN) thin films are directly grown on high in-plane lattice mismatched wide-bandgap silicon carbide (4H-SiC) substrates to explore the twist-dependent structure-property correlations. Concurrently, nanocrystalline h-BN thin film shows strong non-linear second-harmonic generation and ultra-low cross-plane thermal conductivity at room temperature, which are attributed to the twisted domain edges between van der Waals stacked nanocrystals with random in-plane orientations. First-principles calculations based on time-dependent density functional theory manifest strong even-order optical nonlinearity in twisted h-BN layers. This work unveils that directly deposited 2D nanocrystalline thin film on 3D substrates could provide easily accessible twist-interfaces, therefore enabling a simple and scalable approach to utilize the 2D-twistronics integrated in 3D material devices for next-generation nanotechnology.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Relaxation time approximations in PAOFLOW 2.0

Abstract Regardless of its success, the constant relaxation time approximation has limited validity. Temperature and energy dependent effects are important to match experimental trends even in simple situations. We present the implementation of relaxation time approximation models in the calculation of Boltzmann transport in 2.0 and apply those to model band-structures. In addition, using a self-consistent fitting of the model parameters to experimental conductivity data, we provide a flexible tool to extract scattering rates with high accuracy. We illustrate the approximations using simple models and then apply the method to GaAs, Si, $$\hbox {Mg}_3\hbox {Sb}_2$$ Mg 3 Sb 2 , and $$\hbox {CoSb}_3$$ CoSb 3 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory

Direct simulation of the von Neumann dynamics for a general (pure or mixed) quantum state can often be expensive. One prominent example is the real-time time-dependent density functional theory (rt-TDDFT), a widely used framework for the first principle description of many-electron dynamics in chemical and materials systems. Practical rt-TDDFT calculations often avoid the direct simulation of the von Neumann equation, and solve instead a set of Schrödinger equations, of which the dynamics is equivalent to that of the von Neumann equation. However, the time step size employed by the Schrödinger dynamics is often much smaller. Here, in order to improve the time step size and the overall efficiency of the simulation, we generalize a recent work of the parallel transport (PT) dynamics for simulating pure states [An, Lin, Multiscale Model. Simul. 18, 612, 2020] to general quantum states. The PT dynamics provides the optimal gauge choice, and can employ a time step size comparable to that of the von Neumann dynamics. Going beyond the linear and near adiabatic regime in previous studies, we find that the error of the PT dynamics can be bounded by certain commutators between Hamiltonians, density matrices, and their derived quantities. Such a commutator structure is not present in the Schrödinger dynamics. We demonstrate that the parallel transport-implicit midpoint (PT-IM) method is a suitable method for simulating the PT dynamics, especially when the spectral radius of the Hamiltonian is large. The commutator structure of the error bound, and numerical results for model rt-TDDFT calculations in both linear and nonlinear regimes, confirm the advantage of the PT dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Accelerating Quantum Materials Development with Advances in Transmission Electron Microscopy

Quantum materials are driving a technology revolution in sensing, communication, and computing, while simultaneously testing many core theories of the past century. Materials such as topological insulators, complex oxides, superconductors, quantum dots, color center-hosting semiconductors, and other types of strongly correlated materials can exhibit exotic properties such as edge conductivity, multiferroicity, magnetoresistance, superconductivity, single photon emission, and optical-spin locking. These emergent properties arise and depend strongly on the material’s detailed atomic-scale structure, including atomic defects, dopants, and lattice stacking. In this review, we describe how progress in the field of electron microscopy (EM), including in situ and in operando EM, can accelerate advances in quantum materials and quantum excitations. We begin by describing fundamental EM principles and operation modes. We then discuss various EM methods such as (i) EM spectroscopies, including electron energy loss spectroscopy (EELS), cathodoluminescence (CL), and electron energy gain spectroscopy (EEGS); (ii) four-dimensional scanning transmission electron microscopy (4D-STEM); (iii) dynamic and ultrafast EM (UEM); (iv) complementary ultrafast spectroscopies (UED, XFEL); and (v) atomic electron tomography (AET). We describe how these methods could inform structure–function relations in quantum materials down to the picometer scale and femtosecond time resolution, and how they enable precision positioning of atomic defects and high-resolution manipulation of quantum materials. For each method, we also describe existing limitations to solve open quantum mechanical questions, and how they might be addressed to accelerate progress. Among numerous notable results, our review highlights how EM is enabling identification of the 3D structure of quantum defects; measuring reversible and metastable dynamics of quantum excitations; mapping exciton states and single photon emission; measuring nanoscale thermal transport and coupled excitation dynamics; and measuring the internal electric field and charge density distribution of quantum heterointerfaces- all at the quantum materials’ intrinsic atomic and near atomic-length scale. Finally, we conclude by describing open challenges for the future, including achieving stable sample holders for ultralow temperature (below 10K) atomic-scale spatial resolution, stable spectrometers that enable meV energy resolution, and high-resolution, dynamic mapping of magnetic and spin fields. With atomic manipulation and ultrafast characterization enabled by EM, quantum materials will be poised to integrate into many of the sustainable and energy-efficient technologies needed for the 21st century.

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

Predicting magnetic properties of single-molecule magnets from self-interaction-free density-functional theory (Final Report)

In this project we investigated electronic structure of an intermediate-sized copper-based molecule and magnetic and hyperfine properties of several small non-magnetic and magnetic molecules including transition-metal elements by applying self-interaction corrections to density functional theory (DFT). In a sufficient number of cases, DFT-calculated electronic structure and magnetic properties of single-molecule magnets qualitatively differ from corresponding experimental data. This is partly due to self-interacting electrons within the DFT formalism. Recently, an efficient method to correct the self-interactions was proposed, i.e., Fermi-Lowdin prbital (FLO) based self-interaction corrected (SIC) methodology, within DFT. Henceforth, this method is referred to as FLO-SIC method which exists in FLOSIC code. We used this FLO-SIC method for our studies of electronic structure and magnetic and hyperfine properties of small magnetic molecules and non-magnetic molecules. Our study will provide insight into predictions of magnetic properties of single-molecule magnets where self-interaction corrections play a critical role. There are two components of this project. In the first work, we studied the electronic structure of a planar mononuclear Cu-based molecule in two oxidation states, using DFT with the FLO-SIC method. We chose this system because it is small enough and it includes a transition metal element. We found that the standard FLO-SIC method takes too much compute time even for the small transition-metal molecule and so we slightly modified the method in order to expedite the process. In the dianionic state, we found that the FLO-SIC spin density agrees quantitatively with accurate quantum chemistry methods, while DFT spin density without self-interaction corrections are severely deviated from the quantum chemistry methods. We also showed that the energy gap between the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of the dianionic state is larger than that of the monoanionic state. This result is consistent with experimental data. In the second work, we investigated how the interaction between the electron spin and nuclear spin is affected by electron self-interactions within small non-magnetic and magnetic molecules. Such an interaction is called hyperfine interaction. For molecules without significant orbital angular momentum, the hyperfine interaction consists of Fermi contact and dipolar interaction terms. Since the Fermi contact term depends on electron spin density at the nuclear site, it would be highly affected by self-interaction corrections. Therefore, we calculated the hyperfine interaction for the small molecules using DFT with the slightly modified expedited FLO-SIC method which was obtained in the first work, and compared the results to experimental data and DFT calculations without self-interaction corrections. We found significant improvement of the Fermi contact term computed using the FLO-SIC method for small magnetic molecules. Overall, the first and second work provided positive outlook of application of the FLO-SIC method to magnetic molecules and systems including transition-metal elements.

36 MATERIALS SCIENCE↗