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.
Quantum Algorithms for Hamiltonian Simulation
Explore the source record for details and available documents.
Theory of multiple-stellar population synthesis in a non-Hamiltonian setting
Explore the source record for details and available documents.
Third post-Newtonian effective-one-body Hamiltonian in scalar-tensor and Einstein-scalar-Gauss-Bonnet gravity
Explore the source record for details and available documents.
Simultaneous Stoquasticity
Stoquastic Hamiltonians play a role in the computational complexity of the local Hamiltonian problem as well as the study of classical simulability. In particular, stoquastic Hamiltonians can be straightforwardly simulated using Monte Carlo techniques. We address the question of whether two or more Hamiltonians may be made simultaneously stoquastic via a unitary transformation. This question has important implications for the complexity of simulating quantum annealing where quantum advantage is related to the stoquasticity of the Hamiltonians involved in the anneal. We find that for almost all problems no such unitary exists and show that the problem of determining the existence of such a unitary is equivalent to identifying if there is a solution to a system of polynomial (in)equalities in the matrix elements of the initial and transformed Hamiltonians. Solving such a system of equations is NP-hard. We highlight a geometric understanding of this problem in terms of a collection of generalized Bloch vectors.
Simultaneous Stoquasticity
Stoquastic Hamiltonians play a role in the computational complexity of the local Hamiltonian problem as well as the study of classical simulability. In particular, stoquastic Hamiltonians can be straightforwardly simulated using Monte Carlo techniques. We address the question of whether two or more Hamiltonians may be made simultaneously stoquastic via a unitary transformation. This question has important implications for the complexity of simulating quantum annealing where quantum advantage is related to the stoquasticity of the Hamiltonians involved in the anneal. We find that for almost all problems no such unitary exists and show that the problem of determining the existence of such a unitary is equivalent to identifying if there is a solution to a system of polynomial (in)equalities in the matrix elements of the initial and transformed Hamiltonians. Solving such a system of equations is NP-hard. We highlight a geometric understanding of this problem in terms of a collection of generalized Bloch vectors.
Forces acting on free carriers in semiconductors of inhomogeneous composition I, II
An investigation shows that the extended Wannier-Slater theorem with its position-dependent effective mass m cannot be considered a viable description of a semiconductor with compositional changes. The correct description of the behavior of the electron or hole flow in semiconductors in which the chemical composition changes gradually from position to position is more difficult. It is necessary to find a model Hamiltonian. This model Hamiltonian is to some extent arbitry, but it must be representative of a homogeneous sample. The difference between the true and the model Hamiltonian should be small in the sense that this difference may be considered a perturbation. However, a Hamiltonian with the required specifications can always be found provided that the compositional inhomogeneities do not change drastically over the distances of interest.
Alternative descriptions of wave and particle aspects of the harmonic oscillator
The dynamical properties of the wave and particle aspects of the harmonic oscillator can be studied with the help of the time-dependent Schroedinger equation (SE). Especially the time-dependence of maximum and width of Gaussian wave packet solutions allow to show the evolution and connections of those two complementary aspects. The investigation of the relations between the equations describing wave and particle aspects leads to an alternative description of the considered systems. This can be achieved by means of a Newtonian equation for a complex variable in connection with a conservation law for a nonclassical angular momentum-type quantity. With the help of this complex variable, it is also possible to develop a Hamiltonian formalism for the wave aspect contained in the SE, which allows to describe the dynamics of the position and momentum uncertainties. In this case the Hamiltonian function is equivalent to the difference between the mean value of the Hamiltonian operator and the classical Hamiltonian function.
Relativistic and Correlation Effects in CuH, AgH and AuH: Comparison of Various Relativistic Methods
The effects of relativity on the bond lengths, dissociation energies, and harmonic vibrational frequencies of the 1Epsilon(+) electronic ground states of the group IB hydrides CuH, AgH and AuH have been evaluated with a variety of ab initio methods. These properties were investigated with moderately-sized basis sets at the self-consistent field Hartree Fock (SCF HF) level and with second-order Moller-Plesset (MP2) perturbation theory for electron correlation. Comparisons were made between all-electron results using the nonrelativistic Hamiltonian, perturbation theory (PT) at first-order with only the one-electron non-fine structure terms of the Breit-Pauli Hamiltonian, the spin-free Douglas-Kroll (DK) transformed Dirac Hamiltonian and the untransformed Dirac Hamiltonian, and results using two sets of relativistic effective core potentials (RECPs). The expected trends of bond length decrease, dissociation energy increase and harmonic frequency increase with both relativity and correlation are found. Both sets of RECPs are shown to give good results, if accompanied by a reasonable basis set. The DK method is demonstrated to be an inexpensive, reliable approximation to the DHF method.
Charged particle motion in the vicinity of a neutral plane
The motion of a charged particle in a 2-dimensional neutral sheet with linear magnetic field variation is analyzed by Hamiltonian methods. From the basic Hamiltonian formulation, results of Sonnerup are recovered and given intuitive interpretation. The transformed Hamiltonian is used to derive the correct frequency of oscillation and serves as the basis of analytical treatment of perturbed versions of the motion, e.g. with electric field added. For 2-dimensional fields with slightly differing configurations - in particular, with added small field component orthogonal to the sheet, with an X-type null or with a chain of null points - two alternative methods are developed, reducing the problem either to motion in a 2-dimensional potential or to a pair of coupled oscillators.
Redundancy of constraints in the classical and quantum theories of gravitation.
It is shown that in Dirac's version of the quantum theory of gravitation, the Hamiltonian constraints are greatly redundant. If the Hamiltonian constraint condition is satisfied at one point on the underlying, closed three-dimensional manifold, then it is automatically satisfied at every point, provided only that the momentum constraints are everywhere satisfied. This permits one to replace the usual infinity of Hamiltonian constraints by a single condition which may be taken in the form of an integral over the manifold. Analogous theorems are given for the classical Einstein Hamilton-Jacobi equations.
An EPR investigation of the dynamic Jahn-Teller effect in SrCl2:y(2 plus) and SrCl2:Sc(2 plus)
EPR spectra have been observed for SrCl2:Y(2+) and SrCl2:Sc(2+) at liquid helium temperatures. At 1.2 K the spectra were dominated by anisotropic hyperfine patterns whose lineshapes and angular dependences were explained using second order solutions of the effective Hamiltonian for an isolated 2Eg state split by large random internal strains. Pronounced asymmetries in some of the strin produced lineshapes for Srcl2:Sc(2+) are shown to result from second order terms in the solution of the effective Hamiltonian. Coexisting with the anisotropic hyperfine patterns are weak nearly isotropic hyperfine patterns with typical lineshapes. Variations in the apparent intensity of lines in these weak hyperfine patterns as functions of the applied magnetic field direction and temperature imply that these lines result from averaging by vibronic relaxation of a portion of the anisotropic pattern. The effective Hamiltonian parameters for SrCl2:La(2+), SrCl2:y(2+), and SrCl2:SC(2+) are analyzed in terms of crystal field theory modified to include a dynamic Jahn-Teller effect.
Three-electron atoms from a Bethe-Salpeter equation approach
The two-electron Bethe-Salpeter equation is reviewed. The transition to a Hamiltonian form is discussed, and the criterion that the first-order energy shift vanish is used to determine the electron-electron potential that appears in the Hamiltonian. A similar program is carried through for the three-electron case. The fundamental idea of Feldman and Fulton is used for the transition to the Hamiltonian form, although some changes are necessary. The optimum two-electron interaction is the same as that in the two-electron case.
On a Lagrange-Hamilton formalism describing position and momentum uncertainties
According to Heisenberg's uncertainty relation, in quantum mechanics it is not possible to determine, simultaneously, exact values for the position and the momentum of a material system. Calculating the mean value of the Hamiltonian operator with the aid of exact analytic Gaussian wave packet solutions, these uncertainties cause an energy contribution additional to the classical energy of the system. For the harmonic oscillator, e.g., this nonclassical energy represents the ground state energy. It will be shown that this additional energy contribution can be considered as a Hamiltonian function, if it is written in appropriate variables. With the help of the usual Lagrange-Hamilton formalism known from classical particle mechanics, but now considering this new Hamiltonian function, it is possible to obtain the equations of motion for position and momentum uncertainties.
Understanding squeezing of quantum states with the Wigner function
The Wigner function is argued to be the only natural phase space function evolving classically under quadratic Hamiltonians with time-dependent bilinear part. This is used to understand graphically how certain quadratic time-dependent Hamiltonians induce squeezing of quantum states. The Wigner representation is also used to generalize Ehrenfest's theorem to the quantum uncertainties. This makes it possible to deduce features of the quantum evolution, such as squeezing, from the classical evolution, whatever the Hamiltonian.
Research on Spectroscopy, Opacity, and Atmospheres
I propose to continue providing observers with basic data for interpreting spectra from stars, novas, supernovas, clusters, and galaxies. These data will include allowed and forbidden line lists, both laboratory and computed, for the first five to ten ions of all atoms and for all relevant diatomic molecules. I will eventually expand to all ions of the first thirty elements to treat far UV and X-ray spectra, and for envelope opacities. I also include triatomic molecules provided by other researchers. I have also made CDs with Partridge and Schwenke's water data for work on UV stars. The line data also serve as input to my model atmosphere and synthesis programs that generate energy distributions, photometry, limb darkening, and spectra that can be used for planning observations and for fitting observed spectra. The spectrum synthesis programs produce detailed plots with the lines identified. Grids of stellar spectra can be used for radial velocity-, rotation-, or abundance templates and for population synthesis. I am fitting spectra of bright stars to test the data and to produce atlases to guide observers. For each star the whole spectrum is computed from the UV to the far IR. The line data, opacities, models, spectra, and programs are freely distributed on CDs and on my Web site and represent a unique resource for many NASA programs. I am now in full production of new line lists for atoms. I am computing all ions of all elements from H to Zn and the first 5 ions of all the heavier elements, about 800 ions. For each ion I treat as many as 61 even and 61 odd configurations, computing all energy levels and eigenvectors. The Hamiltonian is determined from a scaled-Hartree-Fock starting guess by least squares fitting the observed energy levels. The average energy of each configuration is used in computing scaled-Thomas-Fermi-Dirac wavefunctions for each configuration which in turn are used to compute allowed and forbidden transition integrals. These are multiplied into the LS allowed and forbidden transition arrays. The transition arrays are transformed to the observed coupling to yield the allowed and forbidden line lists. Results are put on the web as they are finished. Provided I get funding,there will be more than 500 million lines. I will then compare ion by ion, to all the laboratory and computed data in the literature and make up a working line list for spectrum synthesis and opacity calculations with the best available data. As the laboratory spectrum analyses are improved, I will redo the calculations with the new energy levels. My original plan when I started the new calculations was to run through all the atoms using my old Cray programs from the 1980's that were limited to 1100 x 1100 arrays in the Hamiltonian for each J. Then I would go back and rerun the more complicated cases with 3000 x 3000 arrays so that I could include many more configurations and more configuration interactions. At present I am limited to 61 even and 61 odd configurations and I try to include everything up through n = 9. The current program runs on Alpha workstations. I decided to test the big program on Fe I and Fe II to see whether there was any great difference in the low configurations compared to those from the Cray program. Besides increasing the number of E1 lines by a factor of 6 to 7.7 million, there was an unexpected result: the electric quadrupole transitions were 10 times stronger than before because the transition integrals are weighted by r(exp 2) ---they become very large for high n, and because there are numerous configuration interactions that mix the low and high configurations. As a check I was able to reproduce Carstang's (1962) lower results by running his three configurations with my program. Since my model atom is still only a subset of a real Fe II ion, the true quadrupole A values are probably larger than mine. The magnetic dipole lines are affected by the mixing but the overall scale does not change. Because of this scovery I decided that there was no point in computing the small array cases. I have been running with as many configurations as I can and with thousands of parameters in the Hamiltonian. The computer runs take much longer to set up and produce than I had expected. I have concentrated on redoing the low iron group spectra, especially to get data for supernova modelers. I have done only Ca I -- Zn I, Ca II -- Zn II, CU I -- Cu XXIX, Zn I - Zn XXX, for practice at high stages of ionization, C I, C II, S I, and CL I and Ag I for people who were working on the laboratory spectra. Check my web site kurucz.harvard.edu for current additions. My latest calculations have been for carbon I and sulphur I, and silicon I is under way using the same elaborate approach as for C I, which took many months to do. These line lists greatly increase the number of lines in the ultraviolet, in the visible, and especially in the infrared. They will increase the opacity in A, F, and G stars. They will account for many unidentified lines in the sun.
Numerical Schemes for the Hamilton-Jacobi and Level Set Equations on Triangulated Domains
Borrowing from techniques developed for conservation law equations, numerical schemes which discretize the Hamilton-Jacobi (H-J), level set, and Eikonal equations on triangulated domains are presented. The first scheme is a provably monotone discretization for certain forms of the H-J equations. Unfortunately, the basic scheme lacks proper Lipschitz continuity of the numerical Hamiltonian. By employing a virtual edge flipping technique, Lipschitz continuity of the numerical flux is restored on acute triangulations. Next, schemes are introduced and developed based on the weaker concept of positive coefficient approximations for homogeneous Hamiltonians. These schemes possess a discrete maximum principle on arbitrary triangulations and naturally exhibit proper Lipschitz continuity of the numerical Hamiltonian. Finally, a class of Petrov-Galerkin approximations are considered. These schemes are stabilized via a least-squares bilinear form. The Petrov-Galerkin schemes do not possess a discrete maximum principle but generalize to high order accuracy.