Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hamiltonian”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Hamilton's Equations with Euler Parameters for Rigid Body Dynamics Modeling

A combination of Euler parameter kinematics and Hamiltonian mechanics provides a rigid body dynamics model well suited for use in strongly nonlinear problems involving arbitrarily large rotations. The model is unconstrained, free of singularities, includes a general potential energy function and a minimum set of momentum variables, and takes an explicit state space form convenient for numerical implementation. The general formulation may be specialized to address particular applications, as illustrated in several three dimensional example problems.

Shivarama, Ravishankar↗
NASA NTRS2006

Dynamic Theory of Relativistic Electrons Stochastic Heating by Whistler Mode Waves with Application to the Earth Magnetosphere

In the Hamiltonian approach an electron motion in a coherent packet of the whistler mode waves propagating along the direction of an ambient magnetic field is studied. The physical processes by which these particles are accelerated to high energy are established. Equations governing a particle motion by group symmetries of the problem were transformed in to a closed pair of nonlinear difference equations. The solutions of these equations have shown there exists the energetic threshold below that the electron motion is regular, and when the initial energy is above the threshold an electron moves stochastically. It is proved that the upper boundary of particle stochastic heating is conditioned by intrinsic property of the particle chaotic motion. Particle energy spectra and pitch angle electron scattering are described by the Fokker-Planck-Kolmogorov equations. It is shown that significant pitch angle diffusion occurs for the Earth radiation belt electrons with energies from a few keV up to a few MeV.

Khazanov, G. V.↗
NASA NTRS1997

Bromine Dioxide, OBrO: The Rotational Spectrum and Molecular Properties

The spectra are well descrie by Hamiltonian which included centrifugal distortion efects for fine and hyperfine terms. The molecular structure and the harmonic force field have been derived, and they as well as fine and hyperfine structure constants, are compared with data of related molecules and ESR data from OBrO isolated in cryogenic salt matrices.

Bromine Dioxide Rotational Spectrum↗
NASA NTRS2004

Spin lifetime tuning in zincblende heterostructures and applications to spin devices

We present analytical expressions for the D'yakonov-Perel' spin relaxation rates under the combined action of bulk and structural inversion asymmetry for zincblende heterostructures when terms up to linear and third order in k are included in the Hamiltonian. We see for heterostructures that, under the right conditions, the lowest-order-in-k component of the spin relaxation tensor can be made to vanish for all spin components at the same time. We study how the inclusion of terms of higher order in k affects these results.

spin filters↗
NASA NTRS2005

Measurement of absolute cross sections for excitation of the 2s(2) S-1 -> 2s2p P-1 degrees transition in O+4

Experimental cross sections are reported for the 1s(2)2s(2) S-1 -> 1s(2)2s2p P-1(o) transition in O+4 located at 19.689 eV. Use is made of the electron energy-loss method, using a merged electron-ion beam geometry. The center-of-mass interaction energies for the measurements in the S-1 -> P-1(o) transition are in the range 18 eV ( below the threshold) to 30 eV. Data are compared with other previous electron energy-loss measurements and with results of a 26 term R-matrix calculation that includes fine structure explicitly via the Breit-Pauli Hamiltonian. Clear resonance enhancement is observed in all experimental and theoretical results near the threshold for this S-1 -> P-1(o) transition.

atomic processes↗
NASA NTRS2006

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗
NASA NTRS2006

Electron Impact Excitation of Forbidden and Allowed Transitions in O(II)

The B-spline R-matrix method is used to investigate the electron impact excitation of forbidden and allowed transitions in singly ionized oxygen. The relativistic effects have been incorporated in the Breit-Pauli Hamiltonian. Flexible non-orthogonal sets of radial functions are used to obtain accurate target description and to represent the scattering functions. The 47 fine-structure levels of the 2s(sup 2)2p(sup 3), 2s2p(sup 4), 2s(sup 2)2p(sup 2)3s, 2s(sup 2)2p(sup 2)3p and 2s(sup 2)2p(sup 2)3d configurations have been included in the scattering calculation. A calculation with 62 levels in the close-coupling expansion using the Breit-Pauli R-matrix (BPRM) method with orthogonal radial functions has also been carried out to check electron correlation, relativistic and channel coupling effects. The present results are in good agreement with the previous 16-level BPRM calculation by Montenegro et a1 (2006 J. Phys. B: At. Mol. Opt. Phys. 39 1863-77) for the forbidden transitions, but differ from the 21-level BPRM calculation by McLaughlin and Bell (1998 J. Phys. B: At. Mol. Opt. Phys. 31 4317-29). Our cross sections for the first forbidden (sup 4)S(sup o)-(sup 2)D(sup o)and resonance (sup 4)S(sup o)-2s(sup 2)p(sup 4) (sup 4)P transitions are in reasonably good agreement with the electron energy-loss and merged-beams experiment.

Tayal, S. S.↗
NASA NTRS2004

Accurate Calculation of Oscillator Strengths for CI II Lines Using Non-orthogonal Wavefunctions

Non-orthogonal orbitals technique in the multiconfiguration Hartree-Fock approach is used to calculate oscillator strengths and transition probabilities for allowed and intercombination lines in Cl II. The relativistic corrections are included through the Breit-Pauli Hamiltonian. The Cl II wave functions show strong term dependence. The non-orthogonal orbitals are used to describe the term dependence of radial functions. Large sets of spectroscopic and correlation functions are chosen to describe adequately strong interactions in the 3s(sup 2)3p(sup 3)nl (sup 3)Po, (sup 1)Po and (sup 3)Do Rydberg series and to properly account for the important correlation and relaxation effects. The length and velocity forms of oscillator strength show good agreement for most transitions. The calculated radiative lifetime for the 3s3p(sup 5) (sup 3)Po state is in good agreement with experiment.

Tayal, S. S.↗
NASA NTRS2007

Strain Dependence of Photoluminescense of Individual Carbon Nanotubes

We have investigated strain dependence of photoluminescense (PL) spectra of single wall carbon nanotubes (SWNT). Nanotubes were sparsely dispersed in a thin PMMA film applied to acrylic bar, and strained in both compression and extension by bending this bar in either direction in a homebuilt four-point bending rig. The average surface strain was measured with high accuracy by a resistive strain gage applied on top of the film. The near infrared imaging and spectroscopy were performed on the inverted microscope equipped with high numerical aperture reflective objective lens and InGaAs CCD cameras. PL was excited with a diode laser at either 658, 730 or 785 nm, linearly polarized in the direction of the strain. We were able to measure (n,m) types and orientation of individual nanotubes with respect to strain direction and strain dependence of their PL maxima. It was found that PL peak shifts with respect to the values measured in SDS micelles are a sum of three components. First, a small environmental shift due to difference in the dielectric constant of the surrounding media, that is constant and independent of the nanotube type. Second, shift due to isotropic compression of the film during drying. Third, shifts produced by the uniaxial loading of the film in the experiment. Second and third shifts follow expression based on the first-order expansion of the TB hamiltonian. Their magnitude is proportional to the nanotube chiral angle and strain, and direction is determined by the nanotube quantum number. PL strain dependence measured for a number of various nanotube types allows to estimate TB carbon-carbon transfer integral.

Nikolaev, Pavel N.↗
NASA NTRS2007

Dynamic Theory of Relativistic Electrons Stochastic Heating by Whistler Mode Waves with Application to the Earth Magnetosphere

In the Hamiltonian approach an electron motion in a coherent packet of the whistler mode waves propagating along the direction of an ambient magnetic field is studied. The physical processes by which these particles are accelerated to high energy are established. Equations governing a particle motion were transformed in to a closed pair of nonlinear difference equations. The solutions of these equations have shown there exists the energetic threshold below that the electron motion is regular, and when the initial energy is above the threshold an electron moves stochastically. Particle energy spectra and pitch angle electron scattering are described by the Fokker-Planck-Kolmogorov equations. Calculating the stochastic diffusion of electrons due to a spectrum of whistler modes is presented. The parametric dependence of the diffusion coefficients on the plasma particle density, magnitude of wave field, and the strength of magnetic field is studies. It is shown that significant pitch angle diffusion occurs for the Earth radiation belt electrons with energies from a few keV up to a few MeV.

Khazanov, G. V.↗
NASA NTRS2006

Energy Levels and Intensity Parameters of Ho3(+) Ions in Y3Al5O12 and Lu3Al5O12

The energy levels of the trivalent lanthanide Ho(sup 3+) in Y3Al5O12 (YAG) and Lu3Al5O12 (LuAG) have been measured. The Stark split levels for the first nine Ho manifolds in these materials have been measured, and the results have been fit to a free ion plus crystal field Hamiltonian to generate a theoretical set of energy levels. Crystal field parameters were varied to determine the best fit between experimental and theoretical energy levels. The energy levels of Ho:LuAG are seen to be very similar to those in Ho:YAG. However, subtle changes resulting from replacing Y(sup 3+) with Lu(sup 3+) in the garnet crystal Y3Al5O12 result in different transition wavelengths in LuAG. This has implications for Ho (sup 5)I7yields (sup 5)I8 lasers operating at approximately 2.1 micrometers. Although the energy levels have been measured previously in Ho:YAG, they have not been measured in Ho:LuAG. A comparison of the energy levels in Ho:YAG measured here show some discrepancies with previous measurements. The consistency of the energy level placement between Ho:LuAG and Ho:YAG indicate that the earlier studies may have some errors in the assignments. Finally, a Judd-Ofelt analysis is performed on Ho:YAG and Ho:LuAG to determine the intensity parameters, and thus, the transition probabilities and branching ratios of the first eight excited manifolds.

Walsh, Brian M.↗
NASA NTRS2009

A Model for Predicting Thermoelectric Properties of Bi2Te3

A parameterized orthogonal tight-binding mathematical model of the quantum electronic structure of the bismuth telluride molecule has been devised for use in conjunction with a semiclassical transport model in predicting the thermoelectric properties of doped bismuth telluride. This model is expected to be useful in designing and analyzing Bi2Te3 thermoelectric devices, including ones that contain such nano - structures as quantum wells and wires. In addition, the understanding gained in the use of this model can be expected to lead to the development of better models that could be useful for developing other thermoelectric materials and devices having enhanced thermoelectric properties. Bi2Te3 is one of the best bulk thermoelectric materials and is widely used in commercial thermoelectric devices. Most prior theoretical studies of the thermoelectric properties of Bi2Te3 have involved either continuum models or ab-initio models. Continuum models are computationally very efficient, but do not account for atomic-level effects. Ab-initio models are atomistic by definition, but do not scale well in that computation times increase excessively with increasing numbers of atoms. The present tight-binding model bridges the gap between the well-scalable but non-atomistic continuum models and the atomistic but poorly scalable ab-initio models: The present tight-binding model is atomistic, yet also computationally efficient because of the reduced (relative to an ab-initio model) number of basis orbitals and flexible parameterization of the Hamiltonian.

Lee, Seungwon↗
NASA NTRS2009

Improved Model of a Mercury Ring Damper

A short document discusses the general problem of mathematical modeling of the three-dimensional rotational dynamics of rigid bodies and of the use of Euler parameters to eliminate the singularities occasioned by the use of Euler angles in such modeling. The document goes on to characterize a Hamiltonian model, developed by the authors, that utilizes the Euler parameters and, hence, is suitable for use in computational simulations that involve arbitrary rotational motion. In this formulation unlike in prior Euler-parameter-based formulations, there are no algebraic constraints. This formulation includes a general potential energy function, incorporates a minimum set of momentum variables, and takes an explicit state-space form convenient for numerical implementation. Practical application of this formulation has been demonstrated by the development of a new and simplified model of the rotational motion of a rigid rotor to which is attached a partially filled mercury ring damper. Models like this one are used in guidance and control of spin-stabilized spacecraft and gyroscope-stabilized seekers in guided missiles.

Fahrenthold, Eric P.↗
NASA NTRS2006

Using Quasiparticle Poisoning To Detect Photons

According to a proposal, a phenomenon associated with excitation of quasiparticles in certain superconducting quantum devices would be exploited as a means of detecting photons with exquisite sensitivity. The phenomenon could also be exploited to perform medium-resolution spectroscopy. The proposal was inspired by the observation that Coulomb blockade devices upon which some quantum logic gates are based are extremely sensitive to quasiparticles excited above the superconducting gaps in their leads. The presence of quasiparticles in the leads can be easily detected via the charge states. If quasiparticles could be generated in the leads by absorption of photons, then the devices could be used as very sensitive detectors of electromagnetic radiation over the spectral range from x-rays to submillimeter waves. The devices in question are single-Cooper-pair boxes (SCBs), which are mesoscopic superconducting devices developed for quantum computing. An SCB consists of a small superconducting island connected to a reservoir via a small tunnel junction and connected to a voltage source through a gate capacitor. An SCB is an artificial two-level quantum system, the Hamiltonian of which can be controlled by the gate voltage. One measures the expected value of the charge of the eigenvectors of this quantum system by use of a radio-frequency single-electron transistor. A plot of this expected value of charge as a function of gate voltage resembles a staircase that, in the ideal case, consists of steps of height 2 e (where e is the charge of one electron). Experiments have shown that depending on the parameters of the device, quasiparticles in the form of "broken" Cooper pairs present in the reservoir can tunnel to the island, giving rise to steps of 1 e. This effect is sometimes called "poisoning." Simulations have shown that an extremely small average number of quasiparticles can generate a 1-e periodic signal. In a device according to the proposal, this poisoning would be turned to advantage. Depending on the wavelength, an antenna or other component would be used to couple radiation into the reservoir, wherein the absorption of photons would break Cooper pairs, thereby creating quasiparticles that, in turn, would tunnel to the island, creating a 1-e signal. On the basis of conservative estimates of device parameters derived from experimental data and computational simulations that fit the data, it has been estimated that the noise equivalent power of a device according to the proposal could be as low as 6 10(exp -22) W/Hz(exp 1/2). It has also been estimated that the spectroscopic resolution (photon energy divided by increment of photon energy) of such a device in visible light would exceed 100.

Echternach, Pierre↗
NASA NTRS2011

Accurate Adaptive Level Set Method and Sharpening Technique for Three Dimensional Deforming Interfaces

In this paper, we demonstrate improved accuracy of the level set method for resolving deforming interfaces by proposing two key elements: (1) accurate level set solutions on adapted Cartesian grids by judiciously choosing interpolation polynomials in regions of different grid levels and (2) enhanced reinitialization by an interface sharpening procedure. The level set equation is solved using a fifth order WENO scheme or a second order central differencing scheme depending on availability of uniform stencils at each grid point. Grid adaptation criteria are determined so that the Hamiltonian functions at nodes adjacent to interfaces are always calculated by the fifth order WENO scheme. This selective usage between the fifth order WENO and second order central differencing schemes is confirmed to give more accurate results compared to those in literature for standard test problems. In order to further improve accuracy especially near thin filaments, we suggest an artificial sharpening method, which is in a similar form with the conventional re-initialization method but utilizes sign of curvature instead of sign of the level set function. Consequently, volume loss due to numerical dissipation on thin filaments is remarkably reduced for the test problems

Kim, Hyoungin↗
NASA NTRS2003

Optimal Feedback Control of Thermal Networks

An improved approach to the mathematical modeling of feedback control of thermal networks has been devised. Heretofore software for feedback control of thermal networks has been developed by time-consuming trial-and-error methods that depend on engineers expertise. In contrast, the present approach is a systematic means of developing algorithms for feedback control that is optimal in the sense that it combines performance with low cost of implementation. An additional advantage of the present approach is that a thermal engineer need not be expert in control theory. Thermal networks are lumped-parameter approximations used to represent complex thermal systems. Thermal networks are closely related to electrical networks commonly represented by lumped-parameter circuit diagrams. Like such electrical circuits, thermal networks are mathematically modeled by systems of differential-algebraic equations (DAEs) that is, ordinary differential equations subject to a set of algebraic constraints. In the present approach, emphasis is placed on applications in which thermal networks are subject to constant disturbances and, therefore, integral control action is necessary to obtain steady-state responses. The mathematical development of the present approach begins with the derivation of optimal integral-control laws via minimization of an appropriate cost functional that involves augmented state vectors. Subsequently, classical variational arguments provide optimality conditions in the form of the Hamiltonian equations for the standard linear-quadratic-regulator (LQR) problem. These equations are reduced to an algebraic Riccati equation (ARE) with respect to the augmented state vector. The solution of the ARE leads to the direct computation of the optimal proportional- and integral-feedback control gains. In cases of very complex networks, large numbers of state variables make it difficult to implement optimal controllers in the manner described in the preceding paragraph.

Papalexandris, Miltiadis↗
NASA NTRS2003

Numerical Modeling of Nanoelectronic Devices

Nanoelectronic Modeling 3-D (NEMO 3-D) is a computer program for numerical modeling of the electronic structure properties of a semiconductor device that is embodied in a crystal containing as many as 16 million atoms in an arbitrary configuration and that has overall dimensions of the order of tens of nanometers. The underlying mathematical model represents the quantummechanical behavior of the device resolved to the atomistic level of granularity. The system of electrons in the device is represented by a sparse Hamiltonian matrix that contains hundreds of millions of terms. NEMO 3-D solves the matrix equation on a Beowulf-class cluster computer, by use of a parallel-processing matrix vector multiplication algorithm coupled to a Lanczos and/or Rayleigh-Ritz algorithm that solves for eigenvalues. In a recent update of NEMO 3-D, a new strain treatment, parameterized for bulk material properties of GaAs and InAs, was developed for two tight-binding submodels. The utility of the NEMO 3-D was demonstrated in an atomistic analysis of the effects of disorder in alloys and, in particular, in bulk In(x)Ga(l-x)As and in In0.6Ga0.4As quantum dots.

Klimeck, Gerhard↗
NASA NTRS2003

Quantum Search in Hilbert Space

A proposed quantum-computing algorithm would perform a search for an item of information in a database stored in a Hilbert-space memory structure. The algorithm is intended to make it possible to search relatively quickly through a large database under conditions in which available computing resources would otherwise be considered inadequate to perform such a task. The algorithm would apply, more specifically, to a relational database in which information would be stored in a set of N complex orthonormal vectors, each of N dimensions (where N can be exponentially large). Each vector would constitute one row of a unitary matrix, from which one would derive the Hamiltonian operator (and hence the evolutionary operator) of a quantum system. In other words, all the stored information would be mapped onto a unitary operator acting on a quantum state that would represent the item of information to be retrieved. Then one could exploit quantum parallelism: one could pose all search queries simultaneously by performing a quantum measurement on the system. In so doing, one would effectively solve the search problem in one computational step. One could exploit the direct- and inner-product decomposability of the unitary matrix to make the dimensionality of the memory space exponentially large by use of only linear resources. However, inasmuch as the necessary preprocessing (the mapping of the stored information into a Hilbert space) could be exponentially expensive, the proposed algorithm would likely be most beneficial in applications in which the resources available for preprocessing were much greater than those available for searching.

Zak, Michail↗