Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hamiltonian structure”

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.

72 records · Page 4

Variational calculations of vibrational properties of ozone

A variational method is used to obtain vibrational-rotational properties for ozone from an experimental quartic force field. Band positions, average structures, matrix elements for calculating infrared intensities, and effective rotational constants are reported for (O-16)3 and its O-18 isotopic species. Also, the degree to which the vibrational energies and properties are converged is investigated as a function of the basis set parameters and basis set size, and of the method of obtaining the Hamiltonian matrix elements. A convenient procedure for assigning the vibrational states is developed for variational vibrational wavefunction expansions.

Carney, G. D.↗

The rotation-inversion spectra and vibration-rotation interaction in NH2D

An additional 86 lines of the NH2D rotation-inversion spectrum were measured between 7.5 and 850 GHz. These data and the existing NH2D and ND2H data are fitted to a Hamiltonian in which the vibration-rotation interaction is included in the form of an off-diagonal inertial tensor term. For comparison with the symmetric ammonias an oblate basis and a Watson-type reduced 'S' set of distortion constants are used. A theory for the interaction term is given which is in excellent agreement with experiment. Structural parameters are derived. The magnitude and relative sign of the two dipole moment components of NH2D are determined and previously measured nitrogen hyperfine structure for both species is reexamined using improved rotational wavefunctions.

Cohen, E. A.↗

Finite element solution of optimal control problems with state-control inequality constraints

It is demonstrated that the weak Hamiltonian finite-element formulation is amenable to the solution of optimal control problems with inequality constraints which are functions of both state and control variables. Difficult problems can be treated on account of the ease with which algebraic equations can be generated before having to specify the problem. These equations yield very accurate solutions. Owing to the sparse structure of the resulting Jacobian, computer solutions can be obtained quickly when the sparsity is exploited.

Bless, Robert R.↗

Composition-property correlations in NiTi-based shape memory alloys from the first principles

Development of shape memory alloys (SMA) has been traditionally accomplished by means of extensive empirical efforts, with a very limited support from phenomenological models and related tools. We present a theoretical framework based on the Hamiltonian formalism, which connects important thermodynamic properties directly to the relative energies (the eigenvalues of the Hamiltonian), computed from the first principles using density functional theory. The estimates based on this formalism allow to establish correlations between composition and physical properties, such as relative energies, phase transition temperatures, hysteresis, thermal expansion, and other thermodynamic information. Importantly, estimates and correlations not only provide a fast answer, but also elucidate underpinning of thermodynamic features in terms of the electronic structure. Using this framework, we consider phase transformations and correlations between the selected properties and composition in the NiTi-based ternary and multicomponent shape memory alloys. Our theoretical guidance facilitates design and development of future alloys.

shape memory alloys↗

Oscillator Strengths for Fine-Structure Transitions in S III

Oscillator strengths and transition probabilities for transitions among the fine-structure levels of the terms belonging to the 3s(sup 2)3p(sup 2), 3s3p(sup 3), 3s(sup 2)3p3d, 3s(sup 2)3p4s, 3s(sup 2)3p4p, and 3s(sup 2)3p4d configurations of S III are calculated using extensive configuration-interaction wave functions. The relativistic effects in intermediate coupling are incorporated by means of the Breit-Pauli Hamiltonian. Small adjustments to the diagonal elements of the Hamiltonian matrices have been made so that the energy splittings are as close as possible to the experimental values. The present results are compared with other available calculations and experiments.

Tayal, S. S.↗

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↗

A general algorithm for solving the algebraic Riccati equation

The generalized eigenvalue problem provides a suitable framework for reliable solutions of many system theoretic, control, and estimation problems. A general algorithm for solving the matrix algebraic Riccati equation (ARE) which utilizes a pencil structure is described here. This algorithm avoids unnecessary inversion of cost or transition matrices, making it a numerically sound way to solve for the gains and/or ARE with singular quadratic costs, for cases satisfying detectability and stabilizability conditions. Examples are solution with discrete dead-beat control, noiseless measurements in Kalman filters and time-delays in discrete-time systems, which cause difficulties in the Hamiltonian standard eigenvalue problem formulation. The ARE algorithm implementatiton and numerical examples are shown.

Walker, R. A.↗

The infrared and microwave spectroscopy of the argon-propyne dimer

Microwave and infrared measurements are reported for the van der Waals complex ArCH3CCH and its isotopomers. The structure is T-shaped with equilibrium center-of-mass separation of 3.73 A and an angle of 82 deg between the molecule symmetry axis and the van der Waals bond. The infrared and microwave spectra are complex due to the effects of a slightly hindered internal rotor. Analysis of the spectral data shows that the dipole moment is almost parallel to the dimer b axis. A modified vibration/rotation Hamiltonian that includes an internal rotor potential is used to show that the barrier to internal rotation is near 10.8/cm.

Blake, Thomas A.↗

Atomic Data and Spectral Line Intensities for Ne III

A number of satellites and rockets have been launched to observe radiation from the Sun and other astrophysical objects. Line radiation is emitted when the electron impact excited levels decay to the lower levels by photon emission. From this radiation, the physical parameters such as electron temperature and density of the astrophysical plasma, elemental abundance, and opacity can be inferred. Ne III lines have been observed in H II regions, Ne-rich filaments in supernovae, and planetary nebulae. The allowed line at 489.50 Angstroms due to the transition 2s(sup 2) 2p(sup 5) (sup 3) P2 (goes to) 2s(sup 2)2p(sup 4)(sup 3)P2 has been identified in the solar spectrum by Vernazza and Reeves using Skylab observations. Other Ne III lines in the solar EUV spectrum have been reported by Thomas and Neupert based on observations from the Solar EUV Rocket Telescope and Spectrograph (SERTS) instrument. Atomic data for Ne III have been calculated by using a set of programs developed at, University College, London. The Superstructure and Distorted Wave (DW) programs have been updated over the years. In the Superstructure program, configuration interaction can be taken into account and radial functions are calculated in a modified Thomas-Fermi-Amaldi potential. This is a statistical potential and depends on parameters lambda 1 which are determined by optimizing the weighted sum of term energies. They are found to be lambda(sub 0)=1.2467, lambda(sub 1)=1.1617, and lambda(sub 2)=1.0663. The relativistic corrections are included by using the Breit-Pauli Hamiltonian as a perturbation to the nonrelativistic Hamiltonian. The same potential is used to calculate reactance matrices in the DW approximation in LS coupling. Collision strengths in intermediate coupling are obtained by using term coupling coefficients obtained from the Superstructure program. In this calculation, the configurations used are 2s(sup 2)2p(sup 4), 2s2p(sup 5), 2s(sup 2)2p(sup 3)3s, 2s(sup 2)p(sup 3)3d giving rise to 57 fine-structure levels in intermediate coupling.

Bhatia, A. K.↗

Symmetry and Error Mitigation in Quantum Alternating Operator Ansatz

Quantum Alternating Operator Ansatz (QAOA) provides a platform for probing quantum advantage in optimization problems on NISQ hardware. By exploring symmetries in the system and the structure of unitaries in the algorithm, we proposed using problem-specific mixing operators to limit the quantum evolution within a subspace where the target state lives. As an example, in case of equality constraint, we studied the use of XY operators that preserve total Pauli Z. This QAOA variant, in the noiseless case, significantly outperforms the vanilla version of QAOA where the constraints are encoded as penalty in the cost Hamiltonian and the search is in the whole Hilbert space. In presence of noise, although the advantage still holds, error mitigation is needed for QAOA. I will overview our work on the design of QAOA in presence of symmetry, the implementation and compilation of XY-QAOA, its performance and error mitigation techniques in presence of noise.

Zhihui Wang↗

Oscillator Strengths of Allowed and Intercombination Transitions in Neutral Sulfur

We have calculated oscillator strengths and transition probabilities of electric-dipole allowed and intercombination transitions from fine-structure levels of the ground 3s(sup 2)3p(sup 4) configuration to the levels belonging to configurations 3s(sup 2)3p(sup 3)4s, 3s(sup 2) 3p(sup 3)5s, 3(sup 2)3p(sup 3)3d, 3s(sup 2)3p(sup 3)4d of neutral sulfur. Extensive configuration-interaction wave functions are used to represent these levels. The relativistic corrections have been included through the Breit-Pauli Hamiltonian. The results are compared with previous theoretical calculations and with measurements.

Tayal, S. S.↗

Approximate constants of motion for classically chaotic vibrational dynamics - Vague tori, semiclassical quantization, and classical intramolecular energy flow

Substantial short time regularity, even in the chaotic regions of phase space, is found for what is seen as a large class of systems. This regularity manifests itself through the behavior of approximate constants of motion calculated by Pade summation of the Birkhoff-Gustavson normal form expansion; it is attributed to remnants of destroyed invariant tori in phase space. The remnant torus-like manifold structures are used to justify Einstein-Brillouin-Keller semiclassical quantization procedures for obtaining quantum energy levels, even in the absence of complete tori. They also provide a theoretical basis for the calculation of rate constants for intramolecular mode-mode energy transfer. These results are illustrated by means of a thorough analysis of the Henon-Heiles oscillator problem. Possible generality of the analysis is demonstrated by brief consideration of classical dynamics for the Barbanis Hamiltonian, Zeeman effect in hydrogen and recent results of Wolf and Hase (1980) for the H-C-C fragment.

Shirts, R. B.↗

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↗

Application of Gauss-Seidel multilevel control to a single axis torsional model

An approach to the application of multilevel control techniques for large space structures is presented. Gauss-Seidel second level controls and an extension of standard linear quadratic regulator techniques are applied to a model consisting of a flexible space vehicle comprising three rigid bodies rotating about a common axis. The method incorporates second order derivatives with respect to time, entailing the use of a two level hierarchy of subsystems. Gauss-Seidel second level control formulations were chosen to avoid the necessity of having as many controls as constraints or using gradient techniques in the choice of the Hamiltonian. The modular nature of the resulting control system allows applications for spacecraft with larger numbers of modules.

Chichester, F. D.↗

Quantum chemical calculation of the equilibrium structures of small metal atom clusters

The application of ab initio quantum mechanical approaches in the study of metal atom clusters requires simplifying techniques that do not compromise the reliability of the calculations. Various aspects of the implementation of the effective core potential (ECP) technique for the removal of the metal atom core electrons from the calculation were examined. The ECP molecular integral formulae were modified to bring out the shell characteristics as a first step towards fulfilling the increasing need to speed up the computation of the ECP integrals. Work on the relationships among the derivatives of the molecular integrals that extends some of the techniques pioneered by Komornicki for the calculation of the gradients of the electronic energy was completed and a formulation of the ECP approach that quite naturally unifies the various state-of-the-art "shape- and Hamiltonian-consistent" techniques was discovered.

Kahn, L. R.↗

On The Computational Capabilities of Physical Systems: Relationship With Conventional Computer Science - Part 2

In the first of this pair of papers, it was proven that there cannot be a physical computer to which one can properly pose any and all computational tasks concerning the physical universe. It was then further proven that no physical computer C can correctly carry out all computational tasks that can be posed to C. As a particular example, this result means that no physical computer that can, for any physical system external to that computer, take the specification of that external system's state as input and then correctly predict its future state before that future state actually occurs; one cannot build a physical computer that can be assured of correctly "processing information faster than the universe does". These results do not rely on systems that are infinite, and/or non-classical, and/or obey chaotic dynamics. They also hold even if one uses an infinitely fast, infinitely dense computer, with computational powers greater than that of a Turing Machine. This generality is a direct consequence of the fact that a novel definition of computation - "physical computation" - is needed to address the issues considered in these papers, which concern real physical computers. While this novel definition does not fit into the traditional Chomsky hierarchy, the mathematical structure and impossibility results associated with it have parallels in the mathematics of the Chomsky hierarchy. This second paper of the pair presents a preliminary exploration of some of this mathematical structure. Analogues of Chomskian results concerning universal Turing Machines and the Halting theorem are derived, as are results concerning the (im)possibility of certain kinds of error-correcting codes. In addition, an analogue of algorithmic information complexity, "prediction complexity", is elaborated. A task-independent bound is derived on how much the prediction complexity of a computational task can differ for two different reference universal physical computers used to solve that task, a bound similar to the "encoding" bound governing how much the algorithm information complexity of a Turing machine calculation can differ for two reference universal Turing machines. Finally, it is proven that either the Hamiltonian of our universe proscribes a certain type of computation, or prediction complexity is unique (unlike algorithmic information complexity), in that there is one and only version of it that can be applicable throughout our universe.

Wolpert, David H.↗

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↗

A semiempirical study of the optimized ground and excited state potential energy surfaces of retinal and its protonated Schiff base

The electronic ground and first excited states of retinal and its Schiff base are optimized for the first time using the semiempirical AM1 Hamiltonian. The barrier for rotation about the C(11)-C(12) double bond is characterized by variation of both the twist angle delta(C(10)-C(11)-C(12)-C(13)) and the bond length d(C(11)-C(12)). The potential energy surface is obtained by varying these two parameters. The calculated ground state rotational barrier is equal to 15.6 kcal/mol for retinal and 20.5 kcal/mol for its Schiff base. The all-trans conformation is more stable by 3.7 kcal/mol than the 11-cis geometry. For the first excited state, S(1,) the 90 degrees twisted geometry represents a saddle point for retinal with the rotational barrier of 14.6 kcal/mol. In contrast, this conformation is an energy minimum for the Schiff base. It can be easily reached at room temperature from the planar minima since it is separated from them by a barrier of only 0.6 kcal/mol. The 90 degrees minimum conformation is more stable than the all-trans by 8.6 kcal/mol. We are thus able to present a reaction path on the S(1) surface of the retinal Schiff base with an almost barrier-less geometrical relaxation into a twisted minimum geometry, as observed experimentally. The character of the ground and first excited singlet states underscores the need for the inclusion of double excitations in the calculations.

Retinaldehyde/chemistry↗