Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional 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 1,009 records · Page 56

Strain-Induced Shifts in Defective Graphite Phonon Modes Predicted by Density Functional Theory

Carbon fiber composites have gained attention as a structural material because of their high strength-to-weight ratio, and understanding the effect of defects on reactivity and mechanical properties is important for the longevity and safety of the composite. Although it is known that strain causes the underlying graphitic vibrational modes to redshift, it is not clear how strain may alter reactivity and defect-induced vibrational changes. To investigate the strain-induced phonon changes of defective carbon fiber composites, density functional theory calculations of graphite are used, including intercalated hydrogen and fluorine defects. By comparing changes in the bond lengths, formation energies, and phonon density of states for uniaxially and biaxially strained graphite, strain was found to generally make defect formation more favorable and the specific behavior changes are dependent on the strain direction and defect identity. Specifically, intercalated fluorine phonons are more sensitive to strain than hydrogen intercalation phonons, and strain applied along the zigzag direction alters the calculated properties more than strain along the armchair direction. Finally, these results highlight the importance of understanding the microstructural effect of deviations from the ideal material because small changes in strain or defect type can significantly alter the behavior of the carbon fiber composite core.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic Structural Optimization of the Air Electrodes for Reversible Protonic Ceramic Electrochemical Cells via IIIA Cation Doping

Reversible protonic ceramic electrochemical cells (R-PCECs) have emerged as a novel technology for clean and efficient energy generation and storage. Improving the oxygen and proton conduction characteristics and stability of air electrodes at intermediate temperatures is crucial for achieving a high performance. Herein, we optimize the electronic structure of a state-of-the-art PrBa 0.8 Ca 0.2 Co 2 O 5+δ (PBCC) air electrode via doping IIIA cations (Al 3+ , Ga 3+ , and In 3+ ). Also, it is shown that PrBa 0.8 Ca 0.2 Co 1.9 Ga 0.10 O 5+δ (PBCCGa 0.10 ) exhibits improved oxygen reaction activity and hydration capability. Density functional theory calculations confirm that Ga doping provides the most favorable electronic structure. An R-PCEC with PBCCGa 0.10 achieves a peak power density of 2.21 W cm -2 and a current density of −4.64 A cm -2 at 1.3 V at 650 °C. Additionally, the PBCCGa 0.10 electrode performs good operational stability in FC mode (for about 100 h), EC mode (for about 100 h), and reversible cyclic testing (over 200 h) at 600 °C.

30 DIRECT ENERGY CONVERSION↗

Excited-state electronic coherence in vinyl bromide ions

We have performed time-dependent density functional theory calculations on the molecule vinyl bromide under ionization by a laser and have found long-lived coherent oscillations in the electron density created by a superposition of orbitals. The superposition of the molecular orbitals is quite pronounced for some intensities and is compared to varied alignments and strengths of the exciting laser. The coherence between the orbitals is also shown to be strongly dependent upon atomic motion.

74 ATOMIC AND MOLECULAR PHYSICS↗

Artificial layered perovskite oxides A(B{sub 0.5}B′{sub 0.5})O{sub 3} as potential solar energy conversion materials

Perovskite oxides with a d{sup 0} electronic configuration are promising photocatalysts and exhibit high electron mobilities. However, their band gaps are too large for efficient solar energy conversion. On the other hand, transition metal cations with partially filled d{sup n} electronic configurations give rise to visible light absorption. In this study, by using hybrid density functional theory calculations, it is demonstrated that the virtues of the two categories of materials can be combined in perovskite oxide A(B{sub 0.5}B′{sub 0.5})O{sub 3} with a layered B-site ordering along the [001] direction. The electronic structures of the four selected perovskite oxide compounds, La(Ti{sub 0.5}Ni{sub 0.5})O{sub 3}, La(Ti{sub 0.5}Zn{sub 0.5})O{sub 3}, Sr(Nb{sub 0.5}Cr{sub 0.5})O{sub 3}, and Sr(Nb{sub 0.5}Fe{sub 0.5})O{sub 3} are calculated and discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Giant Bulk Photovoltaic Effect in Vinylene-Linked Hybrid Heterocyclic Polymer

The bulk photovoltaic effect refers to the generation of a steady photocurrent from a homogeneous noncentrosymmetric material. It offers an alternative to the traditional p–n junction-based photovoltaic mechanism to directly convert sunlight to electricity. In this work, we investigate the bulk photovoltaic effect in low-dimensional polar organic materials with first-principles density functional theory calculations and shift current theory. With a strategy designed to break the inversion symmetry along the polymer chain, we demonstrate that conjugated vinylene-linked hybrid heterocyclic polymers can produce a strong bulk photovoltaic response to light, outperforming benchmark inorganic materials. The high current density results from the delocalized wave functions, composed mainly of carbon p orbitals. The great structural and electronic flexibility of polymers offers a robust paradigm to enhance the shift current response through chemical and physical modifications. The development of polymer blends with polymers of different band gaps potentially enables the utilization of the whole visible light spectrum for energy conversion.

14 SOLAR ENERGY↗

Self-assembled peptide nanotubes as electronic materials: An evaluation from first-principles calculations

In this letter, we report on the evaluation of diphenylalanine (FF), dityrosine (YY), and phenylalanine-tryptophan (FW) self-assembled peptide nanotube structures for electronics and photonics applications. Realistic bulk peptide nanotube material models were used in density functional theory calculations to mimic the well-ordered tubular nanostructures. Importantly, validated functionals were applied, specifically by using a London dispersion correction to model intertube interactions and a range-separated hybrid functional for accurate bandgap calculations. Bandgaps were found consistent with available experimental data for FF, and also corroborate the higher conductance reported for FW in comparison to FF peptide nanotubes. Interestingly, the predicted bandgap for the YY tubular nanostructure was found to be slightly higher than that of FW, suggesting higher conductance as well. In addition, the band structure calculations along the high symmetry line of nanotube axis revealed a direct bandgap for FF. The results enhance our understanding of the electronic properties of these material systems and will pave the way into their application in devices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Microscopic origin of the π states in epitaxial silicene

We investigated the electronic properties of epitaxial silicene on ZrB{sub 2}(0001) thin film grown on Si(111) by means of low-temperature scanning tunneling spectroscopy and density functional theory calculations. The position of silicon atoms and thus, the localization of the valence and conduction states were deducted from the comparison of the spectra and the computed local density of states. We point out the strong contribution of p{sub z} orbitals of specific atoms to those states which indicates the π character of the conduction and valence bands. A clear correlation between hybridization of the orbitals of the Si atoms and the buckling was evidenced.

36 MATERIALS SCIENCE↗

Role of opacity at the 9 keV back lighter energy used in measuring the equation of state of boron at pressures up to a Gbar

Many experiments have been conducted at the National Ignition Facility to measure the equation of state Hugoniot of plastic, boron, and diamond at extreme pressures up to a Gbar. The “Gbar” design employs a strong spherically converging shock launched through a solid ball of material using a hohlraum radiation drive. The shock front conditions are characterized using X-ray radiography, typically at energies near 9 keV. Here we examine how the opacity of boron at 9 keV changes at high pressures and temperatures. Understanding this is vital to unfolding the density in the shock front as pressures exceed 100 Mbar. We compare opacity calculations from a number of methods including the legacy XSN opacity tables, super transition array (STA), average atom (AVAT), and detailed configuration accounting (DCA) methods. We examine how the changing opacity is correlated with the K-shell occupation calculated using the opacity methods as well as other electronic structure methods such as MECCA, Purgatorio, and an all-electron Kohn-Sham density functional theory calculation that uses optimized all-electron norm-conserving Vanderbilt pseudopotentials (ONCV). We also examine the relative contribution of free-free opacity to the overall opacity, which is dominated by bound-free absorption.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comprehensive Computational Study of Partial Lead Substitution in Methylammonium Lead Bromide

Impurities in semiconductors, for example, lead-based hybrid perovskites, have a major influence on their performance as photovoltaic (PV) light absorbers. While impurities could create harmful trap states that lead to nonradiative recombination of charge carriers and adversely affect PV efficiency, they could also potentially increase absorption via midgap energy levels that act as stepping stones for subgap photons or introduce charge carriers via doping. To unearth trends in impurity energy states, we use first principles density functional theory calculations to extensively study partial substitution of Pb in methylammonium lead bromide (MAPbBr 3 ), a representative lead-halide perovskite. Investigation of the density of states and energy levels related to the transition of the substitutional defect from one charge state to another reveals that several elements create midgap energy states in MAPbBr 3 . We machine learned trends and design rules from the computational data and discovered that a few easily computed properties of the bromide compounds of any element can be used to predict the energetics and energy levels of the substitutional defect related to that element. The calculated Fermi level-dependent defect formation energies lead to the observation that substitution by transition metals, Zr, Hf, Nb, and Sc, and group V element Sb can shift the equilibrium Fermi level and change the perovskite conductivity, as determined by the dominant intrinsic point defects. Finally, metal-substituted MAPbBr 3 compounds of Bi, Sc, Ni, and Zr were experimentally investigated, and while there was an improvement in the thin-film morphology and an enhancement in charge carrier lifetimes, no clear evidence of subgap absorption features owing to the substituent being incorporated in the MAPbBr 3 lattice could be seen.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Competing itinerant and local spin interactions in kagome metal FeGe

The combination of a geometrically frustrated lattice, and similar energy scales between degrees of freedom endows two-dimensional Kagome metals with a rich array of quantum phases and renders them ideal for studying strong electron correlations and band topology. The Kagome metal, FeGe is a noted example of this, exhibiting A-type collinear antiferromagnetic (AFM) order at T N ≈ 400 K, then establishes a charge density wave (CDW) phase coupled with AFM ordered moment below T CDW ≈ 110 K, and finally forms a c-axis double cone AFM structure around T Canting ≈ 60 K. Here we use neutron scattering to demonstrate the presence of gapless incommensurate spin excitations associated with the double cone AFM structure of FeGe at temperatures well above T Canting and T CDW that merge into gapped commensurate spin waves from the A-type AFM order. Commensurate spin waves follow the Bose factor and fit the Heisenberg Hamiltonian, while the incommensurate spin excitations, emerging below T N where AFM order is commensurate, start to deviate from the Bose factor around T CDW , and peaks at T Canting . This is consistent with a critical scattering of a second order magnetic phase transition with decreasing temperature. By comparing these results with density functional theory calculations, we conclude that the incommensurate magnetic structure arises from the nested Fermi surfaces of itinerant electrons and the formation of a spin density wave order.

36 MATERIALS SCIENCE↗

Production of neutron-rich heavy nuclei in deep-inelastic 208 Pb + 208 Pb collisions within the stochastic mean-field theory

In deep-inelastic collisions of heavy nuclei, reaction products with a wide range of mass and charge are produced. Such collisions been considered as a possible way to produce superheavy nuclei, as an alternative to fusion reactions. To provide reliable theoretical predictions, it is desired to develop microscopic approaches that correctly and accurately describe nucleon transfer processes in dissipative collisions of heavy nuclei. The purpose of the present work is (1) to investigate the mechanism of nucleon transfers in dissipative collisions of two heavy nuclei, and (2) to explore possible pathways to produce neutron-rich heavy nuclei, through detailed theoretical analyses of fluctuations and correlations in nucleon transfers in 208 Pb + 208 Pb reactions. Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are performed for the collisions of 208 Pb + 208 Pb at 𝐸 c.m. = 832, 936, and 1040 MeV, using the Skyrme SLy4d energy density functional. To calculate fluctuations and correlations in nucleon transfers, we employ the stochastic mean-field (SMF) theory, and the results are compared with another theoretical framework currently available, the time-dependent random phase approximation (TDRPA). Primary and secondary production cross sections are calculated with the SMF theory combined with a statistical model, GEMINI ++ . Using information of nucleon flow across a neck of colliding nuclei in TDHF calculations, we solve quantal diffusion equations for fluctuations and correlations in nucleon transfers based on the SMF theory. From the SMF calculations, we obtain the time evolution of diffusion coefficients as well as fluctuations and correlations in nucleon transfers for a range of initial orbital angular momenta. We compare the results of the SMF calculations with those of TDRPA, showing that TDRPA tends to predict substantially larger fluctuations and correlations in strongly damped collisions of heavy nuclei, which exhibit complex initial angular momentum dependence, while the SMF results provide almost constant (stable) values. Using the obtained fluctuations and correlations, we calculate primary and secondary production cross sections for the 208 Pb + 208 Pb collisions. From the results, we find that both lighter and heavier reaction products as compared to 208 Pb are produced for a wide region in the 𝑁−𝑍 plane as primary products, thanks to the quantal diffusion mechanism in the dissipative collisions. However, we show that cross sections for production of heavy nuclei with 𝑍 ≳ 90 or 𝑁 ≳ 135 are washed out due to secondary particle evaporation and/or fission processes. On the other hand, we find that there remain sizable cross sections for production of neutron-rich nuclei along 𝑁 = 126 with 𝑍< 82, even after secondary disintegration processes. We demonstrate that the secondary production cross sections depend weakly on incident energies, but lower (higher) energy is slightly preferred for production of nuclei with smaller (larger) atomic numbers as compared to 𝑍 = 82. Here, based on the microscopic SMF calculations, it has been shown that deep-inelastic collisions of heavy nuclei, such as 208 Pb + 208 Pb examined in this study, can be a promising means to produce neutron-rich heavy nuclei along 𝑁=126. Discrepancies between the SMF and TDRPA approaches are left unsolved for future investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Structure and electronic properties of mixed (a + c) dislocation cores in GaN

Classical atomistic models and atomic-resolution scanning transmission electron microscopy studies of GaN films reveal that mixed (a + c)-type dislocations have multiple different core structures, including a dissociated structure consisting of a planar fault on one of the (12{sup ¯}10) planes terminated by two different partial dislocations. Density functional theory calculations show that all cores introduce localized states into the band gap, which affects device performance.

36 MATERIALS SCIENCE↗

Breakdown products of gaseous polycyclic aromatic hydrocarbons investigated with infrared ion spectroscopy

We report on a common fragment ion formed during the electron-ionization-induced fragmentation of three different three-ring polycyclic aromatic hydrocarbons (PAHs), fluorene (C{sub 13}H{sub 10}), 9,10-dihydrophenanthrene (C{sub 14}H{sub 12}), and 9,10-dihydroanthracene (C{sub 14}H{sub 12}). The infrared spectra of the mass-isolated product ions with m/z = 165 were obtained in a Fourier transform ion cyclotron resonance mass spectrometer whose cell was placed inside the optical cavity of an infrared free-electron laser, thus providing the high photon fluence required for efficient infrared multiple-photon dissociation. The infrared spectra of the m/z = 165 species generated from the three different precursors were found to be similar, suggesting the formation of a single C{sub 13}H{sub 9}{sup +} isomer. Theoretical calculations using density functional theory (DFT) revealed the fragment's identity as the closed-shell fluorenyl cation. Decomposition pathways from each parent precursor to the fluorenyl ion are proposed on the basis of DFT calculations. The identification of a single fragmentation product from three different PAHs supports the notion of the existence of common decomposition pathways of PAHs in general and can aid in understanding the fragmentation chemistry of astronomical PAH species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DFT study of small fullerene dimer complexes C{sub 20}-N{sub m}@C{sub n} (m = 1-6 and n = 24, 28, 32, 36 and 40)

First principle calculations based on density functional theory were performed to calculate the structural and electronic properties of C{sub 20}-N{sub m}@C{sub n} dimer complexes. The calculated binding energies of the complexes formed are comparable to C{sub 60} dimer which ensures their stability. The bond lengths of these dimer complexes were found to be nearly same as pure complexes C{sub 20}-C{sub n}. Further, nitrogen (N) atoms were encapsulated inside the secondary cage (C{sub n}) of dimer complexes and the number of N atoms depends on diameter of the cage. The HOMO-LUMO gaps of new proposed complexes indicate the increase in gap as compared to pure complexes. Mulliken charge analysis of these complexes has been studied which shows the significant charge transfer from the N atoms to the secondary cage of these complexes. The study propose the formation of the new dimer complexes which are stable and are able to encapsulate atoms which are otherwise reactive in free space.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Fluorescence of PRODAN in Water: a Computational QM/MM MD Study

Fluorescent properties of PRODAN (6-propionyl-2-dimethylaminonaphthalene) in water were studied by means of excited state molecular dynamics simulations employing a combined quantum mechanical and molecular mechanical approach with the time-dependent density functional theory (TD-DFT QM/MM MD). The state of the art coupled cluster method was used to benchmark density functional theory calculations. The influence of the water environment on PRODAN emission was investigated by employing several computational schemes with varying description of the solvent. The issue of the molecular geometry of the excited state PRODAN molecule in water was addressed. The experimental emission spectrum was reproduced for the planar excited state conformer of PRODAN in the extended environment of 300 explicit water molecules. The planar conformer was shown to be predominantly responsible for fluorescence. The twisted isomer is strongly stabilized in water, but rapidly evolve towards a conical intersection, and hence the twisted conformer is fluorescently inactive.

calculations↗

Norm-conserving pseudopotentials and basis sets to explore Lanthanide chemistry in complex environments

A complete set of pseudopotentials and accompanying basis sets for all lanthanide elements are presented based on the relativistic, norm-conserving, separable, dual-space Gaussian-type pseudopotential protocol of Goedecker, Teter and Hutter (GTH) within the generalized gradient approximation (GGA) and the exchange-correlation functional of Perdew, Burke and Ernzerhof (PBE). The corresponding basis sets have been molecularly optimized (MOLOPT) using a contracted form with a single set of Gaussian exponents for s, p and d states. The f states are uncontracted explicitly with Gaussian exponents. Moreover, the Hubbard U values for each lanthanide element, to be used in DFT+U calculations, are also tabulated, allowing for the proper treatment of the strong on-site Coulomb interactions of localized 4f electrons. The accuracy and reliability of our GTH pseudopotentials and companion basis sets optimized for lanthanides is illustrated by a series of test calculations on lanthanide-centered molecules, and solid-state systems, with the most common oxidation states. We anticipate that these pseudopotentials and basis sets will enable larger-scale density functional theory calculations and ab initio molecular dynamics simulations of lanthanide molecules in either gas or condensed phases, as well as of solid state lanthanide-containing materials, allowing to further explore the chemical and physical properties of lanthanide systems.

Lu, Jun-Bo↗

Electronic structure of reduced CeO 2 (111) surfaces interacting with hydrogen as revealed through electron energy loss spectroscopy in comparison with theoretical investigations

Based on both, ab-initio cluster calculations as well as periodic density functional based loss function calculations, we have assigned the origins of the valence electron excitation regime in electron energy loss spectra of well-ordered ceria films in (111) orientation at various states of reduction as well as after exposing the reduced films to hydrogen from the gas phase (Li et al, Angew. Chem. Int. Ed., 58 (2019) 14686–14693). The explicit calculation of intensity distributions using the dipole approximation allow us to draw conclusions about the nature of oxygen vacancies, which occur upon reduction and how those interact with hydrogen to form hydride species. Herein, we also reports a brief discussion of vibrational excitations of hydrogen loaded ceria and corroborates previous interpretations, based on inelastic neutron scattering.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Toward Addressing the Challenge to Predict the Heat Capacities of RDX and HMX Energetic Materials

Availability of heat capacity as function of pressure and temperature is an essential prerequisite for development of a computational multiscale strategy capable to address the evolution of microstructure and energy release in advanced high energy density materials. In the case of 1,3,5-trinitro-1,3,5-triazinane (RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane (HMX) systems as two of the most studied energetic materials, there are substantial gaps in experimental data, with available heat capacities values distributed only in a region close to standard ambient conditions. In this study we demonstrate how these major experimental limitations can be addressed in the case of the RDX and HMX systems based on the combined use of classical and quantum mechanical calculations. We show that by considering ideal gas properties evaluated using quantum mechanical methods, and residual properties obtained from molecular simulations using fully flexible atomistic force field models, excellent agreement can be obtained for the predicted heat capacities to the most recent experimental values. An important advantage of the current computational methodology is that it allows evaluation of both constant-volume and constant-pressure heat capacities for a broad interval of temperatures and pressures, which encompasses solid and liquid phases conditions. In the case of the solid α and γ phases of RDX and the β phase of HMX, the predicted results follow closely both the available experimental data at standard ambient conditions and the results obtained using density functional theory calculations at high pressures, a regime where experimental data are not available. A perspective to expand the current methodology is also discussed.

36 MATERIALS SCIENCE↗