Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio data”

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,099 records · Page 61

A Review on Lignin Liquefaction: Advanced Characterization of Structure and Microkinetic Modeling

Lignin liquefaction microkinetics is a move toward a more first principles (i.e., ab initio) based understanding at the molecular level in reaction engineering. While the microkinetic modelling of reactions to obtain kinetic rate parameters of chemical reactions have been widely used in the field of gas phase combustion and heterogeneous catalysis, this approach has not been as thoroughly developed in the area of biomass thermochemical reactions (e.g., lignin pyrolysis, hydrothermal liquefaction). The difficulties in establishing the structure of complex heterogeneous materials, like lignin, is perhaps the main challenge in developing rational microkinetic descriptions of biomass thermochemical reactions. In this manuscript we review the current state of the art and the challenges to develop microkinetic models for lignin liquefaction technologies (e.g., pyrolysis, hydrothermal liquefaction, solvolysis). A general strategy for the development of microkinetic models for lignin liquefaction technologies is discussed. Here, the first hurdle is to obtain sufficiently rich experimental data of lignin underlying polymeric structure and methodologies to use this data to build realistic lignin structural representations. Some analytical techniques for lignin structural characterization and their associated data, as well as a correlation for calculating the degree of macromolecular lignin branching are discussed. The presence of small lignin oligomeric structures and the role of these structures in lignin pyrolysis is also addressed. The ways in which elementary deconstruction and repolymerization reactions occur within this structure to form a liquid intermediate, and how these deconstruction products continue to interact with each other until they are removed from the liquid intermediate is thoroughly discussed. Further, experimental work with model compounds and the effect of reaction parameters (e.g., temperature, pressure, vapor residence time) are reviewed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accurate Quartic Force Fields and Vibrational Frequencies for HCN and HNC

The quartic force fields of HCN and HNC are determined using atomic natural orbital one-particle basis sets of spdf/spd and spdfg/spdf quality in conjunction with the CCSD(T) electron correlation method (singles and doubles coupled-cluster theory plus a perturbational estimate of the effects of connected triple excitations). The HCN force field is in good agreement with a recent experimentally derived force field and also with the force field recently computed by Wong and Bacskay. On the basis of the good agreement obtained for HCN, it is argued that the ab initio quartic force field for HNC is superior to a prior force field derived from experiment. The harmonic frequencies of HNC are predicted to be 3822 +/- 10,472 +/- 5, and 2051 +/-10/cm for omega1, omega2, and omega3, respectively; the experimentally derived values are above these values and fall outside the estimated uncertainties. Using the quartic force field, spectroscopic constants are predicted for HNC based on a vibrational second-order perturbation theory analysis. It is also asserted that the gas-phase fundamental v(sub 3) for HNC is slightly lower than the matrix isolation value. The range of validity of the quartic force fields is investigated by comparison of variational vibrational energies computed with the quartic force fields to those obtained from our recently reported global HCN/HNC potential energy surface and also to experimental data.

Lee, Timothy J.↗

Solid state {sup 1}H spin-lattice relaxation and isolated-molecule and cluster electronic structure calculations in organic molecular solids: The relationship between structure and methyl group and t-butyl group rotation

We report ab initio density functional theory electronic structure calculations of rotational barriers for t-butyl groups and their constituent methyl groups both in the isolated molecules and in central molecules in clusters built from the X-ray structure in four t-butyl aromatic compounds. The X-ray structures have been reported previously. We also report and interpret the temperature dependence of the solid state {sup 1}H nuclear magnetic resonance spin-lattice relaxation rate at 8.50, 22.5, and 53.0 MHz in one of the four compounds. Such experiments for the other three have been reported previously. We compare the computed barriers for methyl group and t-butyl group rotation in a central target molecule in the cluster with the activation energies determined from fitting the {sup 1}H NMR spin-lattice relaxation data. We formulate a dynamical model for the superposition of t-butyl group rotation and the rotation of the t-butyl group's constituent methyl groups. The four compounds are 2,7-di-t-butylpyrene, 1,4-di-t-butylbenzene, 2,6-di-t-butylnaphthalene, and 3-t-butylchrysene. We comment on the unusual ground state orientation of the t-butyl groups in the crystal of the pyrene and we comment on the unusually high rotational barrier of these t-butyl groups.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio study of the RbSr electronic structure: Potential energy curves, transition dipole moments, and permanent electric dipole moments

Excited states and the ground state of the diatomic molecule RbSr were calculated by post Hartree-Fock molecular orbital theory up to 22 000 cm{sup −1}. We applied a multireference configuration interaction calculation based on multiconfigurational self-consistent field wave functions. Both methods made use of effective core potentials and core polarization potentials. Potential energy curves, transition dipole moments, and permanent electric dipole moments were determined for RbSr and could be compared with other recent calculations. We found a good agreement with experimental spectra, which have been obtained recently by helium nanodroplet isolation spectroscopy. For the lowest two asymptotes (Rb (5s {sup 2}S) + Sr (5s4d {sup 3}P°) and Rb (5p {sup 2}P°) + Sr (5s{sup 2} {sup 1}S)), which exhibit a significant spin-orbit coupling, we included relativistic effects by two approaches, one applying the Breit-Pauli Hamiltonian to the multireference configuration interaction wave functions, the other combining a spin-orbit Hamiltonian and multireference configuration interaction potential energy curves. Using the results for the relativistic potential energy curves that correspond to the Rb (5s {sup 2}S) + Sr (5s4d {sup 3}P°) asymptote, we have simulated dispersed fluorescence spectra as they were recently measured in our lab. The comparison with experimental data allows to benchmark both methods and demonstrate that spin-orbit coupling has to be included for the lowest states of RbSr.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

First-Principles Investigations of the Temperature Dependence of Electronic Structure and Optical Properties of Rutile TiO 2

To gain additional insight into high-temperature functional material properties for applications in optical gas sensing, the temperature effects on the band gap and optical properties of rutile TiO 2 are investigated using ab initio methods. By analyzing the contributions from electron–phonon interaction and lattice thermal expansion, we show that the electron–phonon interaction is the dominant factor for temperature band-gap renormalization. As the temperature increases, the band gap increases until 300 K and then narrows above 300 K. This behavior results from the acoustic phonons, which widen the band gap, dominate below 300 K, while the optical phonons, which narrow the band gap, dominate above 300 K. Our study suggests that the band gap is narrowed by about 138 meV at 1000 K. We also investigated the temperature effects on the dielectric constants, the refractive index, as well as the extinction coefficient. Both the rate of decrease of the refractive index at 650 and 800 nm as well as the experimentally derived band gap agree with experimentally measured data as the temperature increases. Furthermore, our results and computational methods are of interest for developing high-temperature functional materials with applications toward gas sensing.

36 MATERIALS SCIENCE↗

Path integral Monte Carlo approach to the structural properties and collective excitations of liquid $$^3{\text {He}}$$ without fixed nodes

Abstract Due to its nature as a strongly correlated quantum liquid, ultracold helium is characterized by the nontrivial interplay of different physical effects. Bosonic $$^4{\text {He}}$$ 4 He exhibits superfluidity and Bose-Einstein condensation. Its physical properties have been accurately determined on the basis of ab initio path integral Monte Carlo (PIMC) simulations. In contrast, the corresponding theoretical description of fermionic $$^3{\text {He}}$$ 3 He is severely hampered by the notorious fermion sign problem, and previous PIMC results have been derived by introducing the uncontrolled fixed-node approximation. In this work, we present extensive new PIMC simulations of normal liquid $$^3{\text {He}}$$ 3 He without any nodal constraints. This allows us to to unambiguously quantify the impact of Fermi statistics and to study the effects of temperature on different physical properties like the static structure factor $$S({\mathbf {q}})$$ S ( q ) , the momentum distribution $$n({\mathbf {q}})$$ n ( q ) , and the static density response function $$\chi ({\mathbf {q}})$$ χ ( q ) . In addition, the dynamic structure factor $$S({\mathbf {q}},\omega )$$ S ( q , ω ) is rigorously reconstructed from imaginary-time PIMC data. From simulations of $$^3{\text {He}}$$ 3 He , we derived the familiar phonon–maxon–roton dispersion function that is well-known for $$^4{\text {He}}$$ 4 He and has been reported previously for two-dimensional $$^3{\text {He}}$$ 3 He films (Nature 483:576–579 (2012)). The comparison of our new results for both $$S({\mathbf {q}})$$ S ( q ) and $$S({\mathbf {q}},\omega )$$ S ( q , ω ) with neutron scattering measurements reveals an excellent agreement between theory and experiment.

97 MATHEMATICS AND COMPUTING↗

Modeling Intercalation Chemistry with Multiredox Reactions by Sparse Lattice Models in Disordered Rocksalt Cathodes

Modern battery materials can contain many elements with substantial site disorder, and their configurational state has been shown to be critical for their performance. The intercalation voltage profile is a critical parameter to evaluate the performance of energy storage. The application of commonly used cluster expansion techniques to model the intercalation thermodynamics of such systems ab initio is challenged by the combinatorial increase in configurational degrees of freedom as the number of species grows. Such challenges necessitate the efficient generation of lattice models without overfitting and proper sampling of the configurational space under the requirement of charge balance in ionic systems. In this work, we introduce a combined approach that addresses these challenges by (1) constructing a robust cluster expansion Hamiltonian using the sparse regression technique, including -norm regularization and structural hierarchy; and (2) implementing semigrand-canonical Monte Carlo to sample charge-balanced ionic configurations using the table-exchange method and an ensemble average approach. These techniques are applied to a disordered rocksalt oxyfluoride (LMNOF) that is part of a family of promising earth-abundant cathode materials. The simulated voltage profile is found to be in good agreement with experimental data and particularly provides a clear demonstration of the and oxygen contributions to the redox potential as a function of content.

25 ENERGY STORAGE↗

First-principles study on the electronic, optical and thermodynamic properties of ABO 3 (A = La,Sr, B = Fe,Co) perovskites

The electronic, optical and thermodynamic properties of ABO 3 (A = La,Sr, B = Fe,Co) perovskites are investigated using first-principles calculations. The obtained results indicate that SrCoO 3 and SrFeO 3 are metals, while LaCoO 3 and LaFeO 3 are insulators and all of them exhibit strong hybridization of the Fe/Co-3d and O-2p orbitals. By correlating the energy band structures with the peaks of the imaginary part of the dielectric function, we obtained the origin of each electron excitation to provide information about the active bands for the corresponding optical transitions observed in the experiment. Moreover, the Debye temperatures θ D obtained from the phonon frequencies are comparable to the available data. In conclusion, the thermodynamic properties of the Helmholtz free energy F , entropy S , and constant-volume heat capacity C v are investigated based on the phonon spectra.

01 COAL, LIGNITE, AND PEAT↗

Light nuclei with semilocal momentum-space regularized chiral interactions up to third order

In this work, we present a systematic investigation of few-nucleon systems and light nuclei using the current Low Energy Nuclear Physics International Collaboration interactions comprising semilocal momentum-space regularized two- and three-nucleon forces up to third chiral order (N 2 LO). Following our earlier study utilizing the coordinate-space regularized interactions, the two low-energy constants entering the three-body force are determined from the triton binding energy and the differential cross-section minimum in elastic nucleon-deuteron scattering. Predictions are made for selected observables in elastic nucleon-deuteron scattering and in the deuteron breakup reactions, for properties of the A = 3 and A = 4 nuclei, and for spectra of p-shell nuclei up to A = 16. A comprehensive error analysis is performed including an estimation of correlated truncation uncertainties for nuclear spectra. The obtained predictions are generally found to agree with experimental data within errors. Similarly to the coordinate-space regularized chiral interactions at the same order, a systematic overbinding of heavier nuclei is observed, which sets in for A~10 and increases with A.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

First-principles Hubbard U approach for small molecule binding in metal-organic frameworks

We apply first-principles approaches with Hubbard U corrections for calculation of small molecule binding energetics to open-shell transition metal atoms in metal-organic frameworks (MOFs). Using density functional theory with van der Waals dispersion-corrected functionals, we determine Hubbard U values ab initio through an established linear response procedure for M-MOF-74, for a number of different metal centers (M = Ti, V, Cr, Mn, Fe, Co, Ni, and Cu). While our ab initio U values differ from those used in previous work, we show that they result in lattice parameters and electronic contributions to CO{sub 2}-MOF binding energies that lead to excellent agreement with experiments and previous results, yielding lattice parameters within 3%. In addition, U-dependent calculations for an example system, Co-MOF-74, suggest that the CO{sub 2} binding energy grows monotonically with the value of Hubbard U, with the binding energy shifting 4 kJ/mol (or 0.041 eV) over the range of U = 0-5.4 eV. These results provide insight into an approximate but computationally efficient means for calculation of small molecule binding energies to open-shell transition metal atoms in MOFs and suggest that the approach can be predictive with good accuracy, independent of the cations used and the availability of experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Revisiting Selenium Interactions with Pyrite: From Adsorption to Coprecipitation

Interactions of selenium (Se), a trace element bio-essential at low concentrations but highly toxic at high concentrations, with the most abundant sulfide mineral in the Earth's crust, namely pyrite, was investigated over a wide range of time scales. At the nanosecond scale, selenate Se(VI)O 4 2– adsorption onto the net pyrite surface is shown by ab-initio computations to proceed via the formation of a chemical bond between an oxyanion oxygen atom and a surface Fe atom, weakening the other Se-O bonds and reducing Se atom oxidation state. At the hour-to-day scale, adsorption and coprecipitation of selenate Se(VI)O 4 2– and selenite, Se(IV)O 3 2– , were investigated through wet chemical batch experiments at various pH values at different sulfide concentrations. Selenium removal from solution is slower and weaker for selenate than for selenite. After 24h, only 10% of selenate, against 60% of selenite (up to 100% in the presence of sulfide), is removed by the pyrite surface. Independently of its original oxidation state, adsorbed Se is completely reduced to elemental trigonal selenium via adsorption, precipitation or coprecipitation, as shown by XANES spectroscopy. Here, our EXAFS results, compared to published data on Se-rich pyrite, show a Se to S substitution within the pyrite structure. The reductive coprecipitation mechanism of selenium with pyrite represents valuable new insights for improving our understanding of modern and ancient biogeochemical cycles involving Se. In addition, several industries can benefit from direct applications of our findings, such as water treatment, green technologies and sustainable mining.

Adsorption↗

Depolarizing collisions with hydrogen: Neutral and singly ionized alkaline earths

Depolarizing collisions are elastic or quasielastic collisions that equalize the populations and destroy the coherence between the magnetic sublevels of atomic levels. In astrophysical plasmas, the main depolarizing collider is neutral hydrogen. We consider depolarizing rates on the lowest levels of neutral and singly ionized alkali earths Mg I, Sr I, Ba I, Mg II, Ca II, and Ba II, due to collisions with H°. We compute ab initio potential curves of the atom-H° system and solve the quantum mechanical dynamics. From the scattering amplitudes, we calculate the depolarizing rates for Maxwellian distributions of colliders at temperatures T ≤ 10,000 K. A comparative analysis of our results and previous calculations in the literature is completed. We discuss the effect of these rates on the formation of scattering polarization patterns of resonant lines of alkali earths in the solar atmosphere, and their effect on Hanle effect diagnostics of solar magnetic fields.

79 ASTRONOMY AND ASTROPHYSICS↗

Accurate Quartic Force Fields and Vibrational Frequencies for HCN and HNC

The quartic force fields of HCN and HNC are determined using atomic natural orbital one-particle basis sets of spdf/spd and spdfg/spdf quality in conjunction with the CCSD(T) electron correlation method (singles and doubles coupled-cluster theory plus a perturbation estimate of the effects of connected triple excitations). The HCN force field is in good agreement with a recent experimentally derived force field and also with the force field recently computed by Wong and Bacskay. On the basis of the good agreement obtained for HCN, it is argued that the ab initio quartic force field for HNC is superior to a prior force field derived from experiment. The harmonic frequencies of HNC are predicted to be 3822 +/- 10, 472 +/- 5, and 2051 +/- 10 cm(exp -1) for omega(sub 1), omega(sub 2), and omega(sub 3), respectively; the experimentally derived values are above these values and fall outside the estimated uncertainties. Using the quartic force field, spectroscopic constants are predicted for HNC based on a vibrational second-order perturbation theory analysis. It is also asserted that the gas-phase fundamental nu(sub 3) for HNC is slightly lower than the matrix isolation value. The range of validity of the quartic force fields is investigated by comparison of variational vibrational energies computed with the quartic force fields to those obtained from our recently reported global HCN/HNC potential energy surface and also to experimental data.

Lee, Timothy J.↗

Thermophysical properties of liquid Ni around the melting temperature from molecular dynamics simulation

Thermophysical properties of liquid nickel (Ni) around the melting temperature are investigated by means of classical molecular dynamics (MD) simulation, using three different embedded atom method potentials to model the interactions between the Ni atoms. Melting temperature, enthalpy, static structure factor, self-diffusion coefficient, shear viscosity, and thermal diffusivity are compared to recent experimental results. Using ab initio MD simulation, we also determine the static structure factor and the mean-squared displacement at the experimental melting point. For most of the properties, excellent agreement is found between experiment and simulation, provided the comparison relative to the corresponding melting temperature. We discuss the validity of the Hansen-Verlet criterion for the static structure factor as well as the Stokes-Einstein relation between self-diffusion coefficient and shear viscosity. The thermal diffusivity is extracted from the autocorrelation function of a wavenumber-dependent temperature fluctuation variable.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio study of the structure and properties of amorphous silicon hydride from accelerated molecular dynamics simulations

This paper presents a large-scale ab initio simulation study of amorphous silicon hydride (a-Si 1-x H x ) with an emphasis on the structure and properties of the material across a range of hydrogen concentration by combining accelerated molecular dynamics (MD) simulations with first-principles density-functional calculations. The accelerated MD scheme relied on classical metadynamics, which enabled the development of 2500+ high-quality structural models of a-Si 1-x H x , with system sizes ranging from 150 to 6000 atoms and hydrogen concentrations vary from 6 to 20 at. %. The resulting amorphous networks were found to be completely free from any coordination defects and that they all exhibited a pristine band-gap in their electronic spectrum. The microstructural properties of hydrogen distributions were examined with an emphasis on the presence of isolated and clustered environments of hydrogen atoms. The results were compared with experimental data obtained from X-ray diffraction, infrared spectroscopy and nuclear magnetic resonance studies.

36 MATERIALS SCIENCE↗

Mobility of singly-charged lanthanide cations in rare gases: Theoretical assessment of the state specificity

High quality, ab initio calculations are reported for the potential energy curves governing the interactions of four singly-charged lanthanide ions (Yb{sup +}, Eu{sup +}, Lu{sup +}, and Gd{sup +}) with the rare gases (RG = He–Xe). Scalar-relativistic coupled cluster calculations are used for the first three S-state ions, but for Gd{sup +}({sup 10}D°) it is necessary to take the interaction anisotropy into account with the help of the multi-reference technique. The potential energy curves are used to determine the ion mobility and other transport properties describing the motion of the ions through the dilute RG, both as functions of the temperature, T, in the low-field limit, and at fixed T as functions of the ratio of the electrostatic field strength to the gas number density, E/N. The calculated mobilities are in good agreement with the very limited experimental data that have become available recently. The calculations show a pronounced dependence of the transport properties on the electronic configuration of the ion, as well as a significant effect of the spin-orbit coupling on the transport properties of the Gd{sup +} ion, and predict that state-specific mobilities could be detectable in Gd{sup +}–RG experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermoelectric properties of IV–VI-based heterostructures and superlattices

Doping in a manner that introduces anisotropy in order to reduce thermal conductivity is a significant focus in thermoelectric research today. By solving the semiclassical Boltzmann transport equations in the constant scattering time (τ) approximation, in conjunction with ab initio electronic structure calculations, within Density Functional Theory, we compare the Seebeck coefficient (S) and figure of merit (ZT) of bulk PbTe to PbTe/SnTe/PbTe heterostructures and PbTe doping superlattices (SLs) with periodically doped planes. Bismuth and Thallium were used as the n- and p-type impurities, respectively. The effects of carrier concentration are considered via chemical potential variation in a rigid band approximation. The impurity bands near the Fermi level in the electronic structure of PbTe SLs are of Tl s- and Bi p-character, and this feature is independent of the doping concentration or the distance between impurity planes. We observe the impurity bands to have a metallic nature in the directions perpendicular to the doping planes, yet no improvement on the values of ZT is found when compared to bulk PbTe. For the PbTe/SnTe/PbTe heterostructures, the calculated S presents good agreement with recent experimental data, and an anisotropic behavior is observed for low carrier concentrations (n<10{sup 18} cm{sup −3}). A large value of ZT{sub ||} (parallel to the growth direction) of 3.0 is predicted for n=4.7×10{sup 18} cm{sup −3} and T=700 K, whereas ZT{sub p} (perpendicular to the growth direction) is found to peak at 1.5 for n=1.7×10{sup 17} cm{sup −3}. Both electrical conductivity enhancement and thermal conductivity reduction are analyzed. - Graphical abstract: Figure of merit for PbTe/SnTe/PbTe heterostructure along the [0 0 1] direction, P.D. Borges, J.E. Petersen, L. Scolfaro, H.W. Leite Alves, T.H. Myers, Improved thermoelectric properties of IV–VI-based heterostructures and superlattices. - Highlights: • Thermoelectric properties of IV–VI-based heterostructures and superlattices. • High figure of merit is predicted for the PbTe/SnTe/PbTe heterostructure. • Nanotechnology has an important role for the development of thermoelectric devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗