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 541 records · Page 30

An optimum PN code search technique for a given apriori signal location density

A method is presented for optimizing (minimizing) the acquisition time for a PN-type spread spectrum system when the a priori probability density function is not uniform. For an assumed a priori Gaussian density function, calculations show that the acquisition time, when the 0.5 probability acquisition time is used as a measure of acquisition time, is reduced by 40% for a cell detection probability of 0.25 and when three sweeps are used. For the same parameters and with four sweeps the acquisition time is reduced by 41%. When the acquisition time probability is set to 0.9 instead of 0.5, the reduction is 25% of the uniform sweep time.

Holmes, J. K.↗

Radial vibrations of a sodium ion inside icosahedral C60

The very high symmetry of icosahedral C60 suggests that, as a first approximation, an atom trapped inside C60 would be subject to a potential that is radially symmetric about the center. All-electron local-density-functional calculations of the total energy of a sodium ion as a function of radial displacement from the center along the fivefold axis of C60 serve to refine such a radial potential. In particular, the calculations suggest studying potentials that have minima displaced from the center. An analytic functional form for a radial potential having a positive cusp at the origin is proposed, and the s-wave radial solutions of the corresponding Schroedinger equation are examined.

Ballester, J. L.↗

In Situ Poling and Imidization of Amorphous Piezoelectric Polyimides

An amorphous piezoelectric polyimide containing polar functional groups has been developed using a combination of experimental and molecular modeling for potential use in high temperature applications. This amorphous polyimide, (Beta-CN)APB/ODPA, has exhibited good thermal stability and piezoelectric response at temperatures up to 150C. Density functional calculations predicted that a partially cured amic acid (open imide ring) possesses a dipole moment four times larger than the fully imidized closed ring. In situ poling and imidization of the partially cured (Beta-CN)APB/ODPA, was studied in an attempt to maximize the degree of dipolar orientation and the resultant piezoelectric response. A positive corona poling was used to minimize localized arcing during poling and to allow use of higher poling fields without dielectric breakdown. The dielectric relaxation strength, remanent polarization, and piezoelectric response were evaluated as a function of the poling profile. The partially cured, corona poled polymers exhibited higher dielectric relaxation strength (delta varepsilon), remanent polarization (Pr) and piezoelectric strain coefficient (d33) than the fully cured, conventionally poled ones.

Park, Cheol↗

Heats of Segregation of BCC Metals Using Ab Initio and Quantum Approximate Methods

Many multicomponent alloys exhibit surface segregation, in which the composition at or near a surface may be substantially different from that of the bulk. A number of phenomenological explanations for this tendency have been suggested, involving, among other things, differences among the components' surface energies, molar volumes, and heats of solution. From a theoretical standpoint, the complexity of the problem has precluded a simple, unified explanation, thus preventing the development of computational tools that would enable the identification of the driving mechanisms for segregation. In that context, we investigate the problem of surface segregation in a variety of bcc metal alloys by computing dilute-limit heats of segregation using both the quantum-approximate energy method of Bozzolo, Ferrante and Smith (BFS), and all-electron density functional theory. In addition, the composition dependence of the heats of segregation is investigated using a BFS-based Monte Carlo procedure, and, for selected cases of interest, density functional calculations. Results are discussed in the context of a simple picture that describes segregation behavior as the result of a competition between size mismatch and alloying effects

Good, Brian↗

Kinetic Monte Carlo Simulation of Cation Diffusion in Low-K Ceramics

Low thermal conductivity (low-K) ceramic materials are of interest to the aerospace community for use as the thermal barrier component of coating systems for turbine engine components. In particular, zirconia-based materials exhibit both low thermal conductivity and structural stability at high temperature, making them suitable for such applications. Because creep is one of the potential failure modes, and because diffusion is a mechanism by which creep takes place, we have performed computer simulations of cation diffusion in a variety of zirconia-based low-K materials. The kinetic Monte Carlo simulation method is an alternative to the more widely known molecular dynamics (MD) method. It is designed to study "infrequent-event" processes, such as diffusion, for which MD simulation can be highly inefficient. We describe the results of kinetic Monte Carlo computer simulations of cation diffusion in several zirconia-based materials, specifically, zirconia doped with Y, Gd, Nb and Yb. Diffusion paths are identified, and migration energy barriers are obtained from density functional calculations and from the literature. We present results on the temperature dependence of the diffusivity, and on the effects of the presence of oxygen vacancies in cation diffusion barrier complexes as well.

Good, Brian↗

IR Spectra and Properties of Solid Acetone, an Interstellar and Cometary Molecule

Mid-infrared spectra of amorphous and crystalline acetone are presented along with measurements of the refractive index and density for both forms of the compound. Infrared band strengths are reported for the first time for amorphous and crystalline acetone, along with IR optical constants. Vapor pressures and a sublimation enthalpy for crystalline acetone also are reported. Positions of (sup 13) C-labeled acetone are measured. Band strengths are compared to gas-phase values and to the results of a density-functional calculation. A 73 percent error in previous work is identified and corrected.

ethers↗

IR Spectra and Properties of Solid Acetone, an Interstellar and Cometary Molecule

Mid-infrared spectra of amorphous and crystalline acetone are presented along with measurements of the refractive index and density for both forms of the compound. Infrared band strengths are reported for the first time for amorphous and crystalline acetone, along with IR optical constants. Vapor pressures and a sublimation enthalpy for crystalline acetone also are reported. Positions of 13C-labeled acetone are measured. Band strengths are compared to gas-phase values and to the results of a density-functional calculation. A 73% error in previous work is identified and corrected.

IR spectroscopy, band strengths, astrochemistry, a↗

Neutron Skin Measurement of 208Pb and 48Ca Using Parity Violating Electron Scattering

Parity-violating electron scattering experiments (PVES) provides a clean probe of neutron densities that is model independent and free from strong interaction uncertainties in interpretation. The PREX-2 and CREX experiments were run in 2019 and 2020 at Jefferson laboratory measured the nucleon skin thickness, the difference between the r.m.s. neutron radius Rn and the r.m.s. proton radius Rp, of 208Pb and 48Ca via parity violating electroweak asymmetry in the elastic scattering of longitudinally polarized electrons. PREX-2 experiment was performed with 950 MeV electrons scattered at a 5¿ angle with Q2 = 0.00616 ± 0.00004 (GeV/c)2, while CREX used 2182 MeV electrons at the same angle with Q2 = 0.0297 ± 0.0002 (GeV/c)2. For PREX-2 the measured asymmetry was AP V = 550 ± 16 (stat.) ± (syst.) ppb, which corresponds to Rskin = 0.278 ± 0.078 (exp.) ± 0.012 (theo.) fm. The CREX asymmetry was AP V = 2668 ± 106 (stat.) ± 40 (syst.) ppb, which corresponds to Rskin = 0.121 ± 0.026 (exp.) ±0.024 (model) fm. One of the crucial systematic uncertainty that PREX-2 and CREX were sensitive to was the non-parity violating asymmetries that resulted from the helicity correlated false asymmetries in the polarized electron beam. There was a lot of work put towards understanding and suppressing the asymmetries arising from these effects. The parity violating asymmetry measurement required a very precise determination of the electron beam polarization. To accurately determine the beam polarization, a Compton polarimeter was used during both PREX-2 and CREX. A careful alignment of the laser to the Fabry-Perot cavity, data analysis and systematic control was employed to get a precise beam polarization result for the experiments. The PREX-2 measurement has broad implications for increasing our knowledge about neutron star structure and the equation of state of nuclear matter. The combined PREX and CREX results will have implications for future energy density functional calculations and the theory of nuclear structure.

Premathilake, Sachinthani↗

Local and average structure of Mn- and La-substituted BiFeO{sub 3}

The local and average structure of solid solutions of the multiferroic perovskite BiFeO{sub 3} is investigated by synchrotron X-ray diffraction (XRD) and electron density functional theory (DFT) calculations. The average experimental structure is determined by Rietveld refinement and the local structure by total scattering data analyzed in real space with the pair distribution function (PDF) method. With equal concentrations of La on the Bi site or Mn on the Fe site, La causes larger structural distortions than Mn. Structural models based on DFT relaxed geometry give an improved fit to experimental PDFs compared to models constrained by the space group symmetry. Berry phase calculations predict a higher ferroelectric polarization than the experimental literature values, reflecting that structural disorder is not captured in either average structure space group models or DFT calculations with artificial long range order imposed by periodic boundary conditions. Only by including point defects in a supercell, here Bi vacancies, can DFT calculations reproduce the literature results on the structure and ferroelectric polarization of Mn-substituted BiFeO{sub 3}. The combination of local and average structure sensitive experimental methods with DFT calculations is useful for illuminating the structure-property-composition relationships in complex functional oxides with local structural distortions. - Graphical abstract: The experimental and simulated partial pair distribution functions (PDF) for BiFeO{sub 3}, BiFe{sub 0.875}Mn{sub 0.125}O{sub 3}, BiFe{sub 0.75}Mn{sub 0.25}O{sub 3} and Bi{sub 0.9}La{sub 0.1}FeO{sub 3}.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coupling computational thermodynamics with density-function-theory based calculations to design L12 precipitates in Fe Ni based alloys

Achieving a high-volume fraction of thermodynamically stable L12-type precipitates that are resistant to coarsening is of great importance for the development of low-cost FeNi based austenitic steels. With the aid of computational thermodynamics, this work designed two model alloys: Fe-37.4Ni-6.1Al-2.9Ti (FNAT) and Fe-45.2Ni-5.9Al-8.5Si (FNAS). Both alloys were designed to contain a similar amount of L1 2 precipitate in Fe-Ni matrix without forming other precipitates. Density-Function-Theory (DFT) calculation was coupled with computational thermodynamics to predict the critical radius at which the precipitates change shape from spherical to cuboidal. The calculation results suggest that critical radius for the FNAT alloy is about two orders of magnitude larger than that for the FeNiAlSi alloy. Phase stability and morphology of the L1 2 precipitates in these two alloys were experimentally investigated through X-ray diffraction, atom probe tomography, and scanning and transmission electron microscopy. The L1 2 precipitates in the FeNiAlSi system were found to be cuboidal and rod shaped, with much larger size than the spherical ones in the FeNiAlTi system, agreeing with the calculation results. This work suggested that coupling computational thermodynamics with DFT calculations can be reliably used to design L1 2 precipitates in FeNi based alloys.

36 MATERIALS SCIENCE↗

Hybrid functional calculations of electronic structure and carrier densities in rare-earth monopnictides

The structural parameters and electronic structure of rare-earth pnictides are calculated using density functional theory (DFT) with the Heyd, Scuseria, and Ernzerhof (HSE06) screened hybrid functional. Specifically, we focus on RE-V compounds, with RE=La, Gd, Er, and Lu, and V=As, Sb, and Bi, and analyze the effects of spin-orbit coupling and treating the RE4$\textit{f}$ electrons as valence electrons in the projector augmented wave approach. The results of HSE06 calculations are compared with DFT within the generalized gradient approximation (DFT-GGA) and other previous calculations. We find that all these RE-V compounds are semimetals with electron pockets at the $\textit{X}$ point and hole pockets at Γ. Whereas in DFT-GGA the carrier density is significantly overestimated, the computed carrier densities using HSE06 is in good agreement with the available experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Surface Studies of β-1,3,5,7-Tetranitro-1,3,5,7-Tetrazoctane and Pentaerythritol Tetranitrate from Density Functional Tight-Binding Calculations and Implications on Crystal Shape

Here, we use density functional tight-binding (DFTB) theory to calculate the surface energies of two energetic crystals: monoclinic β-1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (β-HMX) and tetragonal pentaerythritol tetranitrate (PETN). The results are then employed to determine crystal shapes using the Bravais–Friedel–Donnay–Harker, attachment energy, and surface energy models. We find that energy-based models yield predictions in good agreement with experimental observations. Additionally, we propose a simple model that reframes surface energy as a measure of the lost intermolecular interactions during the formation of a surface from the bulk. The model accurately captures the results from the DFTB calculations and enables us to explain and predict surface energies as a function of the local molecular environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hybrid functional pseudopotentials

The consistency between the exchange-correlation functional used in pseudopotential construction and in the actual density functional theory calculation is essential for the accurate prediction of fundamental properties of materials. However, routine hybrid density functional calculations at present still rely on generalized gradient approximation pseudopotentials due to the lack of hybrid functional pseudopotentials. Here, we present a scheme for generating hybrid functional pseudopotentials, and we analyze the importance of pseudopotential density functional consistency for hybrid functionals. For the PBE0 hybrid functional, we benchmark our pseudopotentials for structural parameters and fundamental electronic gaps of the Gaussian-2 (G2) molecular dataset and some simple solids. Our results demonstrate that using our PBE0 pseudopotentials in PBE0 calculations improves agreement with respect to all-electron calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

First-principles calculation of the configurational energy density of states for a solid-state ion conductor with a variant of the Wang and Landau algorithm

In this work, a variant of the Wang and Landau algorithm for calculation of the configurational energy density of states is proposed. The algorithm was developed for the purpose of using first-principles simulations, such as density functional theory, to calculate the partition function of disordered sublattices in crystal materials. The expensive calculations of first-principles methods make a parallel algorithm necessary for a practical computation of the configurational energy density of states within a supercell approximation of a solid-state material. The algorithm developed in this work is tested with the two-dimensional (2d) Ising model to bench mark the algorithm and to help provide insight for implementation to a materials science application. Tests with the 2d Ising model revealed that the algorithm has good performance compared to the original Wang and Landau algorithm and the 1/$\textit{t}$ algorithm, in particular the short iteration performance. Further, a proof of convergence is presented within an adiabatic assumption, and the analysis is able to correctly predict the time dependence of the modification factor to the density of states. The algorithm was then applied to the lithium and lanthanum sublattice of the solid-state lithium ion conductor Li 0.5 La 0.5 TiO 3 . This was done to help understand the disordered nature of the lithium and lanthanum. The results find, overall, that the algorithm performs very well for the 2d Ising model and that the results for Li 0.5 La 0.5 TiO 3 are consistent with experiment while providing additional insight into the lithium and lanthanum ordering in the material. The primary result is that the lithium and lanthanum become more mixed between layers along the c axis for increasing temperature. In part, the simulation of the disordered Li 0.5 La 0.5 TiO 3 system serves as a benchmark for what size systems are currently and in the near future practical to calculate with density functional theory methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics

Accurately modeling dense plasmas over wide-ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/$\mathrm{cm^3}$ to 16 g/$\mathrm{cm^3}$ and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature method. Here, we find that the computed EOS and Hugoniot are in good agreement with path integral Monte Carlo results and the sesame database. Additionally, we calculate the ion-ion structure factor and viscosity for selected points. All results presented are at the level of full Kohn-Sham DFT-MD, free of empirical parameters, average-atom, and orbital-free approximations employed previously at such conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Use of density functional theory method to calculate structures of neutral carbon clusters C{sub n} (3 ≤ n ≤ 24) and study their variability of structural forms

In this work, we present modifications to the well-known basin hopping (BH) optimization algorithm [D. J. Wales and J. P. Doye, J. Phys. Chem. A 101, 5111 (1997)] by incorporating in it the unique and specific nature of interactions among valence electrons and ions in carbon atoms through calculating the cluster’s total energy by the density functional tight-binding (DFTB) theory, using it to find the lowest energy structures of carbon clusters and, from these optimized atomic and electronic structures, studying their varied forms of topological transitions, which include a linear chain, a monocyclic to a polycyclic ring, and a fullerene/cage-like geometry. In this modified BH (MBH) algorithm, we define a spatial volume within which the cluster’s lowest energy structure is to be searched, and introduce in addition a cut-and-splice genetic operator to increase the searching performance of the energy minimum than the original BH technique. The present MBH/DFTB algorithm is, therefore, characteristically distinguishable from the original BH technique commonly applied to nonmetallic and metallic clusters, technically more thorough and natural in describing the intricate couplings between valence electrons and ions in a carbon cluster, and thus theoretically sound in putting these two charged components on an equal footing. The proposed modified minimization algorithm should be more appropriate, accurate, and precise in the description of a carbon cluster. We evaluate the present algorithm, its energy-minimum searching in particular, by its optimization robustness. Specifically, we first check the MBH/DFTB technique for two representative carbon clusters of larger size, i.e., C{sub 60} and C{sub 72} against the popular cut-and-splice approach [D. M. Deaven and K. M. Ho, Phys. Rev. Lett. 75, 288 (1995)] that normally is combined with the genetic algorithm method for finding the cluster’s energy minimum, before employing it to investigate carbon clusters in the size range C{sub 3}-C{sub 24} studying their topological transitions. An effort was also made to compare our MBH/DFTB and its re-optimized results carried out by full density functional theory (DFT) calculations with some early DFT-based studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Biological applications of graphene oxide

Graphene as a 2D material has unique chemical and electronic properties. Because of its unique physical, chemical, and electronic properties, its interesting shape and size make it a promising nanomaterial in many biological applications. However, the lower water-solubility and the irreversible aggregation due to the strong π-π stacking hinder the wide application of graphene nanosheets in biomedical field. Thus, graphene oxide (GO), one derivative of graphene, has been used more frequently in the biological system owing to its relatively higher water solubility and biocompatibility. Recently, it has been demonstrated that nanomaterials with different functional groups on the surface can be used to bind the drug molecules with high affinity. GO has different functional groups such as H, OH and O on its surface; it can be a potential candidate as a drug carrier. The interactions of biomolecules and graphene like structures are long-ranged and very weak. Development of new techniques is very desirable for design of bioelectronics sensors and devices. In this work, we present first-principles spin polarized calculations within density functional theory to calculate effects of charging on DNA/RNA nucleobases on graphene oxide. It is shown that how modify structural and electronic properties of nucleobases on graphene oxide by applied charging.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗