Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio molecular dynamics”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Application-specific machine-learned interatomic potentials: exploring the trade-off between DFT convergence, MLIP expressivity, and computational cost

Machine-learned interatomic potentials (MLIPs) are revolutionizing computational materials science and chemistry by offering an efficient alternative to ab initio molecular dynamics (MD) simulations. However, fitting high-quality MLIPs remains a challenging, time-consuming, and computationally intensive task where numerous trade-offs have to be considered, e.g., How much and what kind of atomic configurations should be included in the training set? Which level of ab initio convergence should be used to generate the training set? Which loss function should be used for fitting the MLIP? Which machine learning architecture should be used to train the MLIP? The answers to these questions significantly impact both the computational cost of MLIP training and the accuracy and computational cost of subsequent MLIP MD simulations. In this study, we use a configurationally diverse beryllium dataset and quadratic spectral neighbor analysis potential. We demonstrate that joint optimization of energy versus force weights, training set selection strategies, and convergence settings of the ab initio reference simulations, as well as model complexity can lead to a significant reduction in the overall computational cost associated with training and evaluating MLIPs. This opens the door to computationally efficient generation of high-quality MLIPs for a range of applications which demand different accuracy versus training and evaluation cost trade-offs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine Learning-Accelerated First-Principles Molecular Dynamics Explains Anomalous Lattice Thermal Expansion in BaZr 0.78 Y 0.22 O 3-δ

Fuel cells are a vital clean energy technology that converts chemical energy directly into electricity with high efficiency, making them a cornerstone of a sustainable energy future. Herein we investigate the thermal and chemical lattice expansion behavior of hydrated BaZr 0.78 Y 0.22 O 3-δ using machine learning-accelerated ab initio molecular dynamics simulations. Here, our results reproduce the experimentally observed non-monotonic and anomalous temperature dependence of lattice expansion, which we attribute to the competing effects of thermal expansion and dehydration—two mechanisms that influence the lattice expansion in opposite directions. The importance of this work lies in its detailed demonstration of how advanced computational techniques can accurately capture complex environmental effects, providing a valuable framework for modeling similar phenomena in a variety of material systems and applications.

Proton conducting fuel cell↗
DOE OSTI2024

Exploring Li-Ion Transport Properties of Li 3 TiCl 6 : A Machine Learning Molecular Dynamics Study

We performed large-scale molecular dynamics simulations based on a machine-learning force field (MLFF) to investigate the Li-ion transport mechanism in cation-disordered Li 3 TiCl 6 cathode at six different temperatures, ranging from 25°C to 100°C. In this work, deep neural network method and data generated by ab − initio molecular dynamics (AIMD) simulations were deployed to build a high-fidelity MLFF. Radial distribution functions, Li-ion mean square displacements (MSD), diffusion coefficients, ionic conductivity, activation energy, and crystallographic direction-dependent migration barriers were calculated and compared with corresponding AIMD and experimental data to benchmark the accuracy of the MLFF. From MSD analysis, we captured both the self and distinct parts of Li-ion dynamics. The latter reveals that the Li-ions are involved in anti-correlation motion that was rarely reported for solid-state materials. Similarly, the self and distinct parts of Li-ion dynamics were used to determine Haven’s ratio to describe the Li-ion transport mechanism in Li 3 TiCl 6 . Obtained trajectory from molecular dynamics infers that the Li-ion transportation is mainly through interstitial hopping which was confirmed by intra- and inter-layer Li-ion displacement with respect to simulation time. Ionic conductivity (1.06 mS/cm) and activation energy (0.29eV) calculated by our simulation are highly comparable with that of experimental values. Overall, the combination of machine-learning methods and AIMD simulations explains the intricate electrochemical properties of the Li 3 TiCl 6 cathode with remarkably reduced computational time. Thus, our work strongly suggests that the deep neural network-based MLFF could be a promising method for large-scale complex materials.

Selvaraj, Selva Chandrasekaran (ORCID:000000029023↗
NASA NTRS2016

Ab Initio Simulations of Phase Stability and Martensitic Transitions in NiTi

For NiTi based alloys, the shape memory effect is governed by a transition from a low-temperature martensite phase to a high-temperature austenite phase. Despite considerable experimental and computational work, basic questions regarding the stability of the phases and the martensitic phase transition remain unclear even for the simple case of binary, equiatomic NiTi. We perform ab initio molecular dynamics simulations to describe the temperature-dependent behavior of NiTi and resolve several of these outstanding issues. Structural correlation functions and finite temperature phonon spectra are evaluated to determine phase stability. We show that finite temperature, entropic effects stabilize the experimentally observed martensite (B19') and austenite (B2) phases while destabilizing the theoretically predicted (B33) phase. Free energy computations based on ab initio thermodynamic integration confirm these results and permit estimates of the transition temperature between the phases. In addition to the martensitic phase transition, we predict a new transition between the B33 and B19' phases. The role of defects in suppressing phase transformation temperatures is discussed.

Simulations↗
NASA NTRS2016

Ab Initio Simulations of Temperature Dependent Phase Stability and Martensitic Transitions in NiTi

For NiTi based alloys, the shape memory effect is governed by a transition from a low-temperature martensite phase to a high-temperature austenite phase. Despite considerable experimental and computational work, basic questions regarding the stability of the phases and the martensitic phase transition remain unclear even for the simple case of binary, equiatomic NiTi. We perform ab initio molecular dynamics simulations to describe the temperature-dependent behavior of NiTi and resolve several of these outstanding issues. Structural correlation functions and finite temperature phonon spectra are evaluated to determine phase stability. In particular, we show that finite temperature, entropic effects stabilize the experimentally observed martensite (B19') and austenite (B2) phases while destabilizing the theoretically predicted (B33) phase. Free energy computations based on ab initio thermodynamic integration confirm these results and permit estimates of the transition temperature between the phases. In addition to the martensitic phase transition, we predict a new transition between the B33 and B19' phases. The role of defects in suppressing these phase transformations is discussed.

nickel titanium↗
DOE OSTI2025

Local structure of zinc–indium–tin oxide films via grazing-incidence x-ray pair-distribution functions and theoretical methods

A detailed experimental and theoretical study on the local (r ≤ 4.5 Å) atomic structure of amorphous and crystalline zinc–indium–tin oxide (ZITO) thin films using grazing-incidence x-ray Pair-Distribution Functions (PDFs), ab initio Molecular Dynamics (MD), and Empirical Potential Structure Refinement (EPSR) Monte Carlo simulations is presented. High-energy synchrotron x rays, a two-dimensional detector, and different incident angles were used to probe the depth uniformity of five (ZnO) 0.15 (In 2 O 3 ) 0.70 (SnO 2 ) 0.15 films that were deposited via pulsed-laser deposition at growth temperatures (T G ) ranging from 25 to 300 °C. Films deposited at T G ≤ 150 °C were amorphous. The partially crystalline (T G = 200 °C) and fully crystalline (T G = 300 °C) films were highly textured. Both crystalline and amorphous structures were investigated using ab initio MD and EPSR Monte Carlo simulations. The density of the amorphous films determined from the experimental data agreed with MD calculations. Coordination numbers, bond lengths, and distortion for metal–oxygen and for both the edge- and corner-shared In–metal shells up to 4.5 Å obtained from PDF analysis closely agreed with MD and EPSR simulations. There is a pronounced decrease in the edge- and corner-shared In–Zn distances arising from the shorter Zn–O bond length, Zn–O tetrahedral coordination, and In–O–Zn angle in amorphous ZITO compared to its crystalline counterpart. A maximum in electrical mobility was observed for the amorphous film just before crystallization occurred. While the peak is broad, consistent with nearly unchanged overall cation–oxygen coordination in the amorphous films, ESPR results indicate that the tetrahedral coordination follows the conductivity trend.

Grazing Incidence X-ray↗
DOE OSTI2026

Exploring the Origins of Anti-Ambipolarity in BBL Polymer: Links to Redox Chemistry, Electronic Structure, and Structural Dynamics

We examine the intrinsic physical-chemical properties of the conjugated ladder-type polymer poly(benzimidazobenzophenanthroline) (BBL) in response to electron transfer. We aim at explaining the origin of the anti-ambipolar behavior behind the observed BBL nonlinear response associated with specific device architectures. To elucidate this point, we use theory and computation based on first principles, including density functional theory optimizations, ab initio molecular dynamics, time-dependent DFT, and Marcus-theory analysis. Our results reveal that this redox response is not simply monotonic but follows an alternating odd/even pattern in which gap narrowing and reopening occur sequentially before near-gapless behavior emerges at high charging. Converging theoretical evidence in this work demonstrates that bell shaped conductivity in BBL originates in its fundamental electronic structure and supramolecular organization.

FOS: Physical sciences↗
DOE OSTI2026

Enhanced Catalytic Dechlorination of Polyvinyl Chloride (PVC) and H2 Production Enabled by Synergistic Gaδ+/Ga0 Active Sites in Liquid Metal Particles

Polyvinyl chloride (PVC) is ubiquitous yet challenging to recycle due to its tendency to thermally decompose above 250 °C, releasing toxic, corrosive chlorinated compounds, and its inability to melt. Here, we report a catalytic strategy for PVC upcycling at 160 °C using gallium liquid metal particles (Ga-LMP) featuring a dynamic Ga-GaOOH core–shell architecture. These catalysts enable concurrent dechlorination and hydrogen evolution, yielding up to 7% H2 (based on initial hydrogen atoms in PVC) along with a highly dechlorinated (>95%) carbonaceous solid and aqueous HCl. Mechanistic investigations combining X-ray photoelectron spectroscopy, infrared spectroscopy, solid-state NMR, inelastic neutron scattering, and ab initio molecular dynamics reveal a synergistic interplay between Gaδ+ sites in the GaOOH shell and metallic Ga0 in the core. Cationic Ga initiates C–Cl bond activation and HCl formation, while progressive reduction of the shell exposes Ga0 sites that promote C–H activation and H2 evolution. Control experiments with a Ga salt and bulk Ga liquid metal confirmed that neither oxidation state alone can achieve both transformations efficiently. This work establishes a dynamic dual-site paradigm for liquid metal catalysis, in which the in situ evolution and coexistence of oxidized and metallic species enable sequential and cooperative bond activation pathways. These findings provide a general design principle for novel liquid metal catalysts that target challenging polymer transformations under mild conditions.

Zingg, Benjamin [ORNL] (ORCID:0009000914530153)↗
DOE OSTI2024

Deep-Learning Interatomic Potential Connects Molecular Structural Ordering to the Macroscale Properties of Polyacrylonitrile

Polyacrylonitrile (PAN) is an important commercial polymer, bearing atactic stereochemistry resulting from nonselective radical polymerization. As such, an accurate, fundamental understanding of governing interactions among PAN molecular units is indispensable for advancing the design principles of final products at reduced processability costs. While ab initio molecular dynamics (AIMD) simulations can provide the necessary accuracy for treating key interactions in polar polymers, such as dipole–dipole interactions and hydrogen bonding, and analyzing their influence on the molecular orientation, their implementation is limited to small molecules only. Herein, we show that the neural network interatomic potentials (NNIPs) that are trained on the small-scale AIMD data (acquired for oligomers) can be efficiently employed to examine the structures and properties at large scales (polymers). NNIP provides critical insight into intra- and interchain hydrogen-bonding and dipolar correlations and accurately predicts the amorphous bulk PAN structure validated by modeling the experimental X-ray structure factor. Furthermore, the NNIP-predicted PAN properties, such as density and elastic modulus, are in good agreement with their experimental values. Overall, the trend in the elastic modulus is found to correlate strongly with the PAN structural orientations encoded in the Hermans orientation factor. In conclusion, this study enables the ability to predict the structure–property relations for PAN and analogues with sustainable ab initio accuracy across scales.

36 MATERIALS SCIENCE↗
DOE OSTI2024

UV-Induced Reaction Pathways in Bromoform Probed with Ultrafast Electron Diffraction

For many chemical reactions, it remains notoriously difficult to predict and experimentally determine the rates and branching ratios between different reaction channels. This is particularly the case for reactions involving short-lived intermediates, whose observation requires ultrafast methods. The UV photochemistry of bromoform (CHBr 3 ) is among the most intensely studied photoreactions. Yet, a detailed understanding of the chemical pathways leading to the production of atomic Br and molecular Br 2 fragments has proven challenging. In particular, the role of isomerization and/or roaming and their competition with direct C–Br bond scission has been a matter of continued debate. Here, in this work, gas-phase ultrafast megaelectronvolt electron diffraction (MeV-UED) is used to directly study structural dynamics in bromoform after single 267 nm photon excitation with femtosecond temporal resolution. The results show unambiguously that isomerization contributes significantly to the early stages of the UV photochemistry of bromoform. In addition to direct C–Br bond breaking within <200 fs, formation of iso-CHBr 3 (Br-CH-Br-Br) is observed on the same time scale and with an isomer lifetime of >1.1 ps. The branching ratio between direct dissociation and isomerization is determined to be 0.4 ± 0.2:0.6 ± 0.2, i.e., approximately 60% of molecules undergo isomerization within the first few hundred femtoseconds after UV excitation. The structure and time of formation of iso-CHBr 3 compare favorably with the results of an ab initio molecular dynamics simulation. The lifetime and interatomic distances of the isomer are consistent with the involvement of a roaming reaction mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Predicting Melt Properties Using Atomistic Simulations With A Highly Accurate Physically Informed Neural Network Interatomic Potential

The use of a recently developed machine learning (ML) interatomic potential for molecular dynamics simulations of aluminum melt properties will be presented. Such properties are critical for process modeling in additive manufacturing, including the melt pool size, solidification, and formation of solidification microstructures. Direct first-principles modeling of these processes is computationally prohibitive whereas simulations employing ML potentials combine the high accuracy of quantum-mechanical methods with high computational speeds. The physically-informed neural network (PINN) method used herein, integrates a high-dimensional regression implemented by an artificial neural network with a physics-based bond-order interatomic potential. PINN potentials can accurately reproduce many properties of aluminum in both crystalline-solid and liquid phases. We examine the accuracy of a PINN Al potential in predicting the density, self-diffusivity, viscosity, and the tension of the liquid surface and liquid-solid interfaces. Comparison with experimental data and ab initio molecular dynamics calculations shows very good agreement for all properties tested.

molecular dynamics↗
DOE OSTI2025

Structural and Dynamic Heterogeneity of Deep Eutectic Solvents Composed of Choline Chloride and Ortho-Phenol Derivatives

Structural, thermal, and dynamic properties of four deep eutectic solvents comprising choline chloride paired with ortho-phenolic derivative hydrogen-bond donors were probed using experiments and molecular simulations. The hydrogen-bond donors include phenol, catechol, o-chlorophenol, and o-cresol, in a 3:1 mixture with the hydrogen-bond acceptor choline chloride. Density, viscosity, and pulsed-field gradient NMR diffusivity measurements were conducted over a range of temperatures. Classical and ab initio molecular dynamics simulation results match experimental data reasonably well. Furthermore, the simulation results were then used to perform a more detailed analysis of the local structure and dynamics of these systems.

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

In-Situ Atomic-Scale Revelation of Amorphous Metallic Iron Formation during Hydrogen-Driven Reduction of Iron Oxides

The transition to hydrogen as a green reductant in metal production is critical for decarbonizing the metallurgical industry, yet atomic-scale mechanisms governing reduction pathways and phase evolution remain unresolved. Using in-situ environmental transmission electron microscopy, we identify a hidden pathway that reveals dynamic formation of amorphous metallic iron (Fe) during the hydrogen-driven reduction of ferrous oxides of Fe 3 O 4 and FeO. Real-time imaging uncovers three coexisting transformation routes: (i) Fe 3 O 4 → FeO, (ii) Fe 3 O 4 → amorphous Fe, and (iii) FeO → amorphous Fe. The resulting amorphous Fe exhibits fluid-like mobility, enabling its rapid aggregation and crystallization into core-shell nanostructures, with a crystalline core enveloped by an amorphous shell. Complementary ab initio molecular dynamics simulations trace the amorphous Fe formation to interfacial strain at the metal/oxide interfaces, where large lattice mismatches destabilize the metal lattice during initial metallization. This interplay between thermodynamics and kinetics governs phase evolution: thermodynamics favors a self-limiting amorphous Fe overlayer, while rapid oxide reduction kinetics drives amorphous overgrowth. Our findings demonstrate that amorphous intermediates bypass rate-limiting crystalline steps, providing mechanistic insights to optimize H 2 -based processes for sustainable steelmaking. In conclusion, these insights bridge the gap between macroscopic process engineering and atomic-scale dynamics, with broader implications for catalysis and nanostructured material synthesis, where oxide reduction pathways critically shape functional phases and microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Be₂C as a neutron moderator for molten salt reactors: A first-principles study of structural, electronic, and defect properties

Beryllium carbide (Be 2 C), valued for its high neutron moderation efficiency and low absorption cross section, is a promising high-temperature neutron moderator for molten salt reactors. Its practical adoption, however, demands significant technological maturation, requiring comprehensive theoretical and experimental studies of its response to different conditions, including high temperature and irradiation. Here, we report initial results on the fundamental properties and radiation-induced defects of Be 2 C, focusing on antisites, vacancies, interstitial atoms, and Frenkel pairs in the Be and C sublattices. Using density functional theory (DFT) and ab initio molecular dynamics (AIMD), we calculate the defects formation and binding energies, evaluating their dependence on the supercell size, charge states, and chemical environment. In general, carbon defects exhibit higher formation energies, greater sensitivity to cell size, and stronger impacts on the density of states compared to beryllium defects, with charged state the effects being more pronounced. Static DFT reveals multiple metastable interstitial configurations, while AIMD identifies ground states as C-C <100> dumbbells and octahedral Be interstitials. In conclusion, the diversity of metastable configurations and defect states complicates the diffusion mechanisms, requiring further molecular dynamics analysis to elucidate the mechanisms and rates of radiation-induced atomic transport, as well as the structural stability of Be2C.

Ab initio modeling↗
DOE OSTI2024

The Role of Nonequilibrium Solvent Effects in Enhancing Direct CO 2 Capture at the Air–Aqueous Amino Acid Interface

Direct air capture (DAC) technologies are limited by the poor understanding of the dynamic role of interfaces in modulating the chemisorption of CO 2 from air into solutions. While the reactivity of aqueous amine-based solvents in the bulk environment is strongly inhibited by nonequilibrium solvent effects, promoting DAC at interfaces posits a possibility to reduce the coupling with the solvent and significantly accelerate DAC. Building on an experimentally proven concept to bring an anionic glycine absorbent to the interface through ion-pairing interactions with a positively charged surfactant, we establish the fundamental time scales for key elementary steps involved in DAC with rate theory and enhanced-sampling ab initio molecular dynamics simulations. Here, we elucidate the mechanism by which water influences the free energy barriers and dynamical crossing-recrossing of those barriers, affecting the reaction rates. We find that water reorganizes to partially dehydrate [-NH 2 ], facilitating S N 2-based CO 2 conversion to a zwitterion, which then releases a proton via overhydration of [-NH 2 ]. The low-density interfacial water favors dehydration over overhydration, leading to a comparatively higher barrier (slower kinetics) for proton release. The barrier-recrossing events neutralize this effect, letting both steps occur at the same time scale (sub-microseconds) and making the overall DAC process faster at the interface than in the bulk water. Such an understanding of environment-sensitive solvent effects on the reaction kinetics will help design tailored interfaces for enhanced CO 2 capture kinetics via control of solvation and ion paring.

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

Axial H-Bonding Solvent Controls Inhomogeneous Spectral Broadening, While Peripheral H-Bonding Solvent Controls Vibronic Broadening: Cresyl Violet in Methanol

The dynamics of the nuclei of both a chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore–solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore–solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both force field and ab initio molecular dynamics trajectories along with the inclusion of only certain solvent molecules in the excited-state calculations, we determine that the methanol molecules axial to the chromophore are responsible for the majority of inhomogeneous broadening, with a single methanol molecule that forms an axial hydrogen bond dominating the response. Furthermore, the strong peripheral hydrogen bonds do not contribute to spectral broadening, as they are very stable throughout the dynamics and do not lead to increased energy-gap fluctuations. We also find that treating the strong peripheral hydrogen bonds as molecular mechanical point charges during the molecular dynamics simulation underestimates the vibronic coupling. Including these peripheral hydrogen bonding methanol molecules in the quantum-mechanical region in a geometry optimization increases the vibronic coupling, suggesting that a more advanced treatment of these strongly interacting solvent molecules during the molecular dynamics trajectory may be necessary to capture the full vibronic spectral broadening.

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

Prediction of L⁢i 3 ⁢F⁡e 8 ⁢B 8 compound with rapid one-dimensional ion diffusion channels

Using a computational crystal structure search in the Li-Fe-B ternary system, we predict a stable phase of L⁢i 3 ⁢F⁡e 8 ⁢B 8 , featuring 1D channels that enable rapid Li-ion transport. Ab initio molecular dynamics simulations show that the Li-ion diffusion coefficient in L⁢i 3 ⁢F⁡e 8⁢ B 8 surpasses that of common electrode and conductive additive materials by several orders of magnitude. The high diffusion in L⁢i 3 ⁢F⁡e 8 ⁢B 8 can be explained by the Frenkel–Kontorova model, which describes an incommensurate state between the Li diffusion chain and the periodic potential field caused by the FeB backbone structure. The favorable lithium-ion diffusivity and mechanical properties of L⁢i 3 ⁢F⁡e 8 ⁢B 8 make it a promising conductive additive for battery materials. Furthermore, an external magnetic field can further manipulate the properties of this material due to its predicted itinerant ferromagnetism, which also offers a platform for exploring spin-dependent phenomena.

1-dimensional systems↗
DOE OSTI2025

Electronic Structure of the Au Nanoparticle-TiO 2 Heterojunction: Influence of Nanoparticle Size, Shape, Oxygen Vacancies, and Temperature

The electronic structure of the gold nanoparticle-titanium dioxide (AuNP-TiO 2 ) heterojunction plays a critical role for charge transfer and recombination dynamics that underpin its photocatalytic function. However, building a representative model to capture the key physics remains a significant challenge. Here, we investigate the influence of the AuNP size and shape, as well as oxygen vacancy (V O ) defects at the anatase-phase TiO 2 (101) surface and the temperature of the heterojunction, on its interfacial electronic properties. Using density functional theory (DFT), we compare the closed-shell Au 20 and open-shell Au 19 clusters interfaced with pristine and V O defect TiO 2 surfaces. We find that the presence of a V O defect transforms pure TiO 2 from a p-type to an n-type semiconductor, reversing the interfacial band bending from downward to upward. For the heterosystem, density of states (DOS) analysis shows that V O minimally affects the open-shell Au 19 -TiO 2 system, but it significantly alters the closed-shell Au 20 -TiO 2 system, converting it from p-type to n-type at the Γ-point at 0 K. Furthermore, ab initio molecular dynamics (AIMD) simulations at 300 K reveal significant thermal fluctuations in AuNP positions relative to the TiO 2 surface. These fluctuations result in dynamic variations in the gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) in all systems studied: Au 19 NP-Pristine TiO 2 , Au 20 NP-Pristine TiO 2 , Au 19 NP-V O defect TiO 2 , and Au 20 NP-V O defect TiO 2 . Complementary Bader charge analysis performed at both 0 K and finite-temperature AIMD snapshots supports the emergence of an upward band bending and the formation of a Schottky barrier at the V O containing heterojunctions. Notably, we find that at finite temperature, an Au atom can dynamically passivate the V O , leading to a pronounced widening of the HOMO-LUMO gap in the AuNP-V O defect TiO 2 heterojunctions. Our computational findings underscore the pivotal role of the V O defect and thermal effects in modulating interfacial band alignment, electronic states, and HOMO-LUMO gaps, providing insights for designing Au-TiO 2 heterojunctions with tailored electronic properties.

Mora Perez, Carlos [Lawrence Berkeley National Lab↗