Search NASA⌕ Search

SEARCH · Search NASA

Results for “MOLECULAR DYNAMICS METHOD”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Computational determination of a primary diffusion mode in γ U-10Mo under irradiation

Low enriched uranium (< 20 % 235 U)-molybdenum (U-Mo) monolithic fuel is the primary candidate for high performance research and test reactors, and is in the process of being qualified to replace highly enriched uranium (≥ 20 % 23 5U) fuel. As part of the qualification process, it is critical to understand and predict the behavior of fission gas bubbles under irradiation, which affects fuel swelling and fuel failure. Mechanistic fuel models are being developed that can both reproduce the existing experimental data for fuel swelling, and be further applied to irradiation conditions beyond the experimental scope. Diffusion of species under irradiation conditions is an important parameter in the mechanistic fuel models; however, no temperature-relevant experimental diffusion data exists. In the present work, radiation-enhanced diffusion coefficients of U, Mo, and Xe in λU-10wt.%Mo were calculated in the temperature range between 300 K and 1400 K via rate-theory models and molecular dynamics simulations with an embedded-atom method interatomic potential for the U-Mo-Xe system. Accordingly, total diffusion coefficients under relevant irradiation conditions are determined using previously obtained intrinsic thermal diffusion and radiation-driven diffusion coefficients, as well as the newly calculated radiation-enhanced diffusion coefficients presented herein. Radiation-enhanced diffusion of U and Mo was dominant in the intermediate temperature range, whereas radiation-enhanced diffusion of Xe did not significantly contribute to total diffusion of Xe at the relevant fission rate densities. Radiation-enhanced diffusion of Xe became faster than both intrinsic thermal diffusion and radiation-driven diffusion at a fission rate density of 5 x 10 22 fissions/m 3 /s, which is higher than the typical fission rate density range in research reactors. The temperature regime that radiation-enhanced diffusion of each element dominated was dependent on the fission rate density. Finally, the total diffusion coefficients of U, Mo, and Xe, updated in this work, will be utilized as parameters in the mechanistic fuel models to help predict the behavior of fission gas bubbles under irradiation more accurately.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Computational methods in solution-based plastics purification

Plastic waste can be recycled into resins with near-virgin properties by solution-based purification processes that selectively dissolve polymers, remove contaminants, or detach printing residues. Here, in this review, we examine computational methods for predicting the behavior governing solution-based plastic purification, motivated by the vast polymer–solvent–contaminant compositional space. We discuss thermodynamic and machine learning methods for predicting polymer–solvent and polymer–contaminant interaction and review physics-based molecular dynamics simulations that resolve molecular-scale phenomena within polymer matrices inaccessible to screening methods. We highlight how these methods have informed experimental design for dissolution-based recycling and solvent-based contaminant removal. Finally, we discuss the prospective role of agentic AI in integrating these computational tools with real-time sorting data to adapt purification conditions to the compositional variability of real post-consumer feedstocks. This review charts a path toward computationally guided solution-based purification workflows that can respond to the complexity inherent in plastic waste streams.

Altamimi, Ali [Univ. of Wisconsin, Madison, WI (Un↗
DOE OSTI2023

Simulating dielectric spectra: A demonstration of the direct electric field method and a new model for the nonlinear dielectric response

Here we demonstrate a method to compute the dielectric spectra of fluids in molecular dynamics (MD) by directly applying electric fields to the simulation. We obtain spectra from MD simulations with low magnitude electric fields (≈0.01 V/Å) in agreement with spectra from the fluctuation–dissipation method for water and acetonitrile. We examine this method’s trade-off between noise at low field magnitudes and the nonlinearity of the response at higher field magnitudes. We then apply the Booth equation to describe the nonlinear response of both fluids at low frequency (0.1 GHz) and high field magnitude (up to 0.5 V/Å). We develop a model of the frequency-dependent nonlinear response by combining the Booth description of the static nonlinear dielectric response of fluids with the frequency-dependent linear dielectric response of the Debye model. We find good agreement between our model and the MD simulations of the nonlinear dielectric response for both acetonitrile and water.

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

Solvation at metal/water interfaces: An ab initio molecular dynamics benchmark of common computational approaches

Determining the influence of the solvent on electrochemical reaction energetics is a central challenge in our understanding of electrochemical interfaces. To date, it is unclear how well existing methods predict solvation energies at solid/liquid interfaces, since they cannot be assessed experimentally. Ab initio molecular dynamics (AIMD) simulations present a physically highly accurate, but also a very costly approach. In this work, we employ extensive AIMD simulations to benchmark solvation at charge-neutral metal/water interfaces against commonly applied continuum solvent models. We consider a variety of adsorbates including *CO, *CHO, *COH, *OCCHO, *OH, and *OOH on Cu, Au, and Pt facets solvated by water. The surfaces and adsorbates considered are relevant, among other reactions, to electrochemical CO 2 reduction and the oxygen redox reactions. We determine directional hydrogen bonds and steric water competition to be critical for a correct description of solvation at the metal/water interfaces. As a consequence, we find that the most frequently applied continuum solvation methods, which do not yet capture these properties, do not presently provide more accurate energetics over simulations in vacuum. We find most of the computed benchmark solvation energies to linearly scale with hydrogen bonding or competitive water adsorption, which strongly differ across surfaces. Furthermore, we determine solvation energies of adsorbates to be non-transferable between metal surfaces, in contrast to standard practice.

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

Assessing electrogenetic activation via a network model of biological signal propagation

Introduction: Molecular communication is the transfer of information encoded by molecular structure and activity. We examine molecular communication within bacterial consortia as cells with diverse biosynthetic capabilities can be assembled for enhanced function. Their coordination, both in terms of engineered genetic circuits within individual cells as well as their population-scale functions, is needed to ensure robust performance. We have suggested that “electrogenetics,” the use of electronics to activate specific genetic circuits, is a means by which electronic devices can mediate molecular communication, ultimately enabling programmable control. Methods: Here, we have developed a graphical network model for dynamically assessing electronic and molecular signal propagation schemes wherein nodes represent individual cells, and their edges represent communication channels by which signaling molecules are transferred. We utilize graph properties such as edge dynamics and graph topology to interrogate the signaling dynamics of specific engineered bacterial consortia. Results: We were able to recapitulate previous experimental systems with our model. In addition, we found that networks with more distinct subpopulations (high network modularity) propagated signals more slowly than randomized networks, while strategic arrangement of subpopulations with respect to the inducer source (an electrode) can increase signal output and outperform otherwise homogeneous networks. Discussion: We developed this model to better understand our previous experimental results, but also to enable future designs wherein subpopulation composition, genetic circuits, and spatial configurations can be varied to tune performance. We suggest that this work may provide insight into the signaling which occurs in synthetically assembled systems as well as native microbial communities.

Chun, Kayla↗
DOE OSTI2021

Modeling protein structures from predicted contacts with modern molecular dynamics potentials: accuracy, sensitivity, and refinement

Protein structure prediction has become increasingly popular and successful in recent years. An essential step for fragment-free, template-free methods is the generation of a final three-dimensional protein model from a set of predicted amino acid contacts that are often described by interresidue pairwise atomic distances. Here we explore the use of modern, open-source molecular dynamics (MD) engines, which have been continually developed over the last three decades with all-atom Hamiltonians to model biomolecular structure and dynamics, to generate accurate protein structures starting from a set of inferred pairwise distances. Additionally, the ability of MD empirical physical potentials to correct inaccuracies in the predicted geometries is tested. We rigorously characterize the effect of modeling parameters on results, the effect of different amounts of error in the predicted distances on the final structures, and test the ability of post-processing analysis to sort the best models out of a set of statistical replicas. We find that with exact distances and with noisy distances, the method can produce excellent structural models, and that the molecular dynamics force field seems to help correct errors in distance predictions, resisting the effects of applied noise.

Davidson, Russ↗
DOE OSTI2024

Machine Learning Potentials with the Iterative Boltzmann Inversion: Training to Experiment

Methodologies for training machine learning potentials (MLPs) with quantum-mechanical simulation data have recently seen tremendous progress. Experimental data have a very different character than simulated data, and most MLP training procedures cannot be easily adapted to incorporate both types of data into the training process. Here, we investigate a training procedure based on iterative Boltzmann inversion that produces a pair potential correction to an existing MLP using equilibrium radial distribution function data. By applying these corrections to an MLP for pure aluminum based on density functional theory, we observe that the resulting model largely addresses previous overstructuring in the melt phase. Interestingly, the corrected MLP also exhibits improved performance in predicting experimental diffusion constants, which are not included in the training procedure. The presented method does not require autodifferentiating through a molecular dynamics solver and does not make assumptions about the MLP architecture. Our results suggest a practical framework for incorporating experimental data into machine learning models to improve the accuracy of molecular dynamics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Active species in chloroaluminate ionic liquids catalyzing low-temperature polyolefin deconstruction

Abstract Chloroaluminate ionic liquids selectively transform (waste) polyolefins into gasoline-range alkanes through tandem cracking-alkylation at temperatures below 100 °C. Further improvement of this process necessitates a deep understanding of the nature of the catalytically active species and the correlated performance in the catalyzing critical reactions for the tandem polyolefin deconstruction with isoalkanes at low temperatures. Here, we address this requirement by determining the nuclearity of the chloroaluminate ions and their interactions with reaction intermediates, combining in situ 27 Al magic-angle spinning nuclear magnetic resonance spectroscopy, in situ Raman spectroscopy, Al K-edge X-ray absorption near edge structure spectroscopy, and catalytic activity measurement. Cracking and alkylation are facilitated by carbenium ions initiated by AlCl 3 - tert -butyl chloride (TBC) adducts, which are formed by the dissociation of Al 2 Cl 7 − in the presence of TBC. The carbenium ions activate the alkane polymer strands and advance the alkylation cycle through multiple hydride transfer reactions. In situ 1 H NMR and operando infrared spectroscopy demonstrate that the cracking and alkylation processes occur synchronously; alkenes formed during cracking are rapidly incorporated into the carbenium ion-mediated alkylation cycle. The conclusions are further supported by ab initio molecular dynamics simulations coupled with an enhanced sampling method, and model experiments using n-hexadecane as a feed.

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

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Theoretical Model for Computing Freezing Point Depression of Lithium-Ion Battery Electrolytes

Reliable prediction of freezing point depression in liquid electrolytes will accelerate the development of improved Li-ion batteries which can operate in low temperature environments. In this work we establish a computational methodology to calculate activity coefficients and liquidus lines for battery-relevant liquid electrolytes. Electronic structure methods are used in conjuction with classical molecular dynamics simulations and theoretical expressions for Born solvation energy, ion-atmosphere effects from Debye-Hückel theory and solvent entropic effects. The framework uses no a priori knowledge beyond neat solvent properties and the concentration of salt. LiPF 6 in propylene carbonate (PC), LiPF 6 in dimethyl carbonate (DMC) and LiClO 4 in DMC are investigated up to 1 molal with accuracy better than 3 °C when compared to experimental freezing point measurements. Finally, we find that the difference in freezing point depression between the propylene carbonate-based electrolyte and the dimethyl carbonate electrolytes originates from the difference in the solvent dielectric constant.

25 ENERGY STORAGE↗
DOE OSTI2024

FY24 Progress Report on Viscosity and Thermal Conductivity Measurements of Nuclear Industry Relevant Chloride Salts: An Experimental and Computational Study

As presented in this report, experimental and computational techniques were performed to assess the viscosity and thermal conductivity of key alkali and actinide chloride mixtures for molten salt reactor developers. These mixtures were pure LiCl, NaCl-KCl, LiCl-NaCl, LiCl-KCl, LiCl-NaCl-KCl, and NaCl-UCl 3 . Experimental measurements of viscosity were performed with a rolling ball viscometer, whereas experimental measurements of thermal conductivity were performed with a variable gap apparatus. Additional benchmarking work was performed using both property measurement systems to prepare for x-ray radiography in stainless-steel crucibles for viscosity and to ensure that calibration methods were accurate for thermal conductivity before assessing the NaCl-UCl 3 system. Validation data for the NaCl-UCl 3 in literature are minimal. Details on the calibration methods, salt measurement processes, and sources of error and uncertainty are discussed in detail for both property measurements. The computational methods described herein involved ab-initio molecular dynamics (AIMD) calculations using CP2K. The calculations were performed for the LiCl-KCl-NaCl and NaCl-UCl 3 systems. These calculations not only provided thermophysical property estimations for comparison to experimental data, but they also allowed for the determination of diffusion coefficients, coordination numbers, and radial distribution functions to provide insight into ion mobility and local coordination environments, which is linked to macroscopic property trends.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Mechanistic Studies of a Primitive Homolog of Nitrogenase Involved in Coenzyme F430 Biosynthesis

Methyl-coenzyme M reductase (MCR) is the key enzyme in the biological formation and anaerobic oxidation of methane (AOM). Methane is a potent greenhouse gas and the major component of natural gas. Given the abundance of natural gas reserves in remote areas, there is great current interest in a scalable bio-based process for the conversion of methane to liquid fuel or other high-value commodity chemicals. MCR holds much promise for use in such a methane bioconversion strategy. However, MCR cannot currently be produced in an active form in a heterologous host, due in large part to the lack of genetic and biochemical information about the production of holo MCR. In an effort to overcome this deficiency, our laboratory elucidated the biosynthetic pathway of the unique nickel-containing coenzyme of MCR, F430. The key step in coenzyme F430 biosynthesis (Cfb) was found to involve an unprecedented reductive cyclization reaction that converts Ni-sirohydrochlorin a , c -diamide to 15,17 3 -seco-F430-17 3 -acid. This remarkable transformation, which involves a 6-electron reduction of the isobacteriochlorin ring system, cyclization of the c -acetamide side chain to form a γ-lactam ring, and the formation of 7 stereocenters, is catalyzed by a primitive homolog of nitrogenase (CfbCD). Nitrogenase is a two-component metalloenzyme that catalyzes the ATP-dependent reduction of dinitrogen to ammonia (nitrogen fixation). Homologs of nitrogenase are also involved in the biosynthesis of the photosynthetic pigments chlorophyll and bacteriochlorophyll. Phylogenetic analysis of the CfbCD complex suggests that it is representative of a more ancient lineage of the nitrogenase superfamily, and a thorough investigation of its structure and function is likely to shed light on the mechanisms and evolution of these important metalloenzymes that catalyze multi-electron redox reactions. Moreover, a detailed understanding of the mechanism of the CfbCD complex may aid in the development of specific inhibitors to help reduce natural greenhouse gas emissions and can be exploited for the heterologous production of MCR for methane bioconversion. Towards these goals, the following Specific Aims will be pursued to determine the: 1) Identity of the CfbCD reaction product. The exact reaction catalyzed by CfbCD, including the number of electrons transferred and whether it involves enzymatic cyclization, is unclear. Several approaches, including reaction stoichiometry measurements, spectroelectrochemistry, and magnetic resonance spectroscopy will be applied to elucidate the structure of the reaction product and establish whether CfbCD is a reductase or reductive cyclase. 2) Structure, conformational dynamics, and oligomerization state changes of CfbCD. Significant insight into the mechanism and allosteric regulation of CfbCD can be obtained by assessing changes in the structure and dynamics of the complex during the catalytic cycle. To accomplish this, a combination of size-exclusion chromatography, hydrogen-deuterium exchange mass spectrometry, molecular dynamics simulations, and high-resolution structural methods will be employed. 3) Source, order, and stereochemistry of proton additions during CfbCD catalysis. Details regarding the order and stereochemistry of proton additions during the CfbCD reaction will be uncovered using a combined spectroscopic and computational approach. Complementary mechanistic studies employing site-directed mutagenesis and substrate analogs will establish the identity of active site acid residues and the possible involvement of substrate-assisted catalysis during the CfbCD reaction.

09 BIOMASS FUELS↗
DOE OSTI2020

Uncertainty Quantification in Atomistic Modeling of Metals and Its Effect on Mesoscale and Continuum Modeling: A Review

The design of next-generation alloys through the integrated computational materials engineering (ICME) approach relies on multiscale computer simulations to provide thermodynamic properties when experiments are difficult to conduct. Atomistic methods such as density functional theory (DFT) and molecular dynamics (MD) have been successful in predicting properties of never before studied compounds or phases. However, uncertainty quantification (UQ) of DFT and MD results is rarely reported due to computational and UQ methodology challenges. Over the past decade, studies that mitigate this gap have emerged. These advances are reviewed in the context of thermodynamic modeling and information exchange with mesoscale methods such as the phase-field method (PFM) and calculation of phase diagrams (CALPHAD). The importance of UQ is illustrated using properties of metals, with aluminum as an example, and highlighting deterministic, frequentist, and Bayesian methodologies. Finally, challenges facing routine uncertainty quantification and an outlook on addressing them are also presented.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effect of an oxide layer on high velocity impact of tantalum particles characterized using molecular dynamics

Micro-cold spray, also commonly referred to as the aerosol deposition method or vacuum kinetic spraying, is a process for depositing nanostructured films of metals and ceramics by impacting solid nanoparticles at high velocity onto a substrate at room temperature and at low pressures. Tantalum films deposited by micro-cold spray are desirable because tantalum is non-reactive to certain molten metals and can therefore serve as a non-reactive barrier coating. However, unlike most other metals deposited by micro-cold spray, tantalum quickly forms an oxide layer when exposed to air. In the deposition of films using traditional cold spray with much larger particles, native oxides must fracture for particle–substrate bonding to occur. Furthermore, the fraction of the particle that is oxidized in micro-cold spray and the strain rates experienced by the particles upon impact are both much larger than in cold spray, suggesting that the role of the oxide on film deposition may be distinct. The effect of this oxide layer on particle deformation and adherence to a substrate under high strain rates is studied using molecular dynamics simulations. The effect of oxide layer thickness, particle diameter, and impact velocity are examined.

42 ENGINEERING↗
DOE OSTI2021

The effect of protein mutations on drug binding suggests ensuing personalised drug selection

The advent of personalised medicine promises a deeper understanding of mechanisms and therefore therapies. However, the connection between genomic sequences and clinical treatments is often unclear. We studied 50 breast cancer patients belonging to a population-cohort in the state of Qatar. From Sanger sequencing, we identified several new deleterious mutations in the estrogen receptor 1 gene (ESR1). The effect of these mutations on drug treatment in the protein target encoded by ESR1, namely the estrogen receptor, was achieved via rapid and accurate protein–ligand binding affinity interaction studies which were performed for the selected drugs and the natural ligand estrogen. Four nonsynonymous mutations in the ligand-binding domain were subjected to molecular dynamics simulation using absolute and relative binding free energy methods, leading to the ranking of the efficacy of six selected drugs for patients with the mutations. Our study shows that a personalised clinical decision system can be created by integrating an individual patient’s genomic data at the molecular level within a computational pipeline which ranks the efficacy of binding of particular drugs to variant proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Supervised learning and the finite-temperature string method for computing committor functions and reaction rates

A central object in the computational studies of rare events is the committor function. Though costly to compute, the committor function encodes complete mechanistic information of the processes involving rare events, including reaction rates and transition-state ensembles. Under the framework of transition path theory, Rotskoff et al. [Proceedings of the 2nd Mathematical and Scientific Machine Learning Conference, Proceedings of Machine Learning Research (PLMR, 2022), Vol. 145, pp. 757–780] proposes an algorithm where a feedback loop couples a neural network that models the committor function with importance sampling, mainly umbrella sampling, which collects data needed for adaptive training. Here, in this work, we show additional modifications are needed to improve the accuracy of the algorithm. The first modification adds elements of supervised learning, which allows the neural network to improve its prediction by fitting to sample-mean estimates of committor values obtained from short molecular dynamics trajectories. The second modification replaces the committor-based umbrella sampling with the finite-temperature string (FTS) method, which enables homogeneous sampling in regions where transition pathways are located. We test our modifications on low-dimensional systems with non-convex potential energy where reference solutions can be found via analytical or finite element methods, and show how combining supervised learning and the FTS method yields accurate computation of committor functions and reaction rates. We also provide an error analysis for algorithms that use the FTS method, using which reaction rates can be accurately estimated during training with a small number of samples. The methods are then applied to a molecular system in which no reference solution is known, where accurate computations of committor functions and reaction rates can still be obtained.

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

Evidence for nitrogen binding to surface defects for topological insulator Bi 2 Se 3

Using scanning tunneling spectroscopy and theoretical simulations we have studied the effects of nitrogen gas exposure on the electronic density of states of Bi 2 Se 3 , a well-studied topological insulator. In carefully controlled measurements, Bi 2 Se 3 crystals were initially cleaved in a helium gas environment and then exposed to a 22 SCFH flow of ultra-high purity N 2 gas. We observe a resulting change in the spectral curves, with the exposure effect saturating after approximately 50 min, ultimately bringing the Dirac point about 50 meV closer to the Fermi level. Furthermore, these results are compared to density functional theory calculations, which support a picture of N 2 molecules physisorbing near Se vacancies. Furthermore, ab initio molecular dynamics simulations aided by a Blue Moon ensemble method reveal the dissociative adsorption of N 2 molecules which then bind strongly to Se vacancies at the surface. In this scenario, the binding of the N atom to a Se vacancy site removes the surface defect state created by the vacancy and changes the position of the Fermi energy with respect to the Dirac point.

36 MATERIALS SCIENCE↗