Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional approximation”

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,351 records · Page 75

Variance of phase fluctuations of waves propagating through a random medium

As an electromagnetic wave propagates through a random scattering medium, such as a forest, its energy is attenuated and random phase fluctuations are induced. The magnitude of the random phase fluctuations induced is important in estimating how well a Synthetic Aperture Radar (SAR) can image objects within the scattering medium. The two-layer random medium model, consisting of a scattering layer between free space and ground, is used to calculate the variance of the phase fluctuations induced between a transmitter located above the random medium and a receiver located below the random medium. The scattering properties of the random medium are characterized by a correlation function of the random permittivity fluctuations. The effective permittivity of the random medium is first calculated using the strong fluctuation theory, which accounts for large permittivity fluctuations of the scatterers. The distorted Born approximation is used to calculate the first-order scattered field. A perturbation series for the phase of the received field in the Rytov approximation is then introduced and the variance of the phase fluctuations is also calculated assuming that the transmitter and receiver are in the paraxial limit of the random medium, which allows an analytic solution to be obtained. Results are compared using the paraxial approximation, scalar Green's function formulation, and dyadic Green's function formulation. The effects studied are the dependence of the variance of the phase fluctuations on receiver location in lossy and lossless regions, medium thickness, correlation length and fractional volume of scatterers, depolarization of the incident wave, ground layer permittivity, angle of incidence, and polarization.

Chu, Nelson C.↗

Microstructure of the interplanetary magnetic field near 4 and 5 AU

Interplanetary magnetic field fluctuations near 4 and 5 AU are investigated on the basis of Pioneer 10 and 11 vector magnetic field measurements in order to identify extant wave modes and propagation directions as functions of frequency in the range 0.00004 to 0.009 Hz. Analysis of the approximately 500 spectral matrices as functions of frequency obtained from each spacecraft indicates the distributions of maximum variance directions to be peaked along the normal to the plane of the minimum variance direction and the mean magnetic field, as would be expected for a planar, small-amplitude Alfven wave with a minimum variance direction close to the wave vector. It is thus inferred that single small-amplitude Alfven waves are occasionally present in the temporal and frequency bands analyzed, however waves with wave vectors not well approximated by the minimum variance direction, Alfven wave ensembles, finite-amplitude waves and nonplanar fluctuations may also be present. Observations are also noted to suggest the predominance of the Alfven mode over the fast and slow modes, however do not support theoretical models of outward propagating plane Alfven waves of solar origin.

Parker, G. D.↗

Impact of solvation on the GW quasiparticle spectra of molecules

First-principles calculations for electrochemistry require accurate treatment of both electronic structure and solvation. The perturbative GW approximation starting from density functional theory (DFT) calculations accurately models materials systems with varying dimensionality. Continuum solvation models enable efficient treatment of solvation effects in DFT calculations, but their applications with beyond-DFT electronic structure methods such as GW have been limited. Here, we introduce the frequency-dependent liquid polarizability from a nonlocal continuum solvation model in the screened Coulomb interaction of full-frequency GW calculations with a solvated DFT starting point. In this study, we show that the liquid screening contributions substantially reduce the HOMO–LUMO gap of molecules by 3–5 eV, while solvent effects on the DFT starting point negligibly impact the GW gap. The resulting framework facilitates the simultaneous electronic and solvation accuracy needed for first-principles electrochemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First-principles studies of the atomic, elastic and vibrational properties of BaNbO{sub 3}

We have employed the ab initio pseudopotential method, within a generalized gradient approximation of the density functional theory, a linear response approach, to determine atomic, elastic and lattice dynamical properties of the BaNbO{sub 3}. The obtained structural and elastic parameter are in good agreement with previous experimental and theoretical results. The analysis of B/G ratio indicates that BaNbO{sub 3} is a ductile material. The calculated lattice dynamical properties compared with phonon spectrum of the similar structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First-Principles Integrated Adsorption Modeling for Selective Capture of Uranium from Seawater by Polyamidoxime Sorbent Materials

Nuclear power is a relatively carbon-free energy source that has the capacity to be utilized today in an effort to stem the tides of global warming. The growing demand for nuclear energy, however, could put significant strain on our uranium ore resources, and the mining activities utilized to extract that ore can leave behind long-term environmental damage. A potential solution to enhance the supply of uranium fuel is to recover uranium from seawater using amidoximated adsorbent fibers. This technology has been studied for decades but is currently plagued by the material’s relatively poor selectivity of uranium over its main competitor vanadium. In this work, we investigate the binding schemes between uranium, vanadium, and the amidoxime functional groups on the adsorbent surface. Using quantum chemical methods, binding strengths are approximated for a set of complexation reactions between uranium and vanadium with amidoxime functionalities. Those approximations are then coupled with a comprehensive aqueous adsorption model developed in this work to simulate the adsorption of uranium and vanadium under laboratory conditions. Experimental adsorption studies with uranium and vanadium over a wide pH range are performed, and the data collected are compared against simulation results to validate the model. It was found that coupling ab initio calculations with process level adsorption modeling provides accurate predictions of the adsorption capacity and selectivity of the sorbent materials. Furthermore, this work demonstrates that this multiscale modeling paradigm could be utilized to aid in the selection of superior ligands or ligand compositions for the selective capture of metal ions. Furthermore, this first-principles integrated modeling approach opens the door to the in silico design of next-generation adsorbents with potentially superior efficiency and selectivity for uranium over vanadium in seawater.

36 MATERIALS SCIENCE↗

Half-metallic ferromagnetism in Cr-doped semiconducting Ge-chalcogenide: Density functional approach

A supercell approach has been used to calculate the electronic and magnetic properties of Cr-doped Ge chalcogenide, Ge{sub 1−x}Cr{sub x}Te (x = 0.25 and 0.125). The calculations have been performed using full potential Linear Augmented Plane Wave (FPLAPW) method within generalized gradient approximation (GGA) as exchange-correlation (XC) potential. The calculated results show that the doping of Cr induces the 100% spin polarization at Fermi level (EF) and showed the robust half metallic ferromagnetism in this compound. Thus, the compound at both dopant concentrations behave as dilute magnetic semiconductor (DMS) showing metallic property in majority and semiconducting for minority spin channels which is best suited for spintronic applications. The total magnetic moments of this compound are mainly due to Cr-d states present at E{sup F} with negligible contribution from electronic states of other atoms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Exploring the dependence of gas cooling and heating functions on the incident radiation field with machine learning

ABSTRACT Gas cooling and heating functions play a crucial role in galaxy formation. But, it is computationally expensive to exactly compute these functions in the presence of an incident radiation field. These computations can be greatly sped up by using interpolation tables of pre-computed values, at the expense of making significant and sometimes even unjustified approximations. Here, we explore the capacity of machine learning to approximate cooling and heating functions with a generalized radiation field. Specifically, we use the machine learning algorithm XGBoost to predict cooling and heating functions calculated with the photoionization code cloudy at fixed metallicity, using different combinations of photoionization rates as features. We perform a constrained quadratic fit in metallicity to enable a fair comparison with traditional interpolation methods at arbitrary metallicity. We consider the relative importance of various photoionization rates through both a principal component analysis (PCA) and calculation of SHapley Additive exPlanation (shap) values for our XGBoost models. We use feature importance information to select different subsets of rates to use in model training. Our XGBoost models outperform a traditional interpolation approach at each fixed metallicity, regardless of feature selection. At arbitrary metallicity, we are able to reduce the frequency of the largest cooling and heating function errors compared to an interpolation table. We find that the primary bottleneck to increasing accuracy lies in accurately capturing the metallicity dependence. This study demonstrates the potential of machine learning methods such as XGBoost to capture the non-linear behaviour of cooling and heating functions.

79 ASTRONOMY AND ASTROPHYSICS↗

Push it to the limit: comparing periodic and local approaches to density functional theory for intermolecular interactions

With the aim of systematically comparing two popular approaches to density functional theory–all-electron calculations with local basis sets, and periodic calculations employing plane wave basis sets and norm-conserving pseudopotentials–we have computed complete-basis binding energies across the S22 set of intermolecular interactions, a dataset consisting of noncovalent interactions of small- and medium-sized molecules containing first- and second-row atoms, using the Troullier-Martins norm-conserving pseudopotentials with SPW92, a local spin-density approximation; and PBE, a generalised gradient approximation. Here, we have found that it is challenging to reach the basis set limit with these periodic calculations; for the methods and systems examined, a minimum vacuum distance of 30 Å between a system and its nearest images is necessary–unless some form of dipole correction is employed–as is a kinetic energy cutoff of at least 80 Ry. The trends in convergence with respect to vacuum size and kinetic energy cutoff are largely independent of the level of density functional approximation employed. A sense of the impact of each hyperparameter on basis set error provides a foundation for ensuring quality calculations in future studies and allows us to quantify the basis set errors incurred in existing studies on similar systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

Predicting critical temperatures of iron(II) spin crossover materials: Density functional theory plus U approach

A first-principles study of critical temperatures (T{sub c}) of spin crossover (SCO) materials requires accurate description of the strongly correlated 3d electrons as well as much computational effort. This task is still a challenge for the widely used local density or generalized gradient approximations (LDA/GGA) and hybrid functionals. One remedy, termed density functional theory plus U (DFT+U) approach, introduces a Hubbard U term to deal with the localized electrons at marginal computational cost, while treats the delocalized electrons with LDA/GGA. Here, we employ the DFT+U approach to investigate the T{sub c} of a pair of iron(II) SCO molecular crystals (α and β phase), where identical constituent molecules are packed in different ways. We first calculate the adiabatic high spin-low spin energy splitting ΔE{sub HL} and molecular vibrational frequencies in both spin states, then obtain the temperature dependent enthalpy and entropy changes (ΔH and ΔS), and finally extract T{sub c} by exploiting the ΔH/T − T and ΔS − T relationships. The results are in agreement with experiment. Analysis of geometries and electronic structures shows that the local ligand field in the α phase is slightly weakened by the H-bondings involving the ligand atoms and the specific crystal packing style. We find that this effect is largely responsible for the difference in T{sub c} of the two phases. This study shows the applicability of the DFT+U approach for predicting T{sub c} of SCO materials, and provides a clear insight into the subtle influence of the crystal packing effects on SCO behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of water confined in lyotropic liquid crystals: Molecular dynamics simulations of the dynamic structure factor

The properties of water under confinement are of practical and fundamental interest. Here in this work we study the properties of water in the self-assembled lyotropic phases of gemini surfactants with a focus on testing the standard analysis of quasi-elastic neutron scattering (QENS) experiments. In QENS experiments the dynamic structure factor is measured and fit to models to extract the translational diffusion constant, D T , and rotational relaxation time, τ R . We test this procedure by using simulation results for the dynamic structure factor, extracting the dynamic parameters from the fit as is typically done in experiments, and comparing the values to those directly measured in the simulations. We find that the decoupling approximation, where the intermediate scattering function is assumed to be a product of translational and rotational contributions, is quite accurate. The jump-diffusion and isotropic rotation models, however, are not accurate when the degree of confinement is high. In particular, the exponential approximations for the intermediate scattering function fail for highly confined water and the values of D T and τ R can differ from the measured value by as much as a factor of two. Other models have more fit parameters, however, and with the range of energies and wave-vectors accessible to QENS, the typical analysis appears to be the best choice. In the most confined lamellar phase, the dynamics are sufficiently slow that QENS does not access a large enough time scale and neutron spin echo measurements would be a valuable technique in addition to QENS.

36 MATERIALS SCIENCE↗

Experimentally-Derived Bidirectional Reflectance Distribution Function Data in Support of the Orbital Debris Program Office

The NASA Orbital Debris Program Office (ODPO) has used various optical assets to acquire photometric data of Earth-orbiting objects to define the orbital debris environment. To better characterize and model optical data acquired from ground-based telescopes, the Optical Measurements Center (OMC) at NASA Johnson Space Center emulates illumination conditions seen in space by using equipment and techniques that parallel telescopic observations and source-target-sensor orientations. One of the OMC goals is to improve the size calculation used for optical data by developing an optical-based Size Estimation Model. The current size estimation requires applying a Lambertian phase function, a set albedo value, and range to the observed magnitude. The first step to improving the sampled brightness of laboratory targets is to remove aspect-angle dependencies. Then, the volume of possible object viewing angles is sampled at 21 combinations of azimuth and zenith angles for each solar phase angle. Finally, the acquired images are input into an image processing program that generates approximations for the object’s Bidirectional Reflectance Distribution Function (BRDF) and phase function. The BRDF is a radiometric concept that identifies an object’s material composition by matching a BRDF approximated with photometric data collected by ground-based telescopes with a BRDF generated experimentally from a known object in the laboratory. This paper discusses the validation of experimental BRDF and phase function approximations produced in the OMC and how the findings will be incorporated into ODPO models. A Lambertian sphere is imaged and the subsequent experimental functions are scrutinized to confirm that they correspond to an object that has an isotropic luminance. With the image processing algorithm validated, test objects with varying optical properties are then imaged to confirm that the produced photometric functions are both unique and repeatable. Once the validation is complete, the OMC will be used to evaluate a subset of fragments from a hypervelocity impact test of a mock-up satellite and assess the appropriate phase function and size estimates using BRDF measurements for a large volume of targets composed of various shapes, sizes, and materials.

Orbital Debris Program Office (ODPO)↗

Ab Initio Investigation of Structural Stability and Electronic and Magnetic Properties of the Half-Heusler Alloys: MTiSb (M = Fe, Co, and Ni)

The structural, electronic, and magnetic properties of the half-Heusler MTiSb (M = Fe, Co, Ni) compounds have been examined using the first-principles calculations with the full-potential linear augmented plane wave (FP-LAPW) method within the density functional theory. The electronic structures and magnetic properties of MTiSb (M = Fe, Co, Ni) are studied. It is found that the calculated lattice parameter is in good agreement with the theoretical values. Using the general gradient approximation (GGA), we observe that for the three compounds, the α phase is more stable than β and γ phases. CoTiSb is a semiconductor, NiTiSb is a metal, while the FeTiSb compound is a half-metallic, and the total magnetic moments are 0.00085, − 0.26074, and 0.99353 μ B, respectively, which is in agreement with the Slater–Pauling rule, Furthermore, the origin for the appearance of the half-metallic band gap in the FeTiSb compound was also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Sparse matrix techniques applied to modal analysis of multi-section duct liners

A simplified procedure is presented for analysis of ducts with discretely nonuniform properties. The analysis uses basis functions as the generalized coordinates. The duct eigenfunctions are approximated by finite series of these functions. The emphasis is on solution of the resulting large sparse set of linear equations. Characteristics of sparse matrix algorithms are outlined and some criteria for application are established. Analogies with structural methods are used to illustrate variations which can increase efficiency in generating values for design optimization routines. The effects of basis function selection, number of eigenfunctions and identification and ordering of equations on the sparsity and solution stability are included.

Arnold, W. R.↗

Symmetry-projected wave functions in quantum Monte Carlo calculations

Here, we consider symmetry-projected Hartree-Fock trial wave functions in constrained-path Monte Carlo (CPMC) calculations. Previous CPMC calculations have mostly employed Hartree-Fock (HF) trial wave functions, restricted or unrestricted. The symmetry-projected HF approach results in a hierarchy of wave functions with increasing quality: the more symmetries that are broken and restored in a self-consistent manner, the higher the quality of the trial wave function. This hierarchy is approximately maintained in CPMC calculations: the accuracy in the energy increases and the statistical variance decreases when further symmetries are broken and restored. Significant improvement is achieved in CPMC with the best symmetry-projected trial wave functions over those from simple HF. We analyze and quantify the behavior using the two-dimensional repulsive Hubbard model as an example. In the sign-problem-free region, where CPMC can be made exact but a constraint is deliberately imposed here, spin-projected wave functions remove the constraint bias. Away from half filling, spatial symmetry restoration in addition to that of the spin leads to highly accurate results from CPMC. Since the computational cost of symmetry-projected HF trial wave functions in CPMC can be made to scale algebraically with system size, this provides a potentially general approach for accurate calculations in many-fermion systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗