Search NASA⌕ Search

SEARCH · Search NASA

Results for “MOLECULAR ELECTRONICS”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Femtosecond Core-Level Spectroscopy Reveals Involvement of Triplet States in the Gas-Phase Photodissociation of Fe(CO) 5

Excitation of iron pentacarbonyl [Fe(CO) 5 ], a prototypical photocatalyst, at 266 nm causes the sequential loss of two CO ligands in the gas phase, creating catalytically active, unsaturated iron carbonyls. Despite numerous studies, major aspects of its ultrafast photochemistry remain unresolved because the early excited-state dynamics have so far eluded spectroscopic observation. This has led to the long-held assumption that ultrafast dissociation of gas-phase Fe(CO) 5 proceeds exclusively on the singlet manifold. Herein, we present a combined experimental–theoretical study employing ultrafast extreme ultraviolet transient absorption spectroscopy near the Fe M 2,3 -edge, which features spectral evolution on 100 fs and 3 ps time scales, alongside high-level electronic structure theory, which enables characterization of the molecular geometries and electronic states involved in the ultrafast photodissociation of Fe(CO) 5 . We assign the 100 fs evolution to spectroscopic signatures associated with intertwined structural and electronic dynamics on the singlet metal-centered states during the first CO loss and the 3 ps evolution to the competing dissociation of Fe(CO) 4 along the lowest singlet and triplet surfaces to form Fe(CO) 3 . Calculations of transient spectra in both singlet and triplet states as well as spin–orbit coupling constants along key structural pathways provide evidence for intersystem crossing to the triplet ground state of Fe(CO) 4 . Finally, our work presents the first spectroscopic detection of transient excited states during ultrafast photodissociation of gas-phase Fe(CO) 5 and challenges the long-standing assumption that triplet states do not play a role in the ultrafast dynamics.

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

Structure of Layers Formed by [2-(3,6-Disubstituted-9 H -carbazol-9-yl)ethyl]phosphonic Acids on Metal Oxides

[2-(9 H -Carbazol-9-yl)ethyl]phosphonic acid (2PACz) and its derivatives are being used extensively as hole-transport layers in organic and perovskite solar cells due to their ability to modify electrode work function, surface wettability, and in some cases, to improve active-layer adhesion, while minimizing interfacial energy losses. The orientation and coverage of surface modifiers significantly impact these functional properties; however, the detailed structure and packing in these overlayers are challenging to investigate, leading to a lack of understanding of how this structure and packing influence performance and limiting the ability to rationally design molecular modifiers. Here, we investigate monolayers of 2-(9 H -carbazol-9-yl)ethyl phosphonic acid derivatives (from here on referred to as X-2PACz) on indium tin oxide (ITO) and alpha phase aluminum oxide (α-Al 2 O 3 ) using a combination of X-ray photoelectron spectroscopy (XPS), near-edge X-ray absorption fine structure (NEXAFS) spectroscopy, and X-ray reflectivity (XRR). By correlating elemental ratios, molecular orientation, and electron density profiles, we directly quantify surface coverage, layer thickness, and molecular tilt across a series of chemically related monolayers. We find that 2PACz and the X-2PACz derivatives form dense monolayers on α-Al 2 O 3 and ITO, with similar surface coverages on either substrate that are inversely proportional to molecular steric bulk. These surface coverage results indicate that X-2PACz is sterically limited in its monolayer surface packing density, as opposed to site limited. Despite surface packing density differences between molecules, the NEXAFS data show a constant average molecular orientation of 61° to 65° between the plane of the carbazole and the substrate for all the molecules on both substrates. These results increase our general understanding of 2PACz derivatives as they become increasingly useful in high performance solar cells.

NEXAFS↗
DOE OSTI2026

Path integral molecular dynamics: A high-fidelity approach to quantum dynamics of electrons

We investigate electron transport in the uniform electron gas using ring-polymer molecular dynamics (RPMD). Working in the weakly coupled, non-degenerate regime, we use RPMD to probe how the onset of quantum diffraction effects at high temperature reshapes electron–electron collisions and leads to a classical-to-quantum crossover in macroscopic transport properties. Static thermodynamics obtained with RPMD are consistent with the weak-coupling equation of state, confirming correct quantum Boltzmann sampling. Real-time transport extracted from mean square displacements exhibits the expected ballistic-to-diffusive transition and a systematic reduction of the electronic self-diffusivity as quantum effects strengthen, due to quantum diffraction modifying electron–electron collisions. Direct ring-polymer scattering simulations reveal diffractive “softening” of binary deflections, providing a micro-to-macro link between collision physics and diffusion. The present study establishes RPMD as a quantitative, trajectory-based tool for electron transport across the classical–quantum crossover and furnishes benchmarks for improving Coulomb-log interpolation models. We outline extensions to multi-component plasmas and a path to incorporate Fermi–Dirac statistics within path-integral dynamics.

Electronic transport↗
DOE OSTI2024

Quantification of Dynamic Scattering Effects in Molecular Crystals using Large Angle Rocking Beam Electron Diffraction

Electron crystallography provides a pathway to solve structure of small crystals (< 1um) in size, and thus overcomes difficult synthesis constraints involved in growing large crystals. Generally, electron diffraction data is collected in the form of integrated intensity using continuous rotation or by precession. The measured intensities in 3D are utilized for structure solution. Using this approach, structure solution of difficult crystals, such as small crystals of zeolites, metal-organic frameworks, molecular crystals, and proteins , can be solved by electron diffraction. However, electron structure solutions are regularly reported with higher Rvalues than x-ray or neutron diffraction. While similar structures are found despite the high R-values, the differences in the measure intensity and theory calculated intensity limit information that can be extracted by electron diffraction. Previous work demonstrated that including multiple-scattering effects significantly reduce the R-values . Thus, it is critical to be able to quantify the dynamical diffraction effects in molecular crystals.

Crystallography↗
DOE OSTI2025

From Electronic Structure to Ion Transport: Photoelectron Spectroscopy and Molecular Dynamics Simulations Reveal the Role of Anions in Lithium Battery Electrolytes

Electrolyte anions are pivotal for lithium battery performance, yet their fundamental electronic structural properties are not well understood. In this work, we employ a combination of negative-ion photoelectron spectroscopy (NIPES), ab initio calculations, and molecular dynamics (MD) simulations to investigate the electronic structures of three representative electrolyte anions. This multiscale approach enables us to elucidate how their intrinsic electronic properties govern anion–solvent interactions in gas-phase clusters, as well as lithium-ion (Li + ) solvation structures and ion transport behavior in the condensed phase. NIPES reveals that difluoro(oxalato)borate (DFOB – ), bis(fluorosulfonyl)imide (FSI – ), and bis(oxalato)borate (BOB – ) all exhibit high electron binding energies, with vertical/adiabatic detachment energies increasing from DFOB – (6.09/5.70 eV) to FSI – (6.80/6.10 eV) to BOB – (6.82/6.40 eV), correlating with enhanced oxidation stability. Ab initio calculations reveal that DFOB – /FSI – –solvent complexes bind Li + ∼ 10 kcal/mol stronger than BOB – series, aligning with the strength of a Li + –anion model. DFOB – exhibits pronounced charge localization on both oxygen and fluorine atoms, enabling their involvement in Li + coordination. In contrast, fluorine atoms in FSI – are largely electron-depleted and remain excluded from direct Li + binding. MD simulations further demonstrate that LiDFOB and LiFSI systems exhibit Li + diffusion coefficients three and five times higher than those of LiBOB across four common solvents. Notably, LiFSI salt in acetonitrile (AN) exhibits the fastest Li + diffusion among 12 electrolyte systems, highlighting the synergistic effect of FSI – and AN in promoting ion mobility. In conclusion, these findings provide a molecular-level understanding of the critical roles of anion and its microsolvation in optimizing Li + diffusion dynamics, once again emphasizing the positioning of FSI – and DFOB – as prime candidates for next-generation electrolytes.

25 ENERGY STORAGE↗
DOE OSTI2025

Origin of Stabilization of Ligand-Centered Mixed Valence Ruthenium Azopyridine Complexes: DFT Insights for Neuromorphic Applications

Redox-driven conductance changes are critical processes in molecular- and coordination-complex-based memristive thin films and devices that are envisioned for neuromorphic technologies, but fundamental mechanisms of conductance switching are not fully understood. Here, we explore charge disproportionation (CD) processes in [Ru II L 2 ](PF 6 ) 2 molecular systems that intrinsically involve interfragment charge transfer (IFCT). Using a combination of ab initio molecular dynamics simulation (AIMD), time-dependent density functional theory (TD-DFT), and density functional theory (DFT) calculations, we investigate the electron transfer mechanisms and the roles of temperature and cell volumetric expansion in facilitating the counterion movements and electronic transitions required for low-cost IFCT and charge redistribution. A detailed analysis of the density of states and TD-DFT calculations highlights that unpaired electrons play a crucial role in low-energy transitions, with the azo (N=N) groups of the ligand serving as the primary sites for electronic transport between molecular fragments, further stabilizing the asymmetric state. Localization of added electrons on azo ligands occurs with negligible change at the Ru centers, supported by atomic volume expansions up to +4.74 bohr 3 , and goes along with a progressive reduction of the HOMO−LUMO gap across redox states, suggesting enhanced conductivity. The TD-DFT analysis reveals a dominant IFCT excitation at 2082.76 nm in the doubly reduced (22) state, while a stabilization energy of 1.20 eV of the asymmetric (13) state relative to the symmetric (22) state is predicted by constrained DFT. Periodic DFT and AIMD simulations emulating a molecular film show that the stabilization of the asymmetric state, relative to a symmetric one, translates in net charge separation values (order of ∼0.33 e) that are strongly linked to increased counterion mobility (average counterion displacements exceeding 0.7 Å per atom during CD events) and the involvement of azo groups in electron redistribution. These findings, which align with previously reported experimental and computational data, provide key insights into the IFCT mechanisms and electronic transport facilitated by azo groups, with important implications for redox-driven memristive and neuromorphic technologies.

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

Combining molecular beam epitaxy and low-energy electron microscopy with in situ magnetic susceptibility measurements within an integrated ultrahigh vacuum system

Quantum two-dimensional materials, including ultrathin superconducting films, are of great current research interest. These films are typically fabricated under ultra-high vacuum (UHV) conditions and are sensitive to the environment—prone to oxidation and contamination when exposed to the atmosphere. This hampers the study of their intrinsic properties by standard ex situ techniques. Here, we present a variable-temperature mutual inductance probe system integrated under UHV with molecular beam epitaxy (MBE) synthesis and low-energy electron microscopy, enabling nondestructive in situ characterization of superconducting thin films. The system employs a reflection-type configuration and reaches a low temperature (∼4 K) using a high-cooling-power, vibration-isolated cryocooler. In conclusion, we demonstrate the system performance by measuring the superconducting critical temperature in a copper-oxide thin film.

2D materials↗
DOE OSTI2025

Ab Initio Molecular Dynamics Study of the Reduction of Acetone by the Hydrated Electron

We have investigated the reaction dynamics of the reduction of acetone by the solvated electron in water using ab initio molecular dynamics simulations at 298 and 373 K. The rate constants derived from the simulations are consistent with experimental observations that the reaction has a low activation energy. Detailed analyses carried out to shed light on the mechanism of the reaction show that solvent reorganization plays a key role as a reaction coordinate, as expected from the Marcus electron transfer theory. Furthermore, constrained density functional theory calculations indicate that the electronic coupling is large, placing the reaction in the adiabatic limit. Indeed, the activation energies and rate constants of the simulations are in accordance with the predictions of adiabatic Marcus theory.

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

Tunable Spin Qubit Pairs in Quantum Dot–Molecule Conjugates

Organic molecules and quantum dots (QDs) have both shown promise as materials that can host quantum bits (qubits). This is in part because of their synthetic tunability. The current work employs a combination of both materials to demonstrate a series of tunable quantum dot–organic molecule conjugates that can both host photogenerated spin-based qubit pairs (SQPs) and sensitize molecular triplet states. The photogenerated qubit pairs, composed of a spin-correlated radical pair (SCRP), are particularly intriguing since they can be initialized in well-defined, nonthermally populated, quantum states. Additionally, the radical pair enables charge recombination to a polarized molecular triplet state, also in a well-defined quantum state. The materials underlying this system are an organic molecular chromophore and electron donor, 9,10-bis(phenylethynyl)anthracene, and a quantum dot acceptor composed of ZnO. We prepare a series of quantum dot–molecule conjugates that possess variable quantum dot size and two different linker lengths connecting the two moieties. Optical spectroscopy revealed that the QD–molecule conjugates undergo photoexcited charge separation to generate long-lived charge-separated radical pairs. The resulting spin states are probed using light-induced time-resolved electron paramagnetic resonance (TR-EPR) spectroscopy, revealing the presence of singlet-generated SCRPs and molecular triplet states. Notably, the EPR spectra of the radical pairs are dependent on the geometry of this highly tunable system. The g value of the ZnO QD anion is size tunable, and the line widths are influenced by radical pair separation. Overall, this work demonstrates the power of synthetic tunability in adjusting the spin specific addressability, satisfying a key requirement of functional qubit systems.

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

Current-Driven Symmetry Breaking and Spin–Orbit Polarization in Chiral Wires

The spin dynamics of electrons in chiral molecular systems remains a topic of intense interest, particularly regarding whether geometric chirality inherently induces spin polarization in current-carrying electrons. In this work, we employ ab initio real-time time-dependent density functional theory (rt-TDDFT) to directly simulate the interplay among charge current, spin, and orbital. This realtime tracking extends beyond perturbative treatments, and we analyze how nonequilibrium currents effectively lift the symmetry constraints of screw rotation and time-reversal symmetry. We find that the emergence of spin and orbital angular momenta is dynamically correlated with a concomitant loss of translational (linear) momentum, which we interpret as an intrinsic consequence of current-driven symmetry lowering. The implications of this mechanism for chirality-induced spin selectivity and spintronics device design are discussed.

chiral wire↗
DOE OSTI2024

Interpretable, extensible linear and symbolic regression models for charge density prediction using a hierarchy of many-body correlation descriptors

Here, density functional theory (DFT) is routinely used to make electronic structure predictions for high-throughput screening of materials and molecules for technologically relevant areas, like the identification of better catalysts, electronic materials, and drug discovery. However, the DFT formalism is limited by (a) its poor (quadratic-to-quartic) scaling, and (b) the need to perform repeated eigenvalue computations of the electronic Hamiltonian as part of its self-consistent field (SCF) iteration procedure to obtain the converged ground state electron density, ρ (r). Approaches that directly predict ρ (r) of a structure with high accuracy can accelerate conventional SCF calculations and can also be used in linearly scaling methods such as orbital-free DFT. To this end, we present a procedure to predict the ground state electron density of molecular and periodic three-dimensional systems directly from the atomic structure with a particular emphasis on physical interpretability. In our framework, ρ (r) is modeled using many-body correlation descriptors that accurately capture the effects of local atomic arrangements in the neighborhood of a grid point. Our use of a linear regression scheme to fit to charge density data enables transparent analysis of the relative contributions of various types of local atomic correlations. By systematically including increasingly complex correlations, our model is shown to accurately predict ρ (r) for a variety of chemically and electronically diverse systems — amorphous Ge, Al(001) slab, crystalline Ga 2 O 3 , molecular benzene, and polyethylene. We then demonstrate a symbolic regression-based protocol to construct easily computable, interpretable features from lower-order correlations that significantly improves our electron density predictions with effectively no increase in the computational cost.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Transforming an Ionic Conductor into an Electronic Conductor via Crystallization: In Situ Evolution of Transference Numbers and Structure in (La,Sr)(Ga,Fe)O 3–x Perovskite Thin Films

Mixed-conducting perovskites are workhorse electrochemically active materials, but typical high-temperature processing compromises their catalytic activity and chemo-mechanical integrity. Low-temperature pulsed laser deposition of amorphous films plus mild thermal annealing is an emerging route to form homogeneous mixed conductors with exceptional catalytic activity, but little is known about the evolution of the oxide-ion transport and transference numbers during crystallization. Here the coupled evolution of ionic and electronic transport behavior and structure in room-temperature-grown amorphous (La,Sr)(Ga,Fe)O 3-x films as they crystallize is explored. In situ ac-impedance spectroscopy with and without blocking electrodes, simultaneous capturing-synchrotron-grazing-incidence X-ray diffraction, dc polarization, transmission electron microscopy, and molecular dynamics simulations to evaluate isothermal and non-isothermal crystallization effects and the role of grain boundaries on transference numbers is combined. Ionic conductivity increases by ≈2 orders of magnitude during crystallization, with even larger increases in electronic conductivity. Consequently, as crystallinity increases, LSGF transitions from a predominantly ionic conductor to a predominantly electronic conductor. The roles of evolving lattice structural order, microstructure, and defect chemistry are examined. Grain boundaries appear relatively nonblocking electronically but significantly blocking ionically. The results demonstrate that ionic transference numbers can be tailored over a wide range by tuning crystallinity and microstructure without having to change the cation composition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electronic Structure Distortions in Chromium Chelates Impair Redox Kinetics in Flow Batteries

Aminopolycarboxylate chelates are emerging as a promising class of electrolyte materials for aqueous redox flow batteries, offering tunable redox potentials, solubility, and pH stability through careful selection of ligands and transition metal ions. Despite their potential, the impact of molecular structure modifications on the electronic and electrochemical properties of these chelates remains underexplored. Here, in this study, we examine how introducing a hydroxyl group, often employed for its solubilizing properties, to the backbone of CrPDTA, a reference chelate material, significantly changes the thermodynamics and kinetics of the chelate's redox process. We correlate changes in molecular and electronic structures to different electrochemical responses resulting from the hydroxyl addition and show that the introduction of this functional group leads to a distortion in the octahedral coordination of chromium. Furthermore, increased anisotropic spin density and nonintegral oxidation state changes in the Cr metal center result in a larger barrier for electron transfer in CrPDTA‐OH. It is demonstrated that preserving a hexacoordinate chelate structure across a broad pH range is crucial for efficient flow battery application and it is emphasized that ligand modifications must avoid distorting the octahedral coordination of the transition metal.

25 ENERGY STORAGE↗
DOE OSTI2025

Tetracene Functionalized Si(111) Achieves Enhanced Solar-to-Chemical Energy Conversion via Molecular Acceptor States

The properties of semiconductor|liquid interfaces play a critical role in determining the efficiency of solar-to-hydrogen (STH) conversion. Here, we investigate how molecular functionalization of Si(111) and Si(111)|TiO 2 surfaces impacts photoelectrochemical (PEC) hydrogen production efficiency. We find that functionalization of ∼3% of the atop sites of Si(111) with either 9-anthracene (Anth) or 5-tetracene (Tet), with the remaining sites passivated by methyl groups, provides substrates with high electronic quality and low surface oxide densities, as determined by X-ray photoelectron spectroscopy (XPS) measurements. Surface photovoltage (SPV) spectroscopy shows that surfaces modified with Anth or Tet exhibit an increased photovoltage, with Tet-functionalized surfaces yielding an additional 192 meV relative to methyl-terminated Si(111), indicating improved charge separation for Si-Tet. Further improvement in onset potential was achieved by replacing a nitrogen-containing TiO 2 atomic layer deposition (ALD) precursor (TDMAT) with a precursor lacking nitrogen (TTIP), which eliminates the parasitic defect band in the TiO 2 overlayer (p-Si(111)-Tet|TTIP-TiO 2 |Pt: V OC = +0.283 ± 0.041 V vs RHE). Density functional theory (DFT) analysis demonstrates that compared with Anth-modified Si(111), the Tet-modified surface exhibits more hybridized Si(111)-Tet states closer to the silicon band edges. Mercury contact current–voltage (I–V, dark) measurements quantified the relative interfacial density of states of Si-Tet, Si-Anth and Si-Me surfaces─revealing that the interfacial state density was highest for Si-Tet. This suggests that such hybridized interfaces serve to capture better photoexcited charge, which enables facile electron transfer to molecular acceptors in solution. Altogether, the data indicate that beneficial hybrid molecular LUMO surface states interacting with the Si conduction band edge results in improved hydrogen evolution (HER) performance for p-Si devices.

Group theory↗
DOE OSTI2025

Probing Electronic Coherence between Core-Level Vacancies at Different Atomic Sites

The detailed understanding of electronic coherence in quantum systems requires measurements on the attosecond timescale. Attosecond x-ray pulses enable the study of electronic coherence in core-excited molecular systems. Here we report on the coherent motion of electrons in the 1,1-difluoroethylene ion following ionization of the K shell of the two nonequivalent carbon sites with a subfemtosecond x-ray pulse. Using the angular streaking technique to track the Auger-Meitner decay, we observe temporal modulations of the emission, indicating the electronic coherence of the core-excited ionic states, and extract a 6.5 ± 0.8 fs average lifetime of the core-level vacancies. A quantum-mechanical model is employed to interpret the measurement, and we find the observed temporal modulations are independent of charge density oscillations. This work opens a new regime of coherent electronic motion, beyond charge migration, where electronic coherence manifests in the nonlocal quantum correlation between atomic sites while charge density oscillation is absent. Our results broaden the landscape of electronic coherence in molecular systems. Published by the American Physical Society 2025

Wang, Jun (ORCID:000000031962895X)↗
DOE OSTI2025

Inverse Trans Influence and Uranium-Arene σ-Bonding Drive Molecular Geometry: Ligand Modification from Hard to Soft Flips the Oxide from Axial to Equatorial

A rare example of an equatorially bound terminal uranium(V) oxo complex in a chelating sulfur-based ligand environment, namely [(mes( Me,Ad ArS) 3 )U V (O eq )(THF)] (2), is presented. Octahedrally coordinated 2 is obtained by reaction of the mesitylene-anchored tris-thiophenolate-coordinated uranium(III) complex [U III ((SAr Ad,Me ) 3 mes)] (1) with the oxygen-atom transfer reagent N 2 O. The observed, equatorially bound oxo ligand in 2 is in stark contrast to its known tris-aryloxide analog, [(mes( Me,Ad ArO) 3 )U V (O ax )(THF)] (A), where the oxo ligand occupies the typically observed axial coordination site. Complexes 1 and 2 are characterized by single-crystal X-ray diffraction analyses and spectroscopic and magnetochemical methods, including 1 H NMR, UV/vis/NIR electronic absorption, as well as EPR spectroscopy and SQUID magnetometry, thus confirming the C S symmetry and the pentavalent oxidation state of 2. Encompassing quantum chemical calculations (DFT and CASPT2) on 2 and its tris-phenolate analog A, support and rationalize the structural and electronic differences. The molecular orbital pictures show that a stabilizing σ-bonding interaction arising from the U–O eq inverse trans influence (ITI) is present in 2 but missing in A. In 2, the sulfur 3p orbitals are closer in energy to the uranium 5f manifold than the arene π-system, leading to an ITI, while U–arene σ- or δ-bonding is not observed. Although the arene orbitals remain separated from the uranium 5f orbitals in A, the absence of an ITI allows the arene a 2u orbital to engage in a σ-type interaction with the metal. Thus, incorporating a tris-thiophenolate to an arene anchor introduces a new design concept in molecular f-element chemistry. This approach stabilizes an equatorially bound U(V) oxo center, contrasting with its tris-phenolate counterpart, where oxo coordination is axial. The observed geometric divergence, driven by competing ITI and U–arene interactions, not only tunes electronic structure but also leads to differentiated reactivity: only the phenolate analogs activate H 2 O, while the thiolates do not.

Hydrocarbons↗
DOE OSTI2025

Simulation of 24,000 Electron Dynamics: Real-Time Time-Dependent Density Functional Theory (TDDFT) with the Real-Space Multigrids (RMG)

Here, we present the theory, implementation, and benchmarking of a real-time time-dependent density functional theory (RT-TDDFT) module within the RMG code, designed to simulate the electronic response of molecular systems to external perturbations. Our method offers insights into nonequilibrium dynamics and excited states across a diverse range of systems, from small organic molecules to large metallic nanoparticles. Benchmarking results demonstrate excellent agreement with established TDDFT implementations and showcase the superior stability of our time integration algorithm, enabling long-term simulations with minimal energy drift. The scalability and efficiency of RMG on massively parallel architectures allow for simulations of complex systems, such as plasmonic nanoparticles with thousands of atoms. Future extensions, including nuclear and spin dynamics, will broaden the applicability of this RT-TDDFT implementation, providing a powerful toolset for studies of photoactive materials, nanoscale devices, and other systems where real-time electronic dynamics is essential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗