International Journal of Pure and Applied Sciences and Technology

Similar documents
Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Linear Approximation with Regularization and Moving Least Squares

Lecture 12: Discrete Laplacian

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

where the sums are over the partcle labels. In general H = p2 2m + V s(r ) V j = V nt (jr, r j j) (5) where V s s the sngle-partcle potental and V nt

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Difference Equations

Inductance Calculation for Conductors of Arbitrary Shape

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

A Robust Method for Calculating the Correlation Coefficient

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

14 The Postulates of Quantum mechanics

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

High resolution entropy stable scheme for shallow water equations

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

A new Approach for Solving Linear Ordinary Differential Equations

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Note 10. Modeling and Simulation of Dynamic Systems

Chapter 1. Probability

4DVAR, according to the name, is a four-dimensional variational method.

Introduction to Density Functional Theory. Jeremie Zaffran 2 nd year-msc. (Nanochemistry)

(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate

Global Sensitivity. Tuesday 20 th February, 2018

Numerical Heat and Mass Transfer

Appendix B. The Finite Difference Scheme

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

The Feynman path integral

Chat eld, C. and A.J.Collins, Introduction to multivariate analysis. Chapman & Hall, 1980

One-sided finite-difference approximations suitable for use with Richardson extrapolation

Grover s Algorithm + Quantum Zeno Effect + Vaidman

The Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Time-Varying Systems and Computations Lecture 6

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

A Hybrid Variational Iteration Method for Blasius Equation

The non-negativity of probabilities and the collapse of state

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

LECTURE 9 CANONICAL CORRELATION ANALYSIS

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

Probabilistic method to determine electron correlation energy

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Markov Chain Monte Carlo (MCMC), Gibbs Sampling, Metropolis Algorithms, and Simulated Annealing Bioinformatics Course Supplement

NUMERICAL DIFFERENTIATION

Entropy generation in a chemical reaction

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

THEOREMS OF QUANTUM MECHANICS

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

Identification of Linear Partial Difference Equations with Constant Coefficients

12. The Hamilton-Jacobi Equation Michael Fowler

SL n (F ) Equals its Own Derived Group

CHAPTER 14 GENERAL PERTURBATION THEORY

New Method for Solving Poisson Equation. on Irregular Domains

Linear Regression Analysis: Terminology and Notation

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Hongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)

This column is a continuation of our previous column

On an Extension of Stochastic Approximation EM Algorithm for Incomplete Data Problems. Vahid Tadayon 1

A Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods

6.3.4 Modified Euler s method of integration

The Geometry of Logit and Probit

DUE: WEDS FEB 21ST 2018

PARTICIPATION FACTOR IN MODAL ANALYSIS OF POWER SYSTEMS STABILITY

Supplemental Material: Causal Entropic Forces

De-noising Method Based on Kernel Adaptive Filtering for Telemetry Vibration Signal of the Vehicle Test Kejun ZENG

Color Rendering Uncertainty

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl

Tensor Smooth Length for SPH Modelling of High Speed Impact

A MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS

ESCI 341 Atmospheric Thermodynamics Lesson 10 The Physical Meaning of Entropy

Lecture 16 Statistical Analysis in Biomaterials Research (Part II)

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

An identification algorithm of model kinetic parameters of the interfacial layer growth in fiber composites

Parameter Estimation for Dynamic System using Unscented Kalman filter

AN IMPROVED PARTICLE FILTER ALGORITHM BASED ON NEURAL NETWORK FOR TARGET TRACKING

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

Lecture Notes on Linear Regression

DETERMINATION OF UNCERTAINTY ASSOCIATED WITH QUANTIZATION ERRORS USING THE BAYESIAN APPROACH

Electron-Impact Double Ionization of the H 2

Chapter 13: Multiple Regression

THE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Lecture 14: Forces and Stresses

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

Advanced Quantum Mechanics

Density matrix. c α (t)φ α (q)

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Uncertainty and auto-correlation in. Measurement

Effects of Ignoring Correlations When Computing Sample Chi-Square. John W. Fowler February 26, 2012

Lecture 21: Numerical methods for pricing American type derivatives

Composite Hypotheses testing

The Study of Teaching-learning-based Optimization Algorithm

Formulas for the Determinant

The Order Relation and Trace Inequalities for. Hermitian Operators

DFT with Planewaves pseudopotential accuracy (LDA, PBE) Fast time to solution 1 step in minutes (not hours!!!) to be useful for MD

An Improved multiple fractal algorithm

Transcription:

Int. J. Pure Appl. Sc. Technol., 4() (03), pp. 5-30 Internatonal Journal of Pure and Appled Scences and Technology ISSN 9-607 Avalable onlne at www.jopaasat.n Research Paper Schrödnger State Space Matrx Parametres Estmaton J. J. Medel, * and R. Palma Computng Research Centre, Real-tme and Automatc Control Department. Computng School, Natonal Polytechnc Insttute, Mexco Cty, Mexco * Correspondng author, e-mal: (jjmedelj@yahoo.com.mx) (Receved: 3-9-; Accepted: -0-) Abstract: The Schrödnger equaton as a quantum superposton descrbes the atomc orbtal as a lnear dfference functons bass set. Another alternatve s that ths equaton could be wrtten as a m-dmensonal state space lnear stochastc system wth unnown matrx coeffcents wth smooth movements around the statonary regon. Thus, the matrx s estmated n a probablty sense consderng the dstrbuton functon and ts stochastc moments. The matrx estmaton result s an alternatve soluton wthout losng the general system propertes. The estmaton s based on a gradent stochastc technque wth an nstrumental varable such that the contrbuton matrx s optmal n a probablty sense. Therefore, the advantages wth respect to tradtonal solutons are focused on estmatng the matrx contrbuton on-lne wth a lnear complexty. Keywords: Schrödnger equaton, Mathematcal methods, Fnte dfference methods, Flters, Control theory.. Introducton The Schrödnger equaton H E, s governed by physcs, chemstry and bology prncples and has an approxmate predctve power [9]. However, for over 80 years snce ts ncepton; t was thought to be nsoluble except for a few specal cases. Snce 000, Naatsuj et al., n [9] formulated a general method solvng the Schrödnger equaton appled to atoms and molecules n analytcal expanson form, havng a hgh numercal and computatonal complexty [5]. In a physcal sense, the Schrödnger equaton as a quantum superposton descrbes the atomc orbtal as a bass functon set. Ths s expanded as a lnear functonal combnaton. It s descrbed as a matrces array n fnte dfferences wth unnown coeffcents [6], [7], [8].

Int. J. Pure Appl. Sc. Technol., 4() (03), 5-30 6 An ntal assumpton s that the orbtal number s ncluded n the lnear expanson. In a sense, m- atomc orbtal can be numbered from to m, wth m Z. The expresson for the -th orbtal s gven n (). c χ () r Where s an orbtal represented as the sum of m-atomc orbtals and ts correspondng contrbuton coeffcent c r, and r,m, m Z, where the Hartree-Foc procedure s used obtanng the expanson coeffcents [].. Prelmnares Quantum calculatons are typcally performed wthn a fnte bass functons set. In these cases, the wave functons under consderaton are represented as vectors. The components correspond to coeffcents n the bass set lnear functons combnaton. The operatons are represented as matrces n fnte dfferences. A mnmum bass set s one n whch each atom s represented by a basc Hartree-Foc calculaton functon [4, 7]. The orbtal s expressed as a lnear electron functons combnaton bass, centred on atomc components. The basc atomc orbtal s the hydrogen atom and s nown as Slater-type orbtal [4]. An atomc orbtal s a mathematcal functon that descrbes the wave-le behavour of ether one or a par of electrons n an atom [7]. Atomc orbtal s typcally descrbed as hydrogen scheme (meanng one-electron) wave functons over space. The wave functon for a complex atom electrons cloud may be seen as beng bult up n an electron confguraton that s a sum of smpler hydrogen atomc orbtals. The orbtal approxmaton s symbolcally represented for the atomc wave functon, and s descrbed n ().... () r r 3 Where each, functon wth, m, m Z, s defned as 3-dmensonal ndvdual electron property. The spn s ncluded n the wave functon as a spn-orbtal product. The wave functon has to be antsymmetrc wth respect to electron nterchange wrtten as a dagonal matrx [7]. The goal of most quantum calculatons s the producton of a wave functon. Ths s acheved f all or are nown. consttuent orbtals are descrbed by { } 3. Results and Dscusson, m Quantum super-poston s a fundamental prncple of quantum mechancs. Mathematcally, t refers to Schrödnger equaton solutons property havng a lnear form n state space of a nown m-th-partcle quantum object. Then, the vector contans the poston coordnates assocated wth quantum objects. The stochastc gradent of a wave functon s symbolzed by. Thus, the state space s defned n (3). (,, ) (,, ) W ( ) (3),, Where (,, ) s the Schrödnger equaton, and W s a stochastc system n Hlbert space wth L norm,,, ) s descrbed n a probablty space. ( Wth states as x, x, then x ɺ x, xɺ, where x y z 4π 4π and 4π 4π 0 wth c.

Int. J. Pure Appl. Sc. Technol., 4() (03), 5-30 7 Therefore c consdered the wor Palma et al., n whch the propertes of the m-dmensonal stochastc estmator for a blac-box model were descrbed n [0], the Schrödnger equaton state space n a matrcal form s gven n (4). xɺ xɺ 0 c x 0 x dw dw (4) A lnear relatonshp for a mult-electron system s bult as follows. Let the Schrödnger equaton as (4) n an affne state space gven by (5) t was necessary to change the rows n (4). Where m x, m and[ dw, ] T [ w, ] T are the nnovaton stochastc processes for (4) and (5), respectvely. x dw w mɺ mɺ c 0 0 m m w w (5) From (), each molecular orbtal,, mɺ mɺ n, n Z,, has an assocated affne state space (6). c 0 0 m m w w (6) And the total state space s gven by (7). Mɺ M ɺ C 0 0 I M M W W (7) m Where C, I R are constant vectors, defned by C : { c }, I : {, m m } and, m C 0 M W A,, W. The lnear relatonshp for mult-electron model s (8). 0 I M W ɺ A W (8) After the Schrödnger state space model s formed, ts goal s to estmate the unnown parametres denoted by A ~. Frst, transformng (8) n dfferences havng as a stochastc sample set taen m Z sze from a random vector. Second, the matrx parametres correspondence has the form I A. Thrd, the whte nose has an approxmated dscrete descrpton. Then (8) n A dfferences has the form (9). A The matrx A s ntegrated by the contrbuton weghts set. It s a matrx wth unnown parametres []. Although the Schrödnger equaton was solved by [5, 9], numercal ntegraton methods are usually nsuffcent for accuracy and has a hgh computatonal cost, havng an exponental growth wth respect to the orbtals numbers consdered [9]. In spte of, the state space contrbuton of (9) s numercally reduced wth a lnear computatonal complexty observng n ths model the stablty condtons [3]. Unfortunately, the new system has unnown contrbuton weghts reflected n the matrx A. W W In a probablty sense, t has an optmal descrpton based on the nnovaton process W (9), data nformaton, was developed usng the stochastc gradent, symbolcally descrbed n (0) [6, 8]. ~ A E{ M } E{ W M } (0) Wth W, M are the process nose and the correlaton matrx, respectvely; and the nnovaton () s bult from (9). W A ()

Int. J. Pure Appl. Sc. Technol., 4() (03), 5-30 8 Then () has a dagonal form (). M W 0 C 0 M 0 () 0 M W 0 I 0 M T M 0 ϑ be an nstrumental varable. Multplyng () by ϑ 0 M Let, wth second probablty moment havng the form (3). E{( ) } ~ W A E{ M } ϑ (3) M, ϑ det( M ) 0, andm ϑ M. Fnally, the optmal stochastc estmator s gven by (0). Schrödnger equaton numercally smulated n state space s shown n Fgures and. Fgure.a) depctng the C components, Fgure.b) depctng the I component. Fgures.c) and.d) show the thrd and fourth components. M l 0.5 M.0 0.3 0.5 0. 5 0 5 0 5 0 5 0 a) b) M 3 0.5 M 4 0. 0. 5 0 5 0 5 0 5 0 0. 0. c) d) Fgure. Numercal smulaton: a) Frst, C5 component, b) Second, I component, c) Thrd, zero component, d) Fourth, zero component.

Int. J. Pure Appl. Sc. Technol., 4() (03), 5-30 9 AA 0.0004 0.000 5 0 5 0 0.000 0.0004 a) b) Fgure. Numercal results: a) Matrx parametres: Frst n red, second n green, thrd n blue, fourth n brown, b) Functonal error and ts nset. Fgure.a) shows matrx parametres numercal evoluton and, Fgure.b) shows the functonal error wth an accuracy order of 0-0. 4. Conclusons Due to Schrödnger beng transformed n state space and developed n dagonal form (), the optmal stochastc estmator for a mult-electron has the form (3). Therefore, the numercal complexty s lneal [0]. Fgure 3.a) shows the tradtonal computatonal complexty soluton method, Fgure 3.b) descrbes a lnear state space descrpton wth lnear computatonal complexty. Fgure 3. Computatonal complexty: a) Tradtonal Schrödnger soluton, b) State Space Schrödnger descrpton. The matrx parametres C s estmated usng the dagonalzaton process for m orbtal wthout calculatng the tradtonal potental terms. For the case m havng lowest tme executon as shown n Fgure 3. Ths has mportant mplcatons for data storage and computng when evaluatng very large complex systems. Acnowledgements Ths document was checed by Englsh Teacher Lawrence Whtehll, England.

Int. J. Pure Appl. Sc. Technol., 4() (03), 5-30 30 References [] C.J. Cramer, Essentals of Computatonal Chemstry, John Wley and Sons Ltd., Chchester, England, 00. [] U.K. Dewangan, Structural parameter dentfcaton A new concept, Int. J. Pure Appl. Sc. Technol., 0() (0), -9. [3] B. Ghanbar, M.G. Porshoouh and B. Rahm, A numercal method for solvng systems of nonlnear ODE s, Int. J. Pure Appl. Sc. Technol., 3() (0), 7-34. [4] F. Jensen, Introducton to Computatonal Chemstry, John Wley and Sons Ltd., Chchester, England, 999. [5] Y. Kuroawa, H. Naashma and H. Naatsuj, Free teratve-complement-nteracton calculatons of the hydrogen molecule, Phys. Rev., A7(005), 0650. [6] E.L. Lehmann and G. Casella, Theory of Pont Estmaton, Sprnger-Verlag, New Yor, 998. [7] I.N. Levne, Quantum Chemstry, Prentce Hall, New Jersey, 99. [8] J.J. Mortensen, K. Kaasbjerg, S.L. Fredersen, J.K. Norsov, J.P. Sethna and K.W. Jacobsen, Bayesan error estmaton n densty-functonal theory, Phys. Rev. Lett., 95(005), 640. [9] H. Naatsuj, H. Naashma, Y. Kuroawa and A. Ishawa, Solvng the Schrödnger equaton of atoms and molecules wthout analytcal ntegraton based on the free teratve-complementnteracton wave functon, Phys. Rev. Lett., 99(007), 4040. [0] R.O. Palma, J.J. Medel and G.A. Garrdo, An m-dmensonal stochastc estmator, Rev. Mex. Fís., 58() (0), 0069.