Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT calculation”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Increasing Iridium Oxide Activity for the Oxygen Evolution Reaction with Hafnium Modification

Synthesis and implementation of highly active, stable, and affordable electrocatalysts for the oxygen evolution reaction (OER) is a major challenge in developing energy efficient and economically viable energy conversion devices such as electrolyzers, rechargeable metal-air batteries, and regenerative fuel cells. The current benchmark electrocatalyst for OER is based on iridium oxide (IrO x ) due to its superior performance and excellent stability. However, large scale applications using IrO x are impractical due to its low abundance and high cost. In this work, we report a highly active hafnium-modified iridium oxide (IrHf x O y ) electrocatalyst for OER. The IrHf x O y electrocatalyst demonstrated ten times higher activity in alkaline conditions (pH = 11) and four times higher activity in acid conditions (pH = 1) than a IrO x electrocatalyst. The highest intrinsic mass activity of the IrHf x O y catalyst in acid conditions was calculated as 6950 A gIrO x -1 at an overpotential (η) of 0.3 V. Combined studies utilizing operando surface enhanced Raman spectroscopy (SERS) and DFT calculations revealed that the active sites for OER are the Ir-O species for both IrO x and IrHf x O y catalysts. The presence of Hf sites leads to more negative charge states on nearby O sites, and shortening the bond lengths of Ir-O, and lowering free energies for OER intermediates to accelerate the OER process.

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

Development of molecular dynamics potential for uranium silicide fuels

Use of uranium–silicide (U-Si) in place of uranium dioxide (UO2) is one of the promising concepts being proposed to increase the accident tolerance of nuclear fuels. This is due to a higher thermal conductivity than UO2 that results in lower centerline temperatures. U-Si also has a higher fissile density, which may enable some new cladding concepts that would otherwise require increased enrichment limits to compensate for their neutronic penalty. However, many critical material properties for U-Si have not been determined experimentally. For example, silicide compounds (U3Si2 and U3Si) are known to become amorphous under irradiation. There was clear independent experimental evidence to support a crystalline to amorphous transformation in those compounds. However, it is still not well understood how the amorphous transformation will affect on fuel behavior. It is anticipated that modeling and simulation may deliver guidance on the importance of various properties and help prioritize experimental work. In order to develop knowledge-based models for use at the engineering scale with a minimum of empirical parameters and increase the predictive capabilities of the developed model, inputs from atomistic simulations are essential. First-principles based density functional theory (DFT) calculations will provide the most reliable information. However, it is probably not possible to obtain kinetic information such as amorphization under irradiation directly from DFT simulations due to size and time limitations. Thus, a more feasible way may be to employ molecular dynamics (MD) simulation. Unfortunately, so far no MD potential is available for U-Si to discover the underlying mechanisms. Here, we will present our recent progress in developing a U-Si potential from ab initio data. This work is supported by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program funded by the U.S. Department of Energy, Office of Nuclear Energy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

Insights into the Mechanism of Fe(II) Adsorption and Oxidation at Fe-Clay Mineral Surfaces from First-Principles Calculations

Interfacial reactivity of redox-active iron-bearing mineral surfaces plays a crucial role in many environmental processes including biogeochemical cycling of various elements and contaminants. Herein, we apply density-functional-theory (DFT) calculations to provide atomistic insights into the heterogeneous reaction between aqueous Fe(II) and the Fe-bearing clay mineral nontronite Fe2Si4O10(OH)2 by studying its adsorption mechanism and interfacial Fe(II)-Fe(III) electron transfer (ET) at edge and basal surfaces. We find that edge-bound Fe(II) adsorption complexes at different surface sites (ferrinol, silanol and mixed) may coexist on both (010) and (110) edge facets, with complexes at ferrinol FeO(H) sites being the most energetically favorable and coupled to proton transfer. Calculation of the ET activation energy suggests that interfacial ET into dioctahedral Fe(III) sheets is probable at the clay edges and occurs predominantly but not exclusively through the complexes adsorbed at ferrinol sites and might also involve mixed sites. No clear evidence is found for complexes on basal surface that are compatible with ET through the basal sheet despite this experimentally hypothesized ET interface. This study suggests a strong pH-dependence of Fe(II) surface complexation at basal versus edge facets and highlights the importance of the protonation state of bridging ligands and proton coupled electron transfer to facilitate ET into Fe-rich clay minerals.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2015

Adsorption of carbon monoxide on small aluminum oxide clusters: Role of the local atomic environment and charge state on the oxidation of the CO molecule

We present extensive density functional theory (DFT) calculations dedicated to analyze the adsorption behavior of CO molecules on small Al{sub x}O{sub y}{sup ±} clusters. Following the experimental results of Johnson et al. [J. Phys. Chem. A 112, 4732 (2008)], we consider structures having the bulk composition Al{sub 2}O{sub 3}, as well as smaller Al{sub 2}O{sub 2} and Al{sub 2}O units. Our electron affinity and total energy calculations are consistent with aluminum oxide clusters having two-dimensional rhombus-like structures. In addition, interconversion energy barriers between two- and one-dimensional atomic arrays are of the order of 1 eV, thus clearly defining the preferred isomers. Single CO adsorption on our charged Al{sub x}O{sub y}{sup ±} clusters exhibits, in general, spontaneous oxygen transfer events leading to the production of CO{sub 2} in line with the experimental data. However, CO can also bind to both Al and O atoms of the clusters forming aluminum oxide complexes with a CO{sub 2} subunit. The vibrational spectra of Al{sub x}O{sub y} + CO{sub 2} provides well defined finger prints that may allow the identification of specific isomers. The Al{sub x}O{sub y}{sup +} clusters are more reactive than the anionic species and the final Al{sub 2}O{sup +} + CO reaction can result in the production of atomic Al and carbon dioxide as observed from experiments. We underline the crucial role played by the local atomic environment, charge density distribution, and spin-multiplicity on the oxidation behavior of CO molecules. Finally, we analyze the importance of coadsorption and finite temperature effects by performing DFT Born-Oppenheimer molecular dynamics. Our calculations show that CO oxidation on Al{sub x}O{sub y}{sup +} clusters can be also promoted by the binding of additional CO species at 300 K, revealing the existence of fragmentation processes in line with the ones experimentally inferred.

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

Absolute NMR shielding scales and nuclear spin–rotation constants in {sup 175}LuX and {sup 197}AuX (X = {sup 19}F, {sup 35}Cl, {sup 79}Br and {sup 127}I)

We present nuclear spin–rotation constants, absolute nuclear magnetic resonance (NMR) shielding constants, and shielding spans of all the nuclei in {sup 175}LuX and {sup 197}AuX (X = {sup 19}F, {sup 35}Cl, {sup 79}Br, {sup 127}I), calculated using coupled-cluster singles-and-doubles with a perturbative triples (CCSD(T)) correction theory, four-component relativistic density functional theory (relativistic DFT), and non-relativistic DFT. The total nuclear spin–rotation constants determined by adding the relativistic corrections obtained from DFT calculations to the CCSD(T) values are in general in agreement with available experimental data, indicating that the computational approach followed in this study allows us to predict reliable results for the unknown spin–rotation constants in these molecules. The total NMR absolute shielding constants are determined for all the nuclei following the same approach as that applied for the nuclear spin–rotation constants. In most of the molecules, relativistic effects significantly change the computed shielding constants, demonstrating that straightforward application of the non-relativistic formula relating the electronic contribution to the nuclear spin–rotation constants and the paramagnetic contribution to the shielding constants does not yield correct results. We also analyze the origin of the unusually large absolute shielding constant and its relativistic correction of gold in AuF compared to the other gold monohalides.

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

Enhancing high-fidelity neural network potentials through low-fidelity sampling

The efficacy of neural network potentials (NNPs) critically depends on the quality of the configurational datasets used for training. Prior research using empirical potentials has shown that well-selected liquid–solid transitional configurations of a metallic system can be translated to other metallic systems. This study demonstrates that such validated configurations can be relabeled using density functional theory (DFT) calculations, thereby enhancing the development of high-fidelity NNPs. Training strategies and sampling approaches are efficiently assessed using empirical potentials and subsequently relabeled via DFT in a highly parallelized fashion for high-fidelity NNP training. Our results reveal that relying solely on energy and force for NNP training is inadequate to prevent overfitting, highlighting the necessity of incorporating stress terms into the loss functions. To optimize training involving force and stress terms, we propose employing transfer learning to fine-tune the weights, ensuring that the potential surface is smooth for these quantities composed of energy derivatives. This approach markedly improves the accuracy of elastic constants derived from simulations in both empirical potential-based NNPs and relabeled DFT-based NNPs. Overall, this study offers significant insights into leveraging empirical potentials to expedite the development of reliable and robust NNPs at the DFT level.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Assessment of amide I spectroscopic maps for a gas-phase peptide using IR-UV double-resonance spectroscopy and density functional theory calculations

The spectroscopy of amide I vibrations has become a powerful tool for exploring protein structure and dynamics. To help with spectral interpretation, it is often useful to perform molecular dynamics (MD) simulations. To connect spectroscopic experiments to simulations in an efficient manner, several researchers have proposed “maps,” which relate observables in classical MD simulations to quantum spectroscopic variables. It can be difficult to discern whether errors in the theoretical results (compared to experiment) arise from inaccuracies in the MD trajectories or in the maps themselves. In this work, we evaluate spectroscopic maps independently from MD simulations by comparing experimental and theoretical spectra for a single conformation of the α-helical model peptide Ac-Phe-(Ala){sub 5}-Lys-H{sup +} in the gas phase. Conformation-specific experimental spectra are obtained for the unlabeled peptide and for several singly and doubly {sup 13}C-labeled variants using infrared-ultraviolet double-resonance spectroscopy, and these spectra are found to be well-modeled by density functional theory (DFT) calculations at the B3LYP/6-31G** level. We then compare DFT results for the deuterated and {sup 13}C{sup 18}O-labeled peptide with those from spectroscopic maps developed and used previously by the Skinner group. We find that the maps are typically accurate to within a few cm{sup −1} for both frequencies and couplings, having larger errors only for the frequencies of terminal amides.

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

Linear-scaling calculation of Hartree-Fock exchange energy with non-orthogonal generalised Wannier functions

We present a method for the calculation of four-centre two-electron repulsion integrals in terms of localised non-orthogonal generalised Wannier functions (NGWFs). Our method has been implemented in the ONETEP program and is used to compute the Hartree-Fock exchange energy component of Hartree-Fock and Density Functional Theory (DFT) calculations with hybrid exchange-correlation functionals. As the NGWFs are optimised in situ in terms of a systematically improvable basis set which is equivalent to plane waves, it is possible to achieve large basis set accuracy in routine calculations. The spatial localisation of the NGWFs allows us to exploit the exponential decay of the density matrix in systems with a band gap in order to compute the exchange energy with a computational effort that increases linearly with the number of atoms. We describe the implementation of this approach in the ONETEP program for linear-scaling first principles quantum mechanical calculations. We present extensive numerical validation of all the steps in our method. Furthermore, we find excellent agreement in energies and structures for a wide variety of molecules when comparing with other codes. We use our method to perform calculations with the B3LYP exchange-correlation functional for models of myoglobin systems bound with O{sub 2} and CO ligands and confirm that the same qualitative behaviour is obtained as when the same myoglobin models are studied with the DFT+U approach which is also available in ONETEP. Finally, we confirm the linear-scaling capability of our method by performing calculations on polyethylene and polyacetylene chains of increasing length.

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

Mechanisms and site requirements for NO and NH 3 oxidation on Cu/SSZ-13

Two series of Cu/SSZ-13 catalysts were synthesized via aqueous solution and solid-state ion exchange using SSZ-13 supports of varying Si/Al ratios. The isolated and multinuclear Cu content of these catalysts were determined by H 2 temperature programmed reduction (H 2 -TPR). Multinuclear Cu in these catalysts, including in situ Cu- dimers formed from ZCu II OH coupling and permanent CuO clusters, are active species for dry NO oxidation. NH 3 oxidation on these catalysts follows an internal SCR (i-SCR) mechanism, i.e., a portion of NH 3 is first oxidized to NO, then NO is selectively reduced by the remaining NH 3 to N 2 . NH 3 oxidation displays distinct kinetic behavior below ~300 °C and above ~400 °C. At low temperature the results indicate that NH 3 -solvated mobile Cu-ions are the active centers. CuO clusters, when present, also contribute to the low temperature activity by catalyzing NH 3 oxidation to NO. At high temperature, in situ Cu-dimers and CuO clusters catalyze NH 3 oxidation to NO, and isolated Cu-ions catalyze SCR to realize the cascade turnovers. For both NO and NH 3 oxidation, Cu-dimers balanced by framework charges of close proximity appear to be more active than Cu-dimers balanced by distant framework charges. However, the former Cu-dimers are less stable than the latter and tend to split into monomers in the presence of vicinal Brønsted acid sites. Via density functional theory (DFT) calculations, the i-SCR mechanism for low temperature NH 3 oxidation, i.e., the energetic favorability for the involvement of the NO intermediate, is justified. Furthermore, the DFT results also agree with experimental data that the formation of Cu-dimers from ZCu II OH dimerization is essential for NH 3 oxidation at high temperature.

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

Thorium Monosilicide, ThSi: An Experimental and Theoretical Study

The present theoretical and experimental combination study investigates the ThSi molecule in detail. Computationally, we utilized high-level multireference and coupled-cluster levels of theory conjoined with large correlation consistent basis sets to study a series of electronic and spin–orbit states of ThSi. Here, we report potential energy curves (PECs), electron configurations at equilibrium distances, spectroscopic constants, energetics, and spin–orbit coupling effects for 16 electronic states of ThSi. The studied 16 electronic states are arranged tightly within 0.9 eV, highlighting the complexity of the electronic spectrum of ThSi. The ground electronic state of ThSi is a single-reference 1 1 Σ + state that derives from the 1σ 2 2σ 2 1π 4 electronic configuration. The Ω = 0 + spin–orbit ground state of ThSi is composed of 1 1 Σ + (47%) and 13Π (44%) electronic states. Our measured bond energy (D0) of ThSi, obtained using resonant two-photon ionization (R2PI) spectroscopy is 3.146(4) eV, where the assigned error limit is given in parentheses in units of the last quoted digits. The computed D0 of ThSi (Ω = 0 + ) at the CBS-C-CCSD(T)-δT(Q)-δDK-δSO level (3.181 eV) is in good agreement with the experimental value. Our derived enthalpy of formation for ThSi, Δ f H 0K o (ThSi(g)), is 971.8(6.0) kJ/mol. Finally, we have performed density functional theory (DFT) calculations for ThSi(1 1 Σ + ) using 16 exchange correlation functionals that span multiple rungs of “Jacob’s ladder” of density functional approximation (DFA) to assess the DFT errors on D 0 , r e , and ω e of ThSi with respect to experimental and ab initio coupled-cluster values.

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

Probing the Active Sites of MoS 2 Based Hydrotreating Catalysts Using Modulation Excitation Spectroscopy

The reactive surface sites of MoS 2 hydrotreating catalysts (unpromoted as well as Co- and Ni-promoted) supported on MgAl 2 O 4 spinel were investigated with respect to the substitution of sulfur by oxygen using in situ XAS coupled with modulation excitation spectroscopy (MES). Specifically, MES experiments were carried out by periodically cycling between a H 2 O and H 2 S containing hydrogen gas mixture at 400 °C. Due to the low fraction of S–O exchange, conventional XANES and EXAFS data hardly showed any changes when these catalysts were exposed to increasing ratios of H 2 O to H 2 S in an H 2 atmosphere. XANES and EXAFS data extracted at the Mo K-edge by MES analysis showed that for approximately 1% of the Mo atoms, sulfur atoms are replaced by oxygen atoms when exposed to H 2 O, causing partial oxidation of these active sites. The reaction is reversible and Mo returns to its initial sulfide phase when H 2 O is removed and H 2 S is supplied in the feed. In the case of Co- and Ni-promoted catalysts, the magnitude of S–O exchange was found to be reduced, indicating the beneficial effect of promotion. MES at the Ni K-edge showed that Ni was oxidized during H 2 O exposure, which in turn delayed the Mo oxidation in the Ni-promoted catalyst. The structures of these catalysts under S–O exchange were modeled using density functional theory (DFT) calculations, showing that the edge atoms are affected strongly. For all three catalysts, OH substitution is more favorable, while O substitution could be possible at high H 2 O pressure for unpromoted MoS 2 . Mo K-edge XANES spectra calculated using these simulated structures support the results obtained from the MES experiments. In conclusion, the presented approach using MES in combination with XAS and supported by DFT can be extended in general to catalysts under operando conditions and is thus a useful tool for determination of the active site on an atomic-scale.

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

Enhancing Catalytic Activity of MoS2 Basal Plane S-Vacancy by Co Cluster Addition

The basal plane of molybdenum disulfide (MoS 2 ) was recently activated for hydrogen evolution reaction (HER) by creating sulfur (S) vacancies (MoS 2-x ). However, the HER activity of those S-vacancies depends on the concentration of S-vacancies, imposing a dilemma for either improving activity per site or increasing overall active site density. In this work, we use density functional theory (DFT) calculations and experiments to show that the HER activities of MoS2-x are greatly enhanced by adding cobalt (Co) clusters on the basal plane. Furthermore, our DFT results show that the highest HER activity is achieved when the Co clusters are anchored on the S-vacancies with the interface of Co-Mo as the preferred active site. Our experiments confirm that the addition of Co enhances the activity per unit active site and increases the electrochemical active surface area. These results demonstrate the basal plane activity of MoS 2-x can be enhanced by decorating S-vacancies with transition-metal clusters.

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

Stability and conductivity of cation- and anion-substituted LiBH 4 -based solid-state electrolytes

The high-temperature phase of LiBH 4 (HT-LiBH 4 ) exhibits a promisingly high lithium ion conductivity but is unstable at room temperature. We use density functional theory (DFT) calculations to investigate the stabilization effect of halogen and alkali cation/anion substitutions on HT-LiBH 4 as well the underlying mechanism for the high lithium ion conductivity. We find that increasing the substituent concentration enhances the stabilization of HT-LiBH 4 (i.e., the DFT energy difference between the low- and high-temperature forms of substituted LiBH 4 is reduced). Cation/anion substitution also leads to a higher Li defect (vacancy, interstitial, and Frenkel) formation energy, thereby reducing the Li defect (vacancy, interstitial, and Frenkel) concentrations. Using DFT migration barriers input into kinetic Monte Carlo simulations and the Materials INTerface (MINT) framework, here we calculate the room-temperature lithium ion conductivities for Li(BH 4 ) 1-x I x (x = 0.25 and 0.5) and Li 1-y K y BH 4 (y = 0.25). Our calculations suggest that the lower I concentration leads to a higher lithium ion conductivity of 5.7 x10 -3 S/cm compared with that of Li(BH 4 ) 0.5 I 0.5 (4.2 ×10 -5 S/cm) because of the formation of more Li-related defects. Based on our findings, we suggest that the stabilization of HT-LiBH 4 -based lithium ion conductors can be controlled by carefully tuning the cation/anion substituent concentrations to maximize the lithium ionic conductivities of the specific system.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Energetics of native defects, solute partitioning, and interfacial energy of Q precipitate in Al-Cu-Mg-Si alloys

The compound Al 3 Cu 2 Mg 9 Si 7 , is known as the Q-phase and forms as a thermodynamically stable precipitate during aging in the quaternary Al-Cu-Mg-Si system. We perform atomic-scale density functional theory (DFT) calculations of defect properties, solute partitioning, and interfacial stability of the Al 3 Cu 2 Mg 9 Si 7 (Q) precipitate. We find: (i) simple native point defect (i.e., vacancies and anti-sites) thermodynamics can partially explain the experimentally observed off-stoichiometry, such as the observed variation of compositions, Al 3+δ Cu 2 Mg 9-δ Si 7 (Mg-deficient and Al-rich) in experiment. (ii) Calculated solute-partitioning energies of common solutes allow us to define general rules for site-preference in the Q-phase in terms of electronic structure and atomic radius. To validate our DFT predictions, we perform atom-probe tomography (APT) experiments for six-different elements (Zn, Ni, Mn, Ti, V, and Zr). The results show that the partitioning behavior of solutes Ni, Zn, and Mn are consistent with DFT predictions, but the transition elements (Ti, V, and Zr), which are anomalously slow diffusers in Al, partition to the Q-phase in constrast to DFT partitioning energies. (iii) For the low energy interface ($11\bar{2}$) Q //(510) Al observed in needle shaped Q-precipitate, we survey various terminations and orientations and derive a low-energy interfacial structure. Here we find this low-energy interfacial model has Cu atoms nearest to the interface, which is in agreement with previous literature on Cu interfacial segregation at the Q'//α-Al interface. The computed interfacial energy (0.52 J/m 2 ) and the corresponding structure will be useful input to future multi-scale modeling of microstructural evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Solution-Phase Growth of Cu Nanowires with Aspect Ratios Greater Than 1000: Multiscale Theory

Penta-twinned metal nanowires are finding widespread application in existing and emerging technologies. However, little is known about their growth mechanisms. Here, we probe the origins of chloride- and alkylamine-mediated, solution-phase growth of penta-twinned Cu nanowires from first principles using multi-scale theory. Using quantum density functional theory (DFT) calculations, we characterize the binding and surface diffusion of Cu atoms on chlorine-covered Cu(100) and Cu(111) surfaces. We find stronger binding and slower diffusion of Cu atoms on chlorinated Cu(111) than on chlorinated Cu(100) – a reversal of the trend for bare Cu surfaces. We also probe inter-facet diffusion and find this proceeds faster from Cu(100) to Cu(111) than the reverse. Using the DFT rates for hopping between individual sites at Ångstrom scales, we calculate coarse-grained, inter-facet rates for nanowires of various lengths – up to hundreds of μm – and diameters in the 10-nm range. We predict nanowires with aspect ratios around 100 based on surface diffusion, alone. We also account for the influence of a self-assembled alkylamine layer that covers most of the {100} facets, but is absent or thin and disordered on the {111} facets and in an “end zone” near the {100}-{111} boundary. With an end zone, we predict a wide range of nanowire aspect ratios in the experimental ranges. Our work reveals the mechanisms by which a halide – chloride – promotes the growth of high-aspect-ratio nanowires.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Extending GPU-accelerated Gaussian integrals in the TeraChem software package to f type orbitals: Implementation and applications

Here, the increasing availability of graphics processing units (GPUs) for scientific computing has prompted interest in accelerating quantum chemical calculations through their use. However, the complexity of integral kernels for high angular momentum basis functions often limits the utility of GPU implementations with large basis sets or for metal containing systems. In this work, we report the implementation of f function support in the GPU-accelerated TeraChem software package through the development of efficient kernels for the evaluation of Hamiltonian integrals. The high efficiency of the resulting code is demonstrated through density functional theory (DFT) calculations on increasingly large organic molecules and transition metal complexes, as well as coupled cluster singles and doubles calculations on water clusters. Preliminary investigations into Ni(I) catalysis with DFT and the photochemistry of MnH(CH 3 ) with complete active space self-consistent field are also carried out. Overall, our GPU-accelerated software appears to be well-suited for fast simulation of large transition metal containing systems, as well as organic molecules.

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

Studies on optical properties of Si220 nanoclusters via time-dependent density functional theory calculations

The optical properties of bare and hydrogen passivated Si220 nanoclusters(NCs) in four typical motifs(i.e., bulk-like, onion-like, bucky-diamond and icosahedral motifs) were studied via time-dependent density functional theory(TD-DFT) calculations. The calculation results show that there is a significant blue shift in the optical absorption spectra when the Si NCs are passivated with hydrogen. A strong absorption peak in the visible light region appears for the hydrogenated bulk-like, onion-like and bucky-diamond Si NCs.

Yang, Wenhua↗