Academia.edu no longer supports Internet Explorer.
To browse Academia.edu and the wider internet faster and more securely, please take a few seconds to upgrade your browser.
2004, Physical Review A
We devise a non-Hermitian Rayleigh-Schrödinger perturbation theory for the single-and the multireference case to tackle both the many-body problem and the decay problem encountered, for example, in the study of electronic resonances in molecules. A complex absorbing potential (CAP) is employed to facilitate a treatment of resonance states that is similar to the well-established boundstate techniques. For the perturbative approach, the full CAP-Schrödinger Hamiltonian, in suitable representation, is partitioned according to the Epstein-Nesbet scheme. The equations we derive in the framework of the single-reference perturbation theory turn out to be identical to those obtained by a time-dependent treatment in Wigner-Weisskopf theory. The multireference perturbation theory is studied for a model problem and is shown to be an efficient and accurate method. Algorithmic aspects of the integration of the perturbation theories into existing ab initio programs are discussed, and the simplicity of their implementation is elucidated.
Theoretical Chemistry Accounts, 2016
The partitioning theory provides optical potentials which lead to the resonance energies. The state of the theory in the eighties can be found in the book "Theory of resonances" of Kukulin et al. [3]. The basic principles presented in this book still hold: use of projectors and analytical continuation. Since that time, the theoretical developments have been influenced by the fast development of computers and the theory entered into the quickly expanding field of non-Hermitian quantum mechanics [4]. Jolicard and Austin [5] incorporated optical potentials in computational schemes, and the justification of their method was given by Riss and Meyer [6]. To avoid confusion, we do not use the term optical potential here, but adopt the more suitable expression complex absorbing potential (CAP), as proposed in Ref. [6]. The acronym CAP refers to an energy-independent complex potential added to the Hamiltonian (see review article [7]). For the past two decades, we contributed to the development of approaches combining CAP and perturbation theory [8-12]. In this article, we recognize a resonance as a pole of a Green's function. In Sect. 2, we first recall the definition of the wave function and the energy in the framework of the partitioning technique. Then, a complex absorbing potential written in the form V CAP = −ı ǫ is added to the Hamiltonian. The operator ǫ generalizes the quantity ǫ of collision theory. The parameter in V CAP anticipates that V CAP does not only produce analytic continuation, but that it is also a perturbation operator. The self-energy is expanded in powers of , and the equivalence between this expansion and a Taylor series is demonstrated in Sect. 3 which focuses on our new findings. Finally, the introduction of a convergence operator allows to discuss the convergence properties of these series. Section 4 is devoted to numerical illustrations. First, a discretized N-dimensional Fano model is used to check the accuracy and the convergence properties of the energies. It is shown Abstract We propose a contribution to the theory of quantum resonances that combines complex absorbing potentials (CAP) with standard perturbation theory. We start from resolvents that depend on two variables, the complex energy z and a perturbation parameter. The wave functions and the energies of the resonances are expanded in powers of. It is shown that the zero-order terms correspond to the standard CAP method and that higher-order corrections are significant. The introduction of a convergence operator allows to control the convergence of the perturbation series. Due to the discretization of the continuum, the series are generally asymptotic. Finally, we relate the perturbation series to numerically convenient Taylor series. The theory is illustrated on two model examples.
Proceedings of Symposia in Pure Mathematics, 2007
Quantum theory makes a sharp distinction between bound states and scattering states, the former associated with point spectrum and the latter with continuous spectrum. Resonances associated with quasi-stationary states bridge this distinction, and have posed mathematical challenges since the beginning of the Schrödinger theory. Here the development of the mathematical underpinnings of resonance theory in atomic physics is reviewed, with particular reference to the rôle of the (DC) Stark effect, and time-independent perturbations of bound states in the two-body problem in atomic and molecular physics.
American Journal of Physics, 2020
The second-order perturbative Stark effect on the ground state of hydrogen is a typical example presented in many standard texts on quantum mechanics. Some texts miss the fact that the scattering states are significant contributors to the perturbative energy correction. The inclusion of scattering states has wider applicability than to just the Stark effect. An explicit calculation involving a finite-square well with a perturbation is used to illustrate the importance of including scattering states into the calculation. The second-order correction to the ground-state energy is obtained in three distinct ways. The first involves altering the problem by imposing additional boundary conditions at large distances to make the positive-energy spectrum discrete. The second makes use of the continuum scattering states directly. The third bypasses the use of scattering states by solving a differential equation for the first-order correction to the wave function.
International Journal of Molecular Sciences, 2002
We present in this paper two new versions of Rayleigh-Schrödinger (RS) and the Brillouin-Wigner (BW) state-specific multi-reference perturbative theories (SS-MRPT) which stem from our state-specific multi-reference coupled-cluster formalism (SS-MRCC), developed with a complete active space (CAS). They are manifestly sizeextensive and are designed to avoid intruders. The combining coefficients c µ for the model functions φ µ are completely relaxed and are obtained by diagonalizing an effective operator in the model space, one root of which is the target eigenvalue of interest. By invoking suitable partitioning of the hamiltonian, very convenient perturbative versions of the formalism in both the RS and the BW forms are developed for the second order energy. The unperturbed hamiltonians for these theories can be chosen to be of both Mφller-Plesset (MP) and Epstein-Nesbet (EN) type. However, we choose the corresponding Fock operator f µ for each model function φ µ , whose diagonal elements are used to define the unperturbed hamiltonian in the MP partition. In the EN partition, we additionally include all the diagonal direct and exchange ladders. Our SS-MRPT thus utilizes a multi-partitioning strategy. Illustrative numerical applications are presented for potential energy surfaces (PES) of the ground ( 1 Σ + ) and the first delta ( 1 ∆) states of CH + which possess pronounced multi-reference character. Comparison of the results with the corresponding full CI values indicates the efficacy of our formalisms.
Chemical Physics, 2009
The retaining the excitation degree (RE) partitioning [R.F. Fink, Chem. Phys. Lett. 428 (2006) 461(20 September)] is reformulated and applied to multi-reference cases with complete active space (CAS) reference wave functions. The generalised van Vleck perturbation theory is employed to set up the perturbation equations. It is demonstrated that this leads to a consistent and well defined theory which fulfils all important criteria of a generally applicable ab initio method: The theory is proven numerically and analytically to be size-consistent and invariant with respect to unitary orbital transformations within the inactive, active and virtual orbital spaces. In contrast to most previously proposed multi-reference perturbation theories the necessary condition for a proper perturbation theory to fulfil the zeroth order perturbation equation is exactly satisfied with the RE partitioning itself without additional projectors on configurational spaces. The theory is applied to several excited states of the benchmark systems CH 2 , SiH 2 , and NH 2 , as well as to the lowest states of the carbon, nitrogen and oxygen atoms. In all cases comparisons are made with full configuration interaction results. The multi-reference (MR)-RE method is shown to provide very rapidly converging perturbation series. Energy differences between states of similar configurations converge even faster.
2006
A formalism is presented for the calculation of relativistic corrections to molecular electronic energies and properties. After a discussion of the Dirac and Breit equations and their first-order Foldy-Wouthuysen ͓Phys. Rev. 78, 29 ͑1950͔͒ transformation, we construct a second-quantization electronic Hamiltonian, valid for all values of the fine-structure constant ␣. The resulting ␣-dependent Hamiltonian is then used to set up a perturbation theory in orders of ␣ 2 , using the general framework of time-independent response theory, in the same manner as for geometrical and magnetic perturbations. Explicit expressions are given to second order in ␣ 2 for the Hartree-Fock model. However, since all relativistic considerations are contained in the ␣-dependent Hamiltonian operator rather than in the wave function, the same approach may be used for other wave-function models, following the general procedure of response theory. In particular, by constructing a variational Lagrangian using the ␣-dependent electronic Hamiltonian, relativistic corrections can be calculated for nonvariational methods as well.
Physical Review A, 1978
DiAiculties in the straightfoward application of the complex-coordinate method to the calculation of resonance states in many-electron systems are examined. For the case of shape resonances, it is shown that many of these difficulties can be avoided by using complex coordinates only after reduction of the system to an effective one-electron problem. Further simplifications are achieved by the use of an inner-projection technique to facilitate the computation of the complex Hamiltonian matrix elements. The method is first illustrated by application to a model-potential problem. Its suitability for studying many-electron problems is demonstrated by calculation of the position and width of a low-energy P' shape resonance in Be. We discuss the modifications necessary to study core-excited (Feshbach) resonances.
Journal of Physics B: Atomic, Molecular and Optical Physics, 2009
In molecular reactions at the microscopic level the appearance of resonances has an important influence on the reactivity. To study when a bound state transitions into a resonance and how these transitions depend on various system parameters such as internuclear distances, one needs to look at the poles of the S-matrix. Using numerical continuation methods and bifurcation theory, we develop efficient and robust methods to trace the parameter dependence of the poles of the S-matrix. Using pseudo-arclength continuation, we can minimize the numerical complexity of our algorithm. As a proof-of-concept we have applied our methods on a number of model problems.
Journal of Computational Chemistry, 2003
The second-order multireference perturbation theory using an optimized partitioning, denoted as MROPT(2), is applied to calculations of various molecular properties-excitation energies, spectroscopic parameters, and potential energy curves-for five molecules: ethylene, butadiene, benzene, N 2 , and O 2 . The calculated results are compared with those obtained with second-and third-order multireference perturbation theory using the traditional partitioning techniques. We also give results from computations using the multireference configuration interaction (MRCI) method. The presented results show very close resemblance between the new method and MRCI with renormalized Davidson correction. The accuracy of the new method is good and is comparable to that of second-order multireference perturbation theory using Møller-Plesset partitioning.
Physical Review Letters, 1970
PHYSI CAL REVIEW LETTERS 27 JULY 1970 electrons ordinarily exist only up to microwave frequencies. However, Eq. (8) indicates that electrons modulated by a laser will give rise to fields at discrete frequencies far beyond T It seems plausible that the optical radiation observed by Schwarz and Hora' is due to the presence of these classical fields. The justification of the use of the classical calculation at frequencies»T ' must await a full quantum mechanical treatment of the problem. Further investigation is needed to determine the details of the modulation process, the extent of the modulation, the details of the radiation mechanism, and the resulting intensities.
Physical Review B, 2005
Using the Green's function formalism, an ab initio theory for band structures of crystals is derived starting from the Hartree-Fock approximation. It is based on the algebraic diagrammatic construction scheme for the self-energy which is formulated for crystal orbitals (CO-ADC). In this approach, the poles of the Green's function are determined by solving a suitable Hermitian eigenvalue problem. The method is not only applicable to the outer valence and conduction bands, it is also stable for inner valence bands where strong electron correlations are effective. The key to the proposed scheme is to evaluate the self-energy in terms of Wannier orbitals before transforming it to a crystal momentum representation. Exploiting the fact that electron correlations are mainly local, one can truncate the lattice summations by an appropriate configuration selection scheme. This yields a flat configuration space; i.e., its size scales only linearly with the number of atoms per unit cell for large systems and, under certain conditions, the computational effort to determine band structures also scales linearly. As a first application of the new formalism, a lithium fluoride crystal has been chosen. A minimal basis set description is studied, and a satisfactory agreement with previous theoretical and experimental results for the fundamental band gap and the width of the F 2p valence band complex is obtained.
International Journal of Quantum Chemistry, 1978
Chemical Physics Letters, 1999
Finite-order perturbation corrections are ambiguous since they depend on the partitioning of the Hamiltonian to a zero-order part and perturbation, and any chosen partitioning can be freely modified, e.g. by level shift projectors. To optimize low-order corrections, an approximate variational procedure is proposed to determine level shift parameters from the first-order Ansatz for the wavefunction. The resulting new partitioning scheme provides significantly better second-order results than those obtained by standard partitions like Epstein-Nesbet or Møller-Plesset. We treat the anharmonic oscillator and the atomic electron correlation energy in He, Be and Ne as numerical test cases.
Using a formalism set up by Bloch an effective Hamiltonian is introduced which, when expanded in powers of a perturbation V, generalizes the Rayleigh-Schrödinger perturbation theory. This effective Hamiltonian provides a simple way of handling degeneracies and the orthonormalization of the perturbed wavefunctions. As an illustration the spin Hamiltonian introduced by Pryce in paramagnetic resonance is derived up to third order terms, and the electron-lattice interaction is discussed showing that the Born - Oppenheimer approximation follows from first-order terms in the case of non-degenerate electronic ground state.
The Journal of Chemical Physics, 2008
A nonlinear conjugate gradient optimization scheme is used to obtain excitation energies within the random phase approximation ͑RPA͒. The solutions to the RPA eigenvalue equation are located through a variational characterization using a modified Thouless functional, which is based upon an asymmetric Rayleigh quotient, in an orthogonalized atomic orbital representation. In this way, the computational bottleneck of calculating molecular orbitals is avoided. The variational space is reduced to the physically-relevant transitions by projections. The feasibility of an RPA implementation scaling linearly with system size N is investigated by monitoring convergence behavior with respect to the quality of initial guess and sensitivity to noise under thresholding, both for well-and ill-conditioned problems. The molecular-orbital-free algorithm is found to be robust and computationally efficient, providing a first step toward large-scale, reduced complexity calculations of time-dependent optical properties and linear response. The algorithm is extensible to other forms of time-dependent perturbation theory including, but not limited to, time-dependent density functional theory.
Journal of Chemical Education, 1998
1971
The ground-state energy of a hydrogen atom is treated, according to Wigner, by the RayleighSchrodinger perturbation method, where unperturbed wave functions correspond to a free particle in a spherical "box" of radius R. We have calculated the perturbed energy up to fifth order in e and as the sign of each term from the second approximation on remains negative, we may conjecture that the perturbation series is divergent for R co. It also appears that in this way it will be difficult to obtain a Stieltjes series in the coupling constant e'. We show that it is possible, by using the Rayleigh-Schrodinger-Lbdin variational energies, to obtain good values for the ground-state energy by adjusting conveniently the value of R. The Padé approximants give bad values for R -+ co.
2015
Resonances are temporarily bound states which lie in the continuum part of the Hamiltonian. If the electronic coordinates of the Hamiltonian are scaled ("dilated") by a complex parameter, = e i (, real), then its complex eigenvalues represent the scattering states (resonant and non-resonant) while the eigenvalues corresponding to the bound states and the ionization and the excitation thresholds remain real and unmodified. These make the study of these transient species amenable to the bound state methods. We developed a quadratically convergent multiconfigurational self-consistent field method (MCSCF), a well-established bound-state technique, combined with a dilated Hamiltonian to investigate resonances. This is made possible by the adoption of a second quantization algebra suitable for a set of "complex conjugate biorthonormal" spin orbitals and a modified step-length constraining algorithm to control the walk on the complex energy hypersurface while searching for the stationary point using a multidimensional Newton-Raphson scheme. We present our computational results for the 2 P Be¯ shape resonances using two different computationally efficient methods that utilize complex scaled MCSCF (i.e., CMCSCF). These two methods are to straightforwardly use CMCSCF energy differences and to obtain energy differences using an approximation to the complex multiconfigurational electron propagator. It is found that, differing from previous computational studies by others, there are actually two 2 P Be¯ shape resonances very close in energy. In addition, N 2 resonances are examined using one of these methods.
1997
Starting from previous work, where Rayleigh-Schrödinger perturbation theory has been reformulated in matrix form, a practical algorithm implementation is described using both full matrix and vector alternatives. An heuristic convergence sufficiency criterion based on Gershgorin discs is also presented. Some numerical examples related to atomic CI computations are reported to illustrate the theoretical framework.
The Journal of Chemical Physics, 2006
The introduction of complex absorbing potentials as numerical tools to stabilize or increase the efficiency of calculations based on wave-packet propagation or on eigenvalue problems has the drawback of causing a modification of the Hamilton operator of the problem. In this work the consequences of such a modification are analyzed and the corrections required in order to properly describe the original physical process are derived. As an example, the decay of excited molecular states is considered: it is shown that the standard time-independent expression for the decay spectrum loses its validity when a complex absorbing potential is introduced in the nuclear Hamilton operator of the problem. To remedy the situation, a new, very stable formula is derived and tested on relevant model studies. Numerical examples are discussed.
Loading Preview
Sorry, preview is currently unavailable. You can download the paper by clicking the button above.