Evaluation of Gaussian Molecular Integrals
|
|
- Hannah Warner
- 5 years ago
- Views:
Transcription
1 The Mathematica Journal Evaluation of Gaussian Molecular Integrals II. Kinetic-Energy Integrals Minhhuy Hô Julio Manuel Hernández-Pérez This article carries out the evaluation of kinetic energy integrals using Gaussian-type functions with arbitrary Cartesian angular values. As an example, we calculate the kinetic matrix for the water molecule in the STO-3G basis set. Introduction In this article, the second of a series describing algorithms for evaluating molecular integrals, we detail the evaluation of the kinetic energy integrals. Detailed accounts of molecular integrals can be found in the references of [1]. The electronic kinetic energy in atomic units (m e = e = = 1 ê 4 p e 0 = 1) involves integrals of the type K = - c a Hr; a, A, al c b Hr; b, B, bl dr, (1) in which c a HrL is an unnormalized Cartesian Gaussian primitive centered at the nucleus A = 9A x, A y, A z =: c a Hr; a, A, al = Hx - A x L a x Iy - A y M a y Hz - A z L a z e -a r-a 2 = Hx - A x L a x e -a x-a x 2 Iy - A y M a y e -a y-a y 2 Hz - A z L a z e -a z-a z 2 = c a Hx; a, A x, a x L ÿ c a Iy; a, A y, a y M ÿ c a Hz; a, A z, a z L. Here a is the orbital exponent and the polynomial represents the angular part, in that the sum of the Cartesian angular momenta a x + a y + a z = 0, 1, 2, 3, corresponds to functions of type s, p, d, f,. One notable property of the Gaussian function, which will be used here, is that the derivative of a Gaussian function can be expressed as a sum containing Gaussians of lower and higher Cartesian angular values. In particular, x c ahx; a, A x, a x L = a x c a Hx; a, A x, a x - 1L - 2 a c a Hx; a, A x, a x + 1L. (2) (3)
2 2 Minhhuy Hô and Julio Manuel Hernández-Pérez Similarly, A x c a Hx; a, A x, a x L = - x c ahx; a, A x, a x L. This property is useful in integrals involving differential operators, such as the kinetic energy; in calculations involving the gradient of the energy or Hamiltonian; and in deriving one of the most important algorithms of Gaussian function computation, the recurrence relation. The kinetic operator, expressed in Cartesian coordinates as = J 2 x y z 2 N, enables separation of the kinetic integrals into three components K = K x + K y + K z. For example, K x = c a Hr; a, A, al 2 x 2 c bhr; b, B, bl dr. Furthermore, if we also separate c a into its Cartesian components, K x becomes a product of three integrals K x = c a Hx; a, A x, a x L - (4) 2 x c bhx; b, B x, b x L dx c a Iy; a, A y, a y M ÿ c b Iy; b, B y, b y M 2 - (5) dy - c a Hz; a, A z, a z L ÿ c b Hz; b, B z, b z L dz. Derivative of Gaussian Functions The last two integrals in K x are simply overlap integrals and can be evaluated as outlined in [1]. Integrating by parts, the first integral is Bc ahx; a, A x, a x L - - x= x c bhx; b, B x, b x L x=- x c ahx; a, A x, a x L x c bhx, b, B x, b x L dxf = x c ahx; a, A x, a x L x c bhx, b, B x, b x L dx where, according to the property mentioned above, the first term can be expanded into two Gaussian products, which vanish at the integral limits. (6)
3 Evaluation of Gaussian Molecular Integrals 3 Next, we consider the first term of the integrand. Expressing equation (3) in terms of (1) and rearranging, we get x Hx - A xl a x e -ahx-a xl2 = a x Hx - A x L a x-1 e -ahx-a xl 2 ahx - A x LD Hx - A x L a x e -ahx-a xl 2 = A a x Hx - A x L a x-1-2 ahx - A xl a x+1 E e -ahx-a xl 2. Substituting the result into equation (6), we have 1 2 Aa x Hx - A x L ax-1-2 ahx - A xl ax+1 E - e -ahx-a xl 2 Ab x Hx - B x L b x-1-2 bhx - B xl b x+1 E e -bhx-b xl 2 dx, and expanding the polynomial part yields 1 2 Aa x b x Hx - A x L ax-1 Hx - B x L bx a b x Hx - A xl a x+1 Hx - B x L b x-1-2 a x bhx - A x L a x-1 Hx - B x L b x a bhx - A x L a x+1 Hx - B x L b x+1 E e -ahx-a xl 2 e -bhx-b xl 2 dx. We can now express the first integral of equation (5) in terms of one-dimensional Gaussian primitives: 1 2 Ba x b x - c a Hx; a, A x, a x - 1L ÿ c b Hx; b, B x, b x - 1L dx - (7) (8) (9) 2 a b x - c a Hx; a, A x, a x + 1L ÿ c b Hx; b, B x, b x - 1L dx - 2 a x b - c a Hx; a, A x, a x - 1L ÿ c b Hx; b, B x, b x + 1L dx + (10) 4 a b - c a Hx; a, A x, a x + 1L ÿ c b Hx; b, B x, b x + 1L dxf. Kinetic-Energy Integral We have now shown that the kinetic-energy integral can be written as a product of overlap integrals. The x component, for example, is given by K x = 1 2 Ba x b x - c a Hx; a, A x, a x - 1L ÿ c b Hx; b, B x, b x - 1L dx - 2 a b x - c a Hx; a, A x, a x + 1L ÿ c b Hx; b, B x, b x - 1L dx - 2 a x b - c a Hx; a, A x, a x - 1L ÿ c b Hx; b, B x, b x + 1L dx + (11)
4 4 Minhhuy Hô and Julio Manuel Hernández-Pérez 4 a b - c a Hx; a, A x, a x + 1L ÿ c b Hx; b, B x, b x + 1L dxf - c a Iy; a, A y, a y M ÿ c b Iy; b, B y, b y M dy - c a Hz; a, A z, a z L ÿ c b Hz; b, B z, b z L dz. Recurrence Relations Using the notation of [1], where the overlap integral of two Gaussians is expressed in terms of its orthogonal components, S = c a Hx; a, A x, a x L ÿ c b Hx; b, B x, b x L dx = - p E AB B F3ê2 sx Ha x, b x L s y Ia y, b y M s z Ha z, b z L, the x component of the kinetic-energy integral, equation (11), can be written as K x = 1 2 E p ABB F3ê2 8ax b x s x Ha x - 1, b x - 1L - 2 a b x s x Ha x + 1, b x - 1L - 2 a x bs x Ha x - 1, b x + 1L + 4 a b s x Ha x + 1, b x + 1L< s y Ia y, b y M p s z Ha z, b z L = E AB B F3ê2 kx Ha x, b x L s y Ia y, b y M s z Ha z, b z L, (12) (13) where the factor 1 2 has been absorbed into the definition of the kinetic integral k xha x, b x L. We will derive the recurrence relation for this term k x Ha x, b x L. For a x = b x = 0, the first three terms of equation (13) are zero, and we are left with k x H0, 0L = 2 a b s x H1, 1L. Analogously, for a x = 0 or b x = 0, we have k x Ha x, 0L = -a x b s x Ha x - 1, 1L + 2 a b s x Ha x + 1, 1L, k x H0, b x L = -a b x s x H1, b x - 1L + 2 a b s x H1, b x + 1L, and the recurrence relation for the kinetic integral function with any two Cartesian angular momenta generalizes to (14) k x Ha x, b x L = 1 x b x s x Ha x - 1, b x - 1L - 2 a b x s x Ha x + 1, b x - 1L - 2 a x b s x Ha x - 1, b x + 1L + 4 a b s x Ha x + 1, b x + 1LD. (15)
5 Evaluation of Gaussian Molecular Integrals 5 Implementation The function Kin evaluates the kinetic integral of two Gaussian primitives; here alpha, beta, RA, RB, LA, and LB are a, b, A, B, a, and b as defined earlier. The first step is the evaluation of the overlap integral as described in [1]. Kin@alpha_, beta_, RA_, RB_, LA_, LB_D := Module@8EAB, Kinetic<, Do@ H* Initial Conditions for overlap recurrence relation *L s@i, 0, 0D = 1; s@i, 1, 0D = -HRA@@iDD - HHalpha * RA@@iDD + beta * RB@@iDDL ê HalphetaLLL; H* Recurrence Index *L s@i_, a_, 0D := -HRA@@iDD - Halpha * RA@@iDD + beta * RB@@iDDL ê HalphetaLL * s@i, a - 1, 0D + HHa - 1L ê H2 * HalphetaLLL * s@i, a - 2, 0D; H* Transfer Equation *L s@i_, a_, b_d := s@i, a + 1, b - 1D + HRA@@iDD - RB@@iDDL * s@i, a, b - 1D; H* Initial Conditions for kinetic integral function *L k@i, 0, 0D = 2 * alpha * beta * s@i, 1, 1D; k@i_, a_, 0D := -a * beta * s@i, a - 1, 1D + 2 * alpha * beta * s@i, a + 1, 1D; k@i_, 0, b_d := -alpha * b * s@i, 1, b - 1D + 2 * alpha * beta * s@i, 1, b + 1D; H* Kinetic Energy Integral *L k@i_, a_, b_d := Ha * b * s@i, a - 1, b - 1D - 2 * a * beta * s@i, a - 1, b + 1D - 2 * alpha * b * s@i, a + 1, b - 1D + 4 * alpha * beta * s@i, a + 1, b + 1DL ê 2;, 8i, 1, 3<D; H* Loop over three Cartesian components *L H* Exponential Overlap Factor *L
6 6 Minhhuy Hô and Julio Manuel Hernández-Pérez EAB = Exp@-Halpha * beta ê HalphetaLL * HRA - RBL.HRA - RBLD; H* Kinetic Energy Integral *L D Kinetic = EAB * HHPi ê HalphetaLL^H3 ê 2LL * Hk@1, LA@@1DD, LB@@1DDD * s@2, LA@@2DD, LB@@2DDD * s@3, LA@@3DD, LB@@3DDD + s@1, LA@@1DD, LB@@1DDD * k@2, LA@@2DD, LB@@2DDD * s@3, LA@@3DD, LB@@3DDD + s@1, LA@@1DD, LB@@1DDD * s@2, LA@@2DD, LB@@2DDD * k@3, LA@@3DD, LB@@3DDDL We describe in detail the evaluation of the kinetic-energy matrix for the water molecule (r OH = bohr, HOH = , the geometry optimized at the HF/STO-3G level). The molecule lies in the y-z plane with Cartesian coordinates in atomic units. R = 880., , <, 80., , <, 80., 0., <<; In the STO-3G basis set, each atomic orbital is approximated by a sum of three Gaussians; their unnormalized primitive contraction coefficients and orbital exponents (taken from [2]) are as follows. PrimCoeff = , , <, , , <, , , <, , , <, , , <, , , <, , , <<; OrbCoeff = , , <, , , <, , , <, , , <, , , <, , , <, , , <<;
7 Evaluation of Gaussian Molecular Integrals 7 Here are basis function origins and Cartesian angular values of the orbitals, listed in the order H 1 H1 sl, H 2 H1 sl, O 3 H1 sl, O 3 H2 sl, O 3 H2 p x L, O 3 I2 p y M, and O 3 H2 p z L. FCenter = 8R@@1DD, R@@2DD, R@@3DD, R@@3DD, R@@3DD, R@@3DD, R@@3DD<; CartAng = 880, 0, 0<, 80, 0, 0<, 80, 0, 0<, 80, 0, 0<, 81, 0, 0<, 80, 1, 0<, 80, 0, 1<<; Specifically, for the kinetic-energy integral of the first primitive of the 1 s orbital of hydrogen atom 1, c 11 HrL, and the first primitive of the 2 p z orbital of the oxygen atom, c 71 HrL, - c 11 HrL c 71 HrL dr = e r -R H Hz L µ (16) e r -R O3 2 dx dy dz, where the indices c ij indicate primitive j of basis function i. The kinetic integral in terms of the kinetic and overlap function is then p E AB B F3ê2 9k x H0, 0L s y H0, 0L s z H0, 1L + s x H0, 0L k y H0, 0L s z H0, 1L + s x H0, 0L s y H0, 0L k z H0, 1L=. We start with a = 80, 0, 0< and b = 80, 0, 1<, and from the first four equations of the module, we use the recurrence scheme to build up the overlap values needed for k z H0, 1L. We will need values up to s z H1, 2L, s z H0, 0L = 1, (17) s z H1, 0L = - H1 z - a H1 z + b O z s z H2, 0L = - H1 z - a H1 z + b O z s z H3, 0L = - H1 z - a H1 z + b O z = , s z H1, 0L HL s zh0, 0L = , s z H1, 1L = s z H2, 0L + HH1 z - O z L s z H1, 0L = , s z H2, 1L = s z H3, 0L + HH1 z - O z L s z H2, 0L = , s z H1, 2L = s z H2, 1L + HH1 z - O z L s z H1, 1L = , s z H2, 0L HL s zh1, 0L = , (18) leading to k z H0, 1L = -a ÿ 1 ÿ s z H1, 0L + 2 a ÿ b ÿ s z H1, 2L =
8 8 Minhhuy Hô and Julio Manuel Hernández-Pérez Similarly, for k x H0, 0L, we will need s x H1, 1L which, in turn, comes from s x H1, 0L and s x H2, 0L. We have s x H1, 0L = 0, s x H2, 0L = , s x H1, 1L = , leading to k x H0, 0L = 2 a b s x H1, 1L = For the y component, s y H1, 0L = , s y H2, 0L = , s y H1, 1L = , leading to k y H0, 0L = With E AB = , the kinetic integral of equation (16) is K = , which we can also obtain from the module Kin. Kin@OrbCoeff@@1, 1DD, OrbCoeff@@7, 1DD, FCenter@@1DD, FCenter@@7DD, CartAng@@1DD, CartAng@@7DDD Three Gaussian primitives c a HrL for each atomic orbital f i HrL result in nine integrals of the type that we have just evaluated. For example, the K 17 element of the kinetic-energy matrix, 3 3 K 17 = N 1 i N 7 j d 1 i d 7 j c 1 i HrL c 7 j HrL dr, i=1 j=1 (19) is Ÿ f1 - HrL J N f 7 HrL dr, derived via a contraction scheme that requires the following normalization factor. NormCoeff@a_, 8ax_, ay_, az_<d := 2 a p 3ê4 H4 al Hax+ay+azLê2 ; HH2 ax - 1L!! H2 ay - 1L!! H2 az - 1L!!L 1ê2
9 Evaluation of Gaussian Molecular Integrals 9 This calculates K 17. Sum@ NormCoeff@OrbCoeff@@1, idd, CartAng@@1DDD * NormCoeff@OrbCoeff@@7, jdd, CartAng@@7DDD * PrimCoeff@@1, idd * PrimCoeff@@7, jdd * Kin@OrbCoeff@@1, idd, OrbCoeff@@7, jdd, FCenter@@1DD, FCenter@@7DD, CartAng@@1DD, CartAng@@7DDD, 8i, 3<, 8j, 3<D For larger basis sets, one needs only to replace the summation upper limit 3 in equation (19) with the appropriate number of primitives belonging to a particular basis function and an additional summation for the basis functions. Finally, here is the resulting kinetic energy matrix. Table@Sum@ NormCoeff@OrbCoeff@@p, idd, CartAng@@pDDD * NormCoeff@OrbCoeff@@q, jdd, CartAng@@qDDD * PrimCoeff@@p, idd * PrimCoeff@@q, jdd * Kin@OrbCoeff@@p, idd, OrbCoeff@@q, jdd, FCenter@@pDD, FCenter@@qDD, CartAng@@pDD, CartAng@@qDDD, 8i, 3<, 8j, 3<D, 8p, 1, 7<, 8q, 1, 7<D êê MatrixForm êê Chop Since the kinetic-energy matrix is symmetrical, we need only to calculate the upper elements. The value of the elements K 11 = K 22 = hartree is the electronic kinetic energy of the hydrogen atom described by the STO-3G basis set (compared to the exact value of 0.5 hartree). Analogously, the element K 33 is the corresponding kinetic energy of an electron in the 1 s orbital of the oxygen atom, and so forth.
10 10 Minhhuy Hô and Julio Manuel Hernández-Pérez Conclusion We have provided an introduction to the evaluation of kinetic-energy integrals involving Gaussian-type basis functions both analytically and by use of recurrence relations. The results are sufficiently general so that no modification of the algorithm is needed when larger basis sets with more Gaussian primitives or primitives with larger angular momenta are employed. References [1] M. Hô and J. M. Hernández-Pérez, Evaluation of Gaussian Molecular Integrals I, The Mathematica Journal, doi: /tmj [2] Basis Set Exchange. (April 12, 2012) bse.pnl.gov/bse/portal. M. Hô and J. M. Hernández-Pérez, Evaluation of Gaussian Molecular Integrals II, The Mathematica Journal, dx.doi.org/doi: /tmj About the Authors Minhhuy Hô received his Ph.D. in theoretical chemistry at Queen s University, Kingston, Ontario, Canada in He is currently a professor at the Centro de Investigaciones Químicas at the Universidad Autónoma del Estado de Morelos in Cuernavaca, Morelos, México. Julio-Manuel Hernández-Pérez obtained his Ph.D. at the Universidad Autónoma del Estado de Morelos in He has been a professor of chemistry at the Facultad de Ciencias Químicas at the Benemérita Universidad Autónoma de Puebla since Minhhuy Hô Universidad Autónoma del Estado de Morelos Centro de Investigaciones Químicas Ave. Universidad, No. 1001, Col. Chamilpa Cuernavaca, Morelos, Mexico CP homh@uaem.mx Julio-Manuel Hernández-Pérez Benemérita Universidad Autónoma de Puebla Facultad de Ciencias Químicas Ciudad Universitaria, Col. San Manuel Puebla, Puebla, Mexico CP julio.hernandez@correo.buap.mx
Evaluation of Gaussian Molecular Integrals
The Mathematica Journal Evaluation of Gaussian Molecular Integrals III. Nuclear-Electron Attraction Integrals Minhhuy Hô Julio Manuel Hernández-Pérez This article carries out the evaluation of nuclear-electron
More informationEvaluation of Gaussian Molecular Integrals
The Mthemtic Journl Evlution of Gussin Moleculr Integrls I. verl Integrls Minhhuy Hô Julio Mnuel Hernández-Pérez This rticle discusses the evlution of moleculr overl integrls for Gussin-tye functions with
More informationApproximation with Interpolation Conditions of an Asymptotic Function Using K-algebraic Lagrange Interpolation
Approximation with Interpolation Conditions of an Asymptotic Function Using K-algebraic Lagrange Interpolation G. Martínez-Guzmán, M. Bustillo-Díaz, A. Rangel-Huerta, G. Juárez-Díaz, A. Ata-Pérez and N.
More informationCrystal growth and characterization of solvated organic charge-transfer complexes built on TTF and 9-dicyanomethylenefluorene derivatives
Electronic Supplementary Material (ESI) for CrystEngComm. This journal is The Royal Society of Chemistry 2015 Crystal growth and characterization of solvated organic charge-transfer complexes built on
More informationLecture 9: Molecular integral. Integrals of the Hamiltonian matrix over Gaussian-type orbitals
Lecture 9: Molecular integral evaluation Integrals of the Hamiltonian matrix over Gaussian-type orbitals Gaussian-type orbitals The de-facto standard for electronic-structure calculations is to use Gaussian-type
More informationLinear Equations and Matrix
1/60 Chia-Ping Chen Professor Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Gaussian Elimination 2/60 Alpha Go Linear algebra begins with a system of linear
More informationand ALTO SOLO C H, Rnnciman. Woman's Mistake f*" M>rqut4te\ ol ting about wilh a crutch 00 ac- Gr,,,,d Ri ''*' d5 L o,
BU AK A k A AD DA AB U XXX HA HUDAY U 3 92 3 > k z j - - Y Bk 73( 3 - Q H 2 H 9 k B 3 Bk 29 K 7 k 2 B k k k k Y Y D k Y A Uk x 22 B B x k - B B B 22 A B A 27 -- Ak A - - j B B D B Q A- U A-AK B k AD A
More informationA One-Slide Summary of Quantum Mechanics
A One-Slide Summary of Quantum Mechanics Fundamental Postulate: O! = a! What is!?! is an oracle! operator wave function (scalar) observable Where does! come from?! is refined Variational Process H! = E!
More informationSupporting Information
Supporting Information Macrocyclic Diorganotin Complexes of -Amino Acid Dithiocarbamates as Hosts for Ion-Pair Recognition Jorge Cruz-Huerta, a Manuel Carillo-Morales, a Ericka Santacruz-Juárez, a Irán
More information( ) = 9φ 1, ( ) = 4φ 2.
Chemistry 46 Dr Jean M Standard Homework Problem Set 6 Solutions The Hermitian operator A ˆ is associated with the physical observable A Two of the eigenfunctions of A ˆ are and These eigenfunctions are
More informationLecture 5 Least-squares
EE263 Autumn 2008-09 Stephen Boyd Lecture 5 Least-squares least-squares (approximate) solution of overdetermined equations projection and orthogonality principle least-squares estimation BLUE property
More informationSupplementary Information
Supplementary Information Eco-Friendly Synthesis of 2,3-Dihydroquinazolin-4(1H)-ones Catalyzed by FeCl 3 /Al 2 O 3 and Analysis of Large 1 H NMR Diastereotopic Effect Isabel Monreal, a Mariano Sánchez-Castellanos,
More informationMathematical Notes for E&M Gradient, Divergence, and Curl
Mathematical Notes for E&M Gradient, Divergence, and Curl In these notes I explain the differential operators gradient, divergence, and curl (also known as rotor), the relations between them, the integral
More informationZeta Functions of Burnside Rings for Symmetric and Alternating Groups
International Journal of Algebra, Vol. 6, 2012, no. 25, 1207-1220 Zeta Functions of Burnside Rings for Symmetric and Alternating Groups David Villa-Hernández Benemérita Universidad Autónoma de Puebla Facultad
More informationU AdS d+1 /CF T d. AdS d+1
U 00 1 2 U 0 1 1 2 U 2 0 3 4 AdS d+1 /CF T d U 0 1 1 2 AdS d+1 "i #i ± #i "i p 2 ( #i + "i)( #i + "i) 2 A := tr HA c A A c S A := tr HA A log A H = H A H A c A A c "ih" + #ih# 2 ( #i + "i)( #i + "i)
More informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationUsing Mathematica to model observations of surface magnetic fields
Using Mathematica to model observations of surface magnetic fields version of -Mar-0 (w/ changes marked in red type) Oleg, Thanks for the plots and further info on how Stokes V is computed from B z. While
More informationLeast Squares Regression
Least Squares Regression Chemical Engineering 2450 - Numerical Methods Given N data points x i, y i, i 1 N, and a function that we wish to fit to these data points, fx, we define S as the sum of the squared
More informationC. Show your answer in part B agrees with your answer in part A in the limit that the constant c 0.
Problem #1 A. A projectile of mass m is shot vertically in the gravitational field. Its initial velocity is v o. Assuming there is no air resistance, how high does m go? B. Now assume the projectile is
More informationFaculty of Science Final Examination Physics 457 Honours Quantum Physics II
Physics 457: Honours Quantum Physics II 1 Faculty of Science Final Examination Physics 457 Honours Quantum Physics II Examiner : Professor C. Gale Associate Examiner: Professor S. Jeon Date: Tuesday, April
More informationThe Hydrogen Atom. Dr. Sabry El-Taher 1. e 4. U U r
The Hydrogen Atom Atom is a 3D object, and the electron motion is three-dimensional. We ll start with the simplest case - The hydrogen atom. An electron and a proton (nucleus) are bound by the central-symmetric
More informationwould represent a 1s orbital centered on the H atom and φ 2px )+ 1 r 2 sinθ
Physical Chemistry for Engineers CHEM 4521 Homework: Molecular Structure (1) Consider the cation, HeH +. (a) Write the Hamiltonian for this system (there should be 10 terms). Indicate the physical meaning
More informationQuantum Orbits. Quantum Theory for the Computer Age Unit 9. Diving orbit. Caustic. for KE/PE =R=-3/8. for KE/PE =R=-3/8. p"...
W.G. Harter Coulomb Obits 6-1 Quantum Theory for the Computer Age Unit 9 Caustic for KE/PE =R=-3/8 F p' p g r p"... P F' F P Diving orbit T" T T' Contact Pt. for KE/PE =R=-3/8 Quantum Orbits W.G. Harter
More informationContents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9
MATH 32B-2 (8W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables Contents Multiple Integrals 3 2 Vector Fields 9 3 Line and Surface Integrals 5 4 The Classical Integral Theorems 9 MATH 32B-2 (8W)
More informationPH.D. PRELIMINARY EXAMINATION MATHEMATICS
UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem
More informationH 2 in the minimal basis
H 2 in the minimal basis Alston J. Misquitta Centre for Condensed Matter and Materials Physics Queen Mary, University of London January 27, 2016 Overview H 2 : The 1-electron basis. The two-electron basis
More informationAPPENDIX Z. USEFUL FORMULAS 1. Appendix Z. Useful Formulas. DRAFT 13:41 June 30, 2006 c J.D Callen, Fundamentals of Plasma Physics
APPENDIX Z. USEFUL FORMULAS 1 Appendix Z Useful Formulas APPENDIX Z. USEFUL FORMULAS 2 Key Vector Relations A B = B A, A B = B A, A A = 0, A B C) = A B) C A B C) = B A C) C A B), bac-cab rule A B) C D)
More informationSUPPORTING INFORMATION
SUPPORTING INFORMATION Self-Assembly of Fluorinated Boronic Esters and 4,4 -Bipyridine into 2:1 N B Adducts and Inclusion of Aromatic Guest Molecules in the Solid State: Application for the Separation
More information7.4* General logarithmic and exponential functions
7.4* General logarithmic and exponential functions Mark Woodard Furman U Fall 2010 Mark Woodard (Furman U) 7.4* General logarithmic and exponential functions Fall 2010 1 / 9 Outline 1 General exponential
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationEhrhart polynomial for lattice squares, cubes, and hypercubes
Ehrhart polynomial for lattice squares, cubes, and hypercubes Eugen J. Ionascu UWG, REU, July 10th, 2015 math@ejionascu.ro, www.ejionascu.ro 1 Abstract We are investigating the problem of constructing
More informationYingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2
Purdue University CHM 67300 Computational Quantum Chemistry REVIEW Yingwei Wang October 10, 2013 Review: Prof Slipchenko s class, Fall 2013 Contents 1 Hartree energy 2 2 Many-body system 2 3 Born-Oppenheimer
More informationParts Manual. EPIC II Critical Care Bed REF 2031
EPIC II Critical Care Bed REF 2031 Parts Manual For parts or technical assistance call: USA: 1-800-327-0770 2013/05 B.0 2031-109-006 REV B www.stryker.com Table of Contents English Product Labels... 4
More informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
More informationBasis sets for electron correlation
Basis sets for electron correlation Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 12th Sostrup Summer School Quantum Chemistry
More informationMATRIX ALGEBRA. or x = (x 1,..., x n ) R n. y 1 y 2. x 2. x m. y m. y = cos θ 1 = x 1 L x. sin θ 1 = x 2. cos θ 2 = y 1 L y.
as Basics Vectors MATRIX ALGEBRA An array of n real numbers x, x,, x n is called a vector and it is written x = x x n or x = x,, x n R n prime operation=transposing a column to a row Basic vector operations
More informationOPTI 511R, Spring 2018 Problem Set 1 Prof. R.J. Jones
OPTI 511R, Spring 2018 Problem Set 1 Prof. R.J. Jones Due in class Thursday, Jan 18, 2018 Math Refresher The first 11 questions are just to get you to review some very basic mathematical concepts needed
More informationDirectional Derivative and the Gradient Operator
Chapter 4 Directional Derivative and the Gradient Operator The equation z = f(x, y) defines a surface in 3 dimensions. We can write this as z f(x, y) = 0, or g(x, y, z) = 0, where g(x, y, z) = z f(x, y).
More informationPhysics 228 Today: Ch 41: 1-3: 3D quantum mechanics, hydrogen atom
Physics 228 Today: Ch 41: 1-3: 3D quantum mechanics, hydrogen atom Website: Sakai 01:750:228 or www.physics.rutgers.edu/ugrad/228 Happy April Fools Day Example / Worked Problems What is the ratio of the
More informationSTARK EFFECT IN P-TYPE DELTA-DOPED QUANTUM WELLS
Progress In Electromagnetics Research Letters, Vol. 2, 45 52, 2008 STARK EFFECT IN P-TYPE DELTA-DOPED QUANTUM WELLS A. M. Miteva Space Research Institute Bulgarian Academy of Sciences Moskovska Str. 6,
More informationAbsence of limit cycles for Kukles-type systems
AVANZA. Vol. iv. Fm - Iit, Uacj (014) 5 13. Isbn: 978-607-94-5-3 Absence of limit cycles for Kukles-type systems Osuna, Osvaldo, Rodríguez-Ceballos, Joel, Villaseñor, Gabriel Abstract In this note, we
More informationMathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.
Mathematical Tripos Part IA Lent Term 205 ector Calculus Prof B C Allanach Example Sheet Sketch the curve in the plane given parametrically by r(u) = ( x(u), y(u) ) = ( a cos 3 u, a sin 3 u ) with 0 u
More informationThe Projector Augmented Wave method
The Projector Augmented Wave method Advantages of PAW. The theory. Approximations. Convergence. 1 The PAW method is... What is PAW? A technique for doing DFT calculations efficiently and accurately. An
More informationChemistry 3502/4502. Exam III. All Hallows Eve/Samhain, ) This is a multiple choice exam. Circle the correct answer.
B Chemistry 3502/4502 Exam III All Hallows Eve/Samhain, 2003 1) This is a multiple choice exam. Circle the correct answer. 2) There is one correct answer to every problem. There is no partial credit. 3)
More informationAn improved CA model with anticipation for one-lane traffic flow
An improved CA model with anticipation for one-lane traffic flow MARÍA ELENA. LÁRRAGA JESÚS ANTONIO DEL RÍ0 Facultad de Ciencias, Computer Science Dept. Universidad Autónoma del Estado de Morelos Av. Universidad
More informationChem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 28, December 08, 2014
Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi Lecture 28, December 08, 2014 Solved Homework Water, H 2 O, involves 2 hydrogen atoms and an oxygen
More informationΨ g = N g (1s a + 1s b ) Themes. Take Home Message. Energy Decomposition. Multiple model wave fns possible. Begin describing brain-friendly MO
Themes Multiple model wave fns possible Different explanations of bonding Some you have seen before Errors vary (not consistently better or worse) Convenience varies Brain v. algebra v. computer Where
More informationAddition of Angular Momenta
Addition of Angular Momenta What we have so far considered to be an exact solution for the many electron problem, should really be called exact non-relativistic solution. A relativistic treatment is needed
More informationAlgorithm to compute the electric field gradient tensor in ionic crystals
INVESTIGACIÓN Revista Mexicana de Física 58 01 13 18 FEBRERO 01 Algorithm to compute the electric field gradient tensor in ionic crystals J.J. Hernández-Gómez, V. Marquina and R.W. Gómez Facultad de Ciencias,
More informationAbstract Algebra: Chapters 16 and 17
Study polynomials, their factorization, and the construction of fields. Chapter 16 Polynomial Rings Notation Let R be a commutative ring. The ring of polynomials over R in the indeterminate x is the set
More informationQuantum Mechanics: Vibration and Rotation of Molecules
Quantum Mechanics: Vibration and Rotation of Molecules 8th April 2008 I. 1-Dimensional Classical Harmonic Oscillator The classical picture for motion under a harmonic potential (mass attached to spring
More informationBound states of the hydrogen atom in parabolic coordinates
INVESTIGACIÓN REVISTA MEXICANA DE FÍSICA 54 6 454 458 DICIEMBRE 008 Bound states of the hydrogen atom in parabolic coordinates G.F. Torres del Castillo Departamento de Física Matemática, Instituto de Ciencias,
More information1.6. Quantum mechanical description of the hydrogen atom
29.6. Quantum mechanical description of the hydrogen atom.6.. Hamiltonian for the hydrogen atom Atomic units To avoid dealing with very small numbers, let us introduce the so called atomic units : Quantity
More informationTopics in Algebra and Analysis
Radmila Bulajich Manfrino José Antonio Gómez Ortega Rogelio Valdez Delgado Topics in Algebra and Analysis Preparing for the Mathematical Olympiad Radmila Bulajich Manfrino Facultad de Ciencias Universidad
More informationMODULE 8 Topics: Null space, range, column space, row space and rank of a matrix
MODULE 8 Topics: Null space, range, column space, row space and rank of a matrix Definition: Let L : V 1 V 2 be a linear operator. The null space N (L) of L is the subspace of V 1 defined by N (L) = {x
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationIntroduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule. Vesa Hänninen
Introduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the electronic energy
More informationMath 101 Study Session Spring 2016 Test 4 Chapter 10, Chapter 11 Chapter 12 Section 1, and Chapter 12 Section 2
Math 101 Study Session Spring 2016 Test 4 Chapter 10, Chapter 11 Chapter 12 Section 1, and Chapter 12 Section 2 April 11, 2016 Chapter 10 Section 1: Addition and Subtraction of Polynomials A monomial is
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part IA Tuesday 7 June 2005 1.30 to 4.30 PAPER 3 Before you begin read these instructions carefully. The examination paper is divided into two sections. Each question in Section II
More informationInverses and Elementary Matrices
Inverses and Elementary Matrices 1-12-2013 Matrix inversion gives a method for solving some systems of equations Suppose a 11 x 1 +a 12 x 2 + +a 1n x n = b 1 a 21 x 1 +a 22 x 2 + +a 2n x n = b 2 a n1 x
More informationQMI PRELIM Problem 1. All problems have the same point value. If a problem is divided in parts, each part has equal value. Show all your work.
QMI PRELIM 013 All problems have the same point value. If a problem is divided in parts, each part has equal value. Show all your work. Problem 1 L = r p, p = i h ( ) (a) Show that L z = i h y x ; (cyclic
More informationWorkshop: One-electron-atom orbitals Quantum aspects of physical chemistry
Workshop: One-electron-atom orbitals Quantum aspects of physical chemistry http://quantum.bu.edu/pltl/7.pdf Last updated Tuesday, November 15, 25 12:58:14-5: Copyright 25 Dan Dill (dan@bu.edu) Department
More informationAFM Review Test Review
Name: Class: Date: AFM Review Test Review Multiple Choice Identify the choice that best completes the statement or answers the question. What are the solutions of the inequality?. q + (q ) > 0 q < 3 q
More informationComplex Analysis Homework 4
Complex Analysis Homework 4 Isaac Defrain Steve Clanton David Holz April 9, 2009 Problem 1 Problem. Show that if f = u + iv is analytic then ( v) ( v) = 0. Explain this geometrically. Since f is analytic,
More informationChapter 2 Solutions. Chapter 2 Solutions 1. 4( z t) ( z t) z t
Chapter Solutions 1 Chapter Solutions.1. 1 z t ( z t) z z ( z t) t t It s a twice differentiable function of ( z t ), where is in the negative z direction. 1 y t ( y, t) ( y4 t) ( y 4 t) y y 8( y 4 t)
More informationRegression. Oscar García
Regression Oscar García Regression methods are fundamental in Forest Mensuration For a more concise and general presentation, we shall first review some matrix concepts 1 Matrices An order n m matrix is
More informationJaynes information entropy of small molecules: Numerical evidence of the Collins conjecture
PHYSICAL REVIEW A VOLUME 56, NUMBER 6 DECEMBER 1997 Jaynes information entropy of small molecules: Numerical evidence of the Collins conjecture Juan Carlos Ramírez, Catalina Soriano, Rodolfo O. Esquivel,
More informationSolutions to Chapter Review Questions, Chapter 0
Instructor s Solutions Manual, Chapter 0 Review Question 1 Solutions to Chapter Review Questions, Chapter 0 1. Explain how the points on the real line correspond to the set of real numbers. solution Start
More informationHomework Problem Set 5 Solutions. E e + H corr (a.u.) a.) Determine the bond dissociation enthalpy of ethane in kcal/mol (1 a.u. = kcal/mol).
Chemistry 380.37 Dr. Jean M. Standard Homework Problem Set 5 Solutions 1. Given below are the sum of electronic and thermal enthalpies, E e + H corr, from Hartree-Fock calculations using a 6-31G(d) basis
More informationSheet 06.6: Curvilinear Integration, Matrices I
Fakultät für Physik R: Rechenmethoden für Physiker, WiSe 5/6 Doent: Jan von Delft Übungen: Benedikt Bruognolo, Dennis Schimmel, Frauke Schwar, Lukas Weidinger http://homepages.physik.uni-muenchen.de/~vondelft/lehre/5r/
More informationHighly accurate evaluation of the few-body auxiliary functions and four-body integrals
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICSB: ATOMIC, MOLECULAR AND OPTICAL PHYSICS J. Phys. B: At. Mol. Opt. Phys. 36 (2003) 1857 1867 PII: S0953-4075(03)58918-0 Highly accurate evaluation of the
More informationSection 6.5. Least Squares Problems
Section 6.5 Least Squares Problems Motivation We now are in a position to solve the motivating problem of this third part of the course: Problem Suppose that Ax = b does not have a solution. What is the
More information7.1 Variational Principle
7.1 Variational Principle Suppose that you want to determine the ground-state energy E g for a system described by H, but you are unable to solve the time-independent Schrödinger equation. It is possible
More informationFigures S1- S12. 1D and 2D NMR Spectra for compound Figure S16. X Ray crystallography of compound 1 20
Supporting Information Sulfur-containing Aristoloxazines and Other Constituents of the Roots of Aristolochia orbicularis María Yolanda Rios,*, Víctor Navarro, M. Ángeles Ramírez-Cisneros and Enrique Salazar-Rios
More informationElectron States of Diatomic Molecules
IISER Pune March 2018 Hamiltonian for a Diatomic Molecule The hamiltonian for a diatomic molecule can be considered to be made up of three terms Ĥ = ˆT N + ˆT el + ˆV where ˆT N is the kinetic energy operator
More informationBACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES
No. of Printed Pages : 5 BCS-054 BACHELOR OF COMPUTER APPLICATIONS (BCA) (Revised) Term-End Examination December, 2015 058b9 BCS-054 : COMPUTER ORIENTED NUMERICAL TECHNIQUES Time : 3 hours Maximum Marks
More informationHandbook of Computational Quantum Chemistry. DAVID B. COOK The Department of Chemistry, University of Sheffield
Handbook of Computational Quantum Chemistry DAVID B. COOK The Department of Chemistry, University of Sheffield Oxford New York Tokyo OXFORD UNIVERSITY PRESS 1998 CONTENTS 1 Mechanics and molecules 1 1.1
More informationSummary of time independent electrodynamics
hapter 10 Summary of time independent electrodynamics 10.1 Electrostatics Physical law oulomb s law charges as origin of electric field Superposition principle ector of the electric field E(x) in vacuum
More informationLecture 24 Origins of Magnetization (A number of illustrations in this lecture were generously provided by Prof. Geoffrey Beach)
Lecture 4 Origins of Magnetization (A number of illustrations in this lecture were generously provided by Prof. Geoffrey Beach) Today 1. Magnetic dipoles.. Orbital and spin angular momenta. 3. Non-interacting
More informationQualification Exam: Mathematical Methods
Qualification Exam: Mathematical Methods Name:, QEID#41534189: August, 218 Qualification Exam QEID#41534189 2 1 Mathematical Methods I Problem 1. ID:MM-1-2 Solve the differential equation dy + y = sin
More informationOrigin of the first Hund rule in He-like atoms and 2-electron quantum dots
in He-like atoms and 2-electron quantum dots T Sako 1, A Ichimura 2, J Paldus 3 and GHF Diercksen 4 1 Nihon University, College of Science and Technology, Funabashi, JAPAN 2 Institute of Space and Astronautical
More informationChemistry 532 Practice Final Exam Fall 2012 Solutions
Chemistry 53 Practice Final Exam Fall Solutions x e ax dx π a 3/ ; π sin 3 xdx 4 3 π cos nx dx π; sin θ cos θ + K x n e ax dx n! a n+ ; r r r r ˆL h r ˆL z h i φ ˆL x i hsin φ + cot θ cos φ θ φ ) ˆLy i
More informationSchrödinger Pauli equation for spin-3/2 particles
INVESTIGACIÓN REVISTA MEXICANA DE FÍSICA 50 3) 306 310 JUNIO 004 Schrödinger Pauli euation for spin-3/ particles G.F. Torres del Castillo Departamento de Física Matemática, Instituto de Ciencias Universidad
More informationEXAM INFORMATION. Radial Distribution Function: B is the normalization constant. d dx. p 2 Operator: Heisenberg Uncertainty Principle:
EXAM INFORMATION Radial Distribution Function: P() r RDF() r Br R() r B is the normalization constant., p Operator: p ^ d dx Heisenberg Uncertainty Principle: n ax n! Integrals: xe dx n1 a x p Particle
More informationThe Central Force Problem: Hydrogen Atom
The Central Force Problem: Hydrogen Atom B. Ramachandran Separation of Variables The Schrödinger equation for an atomic system with Z protons in the nucleus and one electron outside is h µ Ze ψ = Eψ, r
More informationChemistry 1A Midterm Exam 1 February 12, Potentially Useful Information. Violet Blue Green Yellow Orange Red Wavelength (nm)
Chemistry 1A Midterm Exam 1 February 12, 2013 Professor Pines 5 pages total Student Name: Student ID#: Potentially Useful Information Violet Blue Green Yellow Orange Red 400 500 600 700 Wavelength (nm)
More informationClassical Mechanics in Hamiltonian Form
Classical Mechanics in Hamiltonian Form We consider a point particle of mass m, position x moving in a potential V (x). It moves according to Newton s law, mẍ + V (x) = 0 (1) This is the usual and simplest
More informationCollection of formulae Quantum mechanics. Basic Formulas Division of Material Science Hans Weber. Operators
Basic Formulas 17-1-1 Division of Material Science Hans Weer The de Broglie wave length λ = h p The Schrödinger equation Hψr,t = i h t ψr,t Stationary states Hψr,t = Eψr,t Collection of formulae Quantum
More informationG : Quantum Mechanics II
G5.666: Quantum Mechanics II Notes for Lecture 5 I. REPRESENTING STATES IN THE FULL HILBERT SPACE Given a representation of the states that span the spin Hilbert space, we now need to consider the problem
More informationChem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components.
Chem 44 Review for Exam Hydrogenic atoms: The Coulomb energy between two point charges Ze and e: V r Ze r Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative
More information2. A die is rolled 3 times, the probability of getting a number larger than the previous number each time is
. If P(A) = x, P = 2x, P(A B) = 2, P ( A B) = 2 3, then the value of x is (A) 5 8 5 36 6 36 36 2. A die is rolled 3 times, the probability of getting a number larger than the previous number each time
More informationThe Schroedinger equation
The Schroedinger equation After Planck, Einstein, Bohr, de Broglie, and many others (but before Born), the time was ripe for a complete theory that could be applied to any problem involving nano-scale
More informationLegendre s Equation. PHYS Southern Illinois University. October 18, 2016
Legendre s Equation PHYS 500 - Southern Illinois University October 18, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 18, 2016 1 / 11 Legendre s Equation Recall We are trying
More informationPHY 396 K. Solutions for problems 1 and 2 of set #5.
PHY 396 K. Solutions for problems 1 and of set #5. Problem 1a: The conjugacy relations  k,  k,, Ê k, Ê k, follow from hermiticity of the Âx and Êx quantum fields and from the third eq. 6 for the polarization
More informationFORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 2017
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Physics 8.07: Electromagnetism II November 5, 207 Prof. Alan Guth FORMULA SHEET FOR QUIZ 2 Exam Date: November 8, 207 A few items below are marked
More informationQuantum Theory. Thornton and Rex, Ch. 6
Quantum Theory Thornton and Rex, Ch. 6 Matter can behave like waves. 1) What is the wave equation? 2) How do we interpret the wave function y(x,t)? Light Waves Plane wave: y(x,t) = A cos(kx-wt) wave (w,k)
More informationMANY ELECTRON ATOMS Chapter 15
MANY ELECTRON ATOMS Chapter 15 Electron-Electron Repulsions (15.5-15.9) The hydrogen atom Schrödinger equation is exactly solvable yielding the wavefunctions and orbitals of chemistry. Howev er, the Schrödinger
More informationy mx 25m 25 4 circle. Then the perpendicular distance of tangent from the centre (0, 0) is the radius. Since tangent
Mathematics. The sides AB, BC and CA of ABC have, 4 and 5 interior points respectively on them as shown in the figure. The number of triangles that can be formed using these interior points is () 80 ()
More informationPart III: The Nuclear Many-Body Problem
Part III: The Nuclear Many-Body Problem To understand the properties of complex nuclei from first principles Microscopic Valence- Space Interactions Model spaces Many-body perturbation theory (MBPT) Calculating
More informationa non-standard generating function for continuous dual q-hahn polynomials una función generatriz no estándar para polinomios q-hahn duales continuos
Revista de Matemática: Teoría y Aplicaciones 2011 181 : 111 120 cimpa ucr issn: 1409-2433 a non-standard generating function for continuous dual -hahn polynomials una función generatriz no estándar para
More information