Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “function”

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 289 records · Page 16

Physically interpretable approximations of many-body spectral functions

The rational function approximation provides a natural and interpretable representation of response functions such as the many-body spectral functions. We apply the vector fitting (VFIT) algorithm to fit a variety of spectral functions calculated from the Holstein model of electron-phonon interactions. We show that the resulting rational functions are highly efficient in their fitting of sharp features in the spectral functions, and could provide a means to infer physically relevant information from a spectral data set. The position of the peaks in the approximated spectral function are determined by the location of poles in the complex plane. Additionally, we developed a variant of VFIT that incorporates regularization to improve the quality of fits. With this procedure, we demonstrate it is possible to achieve accurate spectral function fits that vary smoothly as a function of physical conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Landmark-Warped Emulators for Models with Misaligned Functional Response

Many computer models output functional data, and in some cases, these functional data have similar, but misaligned, shape characteristics. In this paper, we introduce a general approach for building emulators for computer models that output misaligned functional data when key values in the functional response (landmarks) can be easily identified. This approach has two main parts: modeling the aligned (using the landmarks) functional data, and modeling the functions that map the misaligned data to the aligned space (warping functions). As the warping functions are required to be monotonic, we give special attention to modeling monotonic functional response data. We discuss how our approach can be easily applied for a variety of typical emulators, such as Gaussian processes, Bayesian multivariate adaptive regression splines, and Bayesian additive regression trees, and how sensitivity analysis can be performed. We demonstrate our approach by building emulators for two applications: (1) a high-energy-density physics computer model used to simulate inertial confinement fusion ignition experiments, where model outputs are highly misaligned, and (2) a multiphysics continuum hydrocode used to simulate high-velocity impact experiments, where model outputs are only slightly misaligned. In case (1) traditional methods cannot be applied, while in (2) they can be applied, but the proposed method performs significantly better.

97 MATHEMATICS AND COMPUTING↗

Entropy Pair Functional Theory: Direct Entropy Evaluation Spanning Phase Transitions

We prove that, within the class of pair potential Hamiltonians, the excess entropy is a universal, temperature-independent functional of the density and pair correlation function. This result extends Henderson’s theorem, which states that the free energy is a temperature dependent functional of the density and pair correlation. The stationarity and concavity of the excess entropy functional are discussed and related to the Gibbs–Bugoliubov inequality and to the free energy. We apply the Kirkwood approximation, which is commonly used for fluids, to both fluids and solids. Approximate excess entropy functionals are developed and compared to results from thermodynamic integration. The pair functional approach gives the absolute entropy and free energy based on simulation output at a single temperature without thermodynamic integration. We argue that a functional of the type, which is strictly applicable to pair potentials, is also suitable for first principles calculation of free energies from Born–Oppenheimer molecular dynamics performed at a single temperature. In conclusion, this advancement has the potential to reduce the evaluation the free energy to a simple modification to any procedure that evaluates the energy and the pair correlation function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comparative Density Functional Theory Study of Magnetic Exchange Couplings in Dinuclear Transition-Metal Complexes

Multicenter transition-metal complexes (MCTMs) with magnetically interacting ions have been proposed as components for information-processing devices and storage units. For any practical application of MCTMs as magnetic units, it is crucial to characterize their magnetic behavior, and in particular, the isotropic magnetic exchange coupling, J, between its magnetic centers. Due to the large size of typical MCTMs, density functional theory is the only practical electronic structure method for evaluating the J coupling. Here, we assess the accuracy of different density functional approximations for predicting the magnetic couplings of eight dinuclear transition-metal complexes, including five dimanganese, two dicopper, and one divanadium with known reliable experimental J couplings spanning from ferromagnetic to strong antiferromagnetic. The density functionals considered include global hybrid functionals which mix semilocal density functional approximations and exact exchange with a fixed admixing parameter, six local hybrid functionals where the admixing parameters are extended to be spatially dependent, the SCAN and r 2 SCAN meta-generalized gradient approximations (GGAs), and two widely used GGAs. We found that global hybrids tested in this work have a tendency to over-correct the error in magnetic coupling parameters from the Perdew–Burke–Ernzerhof (PBE) GGA as seen for manganese complexes. The performance of local hybrid density functionals shows no improvement in terms of bias and is scattered without a clear trend, suggesting that more efforts are needed for the extension from global to local hybrid density functionals for this particular property. The SCAN and r 2 SCAN meta-GGAs are found to perform as well as benchmark global hybrids on most tested complexes. We further analyze the charge density redistribution of meta-GGAs as well as global and local hybrid density functionals with respect to that of PBE, in connection to the self-interaction error or delocalization error.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Perspectives on van der Waals Density Functionals: The Case of TiS 2

The van der Waals interaction is of foundational importance for a wide variety of physical systems. In particular, van der Waals forces lie at the heart of potential device technologies that may be realized from the functional organization of layered two-dimensional (2D) nanomaterials. Furthermore, for intermediate to large-scale applications modeling, van der Waals density functionals have become the de facto choice for first-principles calculations. In particular, the vdW-DF family of functionals have provided a systematic approach to this theoretically challenging problem. While much progress has been made, there remains room for improvement in the microscopic description of vdW forces from these density functionals. In this work, we compute benchmark results for the binding energy and the electronic density response to binding in TiS2 via accurate diffusion quantum Monte Carlo calculations. We compare these benchmark data to results obtained from local, semilocal, and van der Waals functionals. In particular, we gauge the quality of the original vdW-DF/vdW-DF2 functionals, as well as updated variants such as vdW-DF-C09, vdW-DF-optB88, vdW-DF-optB86b, and vdW-DF2-B86R. We find a close relationship between the accuracy of predicted interlayer separation distances and binding energies for TiS 2 , with the vdW-DF-optB88 functional performing very well in terms of both quantities. In general, the more recently developed functionals are systematic improvements over older ones. However, when considering the response of the electron density to binding, we find that local-density approximation (LDA) and PBEsol generally outperform the vdW-DF functionals in describing the interlayer charge accumulation with vdW-DF-C09 variants performing the best overall.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Minimally entangled state preparation of localized wave functions on quantum computers

Initializing a single site of a lattice scalar field theory into an arbitrary state with support throughout the quantum register requires $O(2^{nQ})$ entangling gates on a quantum computer with $n_Q$ qubits per site. It is conceivable that instead initializing to functions that are good approximations to states may have utility in reducing the number of required entangling gates. In the case of a single site of a noninteracting scalar field theory, initializing to a symmetric exponential wave function requires $n_Q – 1$ entangling gates, the minimal number necessary to create an entangled state of all $n_Q$ qubits. This is compared with the $2^{nQ – 1} + n_Q – 3 + δ_{nQ, 1}$ required for a symmetric Gaussian wave function. In this work, we explore the initialization of one-site $(n_Q = 4)$, two-site $(n_Q = 3)$, and three-site ( n Q = 3 ) noninteracting scalar field theories with symmetric exponential wave functions using IBM's quantum simulators and quantum devices (Poughkeepsie and Tokyo). With the digitizations attainable with $n_Q = 3,4$, these tensor-product wave functions are found to have large overlap with a Gaussian wave function and provide a suitable low-noise initialization for improvement and Somma inflation. Further, in performing these simulations, we have employed a workflow that interleaves calibrations to mitigate systematic errors in production. The calibrations allow tolerance cuts on gate performance including the fidelity of the symmetrizing Hadamard gate, both in vacuum $(|\textbf{0}\gg^{\bigotimes nQ})$ and in medium ($n_Q –1$ qubits initialized to an exponential function). The results obtained in this work are relevant to systems beyond scalar field theories, such as the deuteron radial wave function, two- and three-dimensional Cartesian-space wave functions, and nonrelativistic multinucleon systems built on a localized eigenbasis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Sampling functions for geophysics

A set of spherical sampling functions is defined such that they are related to spherical-harmonic functions in the same way that the sampling functions of information theory are related to sine and cosine functions. An orderly distribution of (N + 1) squared sampling points on a sphere is given, for which the (N + 1) squared spherical sampling functions span the same linear manifold as do the spherical-harmonic functions through degree N. The transformations between the spherical sampling functions and the spherical-harmonic functions are given by recurrence relations. The spherical sampling functions of two arguments are extended to three arguments and to nonspherical reference surfaces. Typical applications of this formalism to geophysical topics are sketched.

Giacaglia, G. E. O.↗

On One-Dimensional Stretching Functions for Finite-Difference Calculations

The class of one dimensional stretching function used in finite difference calculations is studied. For solutions containing a highly localized region of rapid variation, simple criteria for a stretching function are derived using a truncation error analysis. These criteria are used to investigate two types of stretching functions. One is an interior stretching function, for which the location and slope of an interior clustering region are specified. The simplest such function satisfying the criteria is found to be one based on the inverse hyperbolic sine. The other type of function is a two sided stretching function, for which the arbitrary slopes at the two ends of the one dimensional interval are specified. The simplest such general function is found to be one based on the inverse tangent. The general two sided function has many applications in the construction of finite difference grids.

Vinokur, M.↗

A survey of functional programming language principles

Research in the area of functional programming languages has intensified in the 8 years since John Backus' Turing Award Lecture on the topic was published. The purpose of this paper is to present a survey of the ideas of functional programming languages. The paper assumes the reader is comfortable with mathematics and has knowledge of the basic principles of traditional programming languages, but does not assume any prior knowledge of the ideas of functional languages. A simple functional language is defined and used to illustrate the basic ideas. Topics discussed include the reasons for developing functional languages, methods of expressing concurrency, the algebra of functional programming languages, program transformation techniques, and implementations of functional languages. Existing functional languages are also mentioned. The paper concludes with the author's opinions as to the future of functional languages. An annotated bibliography on the subject is also included.

Holloway, C. M.↗

The single scattering phase functions of Jupiter's clouds

The determination of the single scattering phase functions of Jupiter's clouds and a thin upper haze by Tomasko et al. was refined and extended to seven latitudes in blue and red light. The phase function is well-constrained by the Pioneer 10 and 11 photometric data sets. Multiple scattering models were computed to match the limb darkening at each latitude at up to 15 phase angles from 12 deg to 151 deg. Ground-based observations were used for absolute calibration and to extend the data to lower phase angles. The phase functions were parameterized using the double Henyey-Greenstein function. The three Henyey-Greenstein parameters and the single scattering albedo were determined using a non-linear least squares method for the haze and the clouds below. The phase functions derived for the northen zone and belt are remarkably similar to the phase functions of the corresponding regions in the south, with most of the differences in brightness of the northern and southern features resulting from minor differences in single scattering albedo. Analysis of the Equatorial Region is complicated by the presence of numerous small features, but the phase function required is generally similar to that seen in the more homogeneous regions. Details of the phase functions of the haze and clouds are presented, and the differences between the cloud phase functions at low and high latitudes in red and blue light are discussed.

Doose, L. R.↗

The faintest stars - The luminosity and mass functions at the bottom of the main sequence

We present IR K-band photometry of complete samples of VLM candidates constructed from IIIaF and IVN plates in 10 fields taken as part of the POSSII and UKSRC surveys. Using the I-K colors constructed for these stars we estimate a bolometric luminosity function which extends to M(Bol) = 13.75. We find significant evidence for a luminosity function decreasing toward these luminosities. We also find that our results are consistent with those of studies based on the Nearby Star sample, when those data are presented as a bolometric luminosity function. We convert our observed luminosity function into a mass function, which extends with reasonable statistics to 0.08 solar masses - the H-burning minimum mass. We find significant evidence for features in the mass function at these masses. Specifically, the mass function 'turns over' at 0.25 solar mass, goes through a local minimum at about 0.15 solar mass, and may increase again below 0.1 solar mass - none of these features are predicted by any of the current theories of star formation. Lastly, the mass density we observe just above the H-burning minimum mass makes it difficult to envisage brown dwarfs contributing significant quantities of missing mass without invoking either a mass function in this region significantly steeper than that seen for main-sequence stars, or an extremely low cutoff mass to the mass function.

Tinney, Christopher G.↗

Functional and Multifunctional Polymers: Materials for Smart Structures

The ultimate goal of the research in smart structures and smart materials is the development of a new generation of products/devices which will perform better than products/devices built from passive materials. There are a few examples of multilayer polymer systems which function as smart structures, e.g. a synthetic muscle which is a multilayer assembly of a poly(ethylene) layer, a gold layer, and a poly(pyrrole) layer immersed in a liquid electrolyte. Oxidation and reductions of the active pyrrole layer causes the assembly to reversibly deflect and mimic biological muscles. The drawback of such a setup is slow response times and the use of a liquid electrolyte. We have developed multifunctional polymers which will eliminate the use of a liquid electrolyte, and also because the functionalities of the polymers are within a few hundred angstroms, an improved response time to changes in the external field should be possible. Such multifunctional polymers may be classified as the futuristic 'smart materials.' These materials are composed of a number of different functionalities which work in a synergistic fashion to function as a device. The device performs on the application of an external field and such multifunctional polymers may be scientifically labeled as 'field responsive polymers.' Our group has undertaken a systematic approach to develop functional and multifunctional polymers capable of functioning as field responsive polymers. Our approach utilizes multicomponent polymer systems (block copolymers and graft copolymers), the strategy involves the preparation of block or graft copolymers where the functionalities are limited to different phases in a microphase separated system. Depending on the weight (or volume) fractions of each of the components, different microstructures are possible. And, because of the intimate contact between the functional components, an increase in the synergism between the functionalities may be observed. In this presentation, three examples of multifunctional polymers developed in our labs will be reported. The first class of multifunctional polymers are the microphase separated mixed (ionic and electronic) conducting or MIEC block copolymers. The second class being developed in our labs are the biocompatible conductive materials and the conductive fluids. The final class may be considered microwave active smart polymers.

Arnold, S.↗

Electronic and Mechanical Properties of Hydrogen Functionalized Carbon Nanotubes

We examined the electronic and mechanical properties of hydrogen functionalized carbon nanotubes. The functionalization pattern covers two extreme groups. One group has randomly selected functionalization sites including one to twenty percent of the carbon atoms. The other group has regularly patterned functional sites parallel to the tube axis. Metallic, small-gap semiconducting and large-gap semiconducting carbon nanotubes are studied. The results reveal that the electronic properties of the tubes are very sensitive to the degree of functionalization, with even one percent functionalization being enough to render metallic tubes semiconducting. On the other hand, the mechanical properties, like tensile modulus, are much less sensitive to functionalization. For carbon nanotubes functionalized with specific patterns, the electric properties depends strongly on the nature of the functionalization pattern.

Yang, Liu↗

Efficient Evaluation Functions for Multi-Rover Systems

Evolutionary computation can be a powerful tool in cresting a control policy for a single agent receiving local continuous input. This paper extends single-agent evolutionary computation to multi-agent systems, where a collection of agents strives to maximize a global fitness evaluation function that rates the performance of the entire system. This problem is solved in a distributed manner, where each agent evolves its own population of neural networks that are used as the control policies for the agent. Each agent evolves its population using its own agent-specific fitness evaluation function. We propose to create these agent-specific evaluation functions using the theory of collectives to avoid the coordination problem where each agent evolves a population that maximizes its own fitness function, yet the system has a whole achieves low values of the global fitness function. Instead we will ensure that each fitness evaluation function is both "aligned" with the global evaluation function and is "learnable," i.e., the agents can readily see how their behavior affects their evaluation function. We then show how these agent-specific evaluation functions outperform global evaluation methods by up to 600% in a domain where a set of rovers attempt to maximize the amount of information observed while navigating through a simulated environment.

Agogino, Adrian↗

Slit Function Measurement of An Imaging Spectrograph Using Fourier Transform Techniques

Knowledge of a spectrograph slit function is necessary to interpret the unresolved lines in an observed spectrum. A theoretical slit function can be calculated from the sizes of the entrance slit, the detector aperture when it functions as an exit slit, the dispersion characteristic of the disperser, and the point spread function of the spectrograph. A measured slit function is preferred to the theoretical one for the correct interpretation of the spectral data. In a scanning spectrometer with a single exit slit, the slit function is easily measured. In a fixed grating/or disperser spectrograph, illuminating the entrance slit with a near monochromatic light from a pre-monochrmator or a tunable laser and varying the wavelength of the incident light can measure the slit function. Even though the latter technique had been used successfully for the slit function measurements, it had been very laborious and it would be prohibitive to an imaging spectrograph or a multi-object spectrograph that has a large field of view. We explore an alternative technique that is manageable for the measurements. In the proposed technique, the imaging spectrograph is used as a detector of a Fourier transform spectrometer. This method can be applied not only to an IR spectrograph but also has a potential to a visible/UV spectrograph including a wedge filter spectrograph. This technique will require a blackbody source of known temperature and a bolometer to characterize the interferometer part of the Fourier Transform spectrometer. This pa?er will describe the alternative slit function measurement technique using a Fourier transform spectrometer.

Park, Hongwoo↗

Molecular motors and their functions in plants

Molecular motors that hydrolyze ATP and use the derived energy to generate force are involved in a variety of diverse cellular functions. Genetic, biochemical, and cellular localization data have implicated motors in a variety of functions such as vesicle and organelle transport, cytoskeleton dynamics, morphogenesis, polarized growth, cell movements, spindle formation, chromosome movement, nuclear fusion, and signal transduction. In non-plant systems three families of molecular motors (kinesins, dyneins, and myosins) have been well characterized. These motors use microtubules (in the case of kinesines and dyneins) or actin filaments (in the case of myosins) as tracks to transport cargo materials intracellularly. During the last decade tremendous progress has been made in understanding the structure and function of various motors in animals. These studies are yielding interesting insights into the functions of molecular motors and the origin of different families of motors. Furthermore, the paradigm that motors bind cargo and move along cytoskeletal tracks does not explain the functions of some of the motors. Relatively little is known about the molecular motors and their roles in plants. In recent years, by using biochemical, cell biological, molecular, and genetic approaches a few molecular motors have been isolated and characterized from plants. These studies indicate that some of the motors in plants have novel features and regulatory mechanisms. The role of molecular motors in plant cell division, cell expansion, cytoplasmic streaming, cell-to-cell communication, membrane trafficking, and morphogenesis is beginning to be understood. Analyses of the Arabidopsis genome sequence database (51% of genome) with conserved motor domains of kinesin and myosin families indicates the presence of a large number (about 40) of molecular motors and the functions of many of these motors remain to be discovered. It is likely that many more motors with novel regulatory mechanisms that perform plant-specific functions are yet to be discovered. Although the identification of motors in plants, especially in Arabidopsis, is progressing at a rapid pace because of the ongoing plant genome sequencing projects, only a few plant motors have been characterized in any detail. Elucidation of function and regulation of this multitude of motors in a given species is going to be a challenging and exciting area of research in plant cell biology. Structural features of some plant motors suggest calcium, through calmodulin, is likely to play a key role in regulating the function of both microtubule- and actin-based motors in plants.

Non-NASA Center↗

A Study of Penalty Function Methods for Constraint Handling with Genetic Algorithm

COMETBOARDS (Comparative Evaluation Testbed of Optimization and Analysis Routines for Design of Structures) is a design optimization test bed that can evaluate the performance of several different optimization algorithms. A few of these optimization algorithms are the sequence of unconstrained minimization techniques (SUMT), sequential linear programming (SLP) and the sequential quadratic programming techniques (SQP). A genetic algorithm (GA) is a search technique that is based on the principles of natural selection or "survival of the fittest". Instead of using gradient information, the GA uses the objective function directly in the search. The GA searches the solution space by maintaining a population of potential solutions. Then, using evolving operations such as recombination, mutation and selection, the GA creates successive generations of solutions that will evolve and take on the positive characteristics of their parents and thus gradually approach optimal or near-optimal solutions. By using the objective function directly in the search, genetic algorithms can be effectively applied in non-convex, highly nonlinear, complex problems. The genetic algorithm is not guaranteed to find the global optimum, but it is less likely to get trapped at a local optimum than traditional gradient-based search methods when the objective function is not smooth and generally well behaved. The purpose of this research is to assist in the integration of genetic algorithm (GA) into COMETBOARDS. COMETBOARDS cast the design of structures as a constrained nonlinear optimization problem. One method used to solve constrained optimization problem with a GA to convert the constrained optimization problem into an unconstrained optimization problem by developing a penalty function that penalizes infeasible solutions. There have been several suggested penalty function in the literature each with there own strengths and weaknesses. A statistical analysis of some suggested penalty functions is performed in this study. Also, a response surface approach to robust design is used to develop a new penalty function approach. This new penalty function approach is then compared with the other existing penalty functions.

Ortiz, Francisco↗

Quickbird Satellite in-orbit Modulation Transfer Function (MTF) Measurement Using Edge, Pulse and Impulse Methods for Summer 2003

The spatial characteristics of an imaging system cannot be expressed by a single number or simple statement. However, the Modulation Transfer Function (MTF) is one approach to measure the spatial quality of an imaging system. Basically, MTF is the normalized spatial frequency response of an imaging system. The frequency response of the system can be evaluated by applying an impulse input. The resulting impulse response is termed the Point Spread function (PSF). This function is a measure of the amount of blurring present in the imaging system and is itself a useful measure of spatial quality. An underlying assumption is that the imaging system is linear and shift-independent. The Fourier transform of the PSF is called the Optical Transfer Function (OTF) and the normalized magnitude of the OTF is the MTF. In addition to using an impulse input, a knife-edge in technique has also been used in this project. The sharp edge exercises an imaging system at all spatial frequencies. The profile of an edge response from an imaging system is called an Edge Spread Function (ESF). Differentiation of the ESF results in a one-dimensional version of the Point Spread Function (PSF). Finally, MTF can be calculated through use of Fourier transform of the PSF as stated previously. Every image includes noise in some degree which makes MTF of PSF estimation more difficult. To avoid the noise effects, many MTF estimation approaches use smooth numerical models. Historically, Gaussian models and Fermi functions were applied to reduce the random noise in the output profiles. The pulse-input method was used to measure the MTF of the Landsat Thematic Mapper (TM) using 8th order even functions over the San Mateo Bridge in San Francisco, California. Because the bridge width was smaller than the 30-meter ground sample distance (GSD) of the TM, the Nyquist frequency was located before the first zero-crossing point of the sinc function from the Fourier transformation of the bridge pulse. To avoid the zero-crossing points in the frequency domain from a pulse, the pulse width should be less than the width of two pixels (or 2 GSD's), but the short extent of the pulse results in a poor signal-to-noise ratio. Similarly, for a high-resolution satellite imaging system such as Quickbird, the input pulse width was critical because of the zero crossing points and noise present in the background area. It is important, therefore, that the width of the input pulse be appropriately sized. Finally, the MTF was calculated by taking ratio between Fourier transform of output and Fourier transform of input. Regardless of whether the edge, pulse and impulse target method is used, the orientation of the targets is critical in order to obtain uniformly spaced sub-pixel data points. When the orientation is incorrect, sample data points tend to be located in clusters that result in poor reconstruction of the edge or pulse profiles. Thus, a compromise orientation must be selected so that all spectral bands can be accommodated. This report continues by outlining the objectives in Section 2, procedures followed in Section 3, descriptions of the field campaigns in Section 4, results in Section 5, and a brief summary in Section 6.

Helder, Dennis↗