Density Functional Tight Binding Research Articles

Effects of Crystallinity and Branched Chain on Thermal Degradation of Polyethylene: A SCC-DFTB Molecular Dynamics Study

As a widely used plastic, the aging and degradation of polyethylene (PE) are inevitable problems, whether the goal is to prolong the life of PE products or address the issue of white pollution. Molecular simulation is a vital scientific tool in elucidating the mechanisms and processes of chemical reactions. To obtain the distribution and evolution process of PE’s thermal oxidation products, this work employs the self-consistent charge–density functional tight binding (SCC-DFTB) method to perform molecular simulations of the thermal oxidation of PE with different crystallinity and branched structures. We discovered that crystallinity does not affect the thermal oxidation mechanism of PE, but higher crystallinity makes PE more susceptible to cross-linking and carbon chain growth, reducing the degree of PE carbon chain breakage. The branched structure of PE results in differences in free volumes between the carbon chains, with larger pores leading to a concentrated distribution of O2 and chemical defects subsequently formed. The breakdown of PE is slowed down when chemical defects are localized in low-density regions of the carbon chain. The specifics and mechanism of PE’s thermal oxidation are clearly revealed in this paper, which is essential for understanding the process in depth and for the development of anti-aging PE products.

Quantum Simulations of Radiation Damage in a Molecular Polyethylene Analog.

An atomic-level understanding of radiation-induced damage in simple polymers like polyethylene is essential for determining how these chemical changes can alter the physical and mechanical properties of important technological materials such as plastics. Ensembles of quantum simulations of radiation damage in a polyethylene analog are performed using the Density Functional Tight Binding method to help bind its radiolysis and subsequent degradation as a function of radiation dose. Chemical degradation products are categorized with a graph theory approach, and occurrence rates of unsaturated carbon bond formation, crosslinking, cycle formation, chain scission reactions, and out-gassing products are computed. Statistical correlations between product pairs show significant correlations between chain scission reactions, unsaturated carbon bond formation, and out-gassing products, though these correlations decrease with increasing atom recoil energy. The results present relatively simple chemical descriptors as possible indications of network rearrangements in the middle range of excitation energies. Ultimately, the work provides a computational framework for determining the coupling between nonequilibrium chemistry in polymers and potential changes to macro-scale properties that can aid in the interpretation of future radiation damage experiments on plasticmaterials.

Atomically Precise Control of Topological State Hybridization in Conjugated Polymers.

Realization of topological quantum states in carbon nanostructures has recently emerged as a promising platform for hosting highly coherent and controllable quantum dot spin qubits. However, their adjustable manipulation remains elusive. Here, we report the atomically accurate control of the hybridization level of topologically protected quantum edge states emerging from topological interfaces in bottom-up-fabricated π-conjugated polymers. Our investigation employed a combination of low-temperature scanning tunneling microscopy and spectroscopy, along with high-resolution atomic force microscopy, to effectively modify the hybridization level of neighboring edge states by the selective dehydrogenation reaction of molecular units in a pentacene-based polymer and demonstrate their reversible character. Density functional theory, tight binding, and complete active space calculations for the Hubbard model were employed to support our findings, revealing that the extent of orbital overlap between the topological edge states can be finely tuned based on the geometry and electronic bandgap of the interconnecting region. These results demonstrate the utility of topological edge states as components for designing complex quantum arrangements for advanced electronic devices.

Switching of electrochemical selectivity due to plasmonic field-induced dissociation

Electrochemical reactivity is known to be dictated by the structure and composition of the electrocatalyst-electrolyte interface. Here, we show that optically generated electric fields at this interface can influence electrochemical reactivity insofar as to completely switch reaction selectivity. We study an electrocatalyst composed of gold-copper alloy nanoparticles known to be active toward the reduction of CO2 to CO. However, under the action of highly localized electric fields generated by plasmonic excitation of the gold-copper alloy nanoparticles, water splitting becomes favored at the expense of CO2 reduction. Real-time time-dependent density functional tight binding calculations indicate that optically generated electric fields promote transient-hole-transfer-driven dissociation of the O─H bond of water preferentially over transient-electron-driven dissociation of the C─O bond of CO2. These results highlight the potential of optically generated electric fields for modulating pathways, switching reactivity on/off, and even directing outcomes.

How does the integration of amino acids enhance the bonding ability of triazine-based inhibitors with iron surfaces: Insights from SCC-DFTB, COSMO-RS and molecular dynamics simulations

Theoretical investigations into the adsorption of molecules on metal surfaces play a pivotal role in the development of more efficient corrosion inhibitors. This study offers a novel approach by combining Density Functional Tight Binding (DFTB), Molecular Dynamics (MD) simulations, and conductor-like screening model for realistic solvation (COSMO-RS) calculations to comprehensively assess the interaction and diffusion properties of glycine (Gly), 2,2′-azanediyldiacetic acid (IDA), and 5-aminopentanoic acid (5-APA) with two 1,3,5-triazine (Tris) derivatives (Tris-Gly and Tris-IDA) for corrosion protection of carbon steel in acidic medium. The solvation properties of these molecules were evaluated alongside their interaction strength with the Fe(1 1 0) surface, considering both neutral and protonated forms. DFTB simulations revealed a clear trend in interaction energies, following the order: Tris-IDA > Tris-Gly > 5-APA > IDA > Gly. It ranges from −1.191 eV to −1.853 eV, with the highest stability observed for protonated Tris-IDA (−1.853 eV), while neutral Glycine exhibits the lowest interaction energy (−1.263 eV). While these theoretical results diverged slightly from experimental observations, with Tris-Gly outperforming Tris-IDA, this discrepancy was attributed to steric effects and solvent interactions. The DFTB results also indicated strong covalent interactions, particularly involving acetic groups, between the inhibitors’ orbitals and the iron’s 3d orbitals. MD simulations further demonstrated the impact of protonation on inhibitor mobility, with reduced diffusion coefficients due to enhanced electrostatic interactions and hydrogen bonding. COSMO-RS solvation analysis confirmed the strong hydrogen bonding capabilities of these molecules with water. Overall, this study elucidates the mechanisms of corrosion inhibition by integrating multiple simulation techniques, providing key molecular insights that can guide the design of more effective corrosion inhibitors. The findings emphasize the significance of electronic interactions and molecular structure in enhancing inhibitor performance, making this a valuable contribution to corrosion science.

Imidazole derivative for corrosion protection: A comprehensive theoretical and experimental study on mild steel in acidic environments

The inhibition efficiency derivative imidazole, 2-mercapto-5,5-diphenyl-3,5-dihydro-4 H-imidazol-4-one (AM0), against corrosion in a hydrochloric acid (HCl) environment was studied. The research focused on the compound's ability to adsorb onto mild steel surfaces. Various experimental techniques, including potentiodynamic polarization (PDP), electrochemical frequency modulation (EFM), and impedance spectroscopy (EIS), were employed to analyze the inhibitory effects. Results showed that the effectiveness of AM0 increased with concentration, reaching an inhibitory efficiency of 95.1 %. Polarization studies revealed that the inhibitor displayed mixed-type behavior. The adsorption process was modeled using the Langmuir adsorption model. Surface characterization techniques such as scanning electron microscopy (SEM), energy dispersive X-ray analysis (EDS), and X-ray photoelectron spectroscopy (XPS) were utilized to examine the inhibited corrosion surfaces. Theoretical calculations based on density functional theory (DFT) and density functional tight binding (DFTB) suggested that inhibitors with electron-accepting capabilities exhibited strong interactions with the iron surface, shedding light on specific bond formations that contribute to the stability and interaction dynamics at the molecular interface. This study provides valuable information on the corrosion inhibition potential of the new imidazole compound, offering important pointers for the development of effective corrosion inhibitors in acidic environments.

Implications for electronic structures of S-doped graphene nanoribbons from a DFTB algorithm at atomic scale

Self-consistent charge density functional tight binding simulations were used to provide implications for the effect of S atoms’ doping on geometrical structures and electronic states in armchair graphene nanoribbons with different widths. The geometric configurations, energy, charge density differences, Mulliken charges, and molecular orbitals at HOMO and LUMO energy levels are characterized. Besides the changes of bond length and bond angles, the calculation results indicate that the band structures and bandgaps of the graphene nanoribbons is affected by the nanoribbon width as well as S doping positions, where the ribbons exhibit metallic or semiconductor properties. Doping S atoms result significant changes in the charge density differences and Mulliken charges on the atoms near the doped atom. At the same time, the HOMOs and LUMOs also present differences.

Impact of heteroatoms and chemical functionalisation on crystal structure and carrier mobility of organic semiconductors

The substitution of heteroatoms and the functionalisation of molecules are established strategies in chemical synthesis. They target the precise tuning of the electronic properties of hydrocarbon molecules to improve their performance in various applications and increase their versatility. Modifications to the molecular structure often lead to simultaneous changes in the morphology such as different crystal structures. These changes can have a stronger and unpredictable impact on the targeted property. The complex relationships between substitution/functionalization in chemical synthesis and the resulting modifications of properties in thin films or crystals are difficult to predict and remain elusive. Here we address these effects for charge carrier transport in organic crystals by combining simulations of carrier mobilities with crystal structure prediction based on density functional theory and density functional tight binding theory. This enables the prediction of carrier mobilities based solely on the molecular structure and allows for the investigation of chemical modifications prior to synthesis and characterisation. Studying nine specific molecules with tetracene and rubrene as reference compounds along with their combined modifications of the molecular cores and additional functionalisations, we unveil systematic trends for the carrier mobilities of their polymorphs. The positive effect of phenyl groups that is responsible for the marked differences between tetracene and rubrene can be transferred to other small molecules such as NDT and NBT leading to a mobility increase by large factors of about five.

Molecular origin of viscoelasticity and influence of methylation in mesophase pitch

The viscoelastic and thermomechanical properties of pitches are responsible for their melt-spinning behavior, a critical step for manufacturing high-performance pitch-based carbon fibers. Here, we systematically explore the impact of methyl group modifications on the viscoelastic and thermal properties of mesophase pitches. We employ a range of atomistic modeling approaches, including Density Functional Theory (DFT), Density Functional Tight Binding (DFTB), and Classical Molecular Mechanics (MM), to provide detailed insights into the molecular interactions and structural changes. Our results revealed the molecular mechanisms that promote layered structures leading to the anisotropic nature of the viscoelastic behavior of mesophase pitch. Furthermore, we propose a modified molecular representation of naphthalene-based mesophase pitch based on the analysis of x-ray diffraction measurements. This study provides fundamental insights into the molecular structures of mesophase pitch and the role of methyl groups controlling its viscosity, which offer valuable insights into mesophase-based carbon fiber production.

A Mechanism of Resistance Switching in CNT Based Memory Devices

Carbon nanotubes (CNTs) show promise for high energy efficiency electronic devices due to their high electrical conductivity. One possible application is in the fabrication of resistance-change computer memory, where the application of an external electrical bias can use a CNT layer as the switching element in a RAM cell. Memory cells of this type have been extensively demonstrated, where the electrical resistance of a switching layer of CNTs can be reversibly controlled by the application of external bias [1].To better understand these devices, we have developed a nanoscale model of their electrical switching. Towards this end, first-principles calculations of the electrical conductivity of CNT-CNT junctions were performed using the Density Functional Tight Binding (DFTB) in combination with the Non-Equilibrium Greens functions (NEGF) approach, enabling an improved understanding of the tunnelling of electrons across a range of CNT junctions [2]. These results demonstrate that the conductivity of CNT contacts depends on the overlap area between nanotubes and exponentially on the distances between the carbon atoms of the interacting CNTs. Secondly, we employed coarse-grained molecular dynamics [3,4] to produce nanoscale dense CNT film models like those used in the devices. We devise a set of metrics for characterising the structural properties of CNT films, and apply them to a set of model films designed to mimic those used in microelectronics devices. In order to understand how the geometrical properties of CNTs affect the film structure, we analyse films with different CNT chirality and length distributions. Unlike previous studies that focused on relatively low-density systems below ∼0.4 g/cm3, we also consider films with density of ∼0.6 g/cm3, as employed in experimentally fabricated CNT-based NRAM devices. By adding amorphous carbon to the film models, in varying amounts, we are then able to assess its effects. Finally, these two elements were combined to produce electrical simulations of the conductivity of the CNT films connected to TiN electrodes. These electrodes are shown to be partially oxidized during the device fabrication. The electrical conductivity in a device is characterized by formation of conduction paths formed by CNTs between the electrodes and connected via the CNT junctions. The conductivity is modulated by charging of CNTs and amorphous carbon present in the film.Using these simulations, we propose a resistance switching mechanism in TiN/CNT/TiN devices based on the physical movement of CNTs. During the initial operation of the device, oxygen is released from the oxidized bottom electrode and this readily reacts with individual CNTs to form charge trapping sites. Once these sites have become charged, they enable controllable movement of CNTs at the bottom electrode under the application of electric field. By reversibly connecting and disconnecting such CNTs from the bottom electrode, the electrical conductivity through the CNT film as a whole can be controlled.

The Impact of Carbon Corrosion Towards Oxygen Reduction Reaction Activity on Atomically Dispersed M-N-C Catalysts

Atomically dispersed Pt group metal free (PGM free) M-N-C catalysts offer a promising pathway towards a clean energy future.1 This class of materials has proven to be a viable substitute for commercial Pt/C to catalyze oxygen reduction reaction (ORR) with similar activity at only a fraction of the raw material cost.2 To gain a wider market adoption, however, M-N-C catalysts need to significantly improve its durability compared to Pt and its alloys in a real-world operating condition. Towards this end, substantial efforts have focused on making more durable M-N-C catalysts while simultaneously maintaining and improving its ORR activity.3,4 As part of the effort to understand the fundamental aspects that control M-N-C catalysts’ durability during ORR condition, we employ a combined approach using first principles density functional theory (DFT) and density functional tight binding theory (DFTB) to identify the impact of carbon corrosion occurring on M-N-C catalyst towards ORR activity. Our starting atomic representation of the ORR active site is the “C10” type local structure composed of FeN4C10 local defect5 embedded in a monolayer graphene support of varying cell size. As carbon corrosion occurs and consequently creating local carbon vacancies surrounding the active site, we then obtain a family of FeN4C11 and FeN4C12 type of local structures. Our DFT calculations of OH binding energy on these corrosion derived structures suggest that the presence of C vacancies and the subsequent H passivation on dangling C introduces elastic structural deformation to various degree that decreases the theoretical ORR limiting potential. This decrease in DFT calculated ORR limiting potential across a family of FeN4C10, FeN4C11, and FeN4C12 local structures can be attributed to a mechanical energy arising from the surface reconstructions. Reference D. A. Cullen et al., Nat Energy, 6, 462–474 (2021).G. Wu, K. L. More, C. M. Johnston, and P. Zelenay, Science, 332, 443–447 (2011).E. F. Holby, G. Wang, and P. Zelenay, ACS Catal., 10, 14527–14539 (2020).T. Patniboon and H. A. Hansen, ACS Catal., 11, 13102–13118 (2021).H. T. Chung et al., Science, 357, 479–484 (2017).

Density Functional Tight-Binding Models for Band Structures of Transition-Metal Alloys and Surfaces across the d-Block.

First-principles electronic structure simulations are an invaluable tool for understanding chemical bonding and reactions. While machine-learning models such as interatomic potentials significantly accelerate the exploration of potential energy surfaces, electronic structure information is generally lost. Particularly in the field of heterogeneous catalysis, simulated electron band structures provide fundamental insights into catalytic reactivity. This ab initio knowledge is preserved in semiempirical methods such as density functional tight binding (DFTB), which extend the accessible computational length and time scales beyond first-principles approaches. In this paper we present Shell-Optimized Atomic Confinement (SOAC) DFTB electronic-part-only parametrizations for bulk and surface band structures of all d-block transition metals that enable efficient predictions of electronic descriptors for large structures or high-throughput studies on complex systems outside the computational reach of density functional theory.

Theoretical and Molecular Dynamics Simulation Approaches to Praseodymium Content Effect on Radiation Shielding and Time Dependent Volume‐Collapse for a Newly Produced Borate Glass

AbstractGamma‐ray, fast neutron and charged particle shielding characteristics of Praseodymium (Pr3+) activated glasses with composition 10ZnF2‐(10‐x)SrO‐15Li2O‐65H3BO3‐xPr2O3(where x = 0.1, 0.5, 1.0, 1.5, and 2.0) are investigated. For this, the radiation shielding parameters are estimated by using Photon Shielding and Dosimetry (Phy‐X/PSD), Interaction of Proton, Alpha, Gamma rays, Electrons, and X‐rays with matter (PAGEX), Stopping Powers and Ranges for Electrons (ESTAR), and The Stopping and Range of Ions in Matter (SRIM) codes. Furthermore, the results are compared comprehensively and comparatively. The glasses with increasing Pr3+ content exhibited higher shielding potential. It can be said that the ZFLHBPr2.0 sample could be used as a shielding material in various applied fields. Additionally, the effect of Pr composition on structural transitions and time‐dependent volume‐collapse of this new family of borate glass doped with Pr is analyzed by molecular dynamics study (MD) based density functional tight binding method. The simulation results show that the addition of Pr led to the cubic‐orthorhombic transformation of simulation cell at determined temperature and 1 atm pressure. Also, a sudden decrease in volume are observed within a short time at 8% Pr composition (x = 2.0) atmospheric pressure. This result provides a strong connection between Pr composition and volume change with time and at constant low pressure.

Corrosion inhibition performance and adsorption mechanism of new synthesized symmetrical diarylidenacetone-based inhibitor on C38 steel in 15 % HCl medium: Theoretical insight and experimental validation

Despite the engineering potential of using organic substances to protect vulnerable metallic materials from corrosive environments, the interaction mechanisms and adsorption processes that lead to the formation of hybrid materials remain poorly understood. Further research is necessary to elucidate how these organic compounds interact with metallic surfaces and the specific pathways through which they adsorb to form protective layers. Understanding these mechanisms is crucial for developing more effective corrosion-resistant hybrid materials, which could significantly enhance the properties of metallic components and provide excellent electrochemical stability in various industrial applications. This work aims to explore the corrosion performance and interfacial mechanism of two symmetrical diarylidenacetone-based compounds, namely, (1E,4E)-1,5-bis(4-methoxyphenyl) penta-1,4-dien-3-one (MPD) and 2,6-di((E)-benzylidene) cyclohexan-1-one (BZC), in a bid to control the corrosion kinetics of C38 steel in aggressive environment. This assessment integrated experimental, surface, and theoretical explorations. The results of our study revealed the remarkable ability of MPD and BZC molecules to generate a dense, resistant protective film on the C38 steel surface. This protective film acted as a barrier, effectively blocking the penetration of corrosive ions and their interaction with the C38 substrate. The inhibition efficiencies of the investigated products exhibited a consistent rise with increased amounts of additives. The inhibition efficiency reached 95 % for MPD and 90 % for BZC at 10−3 M, as calculated from the fitting Nyquist diagrams. The interaction between the inhibitor and Fe(110) was studied utilizing density functional tight binding (DFTB) and Molecular dynamics simulations (MDS). The results demonstrate the formation of covalent bonds, with the molecules aligning in a parallel orientation of MPD and BZC. This alignment with experimental findings provides strong validation for the efficacy of the inhibitors in corrosion inhibition.

Laser pulse induced second- and third-harmonic generation of gold nanorods with real-time time-dependent density functional tight binding (RT-TDDFTB) method.

In this study, we investigate second- and third-harmonic generation processes in Au nanorod systems using the real-time time-dependent density functional tight binding method. Our study focuses on the computation of nonlinear signals based on the time dependent dipole response induced by linearly polarized laser pulses interacting with nanoparticles. We systematically explore the influence of various laser parameters, including pump intensity, duration, frequency, and polarization directions, on harmonic generation. We demonstrate all the results using Au nanorod dimer systems arranged in end-to-end configurations, and disrupting the spatial symmetry of regular single nanorod systems is crucial for second-harmonic generation processes. Furthermore, we study the impact of nanorod lengths, which lead to variable plasmon energies, on harmonic generation, and estimates of polarizabilities and hyper-polarizabilities are provided.

Efficient Parameterization of Density Functional Tight-Binding for 5f-Elements: A Th-O Case Study.

Density functional tight binding (DFTB) models for f-element species are challenging to parametrize owing to the large number of adjustable parameters. The explicit optimization of the terms entering the semiempirical DFTB Hamiltonian related to f orbitals is crucial to generating a reliable parametrization for f-block elements, because they play import roles in bonding interactions. However, since the number of parameters grows quadratically with the number of orbitals, the computational cost for parameter optimization is much more expensive for the f-elements than for the main group elements. In this work we present a set of efficient approaches for mitigating the hurdle imposed by the large size of the parameter space. A novel group-by-orbital correction functions for two-center bond integrals was developed. With this approach the number of parameters is reduced, and it grows linearly with the number of elements, maintaining the accuracy and the number of parameters, in the case of f elements, by more than 40%. The parameter optimization step was accelerated by means of the mini-batch BFGS method. This method allows parameter optimizations with much larger training sets than other single batch methods. A stochastic optimizer was employed that helped overcome shallow local minima in the objective function. The proposed algorithm was used to parametrize the DFTB Hamiltonian for the Th-O system, which was subsequently applied to the study of ThO2 nanoparticles. The training set consisted of 6322 unique structures, which is barely feasible with conventional optimization methods. The optimized parameter set, LANL-ThO, displays good agreement with DFT-calculated properties such as energies, forces, and structures for both clusters and bulk ThO2. Benefiting from the fewer number of parameters and lower computational costs for objective function evaluations, this new approach shows its potential applications in DFTB parametrization for elements with high angular momentum, which present a challenge to conventional methods.

Implementation of Nonadiabatic Molecular Dynamics for Intersystem Crossing Based on a Time-Dependent Density-Functional Tight-Binding Method.

Intersystem crossing (ISC) and internal conversion (IC) are types of nonadiabatic transitions that play important roles in a wide range of fields, including photochemistry, photophysics, and photobiology. The nonadiabatic molecular dynamics (NA-MD) method is a powerful tool for computational simulations of dynamic phenomena involving nonadiabatic transitions. In this study, we implemented the NA-MD method, which treats ISC and IC on an equal footing, where the electronic structure is treated at the level of the time-dependent (TD) density-functional tight-binding (DFTB) method, a low-cost semiempirical analog of TD density functional theory (DFT). In particular, the spin-orbit coupling calculation algorithm was implemented in the TD-DFTB framework, and the results showed trends similar to those obtained using TD-DFT. In addition, the NA-MD method successfully reproduced ultrafast ISC of 2-nitronaphthalene.

Computational predictions of cocrystal formation: A benchmark study of 28 assemblies comparing five methods from high-throughput to advanced models.

Cocrystals are assemblies of more than one type of molecule stabilized through noncovalent interactions. They are promising materials for improved drug formulation in which the stability, solubility, or biocompatibility of the active pharmaceutical ingredient (API) is improved by including a coformer. In this work, a range of density functional theory (DFT) and density functional tight binding (DFTB) models are systematically compared for their ability to predict the lattice enthalpy of a broad range of existing pharmaceutically relevant cocrystals. These range from cocrystals containing model compounds 4,4'-bipyridine and oxalic acid to those with the well benchmarked APIs of aspirin and paracetamol, all tested with a large set of alternative coformers. For simple cocrystals, there is a general consensus in lattice enthalpy calculated by the different DFT models. For the cocrystals with API coformers the cocrystals, enthalpy predictions depend strongly on the DFT model. The significantly lighter DFTB models predict unrealistic values of lattice enthalpy even for simple cocrystals.

Electrooxidation of ethylene glycol using Ag-Pt nanotubes supported on silica: Correlating the unexpected O-vacancies creation with catalytic performance

Electrooxidation of small organic compounds is crucial in achieving clean and efficient energy production. In this context, ethylene glycol (EG) has attracted considerable interest due to its notable energy density, which is suitable for applications in direct alcohol fuel cells. In light of this, we prepared Ag-Pt nanotubes (Ag-PtNTs) supported on silica (SiO2) for the electrooxidation of this alcohol. However, our goal was not only to create an electrocatalyst that enhances electrooxidation efficiency but also to explore the interaction between the SiO2 and the nanotubes. Such an exploration aims to contest the prevailing perception of SiO2 as an inert and non-conductive material. Thus, comprehensive physicochemical characterization of the electrocatalyst was conducted, supported by density functional tight-binding (DFTB) calculations, to elucidate the correlation between composition, electronic effects, and catalytic performance, showing that Ag interacts more with the oxygen of SiO2 that Pt. The interaction between SiO2 and Ag-Pt NTs was explored through Mott-Schottky (M-S) analysis and electron paramagnetic resonance (EPR) spectroscopy, revealing the creation of oxygen vacancies (o-vacancies) and heterojunction formation, which enhance charge transfer and catalytic activity, both supported by theoretical calculations. In addition, mechanistic insights into the electrooxidation process were obtained through potential energy surface (PES) assessments, revealing a concerted pathway for glycolaldehyde formation with low barrier energy and favorable thermodynamics. Finally, the synthesized Ag-PtNTs/Si electrocatalyst exhibited an onset potential/current density of 0.45 V/0.09 mV mA cm−2 at 25 °C, demonstrating efficient electrooxidation capabilities. Therefore, this study provides valuable insights into designing efficient and stable electrocatalysts for alcohol electrooxidation, contributing to improving cleaner energy technologies.

Obtaining Robust Density Functional Tight-Binding Parameters for Solids across the Periodic Table.

The density functional tight-binding (DFTB) approach allows electronic structure-based simulations at length and time scales far beyond what is possible with first-principles methods. This is achieved by using minimal basis sets and empirical approximations. Unfortunately, the sparse availability of parameters across the periodic table is a significant barrier to the use of DFTB in many cases. We therefore propose a workflow that allows the robust and consistent parametrization of DFTB across the periodic table. Importantly, our approach requires no element-pairwise parameters so that the parameters can be used for all element combinations and are readily extendable. This is achieved by parametrizing all elements on a consistent set of artificial homoelemental crystals, spanning a wide range of coordination environments. The transferability of the resulting periodic table baseline parameters to multielement systems and unknown structures is explored and the model is extensively benchmarked against previous specialized and general DFTB parametrizations.