Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCalculating eigenvalues and eigenvectors of parameter-dependent hamiltonians using an adaptative wave operator method
272
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We propose a wave operator method to calculate eigenvalues and eigenvectors of large parameter-dependent matrices, using an adaptative active subspace. We consider a hamiltonian which depends on external adjustable or adiabatic parameters, using adaptative projectors which follow the successive eigenspaces when the adjustable parameters are modified. The method can also handle non-hermitian hamiltonians. An iterative algorithm is derived and tested through comparisons with a standard wave operator algorithm using a fixed active space and with a standard block-Davidson method. The proposed approach is competitive, it converges within a few dozen iterations at constant memory cost. We first illustrate the abilities of the method on a 4-D coupled oscillator model hamiltonian. A more realistic application to molecular photodissociation under intense laser fields with varying intensity or frequency is also presented. Maps of photodissociation resonances of H${}_2^+$ in the vicinity of exceptional points are calculated as an illustrative example.
rate research
Read More
Suppose we want to find the eigenvalues and eigenvectors for the linear operator L, and suppose that we have solved this problem for some other nearby operator K. In this paper we show how to represent the eigenvalues and eigenvectors of L in terms of the corresponding properties of K.
The spectral renormalization method was introduced in 2005 as an effective way to compute ground states of nonlinear Schrodinger and Gross-Pitaevskii type equations. In this paper, we introduce an orthogonal spectral renormalization (OSR) method to compute ground and excited states (and their respective eigenvalues) of linear and nonlinear eigenvalue problems. The implementation of the algorithm follows four simple steps: (i) reformulate the underlying eigenvalue problem as a fixed point equation, (ii) introduce a renormalization factor that controls the convergence properties of the iteration, (iii) perform a Gram-Schmidt orthogonalization process in order to prevent the iteration from converging to an unwanted mode; and (iv) compute the solution sought using a fixed-point iteration. The advantages of the OSR scheme over other known methods (such as Newtons and self-consistency) are: (i) it allows the flexibility to choose large varieties of initial guesses without diverging, (ii) easy to implement especially at higher dimensions and (iii) it can easily handle problems with complex and random potentials. The OSR method is implemented on benchmark Hermitian linear and nonlinear eigenvalue problems as well as linear and nonlinear non-Hermitian $mathcal{PT}$-symmetric models.
In previous works it was shown that protein 3D-conformations could be encoded into discrete sequences called dominance partition sequences (DPS), that generated a linear partition of molecular conformational space into regions of molecular conformations that have the same DPS. In this work we describe procedures for building in a cubic lattice the set of 3D-conformations that are compatible with a given DPS. Furthermore, this set can be structured as a graph upon which a combinatorial algorithm can be applied for computing the mean energy of the conformations in a cell.
A method is developed that allows analysis of quantum Monte Carlo simulations to identify errors in trial wave functions. The purpose of this method is to allow for the systematic improvement of variational wave functions by identifying degrees of freedom that are not well-described by an initial trial state. We provide proof of concept implementations of this method by identifying the need for a Jastrow correlation factor, and implementing a selected multi-determinant wave function algorithm for small dimers that systematically decreases the variational energy. Selection of the two-particle excitations is done using quantum Monte Carlo within the presence of a Jastrow correlation factor, and without the need to explicitly construct the determinants. We also show how this technique can be used to design compact wave functions for transition metal systems. This method may provide a route to analyze and systematically improve descriptions of complex quantum systems in a scalable way.
A global solution of the Schrodinger equation for explicitly time-dependent Hamiltonians is derived by integrating the non-linear differential equation associated with the time-dependent wave operator. A fast iterative solution method is proposed in which, however, numerous integrals over time have to be evaluated. This internal work is done using a numerical integrator based on Fast Fourier Transforms (FFT). The case of a transition between two potential wells of a model molecule driven by intense laser pulses is used as an illustrative example. This application reveals some interesting features of the integration technique. Each iteration provides a global approximate solution on grid points regularly distributed over the full time propagation interval. Inside the convergence radius, the complete integration is competitive with standard algorithms, especially when high accuracy is required.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
