RESUMO
Accurate description of nonadiabatic dynamics of molecules at metal surfaces involving electron transfer has been a long-standing challenge for theory. Here, we tackle this problem by first constructing high-dimensional neural network diabatic potentials including state crossings determined by constrained density functional theory, then applying mixed quantum-classical surface hopping simulations to evolve coupled electron-nuclear motion. Our approach accurately describes the nonadiabatic effects in CO scattering from Au(111) without empirical parameters and yields results agreeing well with experiments under various conditions for this benchmark system. We find that both adiabatic and nonadiabatic energy loss channels have important contributions to the vibrational relaxation of highly vibrationally excited CO(v_{i}=17), whereas relaxation of low vibrationally excited states of CO(v_{i}=2) is weak and dominated by nonadiabatic energy loss. The presented approach paves the way for accurate first-principles simulations of electron transfer mediated nonadiabatic dynamics at metal surfaces.
RESUMO
BACKGROUND: Modular acetabular components for total hip arthroplasty (THA) provide intraoperative flexibility; however, polyethylene liner dissociation may occur. This study aimed to examine the incidence and causes of liner dissociation associated with a specific acetabular component design at a single centre. MATERIALS AND METHODS: A retrospective analysis of 7027 patients who underwent primary THA was performed to identify isolated liner dislocations. Patient demographics, clinical presentations, surgical and implant details, and both radiographic and computed tomography (CT) findings were analysed. Patients with liner dislocation were matched to a control group via 2:1 propensity score matching, and a logistic regression analysis was employed to identify associated risk factors. RESULTS: A total of 32 patients (0.45%) experienced liner dislocation at a mean 71.47 ± 60.10 months post surgery. Significant factors contributing to dislocations included the use of a conventional compared with a highly crosslinked polyethylene component (p = 0.049) and screw fixation (p = 0.028). Radiographic and CT analysis highlighted the importance of proper component orientation, revealing that patients experiencing dislocations demonstrated significantly lower acetabular cup anteversion angles (p = 0.001) compared with the control group. Impingement and malposition, identified in 41% and 47% of the cases, respectively, further underscored the multifactorial nature of dislocation risks. CONCLUSIONS: While the overall rate of polyethylene liner dislocation was low, the findings of this study highlight the importance of appropriate cup placement to decrease the risk of dissociation. It further substantiates the influence of impingement and malposition in liner displacement, with increased mechanical stress exerted on the locking mechanism under adverse conditions and the potential risk increase due to screw placement.
Assuntos
Artroplastia de Quadril , Prótese de Quadril , Polietileno , Desenho de Prótese , Falha de Prótese , Humanos , Artroplastia de Quadril/efeitos adversos , Artroplastia de Quadril/instrumentação , Artroplastia de Quadril/métodos , Estudos Retrospectivos , Masculino , Feminino , Prótese de Quadril/efeitos adversos , Estudos de Casos e Controles , Pessoa de Meia-Idade , Idoso , Fatores de Risco , Acetábulo/cirurgia , Tomografia Computadorizada por Raios X , Pontuação de PropensãoRESUMO
Modern software engineering of electronic structure codes has seen a paradigm shift from monolithic workflows toward object-based modularity. Software objectivity allows for greater flexibility in the application of electronic structure calculations, with particular benefits when integrated with approaches for data-driven analysis. Here, we discuss different approaches to create deep modular interfaces that connect big-data workflows and electronic structure codes and explore the diversity of use cases that they can enable. We present two such interface approaches for the semi-empirical electronic structure package, DFTB+. In one case, DFTB+ is applied as a library and provides data to an external workflow; in another, DFTB+receives data via external bindings and processes the information subsequently within an internal workflow. We provide a general framework to enable data exchange workflows for embedding new machine-learning-based Hamiltonians within DFTB+ or enabling deep integration of DFTB+ in multiscale embedding workflows. These modular interfaces demonstrate opportunities in emergent software and workflows to accelerate scientific discovery by harnessing existing software capabilities.
RESUMO
Predicting the degradation processes of molecules over long time scales is a key aspect of industrial materials design. However, it is made computationally challenging by the need to construct large networks of chemical reactions that are relevant to the experimental conditions that kinetic models must mirror, with every reaction requiring accurate kinetic data. Here, we showcase Kinetica.jl, a new software package for constructing large-scale chemical reaction networks in a fully automated fashion by exploring chemical reaction space with a kinetics-driven algorithm; coupled to efficient machine-learning models of activation energies for sampled elementary reactions, we show how this approach readily enables generation and kinetic characterization of networks containing â¼103 chemical species and ≃104-105 reactions. Symbolic-numeric modeling of the generated reaction networks is used to allow for flexible, efficient computation of kinetic profiles under experimentally realizable conditions such as continuously variable temperature regimes, enabling direct connection between bottom-up reaction networks and experimental observations. Highly efficient propagation of long-time-scale kinetic profiles is required for automated reaction network refinement and is enabled here by a new discrete kinetic approximation. The resulting Kinetica.jl simulation package therefore enables automated generation, characterization, and long-time-scale modeling of complex chemical reaction systems. We demonstrate this for hydrocarbon pyrolysis simulated over time scales of seconds, using transient temperature profiles representing those of tubular flow reactor experiments.
RESUMO
We report the quantitative adsorption structure of pristine graphene on Cu(111) determined using the normal incidence x-ray standing wave technique. The experiments constitute an important benchmark reference for the development of density functional theory approximations able to capture long-range dispersion interactions. Electronic structure calculations based on many-body dispersion-inclusive density functional theory are able to accurately predict the absolute measure and variation of adsorption height when the coexistence of multiple moiré superstructures is considered. This provides a structural model consistent with scanning probe microscopy results.
RESUMO
The role of the inorganic substrate termination, within the organic-inorganic interface, has been well studied for systems that contain strong localised bonding. However, how varying the substrate termination affects coordination to delocalised electronic states, like that found in aromatic molecules, is an open question. Azupyrene, a non-alternant polycyclic aromatic hydrocarbon, is known to bind strongly to metal surfaces through its delocalised π orbitals, thus yielding an ideal probe into delocalised surface-adsorbate interactions. Normal incidence X-ray standing wave (NIXSW) measurements and density functional theory calculations are reported for the adsorption of azupyrene on the (111), (110) and (100) surface facets of copper to investigate the dependence of the adsorption structure on the substrate termination. Structural models based on hybrid density functional theory calculations with non-local many-body dispersion yield excellent agreement with the experimental NIXSW results. No statistically significant difference in the azupyrene adsorption height was observed between the (111) and (100) surfaces. On the Cu(110) surface, the molecule was found to adsorb 0.06 ± 0.04 Å closer to the substrate than on the other surface facets. The most energetically favoured adsorption site on each surface, as determined by DFT, is subtly different, but in each case involved a configuration where the aromatic rings were centred above a hollow site, consistent with previous reports for the adsorption of small aromatic molecules on metal surfaces.
RESUMO
The discrete and charge-separated nature of matter - electrons and nuclei - results in local electrostatic fields that are ubiquitous in nanoscale structures and relevant in catalysis, nanoelectronics and quantum nanoscience. Surface-averaging techniques provide only limited experimental access to these potentials, which are determined by the shape, material, and environment of the nanostructure. Here, we image the potential over adatoms, chains, and clusters of Ag and Au atoms assembled on Ag(111) and quantify their surface dipole moments. By focusing on the total charge density, these data establish a benchmark for theory. Our density functional theory calculations show a very good agreement with experiment and allow a deeper analysis of the dipole formation mechanisms, their dependence on fundamental atomic properties and on the shape of the nanostructures. We formulate an intuitive picture of the basic mechanisms behind dipole formation, allowing better design choices for future nanoscale systems such as single-atom catalysts.
RESUMO
The reactive chemistry of molecular hydrogen at surfaces, notably dissociative sticking and hydrogen evolution, plays a crucial role in energy storage and fuel cells. Theoretical studies can help to decipher underlying mechanisms and reaction design, but studying dynamics at surfaces is computationally challenging due to the complex electronic structure at interfaces and the high sensitivity of dynamics to reaction barriers. In addition, ab initio molecular dynamics, based on density functional theory, is too computationally demanding to accurately predict reactive sticking or desorption probabilities, as it requires averaging over tens of thousands of initial conditions. High-dimensional machine learning-based interatomic potentials are starting to be more commonly used in gas-surface dynamics, yet robust approaches to generate reliable training data and assess how model uncertainty affects the prediction of dynamic observables are not well established. Here, we employ ensemble learning to adaptively generate training data while assessing model performance with full uncertainty quantification (UQ) for reaction probabilities of hydrogen scattering on different copper facets. We use this approach to investigate the performance of two message-passing neural networks, SchNet and PaiNN. Ensemble-based UQ and iterative refinement allow us to expose the shortcomings of the invariant pairwise-distance-based feature representation in the SchNet model for gas-surface dynamics.
RESUMO
Many-body dispersion (MBD) is a powerful framework to treat van der Waals (vdW) dispersion interactions in density-functional theory and related atomistic modeling methods. Several independent implementations of MBD with varying degree of functionality exist across a number of electronic structure codes, which both limits the current users of those codes and complicates dissemination of new variants of MBD. Here, we develop and document libMBD, a library implementation of MBD that is functionally complete, efficient, easy to integrate with any electronic structure code, and already integrated in FHI-aims, DFTB+, VASP, Q-Chem, CASTEP, and Quantum ESPRESSO. libMBD is written in modern Fortran with bindings to C and Python, uses MPI/ScaLAPACK for parallelization, and implements MBD for both finite and periodic systems, with analytical gradients with respect to all input parameters. The computational cost has asymptotic cubic scaling with system size, and evaluation of gradients only changes the prefactor of the scaling law, with libMBD exhibiting strong scaling up to 256 processor cores. Other MBD properties beyond energy and gradients can be calculated with libMBD, such as the charge-density polarization, first-order Coulomb correction, the dielectric function, or the order-by-order expansion of the energy in the dipole interaction. Calculations on supramolecular complexes with MBD-corrected electronic structure methods and a meta-review of previous applications of MBD demonstrate the broad applicability of the libMBD package to treat vdW interactions.
RESUMO
Polycrystalline boron-doped diamond (BDD) is widely used as a working electrode material in electrochemistry, and its properties, such as its stability, make it an appealing support material for nanostructures in electrocatalytic applications. Recent experiments have shown that electrodeposition can lead to the creation of stable small nanoclusters and even single gold adatoms on the BDD surfaces. We investigate the adsorption energy and kinetic stability of single gold atoms adsorbed onto an atomistic model of BDD surfaces by using density functional theory. The surface model is constructed using hybrid quantum mechanics/molecular mechanics embedding techniques and is based on an oxygen-terminated diamond (110) surface. We use the hybrid quantum mechanics/molecular mechanics method to assess the ability of different density functional approximations to predict the adsorption structure, energy, and barrier for diffusion on pristine and defective surfaces. We find that surface defects (vacancies and surface dopants) strongly anchor adatoms on vacancy sites. We further investigated the thermal stability of gold adatoms, which reveals high barriers associated with lateral diffusion away from the vacancy site. The result provides an explanation for the high stability of experimentally imaged single gold adatoms on BDD and a starting point to investigate the early stages of nucleation during metal surface deposition.
RESUMO
Mixed quantum-classical (MQC) methods for simulating the dynamics of molecules at metal surfaces have the potential to accurately and efficiently provide mechanistic insight into reactive processes. Here, we introduce simple two-dimensional models for the scattering of diatomic molecules at metal surfaces based on recently published electronic structure data. We apply several MQC methods to investigate their ability to capture how nonadiabatic effects influence molecule-metal energy transfer during the scattering process. Specifically, we compare molecular dynamics with electronic friction, Ehrenfest dynamics, independent electron surface hopping, and the broadened classical master equation approach. In the case of independent electron surface hopping, we implement a simple decoherence correction approach and assess its impact on vibrationally inelastic scattering. Our results show that simple, low-dimensional models can be used to qualitatively capture experimentally observed vibrational energy transfer and provide insight into the relative performance of different MQC schemes. We observe that all approaches predict similar kinetic energy dependence but return different vibrational energy distributions. Finally, by varying the molecule-metal coupling, we can assess the coupling regime in which some MQC methods become unsuitable.
RESUMO
X-ray photoemission and X-ray absorption spectroscopy are important techniques to characterize chemical bonding at surfaces and are often used to identify the strength and nature of adsorbate-substrate interactions. In this study, we judge the ability of X-ray spectroscopic techniques to identify different regimes of chemical bonding at metal-organic interfaces. To achieve this, we sample different interaction strength regimes in a comprehensive and systematic way by comparing two topological isomers, azulene and naphthalene, adsorbed on three metal substrates with varying reactivity, namely the (111) facets of Ag, Cu, and Pt. Using density functional theory, we simulate core-level binding energies and X-ray absorption spectra of the molecular carbon species. The simulated spectra reveal three distinct characteristics based on the molecule-specific spectral features which we attribute to types of surface chemical bonding with varying strength. We find that weak physisorption only leads to minor changes compared to the gas-phase spectra, weak chemisorption leads to charge transfer and significant spectral changes, and strong chemisorption leads to a loss of the molecule-specific features in the spectra. The classification we provide is aimed at assisting interpretation of experimental X-ray spectra for complex metal-organic interfaces.
RESUMO
Independent electron surface hopping (IESH) is a computational algorithm for simulating the mixed quantum-classical molecular dynamics of adsorbate atoms and molecules interacting with metal surfaces. It is capable of modeling the nonadiabatic effects of electron-hole pair excitations on molecular dynamics. Here, we present a transparent, reliable, and efficient implementation of IESH, demonstrating its ability to predict scattering and desorption probabilities across a variety of systems, ranging from model Hamiltonians to full dimensional atomistic systems. We further show how the algorithm can be modified to account for the application of an external bias potential, comparing its accuracy to results obtained using the hierarchical quantum master equation. Our results show that IESH is a practical method for modeling coupled electron-nuclear dynamics at metal surfaces, especially for highly energetic scattering events.
RESUMO
The design of molecules and materials with tailored properties is challenging, as candidate molecules must satisfy multiple competing requirements that are often difficult to measure or compute. While molecular structures produced through generative deep learning will satisfy these patterns, they often only possess specific target properties by chance and not by design, which makes molecular discovery via this route inefficient. In this work, we predict molecules with (Pareto-)optimal properties by combining a generative deep learning model that predicts three-dimensional conformations of molecules with a supervised deep learning model that takes these as inputs and predicts their electronic structure. Optimization of (multiple) molecular properties is achieved by screening newly generated molecules for desirable electronic properties and reusing hit molecules to retrain the generative model with a bias. The approach is demonstrated to find optimal molecules for organic electronics applications. Our method is generally applicable and eliminates the need for quantum chemical calculations during predictions, making it suitable for high-throughput screening in materials and catalyst design.
Assuntos
Eletrônica , Ensaios de Triagem em Larga EscalaRESUMO
The computational prediction of the structure and stability of hybrid organic-inorganic interfaces provides important insights into the measurable properties of electronic thin film devices, coatings, and catalyst surfaces and plays an important role in their rational design. However, the rich diversity of molecular configurations and the important role of long-range interactions in such systems make it difficult to use machine learning (ML) potentials to facilitate structure exploration that otherwise requires computationally expensive electronic structure calculations. We present an ML approach that enables fast, yet accurate, structure optimizations by combining two different types of deep neural networks trained on high-level electronic structure data. The first model is a short-ranged interatomic ML potential trained on local energies and forces, while the second is an ML model of effective atomic volumes derived from atoms-in-molecules partitioning. The latter can be used to connect short-range potentials to well-established density-dependent long-range dispersion correction methods. For two systems, specifically gold nanoclusters on diamond (110) surfaces and organic π-conjugated molecules on silver (111) surfaces, we train models on sparse structure relaxation data from density functional theory and show the ability of the models to deliver highly efficient structure optimizations and semi-quantitative energy predictions of adsorption structures.
RESUMO
Different types of immune cells are involved in atherogenesis and may act atheroprotective or atheroprogressive. Here, we describe an in vitro approach to analyze CD11c+ cells and CD11c+-derived ApoE in atherosclerosis. The major steps include harvesting mouse bone marrow, plating cells in culture dishes, treating them with differentiation factors, and collecting cells after removal of undesirable populations. This protocol can be adapted for CD11c+ cells in different contexts, thus, serving as models for different diseases and to analyze cell-specific molecules. For complete details on the use and execution of this protocol, please refer to Sauter et al. (2021).
Assuntos
Medula Óssea , Células Dendríticas , Animais , Apolipoproteínas E , Células da Medula Óssea , Antígeno CD11c , CamundongosRESUMO
Molecular energy transfer and reactions at solid surfaces depend on the molecular orientation relative to the surface. While such steric effects have been largely understood in electronically adiabatic processes, the orientation-dependent energy transfer in NO scattering from Au(111) was complicated by electron-mediated nonadiabatic effects, thus lacking a clear interpretation and posing a great challenge for theories. Herein, we investigate the stereodynamics of adiabatic and nonadiabatic energy transfer via molecular dynamics simulations of NO(v = 3) scattering from Au(111) using realistic initial orientation distributions based on accurate neural network fitted adiabatic potential energy surface and electronic friction tensor. Our results reproduce the observed stronger vibrational relaxation for N-first orientation and enhanced rotational rainbow for O-first orientation, and demonstrate how adiabatic anisotropic interactions steer molecules into the more attractive N-first orientation to experience more significant energy transfer. Remaining disagreements with experiment suggest the direction for further developments of nonadiabatic theories for gas-surface scattering.
RESUMO
Here, we describe an in vivo approach to visualize CD11c+ cells in atherosclerosis. In particular, we use a protocol for X-Gal staining of immune cells within atherosclerotic plaques, which can be used as an alternative to analyze plaque composition and cell-specific molecules in atherogenesis. LacZ knockin mice have to be bred to mice carrying the CD11ccre recombinase-both brought onto an ApoE-/- background-to be able to visualize this cell type of interest in the plaques by X-Gal staining. With this approach, different immune cells in atherogenesis can be examined. For complete details on the use and execution of this protocol, please refer to Sauter et al. (2021).
Assuntos
Aterosclerose , Placa Aterosclerótica , Animais , Aterosclerose/genética , Antígeno CD11c/genética , Óperon Lac/genética , Camundongos , Camundongos Knockout , Placa Aterosclerótica/genéticaRESUMO
Defects play a critical role for the functionality and performance of materials, but the understanding of the related effects is often lacking, because the typically low concentrations of defects make them difficult to study. A prominent case is the topological defects in two-dimensional materials such as graphene. The performance of graphene-based (opto-)electronic devices depends critically on the properties of the graphene/metal interfaces at the contacting electrodes. The question of how these interface properties depend on the ubiquitous topological defects in graphene is of high practical relevance, but could not be answered so far. Here, we focus on the prototypical Stone-Wales (S-W) topological defect and combine theoretical analysis with experimental investigations of molecular model systems. We show that the embedded defects undergo enhanced bonding and electron transfer with a copper surface, compared to regular graphene. These findings are experimentally corroborated using molecular models, where azupyrene mimics the S-W defect, while its isomer pyrene represents the ideal graphene structure. Experimental interaction energies, electronic-structure analysis, and adsorption distance differences confirm the defect-controlled bonding quantitatively. Our study reveals the important role of defects for the electronic coupling at graphene/metal interfaces and suggests that topological defect engineering can be used for performance control.
RESUMO
Accurate and efficient methods to simulate nonadiabatic and quantum nuclear effects in high-dimensional and dissipative systems are crucial for the prediction of chemical dynamics in the condensed phase. To facilitate effective development, code sharing, and uptake of newly developed dynamics methods, it is important that software implementations can be easily accessed and built upon. Using the Julia programming language, we have developed the NQCDynamics.jl package, which provides a framework for established and emerging methods for performing semiclassical and mixed quantum-classical dynamics in the condensed phase. The code provides several interfaces to existing atomistic simulation frameworks, electronic structure codes, and machine learning representations. In addition to the existing methods, the package provides infrastructure for developing and deploying new dynamics methods, which we hope will benefit reproducibility and code sharing in the field of condensed phase quantum dynamics. Herein, we present our code design choices and the specific Julia programming features from which they benefit. We further demonstrate the capabilities of the package on two examples of chemical dynamics in the condensed phase: the population dynamics of the spin-boson model as described by a wide variety of semiclassical and mixed quantum-classical nonadiabatic methods and the reactive scattering of H2 on Ag(111) using the molecular dynamics with electronic friction method. Together, they exemplify the broad scope of the package to study effective model Hamiltonians and realistic atomistic systems.