We explore the character angle dependence of dislocation-solute interactions in a face-centered cubic random Fe0.70Ni0.11Cr0.19 alloy through molecular dynamics (MD) simulations of dislocation mobility. Using the MD mobility data, we determine the phonon and thermally activated solute drag parameters which govern mobility for each dislocation character angle. The resulting parameter set indicates that, surprisingly, the solute energy barrier does not depend on character angle. Instead, only the zero-temperature flow stress—which is dictated by the activation area for thermal activation—is dependent on character angle. By analyzing the line roughness from MD simulations and the geometry of a bowing dislocation line undergoing thermal activation, we conclude that the character angle dependence of the activation area in this alloy is governed by the dislocation line tension, rather than the dislocation-solute interaction itself. Our findings motivate further investigation into the line geometry of dislocations in solid solutions.
The fundamental interactions between an edge dislocation and a random solid solution are studied by analyzing dislocation line roughness profiles obtained from molecular dynamics simulations of Fe0.70Ni0.11Cr0.19 over a range of stresses and temperatures. These roughness profiles reveal the hallmark features of a depinning transition. Namely, below a temperature-dependent critical stress, the dislocation line exhibits roughness in two different length scale regimes which are divided by a so-called correlation length. This correlation length increases with applied stress and at the critical stress (depinning transition or yield stress) formally goes to infinity. Above the critical stress, the line roughness profile converges to that of a random noise field. Motivated by these results, a physical model is developed based on the notion of coherent line bowing over all length scales below the correlation length. Above the correlation length, the solute field prohibits such coherent line bow outs. Using this model, we identify potential gaps in existing theories of solid solution strengthening and show that recent observations of length-dependent dislocation mobilities can be rationalized.
There is substantial demand for theoretical/computational tools that can produce correct predictions of the geometric structure and band gap to accelerate the design and screening of new materials with desirable electronic properties. DFT-based methods exist that reliably predict electronic structure given the correct geometry. Similarly, when good spectroscopic data are available, these same methods may, in principle, be used as input to the inverse problem of generating a good structural model. The same is generally true for gas-phase systems, for which the choice of method is different, but factors that guide its selection are known. Despite these successes, there are shortcomings associated with DFT for the prediction of materials' electronic structure. The present paper offers a perspective on these shortcomings. Fundamentally, the shortcomings associated with DFT stem from a lack of knowledge of the exact functional form of the exchange–correlation functional. Inaccuracies therefore arise from using an approximate functional. These inaccuracies can be reduced by judicious selection of the approximate functional. Other apparent shortcomings present due to misuse or improper application of the method. One of the most significant difficulties is the lack of a robust method for predicting electronic and geometric structure when only qualitative (connectivity) information is available about the system/material. Herein, some actual shortcomings of DFT are distinguished from merely common improper applications of the method. The role of the exchange functional in the predicted relationship between geometric structure and band gap is then explored, using fullerene, 2D polymorphs of elemental phosphorus and polyacetylene as case studies. The results suggest a potentially fruitful avenue of investigation by which some of the true shortcomings might be overcome, and serve as the basis for an appeal for high-accuracy experimental structure data to drive advances in theory.
In order to study the effects of Ni oxidation barriers on H diffusion in Zr, a Ni-Zr-H potential was developed based on an existing Ni-Zr potential. Using this and existing binary potentials H diffusion characteristics were calculated and some limited findings for the performance of Ni on Zr coatings are made.
Austenitic stainless steels (Fe-Cr-Ni) are resistant to hydrogen embrittlement but have not been studied using molecular dynamics simulations due to the lack of an Fe-Cr-Ni-H interatomic potential. Herein we describe our recent progress towards molecular dynamics studies of hydrogen effects in Fe-Cr-Ni stainless steels. We first describe our Fe-Cr-Ni-H interatomic potential and demonstrate its characteristics relevant to mechanical properties. We then demonstrate that our potential can be used in molecular dynamics simulations to derive Arrhenius equation of hydrogen diffusion and to reveal twinning and phase transformation deformation mechanisms in stainless steels.
Ionic liquids are a unique class of materials with several potential applications in electrochemical energy storage. When used in electrolytes, these highly coordinating solvents can influence device performance through their high viscosities and strong solvation behaviors. In this work, we explore the effects of pyrrolidinium cation structure and Li+ concentration on transport processes in ionic liquid electrolytes. We present correlated experimental measurements and molecular simulations of Li+ mobility and O2 diffusivity, and connect these results to dynamic molecular structural information and device performance. In the context of Li-O2/Li-air battery chemistries, we find that Li+ mobility is largely influenced by Li+-anion coordination, but that both Li+ and O2 diffusion may be affected by variations of the pyrrolidinium cation and Li+ concentration.
Fe-Ni-Cr stainless-steels are important structural materials because of their superior strength and corrosion resistance. Atomistic studies of mechanical properties of stainless-steels, however, have been limited by the lack of high-fidelity interatomic potentials. Here using density functional theory as a guide, we have developed a new Fe-Ni-Cr embedded atom method potential. We demonstrate that our potential enables stable molecular dynamics simulations of stainless-steel alloys at high temperatures, accurately reproduces the stacking fault energy—known to strongly influence the mode of plastic deformation (e.g., twinning vs. dislocation glide vs. cross-slip)—of these alloys over a range of compositions, and gives reasonable elastic constants, energies, and volumes for various compositions. The latter are pertinent for determining short-range order and solute strengthening effects. Our results suggest that our potential is suitable for studying mechanical properties of austenitic and ferritic stainless-steels which have vast implementation in the scientific and industrial communities. Published 2018. This article is a U.S. Government work and is in the public domain in the USA.
The efficiency of dye-sensitized solar cells (DSSCs) is strongly influenced by dye molecule orientation and interactions with the substrate. Understanding the factors controlling the surface orientation of sensitizing organic molecules will aid in the improvement of both traditional DSSCs and other devices that integrate molecular linkers at interfaces. Here, we describe a general approach to understand relative dye-substrate orientation and provide analytical expressions predicting orientation. We consider the effects of substrate, solvent, and protonation state on dye molecule orientation. In the absence of solvent, our model predicts that most carboxylic acid-functionalized molecules prefer to lie flat (parallel) on the surface, due to van der Waals interactions, as opposed to a tilted orientation with respect to the surface that is favored by covalent bonding of the carboxylic acid group to the substrate. When solvation effects are considered, however, the molecules are predicted to orient perpendicular to the surface. We extend this approach to help understand and guide the orientation of metal-organic framework (MOF) thin-film growth on various metal-oxide substrates. A two-part analytical model is developed on the basis of the results of DFT calculations and ab initio MD simulations that predicts the binding energy of a molecule by chemical and dispersion forces on rutile and anatase TiO2 surfaces, and quantifies the dye solvation energy for two solvents. The model is in good agreement with the DFT calculations and enables rapid prediction of dye molecule and MOF linker binding preference on the basis of the size of the adsorbing molecule, identity of the surface, and the solvent environment. We establish the threshold molecular size, governing dye molecule orientation, for each condition.
Ni3(2,3,6,7,10,11-hexaiminotriphenylene)2 is a π-stacked layered metal-organic framework material with extended π-conjugation that is analogous to graphene. Published experimental results indicate that the material is semiconducting, but all theoretical studies to date predict the bulk material to be metallic. Given that previous experimental work was carried out on specimens containing complex nanocrystalline microstructures and the tendency for internal interfaces to introduce transport barriers, we apply DFT to investigate the influence of internal interface defects on the electronic structure of Ni3(HITP)2. The results show that interface defects can introduce a transport barrier by breaking the π-conjugation and/or decreasing the dispersion of the electronic bands near the Fermi level. We demonstrate that the presence of defects can open a small gap, in the range of 15-200 meV, which is consistent with the experimentally inferred hopping barrier.
Metal organic frameworks (MOFs) are extended, nanoporous crystalline compounds consisting of metal ions interconnected by organic ligands. Their synthetic versatility suggest a disruptive class of opto - electronic materials with a high degree of electrical tunability and without the property - degrading disorder of organic conductors. In this project we determined the factors controlling charge and energy transport in MOFs and evaluated their potential for thermoelectric energy conversion. Two strategies for a chieving electronic conductivity in MOFs were explored: 1) using redox active 'guest' molecules introduced into the pores to dope the framework via charge - transfer coupling (Guest@MOF), 2) metal organic graphene analogs (MOGs) with dispersive band structur es arising from strong electronic overlap between the MOG metal ions and its coordinating linker groups. Inkjet deposition methods were developed to facilitate integration of the guest@MOF and MOG materials into practical devices.
Metal-organic frameworks (MOFs) are highly ordered, functionally tunable supramolecular materials with the potential to improve dye-sensitized solar cells (DSSCs). Several recent reports have indicated that photocurrent can be generated in Grätzel-type DSSC devices when MOFs are used as the sensitizer. However, the specific role(s) of the incorporated MOFs and the potential influence of residual MOF precursor species on device performance are unclear. Herein, we describe the assembly and characterization of a simplified DSSC platform in which isolated MOF crystals are used as the sensitizer in a planar device architecture. We selected a pillared porphyrin framework (PPF) as the MOF sensitizer, taking particular care to avoid contamination from light-absorbing MOF precursors. Photovoltaic and electrochemical characterization under simulated 1-sun and wavelength-selective illumination revealed photocurrent generation that is clearly ascribable to the PPF MOF. Continued refinement of highly versatile MOF structure and chemistry holds promise for dramatic improvements in emerging photovoltaic technologies. (Figure Presented).
Li-air or Li-oxygen batteries promise significantly higher energies than existing commercial battery technologies, yet their development has been hindered by a lack of suitable electrolytes. In this article, we evaluate the physical properties of varied electrolyte compositions to form generalized criteria for electrolyte design. We show that oxygen transport through non-aqueous electrolytes has a critical impact on the discharge rate and capacity of Li-air batteries. Through experiments and molecular dynamics simulations, we highlight that the choice of salt species and concentration have an outsized influence on oxygen solubility, while solvent choice is the major influence on oxygen diffusivity. The stability of superoxide reaction intermediates, key to the oxygen reduction mechanism, is also affected by variations in salt concentration and the choice of solvent. The importance of reactant transport is confirmed through Li-air cell discharge, which demonstrates good agreement between the observed and calculated mass transport-limited currents. These results showcase the impact of electrolyte composition on transport in metal-air batteries and provide guiding principles and simulation-based tools for future electrolyte design.
TlBr has the properties to become the leading radiation detection semiconductor. It has not yet been deployed due to a short lifetime of only hours to weeks. While the rapid structural deteriorations must come from ionic conduction under operating electrical fields, detailed aging mechanisms have not been understood. As a result, progress to extend lifetime has been limited despite extensive studies in the past. We have developed new atomistic simulation capabilities to enable study of ionic conduction under electrical fields. Our combined simulations and experiments indicate that dislocations in TlBr climb under electrical fields. This climb is the root cause for structural deterioration. Hence, we discovered new strengthening methods to reduce aging. Our new atomistic simulation approach can have broader impact on other Sandia programs including battery research. Our project results in 4 publications, a new invention, new LAMMPS capabilities, solution to mission relevant materials, and numerous presentations.
As the world transitions from fossil fuels to clean energy sources in the coming decades, many technological challenges will require chemists and material scientists to develop new materials for applications related to energy conversion, storage, and efficiency. Because of their unprecedented adaptability, metal-organic frameworks (MOFs) will factor strongly in this portfolio. By utilizing the broad synthetic toolkit provided by the fields of organic and inorganic chemistry, MOF pores can be customized to suit a particular application. Of particular importance is the ability to tune the strength of the interaction between the MOF pores and guest molecules. By cleverly controlling these MOF-guest interactions, the chemist may impart new function into the Guest@MOF materials otherwise lacking in vacant MOF. Herein, we highlight the concept of the Guest@MOF as it relates to our efforts to develop these materials for energy-related applicatons. Our work in the areas of H2 and noble gas storage, hydrogenolysis of biomass, light-harvesting, and conductive materials will be discussed. Of relevance to light-harvesting applications, we report for the first time a postsynthetic modification strategy for increasing the loading of a light-sensitive electron-donor molecule in the pores of a functionalized MIL-101 structure. Through the demonstrated versatility of these approaches, we show that, by treating guest molecules as integral design elements for new MOF constructs, MOF science can have a significant impact on the advancement of clean energy technologies.
The metal organic framework material Ni3(2,3,6,7,10,11 - hexaiminotriphenylene)2, (Ni3(HITP)2) is composed of layers of extended conjugated planes analogous to graphene. We carried out Density functional theory (DFT) calculations to model the electronic structure of bulk and monolayer Ni3(HITP)2. The layered 3D material is metallic, similar to graphene. Our calculations predict that there is appreciable band dispersion not only in-plane, but perpendicular to the stacking planes as well, suggesting that, unlike graphene, the conductivity may be nearly isotropic. In contrast, a 2D monolayer of the material exhibits a band gap, consistent with previously published results. Insight obtained from studies of the evolution of the material from semiconducting to metallic as the material is transitioned from 2D to 3D suggests the possibility of modifying the material to render it semiconducting by changing the metal center and inserting spacer moieties between the layers. Furthermore, the DFT calculations predict that the modified material will be structurally stable and exhibit a band gap.
Carbon nanostructures, such as nanotubes and graphene, are of considerable interest due to their unique mechanical and electrical properties. The materials exhibit extremely high strength and conductivity when defects created during synthesis are minimized. Atomistic modeling is one technique for high resolution studies of defect formation and mitigation. To enable simulations of the mechanical behavior and growth mechanisms of C nanostructures, a high-fidelity analytical bond-order potential for the C is needed. To generate inputs for developing such a potential, we performed quantum mechanical calculations of various C structures.
Current austenitic stainless steel storage reservoirs for hydrogen isotopes (e.g. deuterium and tritium) have performance and operational life-limiting interactions (e.g. embrittlement) with H-isotopes. Aluminum alloys (e.g.AA2219), alternatively, have very low H-isotope solubilities, suggesting high resistance towards aging vulnerabilities. This report summarizes the work performed during the life of the Lab Directed Research and Development in the Nuclear Weapons investment area (165724), and provides invaluable modeling and experimental insights into the interactions of H isotopes with surfaces and bulk AlCu-alloys. The modeling work establishes and builds a multi-scale framework which includes: a density functional theory informed bond-order potential for classical molecular dynamics (MD), and subsequent use of MD simulations to inform defect level dislocation dynamics models. Furthermore, low energy ion scattering and thermal desorption spectroscopy experiments are performed to validate these models and add greater physical understanding to them.
In this project we developed t he atomistic models needed to predict how graphene grows when carbon is deposited on metal and semiconductor surfaces. We first calculated energies of many carbon configurations using first principles electronic structure calculations and then used these energies to construct an empirical bond order potentials that enable s comprehensive molecular dynamics simulation of growth. We validated our approach by comparing our predictions to experiments of graphene growth on Ir, Cu and Ge. The robustness of ou r understanding of graphene growth will enable high quality graphene to be grown on novel substrates which will expand the number of potential types of graphene electronic devices.
TlBr is promising for g- and x- radiation detection, but suffers from rapid performance degradation under the operating external electric fields. To enable molecular dynamics (MD) studies of this degradation, we have developed a Stillinger-Weber type of TlBr interatomic potential. During this process, we have also addressed two problems of wider interests. First, the conventional Stillinger-Weber potential format is only applicable for tetrahedral structures (e.g., diamond-cubic, zinc-blende, or wurtzite). Here we have modified the analytical functions of the Stillinger-Weber potential so that it can now be used for other crystal structures. Second, past modifications of interatomic potentials cannot always be applied by a broad community because any new analytical functions of the potential would require corresponding changes in the molecular dynamics codes. Here we have developed a polymorphic potential model that simultaneously incorporates Stillinger-Weber, Tersoff, embedded-atom method, and any variations (i.e., modified functions) of these potentials. As a result, we have implemented this polymorphic model in MD code LAMMPS, and demonstrated that our TlBr potential enables stable MD simulations under external electric fields.
Graphene films grown by vapour deposition tend to be polycrystalline due to the nucleation and growth of islands with different in-plane orientations. Here, using low-energy electron microscopy, we find that micron-sized graphene islands on Ir(111) rotate to a preferred orientation during thermal annealing. We observe three alignment mechanisms: the simultaneous growth of aligned domains and dissolution of rotated domains, that is, â €- ripeningâ €™; domain boundary motion within islands; and continuous lattice rotation of entire domains. By measuring the relative growth velocity of domains during ripening, we estimate that the driving force for alignment is on the order of 0.1â €‰meV per C atom and increases with rotation angle. A simple model of the orientation-dependent energy associated with the moiré corrugation of the graphene sheet due to local variations in the graphene-substrate interaction reproduces the results. This work suggests new strategies for improving the van der Waals epitaxy of 2D materials.
Metal-organic frameworks (MOFs) are crystalline nanoporous materials comprised of organic electron donors linked to metal ions by strong coordination bonds. Applications such as gas storage and separations are currently receiving considerable attention, but if the unique properties of MOFs could be extended to electronics, magnetics, and photonics, the impact on material science would greatly increase. Recently, we obtained "emergent properties," such as electronic conductivity and energy transfer, by infiltrating MOF pores with "guest" molecules that interact with the framework electronic structure. In this Perspective, we define a path to emergent properties based on the Guest@MOF concept, using zinc-carboxylate and copper-paddlewheel MOFs for illustration. Energy transfer and light harvesting are discussed for zinc carboxylate frameworks infiltrated with triplet-scavenging organometallic compounds and thiophene- and fullerene-infiltrated MOF-177. In addition, we discuss the mechanism of charge transport in TCNQ-infiltrated HKUST-1, the first MOF with electrical conductivity approaching conducting organic polymers. These examples show that guest molecules in MOF pores should be considered not merely as impurities or analytes to be sensed but also as an important aspect of rational design.
Electro-optical organic materials hold great promise for the development of high-efficiency devices based on exciton formation and dissociation, such as organic photovoltaics (OPV) and organic light-emitting devices (OLEDs). However, the external quantum efficiency (EQE) of both OPV and OLEDs must be improved to make these technologies economical. Efficiency rolloff in OLEDs and inability to control morphology at key OPV interfaces both reduce EQE. Only by creating materials that allow manipulation and control of the intimate assembly and communication between various nanoscale excitonic components can we hope to first understand and then engineer the system to allow these materials to reach their potential. The aims of this proposal are to: 1) develop a paradigm-changing platform for probing excitonic processes composed of Crystalline Nanoporous Frameworks (CNFs) infiltrated with secondary materials (such as a complimentary semiconductor); 2) use them to probe fundamental aspects of excitonic processes; and 3) create prototype OPVs and OLEDs using infiltrated CNF as active device components. These functional platforms will allow detailed control of key interactions at the nanoscale, overcoming the disorder and limited synthetic control inherent in conventional organic materials. CNFs are revolutionary inorganic-organic hybrid materials boasting unmatched synthetic flexibility that allow tuning of chemical, geometric, electrical, and light absorption/generation properties. For example, bandgap engineering is feasible and polyaromatic linkers provide tunable photon antennae; rigid 1-5 nm pores provide an oriented, intimate host for triplet emitters (to improve light emission in OLEDs) or secondary semiconducting polymers (creating a charge-separation interface in OPV). These atomically engineered, ordered structures will enable critical fundamental questions to be answered concerning charge transport, nanoscale interfaces, and exciton behavior that are inaccessible in disordered systems. Implementing this concept also creates entirely new dimensions for device fabrication that could both improve performance, increase durability, and reduce costs with unprecedented control of over properties. This report summarizes the key results of this project and is divided into sections based on publications that resulted from the work. We begin in Section 2 with an investigation of light harvesting and energy transfer in a MOF infiltrated with donor and acceptor molecules of the type typically used in OPV devices (thiophenes and fullerenes, respectively). The results show that MOFs can provide multiple functions: as a light harvester, as a stabilizer and organizer or the infiltrated molecules, and as a facilitator of energy transfer. Section 3 describes computational design of MOF linker groups to accomplish light harvesting in the visible and facilitate charge separation and transport. The predictions were validated by UV-visible absorption spectroscopy, demonstrating that rational design of MOFs for light-harvesting purposes is feasible. Section 4 extends the infiltration concept discussed in Section to, which we now designate as "Molecule@MOF" to create an electrically conducting framework. The tailorability and high conductivity of this material are unprecedented, meriting publication in the journal Science and spawning several Technical Advances. Section 5 discusses processes we developed for depositing MOFs as thin films on substrates, a critical enabling technology for fabricating MOF-based electronic devices. Finally, in Section 6 we summarize results showing that a MOF thin film can be used as a sensitizer in a DSSC, demonstrating that MOFs can serve as active layers in excitonic devices. Overall, this project provides several crucial proofs-of- concept that the potential of MOFs for use in optoelectronic devices that we predicted several years ago [ 3 ] can be realized in practice.
Biological imaging and assay technologies rely on fluorescent organic dyes as reporters for a number of interesting targets and processes. However, limitations of organic dyes such as small Stokes shifts, spectral overlap of emission signals with native biological fluorescence background, and photobleaching have all inhibited the development of highly sensitive assays. To overcome the limitations of organic dyes for bioassays, we propose to develop lanthanide-based luminescent dyes and demonstrate them for molecular reporting applications. This relatively new family of dyes was selected for their attractive spectral and chemical properties. Luminescence is imparted by the lanthanide atom and allows for relatively simple chemical structures that can be tailored to the application. The photophysical properties offer unique features such as narrow and non-overlapping emission bands, long luminescent lifetimes, and long wavelength emission, which enable significant sensitivity improvements over organic dyes through spectral and temporal gating of the luminescent signal.Growth in this field has been hindered due to the necessary advanced synthetic chemistry techniques and access to experts in biological assay development. Our strategy for the development of a new lanthanide-based fluorescent reporter system is based on chelation of the lanthanide metal center using absorbing chromophores. Our first strategy involves "Click" chemistry to develop 3-fold symmetric chelators and the other involves use of a new class of tetrapyrrole ligands called corroles. This two-pronged approach is geared towards the optimization of chromophores to enhance light output.