Search NASA⌕ Search

SEARCH · Search NASA

Results for “EXTRAPOLATION”

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 1,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Reliable Potential Energy Surfaces for the Reactions of H 2 O with ThO 2 , PaO 2 + , UO 2 2+ , and UO 2 +

The potential energy surfaces for the reactions of H 2 O with ThO 2 , PaO 2 + , UO 2 2+ , and UO 2 + have been calculated at the coupled cluster CCSD(T) level extrapolated to the complete basis set limit with additional corrections including scalar relativistic and spin-orbit. The reactions proceed by the formation of an initial Lewis acid-base adduct (H 2 O)AnO 2 0/+/2+ followed by a proton transfer to generate the dihydroxide AnO(OH) 2 0/+/2+ . The results are in excellent agreement with mass spectrometry experiments and prior calculations of hydrolysis reactions of the group 4 transition metal dioxides MO 2 . The differences in the energies of the stationary points on the potential energy surface are explained in terms of the charges on the system and the populations on the metal center. The use of an improved starting point for the coupled cluster CCSD(T) calculations based on density functional theory with the PW91 exchange-correlation functional or Brueckner orbitals is described. The importance of including second-order spin-orbit corrections for closed-shell molecules is also described. These improvements in the calculations are then correlated with the 5f populations on the actinide.

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

Extension of structure-reactivity correlations for the hydrogen abstraction reaction by bromine atom and comparison to chlorine atom and hydroxyl radical

Recently we presented structure-reactivity correlations for the gas-phase ambient-temperature rate constants for hydrogen abstraction from sp 3 -hybridized carbon by chlorine atom and hydroxyl radical (Cl•/HO• + HCR 3 → HCl/HOH + •CR 3 ); the reaction enthalpy effect was represented by the independent variable Δ r H and the polar effect by the independent variables F and R, the Hammett constants for field/inductive and resonance effects. Both these reactions are predominantly exothermic and have early transition states. Here we present a parallel treatment for Br• whose reaction is significantly endothermic with a correspondingly late transition state. In spite of lower expectations because the available data base is less extensive and much more scattered and because long temperature extrapolations are often required, the resulting least-squares fit (log k 298,Br = –0.147 Δ r H –4.32 ΣF –4.28 ΣR –12.38 with r 2 = 0.92) was modestly successful and useful for initial predictions. The coefficient of Δ r H was ~4-fold greater, indicative of the change from an early to a late transition state; meanwhile the sizable coefficients of ΣF and ΣR indicate the persistence of the polar effect. Although the mean unsigned deviation of 0.79 log k 298 units is rather large, it must be considered in the context of a total span of over 15 log units in the data set. Lastly, the major outliers are briefly discussed.

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

The Reaction Rates of Amidogen and Ammonia with Nitrous Oxide: Implications for Combustion Mechanisms

Pulsed laser photolysis experiments with laser-induced fluorescence detection of NH 2 set an upper limit to the rate constant for reaction with N 2 O of k < 1 × 10 −15 cm 3 molecule −1 s −1 at 513 K. Computations were based on geometries and anharmonic frequency analysis (B2PLYP-D3/cc-pVTZ) followed by coupled cluster calculations extrapolated to the infinite basis set limit, with corrections for core−valence electron correlation, scalar relativistic effects, and correlation up to CCSDT(Q). Species that showed multireference character were quantified with MRCI(7,7)+Q/cc-pVTZ theory. Rate constants were obtained for the dominant product channel H 2 NN + NO, along with HN 3 + OH, H 2 NO + N 2 , NNH + HNO, ON(NH) 2 and HNNH + NO. The last channel is slow even at 2500 K, contrary to an early empirical estimate and confirming recent suggestions. Modeling of literature experiments on oxidation of NH 3 by N 2 O shows that all channels are too slow to make a significant impact on the loss of N 2 O in ammonia flames. Similarly, the direct NH 3 + N 2 O reaction is found to be negligibly slow.

Ammonia↗
DOE OSTI2025

The Remarkable I 2 O 3 Molecule: A New View from Theory

Atmospheric iodine chemistry has garnered increasing attention as a result of increased iodine emissions. A key subset of this chemistry involves iodine oxides (I 2 O 2–5 ), which serve as precursors to particle formation. Among these, I 2 O 3 is the simplest iodine oxide involved in particle formation, but it has remained undetected in the atmosphere. Previous theoretical studies have characterized this peculiar molecule, primarily using energies to refine geometries obtained at low levels of theory. Due to the reemerging interest in I 2 O 3 , this study presents geometries optimized at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory─marking the first instance, to the best of our knowledge, where this system has been studied exclusively with CCSD(T). Harmonic vibrational frequencies were computed at the same level of theory. Final energetics were obtained using the very high level CCSDT(Q) method with basis sets up to quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) and extrapolated to the CBS limit to yield CCSDT(Q)/CBS//CCSD(T)/aug-cc-pwCVTZ-PP energies. These energies include harmonic zero-point vibrational energy corrections and scalar relativistic energy corrections. Additionally, this study discovers new isomers along the I 2 O 3 potential energy surface, a novel contribution to the field. The performance of different computational methods and DFT functionals commonly used in atmospheric chemistry is also assessed relative to high-level theoretical methods.

basis sets↗
DOE OSTI2025

Vibrational Signatures of Unrealized Phosphorus Suboxide Intermediates in White Phosphorus Oxidation Reactions

White phosphorus ignition notoriously produces the phosphoric acid anhydride P 4 O 10 , yet the intermediate oxidation steps remain undetermined. Here, we report the first geometric and vibrational characterization of two P 4 O 2 isomers, P 3 OPO and P 3 PO 2 , and substantiate a previously proposed cyclic P 4 O 2 isomer. We formally assign the infrared bands observed by Andrews and Mielke at 898 and 891 cm −1 to the antisymmetric P− O−P vibrations of P 3 OPO species. Additional bands corresponding to terminal P=O and −PO 2 stretches of P 3 OPO and P 3 PO 2 discussed herein also went unrecognized due to the peculiar bonding of P 4 O 2 species compared to oxo-bridged P 4 O x (x = 3−6) species. Sequential addition of oxygen atoms to the P 4 tetrahedron appears to form P 3 OPO and P 3 PO 2 , while cyclic P 4 O 2 is formed from P 2 O dimerization. CCSD(T) geometries, CCSD(T) + MP2[δVPT2] fundamental frequencies, and enthalpies of formation extrapolated using focal-point analysis are reported. The predicted enthalpies of formation relative to tetrahedral P 4 plus 3 Σ g − O 2 for bent-P 3 OPO, extended-P 3 OPO, P 3 PO 2 , and cyclic-P 4 O 2 are −93, −88, −85, and −63 kcal mol −1 , respectively.

Molecular structure↗
DOE OSTI2026

Redox Potentials for Hydrogen Species in Acetonitrile and Tetrahydrofuran

The properties of the fundamental hydrogen species, proton (H+), hydrogen atom (H•), and hydride anion (H-) are critical to a vast range of chemical processes, yet their thermodynamic properties in nonaqueous solvents are not well established. A hybrid supermolecule-continuum approach is used to predict the Gibbs free energies of solvation (?G°solv) and standard redox potentials (E°) for the 2H+/H2 and H•/H- couples in acetonitrile (MeCN) and tetrahydrofuran (THF) following the approach previously used for water. Gas phase geometry optimizations were done using density functional theory (DFT) with the ?B97XD functional, Møller-Plesset second-order perturbation theory (MP2) with augmented correlation-consistent basis sets and composite-correlated molecular orbital theory (G3MP2 and G3MP2B3). Additional single point calculations in gas phase were performed for H- affinity for MeCN and THF at the coupled cluster CCSD(T) level using MP2 geometries. Solvation was included using the self-consistent reaction field SMD model. For H+, the ?G°solv values are predicted to be -252.2 kcal/mol in MeCN and -261.0 kcal/mol in THF. The redox potentials for the 2H+/H2 couple are +0.49 V in MeCN and +0.11 V in THF relative to the aqueous SHE, consistent with the available experimental data. The ?G°solv for H- is predicted to be -80 kcal/mol in MeCN and -68 in THF, and for H•, ?G°solv is predicted to be -1.3 kcal/mol in THF and -2.0 kcal/mol in MeCN. These solvation energies yield respective calculated redox potentials for the H•/H- couple of -0.14 V in MeCN and -0.68 V in THF relative to the aqueous SHE. H2 solvation is slightly exergonic in MeCN (-3.3 kcal/mol) and near thermoneutral in THF (+1.3 kcal/mol). The pKa(H2) is predicted to be 45.4 in MeCN and 47.8 in THF. The proton-coupled electron transfer (PCET) accounting term CG was evaluated giving values on the ferrocene scale within approximately 2 kcal/mol of literature in both solvents. These results provide benchmarks for hydrogen redox chemistry in MeCN and THF and improve consistency across E°, ?G°solv, pKa, and CG. These results provide benchmarks for hydrogen redox chemistry in nonaqueous media challenging previous estimates based on indirect extrapolations and assumptions.

Duda, Damian P.↗
DOE OSTI2020

Thermodynamics of Metal Carbonates and Bicarbonates and Their Hydrates for Mg, Ca, Fe, and Cd Relevant to Mineral Energetics

The heats of formation of carbonate, bicarbonate and bicarbonate/hydroxide metal complexes and hydrated versions these complexes of Mg 2+ , Ca 2+ , Fe 2+ , and Cd 2+ are predicted from atomization energies using correlated molecular orbital theory at the CCSD(T) level extrapolated to the complete basis set limit following the Feller-Peterson-Dixon (FPD) approach. Using the calculated gas phase values and the available experimental solid state values, the cohesive energies of selected minerals were calculated. The gas phase decomposition energies into MO, CO 2 and H 2 O follow the order Mg ~ Ca > Cd ~ Fe, and correlate with the hardness of the metal +2 ions. Gas phase hydration energies show that the order is Mg > Fe > Ca ~ Cd. There are a number of bulk hydrated Mg and Ca complexes that occur as minerals but there are few if any for Fe and Cd, suggesting that a number of factors are important in determining the stability of the bulk mineral hydrates. The FPD heats of formation were used to benchmark a range of density functional theory exchange-correlation functionals, including those commonly used in solid state mineral calculations. None of the functionals provided chemical accuracy agreement (± 1 kcal/mol) with the FPD results The best functionals are ωB97X and ωB97X-D with mean average unsigned errors of 10 kcal/mol.

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

Molecular theory and the effects of solute attractive forces on hydrophobic interactions

The role of solute attractive forces on hydrophobic interactions is studied by coordinated development of theory and simulation results for Ar atoms in water. In this paper, we present a concise derivation of the local molecular field (LMF) theory for the effects of solute attractive forces on hydrophobic interactions, a derivation that clarifies the close relation of LMF theory to the EXP approximation applied to this problem long ago. The simulation results show that change from purely repulsive atomic solute interactions to include realistic attractive interactions diminishes the strength of hydrophobic bonds. For the Ar–Ar rdfs considered pointwise, the numerical results for the effects of solute attractive forces on hydrophobic interactions are opposite in sign and larger in magnitude than predicted by LMF theory. That comparison is discussed from the point of view of quasichemical theory, and it is suggested that the first reason for this difference is the incomplete evaluation within LMF theory of the hydration energy of the Ar pair. With a recent suggestion for the system-size extrapolation of the required correlation function integrals, the Ar–Ar rdfs permit evaluation of osmotic second virial coefficients B 2 . Those B 2 ’s also show that incorporation of attractive interactions leads to more positive (repulsive) values. With attractive interactions in play, B 2 can change from positive to negative values with increasing temperatures. Furthermore, this is consistent with the puzzling suggestions of decades ago that B 2 ≈ 0 for intermediate cases of temperature or solute size. In all cases here, B 2 becomes more attractive with increasing temperature.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Self-Assembly and Critical Aggregation Concentration Measurements of ABA Triblock Copolymers with Varying B Block Types: Model Development, Prediction, and Validation

The dissipative particle dynamics (DPD) simulation technique is a coarse-grained (CG) molecular dynamics-based approach that can effectively capture the hydrodynamics of complex systems while retaining essential information about the structural properties of the molecular species. An advantageous feature of DPD is that it utilizes soft repulsive interactions between the beads, which are CG representation of groups of atoms or molecules. In this study, we used the DPD simulation technique to study the aggregation characteristics of ABA triblock copolymers in aqueous medium. Pluronic polymers (PEG-PPO-PEG) were modeled as two segments of hydrophilic beads and one segment of hydrophobic beads. Tyrosine-derived PEG 5K -b-oligo(desaminotyrosyl tyrosine octyl ester-suberate)-b-PEG 5K (PEG5K-oligo(DTO-SA)-PEG 5K ) block copolymers possess alternate rigid and flexible components along the hydrophobic oligo(DTO-SA) chain, and were modeled as two segments of hydrophilic beads and one segment of hydrophobic, alternate soft and hard beads. The formation, structure, and morphology of the initial aggregation of the polymer molecules in aqueous medium were investigated by following the aggregation dynamics. The dimensions of the aggregates predicted by the computational approach were in good agreement with corresponding results from experiments, for the Pluronic and PEG 5K -oligo(DTO-SA)-PEG 5K block copolymers. In addition, DPD simulations were utilized to determine the critical aggregation concentration (CAC), which was compared with corresponding results from an experimental approach. For Pluronic polymers F68, F88, F108, and F127, the computational results agreed well with experimental measurements of the CAC measurements. For PEG 5K -b-oligo(DTO-SA)-b-PEG 5K block polymers, the complexity in polymer structure made it difficult to directly determine their CAC values via the CG scheme. Therefore, we determined CAC values of a series of triblock copolymers with 3–8 DTO-SA units using DPD simulations, and used these results to predict the CAC values of triblock copolymers with higher molecular weights by extrapolation. In parallel, a PEG 5K -b-oligo(DTO-SA)-b-PEG 5K block copolymer was synthesized, and the CAC value was determined experimentally using the pyrene method. The experimental CAC value agreed well with the CAC value predicted by simulation. Furthermore, these results validate our CG models, and demonstrate an avenue to simulate and predict aggregation characteristics of ABA amphiphilic triblock copolymers with complex structures.

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

Coarse-Grained Molecular Simulation of the Hierarchical Self-Assembly of π-Conjugated Optoelectronic Peptides

Self-assembled aggregates of peptides containing aromatic groups possess optoelectronic properties that make them attractive targets for the fabrication of biocompatible electronics. Molecular-level understanding of how the microscopic peptide chemistry influences the properties of the aggregates is vital for rational peptide design. We construct a coarse-grained model of Asp-Phe-Ala-Gly-OPV3-Gly-Ala-Phe-Asp (DFAG-OPV3-GAFD) peptides containing OPV3 (distyrylbenzene) π-conjugated cores explicitly parameterized against all-atom calculations and perform molecular dynamics simulations of the self-assembly of hundreds of molecules over hundreds of nanoseconds. We observe a hierarchical assembly mechanism wherein ~2-8 peptides assemble into stacks with aligned aromatic cores that subsequently form elliptical aggregates and ultimately a branched network with a fractal dimensionality of ~1.5. The assembly dynamics are well described by a Smoluchowski coagulation process for which we extract rate constants from the molecular simulations to both furnish insight into the microscopic assembly kinetics and extrapolate our aggregation predictions to time and length scales beyond the reach of molecular simulation. Lastly, this study presents new molecular-level understanding of the morphology and dynamics of the spontaneous self-assembly of DFAG-OPV3-GAFD peptides and establishes a systematic protocol to develop coarse-grained models of optoelectronic peptides for the exploration and design of π-conjugated peptides with tunable optoelectronic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Disruption of Hydrogen-Bonding Network Eliminates Water Anomalies Normally Observed on Cooling to Its Calorimetric Glass Transition

Here, the calorimetric glass-transition temperature of water is 136 K, but extrapolation of thermodynamic and relaxation properties of water from ambient temperature to below its homogeneous nucleation temperature T H = 235 K predicts divergence at T S = 228 K. The “no-man’s land” between the T H and glassy water crystallization temperature of 150 K, which is encountered on warming up from the vitrified state, precludes a straightforward reconciliation of the two incompatible temperature dependences of water properties, above 235 K and below 150 K. The addition of lithium chloride to water allows bypassing both T H and T S on cooling, resulting in the dynamics with no features except the calorimetric glass transition, still at 136 K. We show that lithium chloride prevents hydrogen-bonding network completion in water on cooling, as manifested, in particular, in changing microscopic diffusion mechanism of the water molecules. Thus thermodynamic and relaxation peculiarities exhibited by pure water on cooling to its glass transition, such as the existence of the T H and T S , must be associated specifically with the hydrogen-bonding network.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Finding Order in the Disordered Hydration Shell of Rapidly Exchanging Water Molecules around the Heaviest Alkali Cs + and Fr +

Here, we report the structural and dynamical characterization of the intrinsically disordered hydration shells of the heaviest alkali ions, Cs + and Fr + , obtained in ab initio molecular dynamics simulations. The knowledge of solvation and complexation properties of short-lived Fr + is very limited and mostly based on extrapolations from the smaller alkali metal ions. To this end, we provide a critical insight into Fr + solvation, demonstrating an extreme example of disordered solvation with no distinction between the ion-bound and solvent-bound states of water based on the ion–water distance. However, these two states are distinguished through distance-solvent rearrangement correlation, where either coordination number or electric field is employed to treat solvent rearrangement. Utilizing reaction rate theory, we find that the water exchange time scale for Fr + (2.1–2.3 ps) is unexpectedly slower than for Cs + (0.5–1.2 ps), because Fr + experiences stronger nonequilibrium solvent effects. This study provides a new perspective on weak and hydrophobic solvation.

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

Ion–Molecule Complex Dissociation and Formation Dynamics in LiCl Aqueous Solutions from 2D IR Spectroscopy

Ion–molecule complex dynamics as well as water dynamics in concentrated lithium chloride (LiCl) solutions are examined using ultrafast two-dimensional infrared (2D IR) spectroscopy with the CN stretching mode of methyl thiocyanate (MeSCN) as the vibrational probe. In pure water, MeSCN has a narrow symmetric absorption line shape. 2D IR spectral diffusion measurements of the CN stretch give the identical time dependence of water dynamics, as previously observed using the OD stretch of HOD in H 2 O. In concentrated LiCl solutions, the IR absorption spectrum of MeSCN displays two distinct peaks, one corresponding to water H-bonded to the N lone pair of MeSCN (W) and the other corresponding to Li + associated with the N (L). These two species are in equilibrium, and switching of the CN bonding partner from Li + to H 2 O and vice versa was observed and explicated with 2D IR chemical exchange spectroscopy. The MeSCN·Li + complex dissociation time constant, τ LW , and the MeSCN·H 2 O dissociation time constant, τ WL , were determined. The observed τ LW chemical exchange dissociation time constant changes from 60 to 40 ps as the LiCl concentration decreases from ~10.7 to ~7.7 M, mainly due to the increase of the water concentration as the LiCl concentration is reduced. The observed time constants are independent of the model for the chemical reaction. With the assumption of a simple chemical equation, MeSCN·Li + + H 2 O ⇌ MeSCN·H 2 O + Li + , the equilibrium equation rate constants were obtained from the observed chemical exchange time constants. It was determined that the equilibrium rate constants barely change even though the viscosity changes by a factor of 2 and the ionic strength changes by a factor of 1.4. Extrapolation to dilute LiCl solution estimates the τ LW to be ~30 ps. The orientational relaxation (anisotropy decay) of both the W and L complexes was measured using polarization selective 2D IR experiments. The lithium-bonded species undergoes orientational relaxation ~3 times slower than the water-bonded species in each LiCl solution studied. In conclusion, the difference demonstrates the distinct interactions with the medium experienced by the neutral and charged species in the concentrated salt solutions.

14 SOLAR ENERGY↗
DOE OSTI2020

Predicting Elastic Properties of Materials from Electronic Charge Density Using 3D Deep Convolutional Neural Networks

Materials representation plays a key role in machine learning-based prediction of materials properties and new materials discovery. Currently both graph and three-dimensional (3D) voxel representation methods are based on the heterogeneous elements of the crystal structures. Here, we propose to use electronic charge density (ECD) as a generic unified 3D descriptor for materials property prediction with the advantage of possessing close relation with the physical and chemical properties of materials. We developed an ECD-based 3D convolutional neural networks (CNNs) for predicting the elastic properties of materials, in which CNNs can learn effective hierarchical features with multiple convolving and pooling operations. Extensive benchmark experiments over 2170 $Fm\bar3m$ face-centered-cubic materials show that our ECD-based CNNs can achieve good performance for elasticity prediction. Especially, our CNN models based on the fusion of elemental Materials-Agnostic Platform for Informatics and Exploration features and ECD descriptors achieved the best fivefold cross-validation performance. More importantly, we showed that our ECD-based CNN models can achieve significantly better extrapolation performance when evaluated over nonredundant data sets, where there are few neighbor-training samples around test samples. As an additional validation, we evaluated the predictive performance of our models on 329 materials of space group $Fm\bar3m$ by comparing to density functional theory calculated values, which shows a better prediction power of our model for bulk modulus than shear modulus. Because of the unified representation power of ECD, it is expected that our ECD-based CNN approach can also be applied to predict other physical and chemical properties of crystalline materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling Hydrogen Evolution Reaction Kinetics through Explicit Water–Metal Interfaces

Despite the apparent simplicity of the hydrogen evolution reaction (HER) and the decades of research into it, controversy remains in the literature regarding the identity of the active site and the competition between the Heyrovsky and Tafel steps. In this work, we use charge-extrapolated ab initio simulations with explicit water in conjunction with mean-field microkinetic modeling to explore the mechanism for HER on both close-packed (111) and stepped (211) transition metals. First, we show that atop H*, beyond a monolayer of hollow H*, is unlikely to play a role in the HER mechanism, given its very positive adsorption energies. The energetics suggests the Volmer–Heyrovsky mechanism to predominate on fcc transition metals under typical operating conditions. Here, we evaluate our theoretical results vs several experimental observations. We show that the Volmer–Heyrovsky mechanism predicts an activity volcano with its peak at a H* binding ΔG H* ≈ 0 eV, consistent with experiment. In contrast, the Volmer–Tafel volcano shows a broad rate plateau between ΔG H* ≈ 0 eV and ΔG H* ≈ – 0.4 eV. We find our theoretical Tafel slopes to be consistent with experimental ones on a range of transition metals. We show that, in line with experimental observations, the introduction of a CO(g) atmosphere shifts the strong binding metals toward the weak binding leg. Our study suggests that the simple thermodynamic approach to HER activity still holds, even when a detailed kinetic picture is considered.

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

Radiolysis of Thin Water Ice in Electron Microscopy

Little is known about the radiolysis of water ice, especially as compared to the radiolysis of liquid water. In this study, the radiolytic decomposition of thin water ice films is probed for an initial 80 keV electron beam using electron energy loss spectroscopy, EELS. Pre-peaks in the oxygen K-edge spectra give the relative intensities of most of the oxygen containing species produced by radiolytic water decomposition. Contrary to expectations from the extrapolation of liquid water radiolysis data to high dose rates where significant H 2 O 2 production is expected, the main molecular product observed in the oxygen K-edge EELS spectra of water ice is O 2 . Significant mass loss of water and decomposition of most of its radiolytic products are observed for higher exposures, and the O atom seems to have a major role in the subsequent chemistry. An inverse relationship between the formation of O 2 and of ·OH is observed. A new high-dose reaction scheme is proposed. Furthermore, the significance of our results with respect to the radiolysis of water ice and for electron microscopy studies is discussed.

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

Effect of Surface Stoichiometry on Blinking and Hole Trapping Dynamics in CdSe Nanocrystals

Here, we measure the photoinduced carrier dynamics as the surface composition of CdSe nanocrystals is systematically varied from metal rich (~80% surface Cd) to nearly stoichiometric (~50% surface Cd). Using time-resolved optical spectroscopy, we determine that the luminescence lifetime is controlled by the rate of hole trapping at the newly exposed surface selenium atoms. However, the increased rate of the photoluminescence decay is not sufficient to explain the decreased photoluminescence quantum yield, and requires a growing proportion of nanocrystals in a dark or “OFF” state to explain the data. A global kinetic model is proposed that relates the fraction of selenium sites to the rate of hole trapping. A linear relationship between the rate of hole trapping and the fraction of exposed Se sites (x Se ) is observed within the range of accessible stoichiometries (x Se = 0.5–0.2). Extrapolation to higher surface cadmium fractions suggests that not all Se sites serve as effective hole traps. These findings confirm the strong nonlinear dependence of the fluorescence yield on the nanocrystal stoichiometry.

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

Coupled Cluster and Density Functional Studies of Atomic Fluorine Chemisorption on Coronene as Model Systems for Graphene Fluorination

Potential energy curves (PECs) of atomic fluorine adsorption on coronene as a model for graphene or nanocarbon surfaces have been computed. The PECs were obtained by scanning the fluorine atom distance to one of the center carbon atoms of the coronene molecule as model system in a “top” position from 4.0 to1.0 Å in intervals of 0.1 Å using a variety of quantum chemical methods. Various density functional theory (DFT) functionals, such as B3LYP, PBE, PBE0, CAM-B3LYP, and LC-ωPBE; approximate DFT methods such as several levels of the density-functional tight-binding (DFTB) method, as well as ab initio wave function theory methods such as MP2, CCSD, CCSD(T), and G2MS extrapolations, were used to evaluate energies for B3LYP/cc-pVDZ PEC geometries. G2MS is an approximation to the highly accurate CCSD(T)/cc-pVTZ level of theory in our work. We found that fluorine is chemically adsorbed on coronene with a binding energy of 22.9 and 23.3 kcal/mol at the B3LYP/cc-pVDZ and G2MS levels of theory, respectively, and 18.3 and 19.3 kcal/mol after counterpoise correction. Here, we found that pure DFT functionals and their DFTB approximations fail to predict the correct dissociation limit due to the DFT-inherent self-interaction error and various limitations in the DFTB approximation itself.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗