|
|
| Line 1: |
Line 1: |
| {{distinguish2|the [[Born approximation]]}}
| | These days, it's known by everyone that the prices of the merchandise are getting on improving. This is the big reason due to which existing standard of several folks are really low. Many times people experience it tough to buy their preferred things due to not enough financing.<br><br>But, every dilemma has some option and the same expression applies on the problem that's being discussed. One of many most distinguished methods for getting reduce this issue will be to avail discount discount voucher codes. These are the codes which are the online type of the discount promotion that people usually observe inside our daily newspapers. By giving these deals to the shopkeepers, you can get superior concessions about the product that you want. Whenever you arrive at onlineshopping, the aforesaid voucher codes attain exactly the same purpose. With all the aid of these codes you will get good-looking abatement around the goods online websites. However, the technique of using these codes is really not the same as that of discounts. Nevertheless, these have become user friendly. You just must enter the code in the pack that you could see to the holder page of the buying site. Subsequently some volume can automatically be found from your unique level of the merchandise.<br><br>These might be benefited from the variety of web sites. There are always a large numbers of sites which are specifically designed for delivering such valuable vouchers. Aside, some store and affiliate websites can also be serving precisely the same function, however selling and purchase of these products can be completed on these websites. If you're uninformed about these websites which are providing the codes subsequently it's super easy to find them. Search engines aid a great deal to get the details about such websites. On searching regarding the coupon codes you obtain a long list of websites which provide these codes. These codes are often provided by the suppliers for growing their customers.<br><br>Once you reach these website which is offering these codes, then you certainly observe several categories of these deals. These categories may include, publications, music, gifts, flowers etc. if you pick any of these classes, you will get the deals on the same. You're able to choose any of these depending on your requirements. On the other hand, on a number of the affiliate sites you obtain many amenities which enhance your understanding of these codes. These websites aid the people to add their websites in addition to evaluations. Using the help of this, it is possible to present in addition to get the advice. For those who have the sound information about the coupons then you can certainly add your write-ups to these sites. To the other-hand, you're able to get the fair information by examining the opinions of other consumers. Aside, many external sites give you the articles around the voucher codes whereby you're able to consider the comprehensive info.<br><br>There is without doubt in the undeniable fact that these discount voucher codes show to be very useful for that people in addition to suppliers. These deals on one area boost the living standard of the folks by minimizing the costs of the products, to the other-hand they boost the purchaser foot of the shops and producers. More on our site [http://www.houstonbattalion.com/activity/p/16556/ [http://www.houstonbattalion.com/activity/p/16556/ read more here]]. |
| {{multiple issues|
| |
| {{more footnotes|date=March 2013}}
| |
| {{Technical|date=September 2010}}
| |
| }}
| |
| | |
| In [[quantum chemistry]], the computation of the energy and the [[wavefunction]] of an average-size [[molecule]] is a formidable task that is alleviated by the '''Born–Oppenheimer''' ('''BO''') '''approximation''', named after [[Max Born]] and [[J. Robert Oppenheimer]]. For instance the [[benzene]] molecule consists of 12 nuclei and 42 electrons. The [[Schrödinger equation|time independent Schrödinger equation]], which must be solved to obtain the energy and wavefunction of this molecule, is a [[partial differential equation|partial differential eigenvalue equation]] in 162 variables—the spatial coordinates of ''the electrons and the nuclei''. The BO approximation makes it possible to compute the wavefunction in two less complicated consecutive steps. This approximation was proposed in 1927, in the early period of quantum mechanics, by Born and Oppenheimer and is still indispensable in quantum chemistry.
| |
| | |
| In basic terms, it allows the wavefunction of a molecule to be broken into its electronic and nuclear (vibrational, rotational) components.
| |
| | |
| :: <math> \Psi_\mathrm{total} = \psi_\mathrm{electronic} \times \psi_\mathrm{nuclear} </math>
| |
| | |
| In the first step of the BO approximation the ''electronic'' Schrödinger equation is solved, yielding the wavefunction <math>\psi_{\text{electronic}}</math> depending on electrons only. For benzene this wavefunction depends on 126 electronic coordinates. During this solution the nuclei are fixed in a certain configuration, very often the equilibrium configuration. If the effects of the quantum mechanical nuclear motion are to be studied, for instance because a [[molecular vibration|vibrational spectrum]] is required, this electronic computation must be in nuclear coordinates. In the second step of the BO approximation this function serves as a potential in a Schrödinger equation ''containing only the nuclei''—for benzene an equation in 36 variables.
| |
| | |
| The success of the BO approximation is due to the high ratio between nuclear and electronic masses. The approximation is an important tool of quantum chemistry; without it only the lightest molecule, H<sub>2</sub>, could be handled, and all computations of molecular wavefunctions for larger molecules make use of it. Even in the cases where the BO approximation breaks down, it is used as a point of departure for the computations.
| |
| | |
| The electronic energies, constituting the nuclear potential, consist of kinetic energies, interelectronic repulsions and electron–nuclear attractions. In a [[handwaving]] manner the nuclear potential is taken to be an averaged electron–nuclear attraction.
| |
| The BO approximation follows from the inertia of electrons being considered to be negligible in comparison to the atom to which they are bound.
| |
| | |
| In molecular [[spectroscopy]], because the ratios of the periods of the electronic, vibrational and rotational energies are each related to each other on scales in the order of a thousand, the Born–Oppenheimer name has also been attached to the approximation where the energy components are treated separately.
| |
| | |
| :: <math> E_\mathrm{total} = E_\mathrm{electronic} + E_\mathrm{vibrational} + E_\mathrm{rotational}+ E_\mathrm{nuclear}</math>
| |
| | |
| The nuclear spin energy is so small that it is normally omitted.
| |
| | |
| ==Short description==
| |
| The [[Max Born|Born]]–[[Robert Oppenheimer|Oppenheimer]] (BO) [[approximation]] is ubiquitous in [[Quantum chemistry|quantum chemical]] calculations of molecular [[wavefunction]]s. It consists of two steps.
| |
| | |
| In the first step the nuclear [[kinetic energy]] is neglected,{{ref|fn1|a}} that is, the corresponding operator ''T''<sub>n</sub> is subtracted from the total [[molecular Hamiltonian]]. In the remaining electronic Hamiltonian ''H''<sub>e</sub> the nuclear positions enter as parameters. The electron–nucleus interactions are ''not'' removed and the electrons still "feel" the [[Coulomb potential]] of the nuclei clamped at certain positions in space. (This first step of the BO approximation is therefore often referred to as the ''clamped nuclei'' approximation.)
| |
|
| |
| The electronic [[Schrödinger equation]]
| |
| | |
| ::<math> H_\mathrm{e}(\mathbf{r,R} )\; \chi(\mathbf{r,R}) = E_\mathrm{e} \; \chi(\mathbf{r,R}) </math>
| |
| | |
| is solved (out of necessity, approximately). The quantity '''r''' stands for all electronic coordinates and '''R''' for all nuclear coordinates. The electronic energy [[eigenvalue]] ''E''<sub>e</sub> depends on the chosen positions '''R''' of the nuclei. Varying these positions '''R''' in small steps and repeatedly solving the electronic [[Schrödinger equation]], one obtains ''E''<sub>e</sub> as a function of '''R'''. This is the [[potential energy surface]] (PES): ''E''<sub>e</sub>('''R''') . Because this procedure of recomputing the electronic wave functions as a function of an infinitesimally changing nuclear geometry is reminiscent of the conditions for the [[adiabatic theorem]], this manner of obtaining a PES is often referred to as the ''adiabatic approximation'' and the PES itself is called an ''adiabatic surface''.{{ref|fn2|b}}
| |
| | |
| In the second step of the BO approximation the nuclear kinetic energy ''T''<sub>n</sub> (containing partial derivatives with respect to the components of '''R''') is reintroduced and the Schrödinger equation for the nuclear motion{{ref|fn2|c}}
| |
| | |
| :<math>
| |
| \left[ T_\mathrm{n} + E_\mathrm{e}(\mathbf{R})\right] \phi(\mathbf{R}) = E \phi(\mathbf{R}) </math>
| |
| is solved. This second step of the BO approximation involves separation of vibrational, translational, and rotational motions. This can be achieved by application of the [[Eckart conditions]]. The eigenvalue ''E'' is the total energy of the molecule, including contributions from electrons, nuclear vibrations, and overall rotation and translation of the molecule.
| |
| | |
| ===Footnotes===
| |
| *{{note|fn1|a}} This step is often justified by stating that "the heavy nuclei move more slowly than the light [[electron]]s." Classically this statement makes sense only if the [[momentum]] ''p'' of electrons and nuclei is of the same order of magnitude. In that case ''m''<sub>nuc</sub> >> ''m''<sub>elec</sub> implies ''p''<sup>2</sup>/(2''m''<sub>nuc</sub>) << ''p''<sup>2</sup>/(2''m''<sub>elec</sub>). It is easy to show that for two bodies in circular orbits around their center of mass (regardless of individual masses), the momentum of the two bodies is equal and opposite, and that for any collection of particles in the center of mass frame, the net momentum is zero. Given that the center of mass frame is the lab frame (where the molecule is stationary), the momentum of the nuclei must be equal and opposite to that of the electrons. A hand waving justification can be derived from quantum mechanics as well. Recall that the corresponding operators do not contain mass and think of the molecule as a box containing the electrons and nuclei and see [[particle in a box]]. Since the kinetic energy is ''p''<sup>2</sup>/(2''m''), it follows that, indeed, the kinetic energy of the nuclei in a molecule is usually much smaller than the kinetic energy of the electrons, the mass ratio being on the order of 10<sup>4</sup>).
| |
| | |
| *{{note|fn2|b}} It is assumed, in accordance with the [[adiabatic theorem]], that the same electronic state (for instance the electronic ground state) is obtained upon small changes of the nuclear geometry. The method would give a discontinuity (jump) in the PES if electronic state-switching would occur.
| |
| | |
| *{{note|fn3|c}} This equation is time-independent and stationary wavefunctions for the nuclei are obtained, nevertheless it is traditional to use the word "motion" in this context, although classically motion implies time-dependence.
| |
| | |
| ==Derivation of the Born–Oppenheimer approximation==
| |
| It will be discussed how the BO approximation may be derived and under which conditions it is applicable. At the same time we will show how the BO approximation may be improved by including [[vibronic coupling]]. To that end the second step of the BO approximation is generalized to a set of coupled eigenvalue equations depending on nuclear coordinates only. Off-diagonal elements in these equations are shown to be nuclear kinetic energy terms.
| |
|
| |
| It will be shown that the BO approximation can be trusted whenever the PESs, obtained from the solution of the electronic Schrödinger equation, are well separated:
| |
| | |
| :<math> E_0(\mathbf{R}) \ll E_1(\mathbf{R}) \ll E_2(\mathbf{R}) \ll \cdots \text{ for all }\mathbf{R}</math>.
| |
| | |
| We start from the ''exact'' non-relativistic, time-independent molecular Hamiltonian:
| |
| | |
| :<math>
| |
| H= H_\mathrm{e} + T_\mathrm{n} \,
| |
| </math>
| |
| with
| |
| :<math>
| |
| H_\mathrm{e}=
| |
| -\sum_{i}{\frac{1}{2}\nabla_i^2}-
| |
| \sum_{i,A}{\frac{Z_A}{r_{iA}}} + \sum_{i>j}{\frac{1}{r_{ij}}}+ \sum_{A > B}{\frac{Z_A Z_B}{R_{AB}}}
| |
| \quad\mathrm{and}\quad T_\mathrm{n}=-\sum_{A}{\frac{1}{2M_A}\nabla_A^2}.
| |
| </math>
| |
| | |
| The position vectors <math>\mathbf{r}\equiv \{\mathbf{r}_i\}</math>
| |
| of the electrons and the position vectors <math>\mathbf{R}\equiv \{\mathbf{R}_A = (R_{Ax},\,R_{Ay},\,R_{Az})\}</math> of the nuclei are with respect to a Cartesian [[inertial frame]]. Distances between particles are written as <math> r_{iA} \equiv |\mathbf{r}_i - \mathbf{R}_A|</math> (distance between electron ''i'' and nucleus ''A'') and similar definitions hold for <math>r_{ij}\;</math> and <math> R_{AB}\,</math>.
| |
| We assume that the molecule is in a homogeneous
| |
| (no external force) and isotropic (no external torque) space. The only interactions are the
| |
| Coulomb interactions between the electrons and nuclei. The Hamiltonian is expressed in [[atomic units]], so that we do not see Planck's constant, the dielectric constant of the vacuum, electronic charge, or electronic mass in this formula. The only constants explicitly entering the formula are ''Z''<sub>A</sub> and ''M''<sub>A</sub>—the atomic number and mass of nucleus ''A''.
| |
| | |
| It is useful to introduce the total nuclear momentum and to rewrite the nuclear kinetic energy operator as follows:
| |
| | |
| :<math> T_\mathrm{n} = \sum_{A} \sum_{\alpha=x,y,z} \frac{P_{A\alpha} P_{A\alpha}}{2M_A}
| |
| \quad\mathrm{with}\quad
| |
| P_{A\alpha} = -i {\partial \over \partial R_{A\alpha}}. </math>
| |
| | |
| Suppose we have ''K'' electronic eigenfunctions <math>\chi_k (\mathbf{r}; \mathbf{R})</math> of <math>H_\mathrm{e}\,</math>, that is, we have solved
| |
| | |
| :<math>
| |
| H_\mathrm{e}\;\chi_k (\mathbf{r};\mathbf{R}) = E_k(\mathbf{R})\;\chi_k (\mathbf{r};\mathbf{R}) \quad\mathrm{for}\quad k=1,\ldots, K.
| |
| </math>
| |
| | |
| The electronic wave functions <math>\chi_k\,</math> will be taken to be real, which is possible when there are no magnetic or spin interactions. The ''parametric dependence'' of the functions <math>\chi_k\,</math> on the nuclear coordinates is indicated by the symbol after the semicolon. This indicates that, although <math>\chi_k\,</math> is a real-valued function of <math>\mathbf{r}</math>, its functional form depends on <math>\mathbf{R}</math>.
| |
| | |
| For example, in the molecular-orbital-linear-combination-of-atomic-orbitals [[Molecular orbital#Qualitative discussion|(LCAO-MO)]] approximation, <math>\chi_k\,</math> is a molecular orbital (MO) given as a linear expansion of atomic orbitals (AOs). An AO depends visibly on the coordinates of an electron, but the nuclear coordinates are not explicit in the MO. However, upon change of geometry, i.e., change of <math>\mathbf{R}</math>, the LCAO coefficients obtain different values and we see corresponding changes in the functional form of the MO <math>\chi_k\,</math>.
| |
| | |
| We will assume that the parametric dependence is continuous and differentiable, so that it is meaningful to consider
| |
| | |
| :<math>
| |
| P_{A\alpha}\chi_k (\mathbf{r};\mathbf{R}) = - i \frac{\partial\chi_k (\mathbf{r};\mathbf{R})}{\partial R_{A\alpha}} \quad \mathrm{for}\quad \alpha=x,y,z,
| |
| </math>
| |
| which in general will not be zero. | |
| | |
| The total wave function <math>\Psi(\mathbf{R},\mathbf{r})</math> is expanded in terms of <math>\chi_k (\mathbf{r}; \mathbf{R})</math>:
| |
| :<math>
| |
| \Psi(\mathbf{R}, \mathbf{r}) = \sum_{k=1}^K \chi_k(\mathbf{r};\mathbf{R}) \phi_k(\mathbf{R}),
| |
| </math>
| |
| with
| |
| :<math>
| |
| \langle\,\chi_{k'}(\mathbf{r};\mathbf{R})\,|\, \chi_k(\mathbf{r};\mathbf{R})\rangle_{(\mathbf{r})} = \delta_{k' k}
| |
| </math>
| |
| and where the subscript <math>(\mathbf{r})</math> indicates that the integration, implied by the [[bra–ket notation]], is over electronic coordinates only. By definition, the matrix with general element
| |
| :<math> \big(\mathbb{H}_\mathrm{e}(\mathbf{R})\big)_{k'k} \equiv \langle \chi_{k'}(\mathbf{r};\mathbf{R})
| |
| | H_\mathrm{e} |
| |
| \chi_k(\mathbf{r};\mathbf{R})\rangle_{(\mathbf{r})} = \delta_{k'k} E_k(\mathbf{R})
| |
| </math>
| |
| is diagonal. After multiplication by the real function <math>\chi_{k'}(\mathbf{r};\mathbf{R})</math> from the left and integration over the electronic coordinates <math>\mathbf{r}</math> the total Schrödinger equation
| |
| :<math>
| |
| H\;\Psi(\mathbf{R},\mathbf{r}) = E \; \Psi(\mathbf{R},\mathbf{r})
| |
| </math>
| |
| is turned into a set of ''K'' coupled eigenvalue equations depending on nuclear coordinates only
| |
| | |
| :<math> \left[ \mathbb{H}_\mathrm{n}(\mathbf{R}) + \mathbb{H}_\mathrm{e}(\mathbf{R}) \right]
| |
| \; \boldsymbol{\phi}(\mathbf{R}) = E\; \boldsymbol{\phi}(\mathbf{R}).
| |
| </math>
| |
| | |
| The column vector <math>\boldsymbol{\phi}(\mathbf{R})</math> has elements <math>\phi_k(\mathbf{R}),\; k=1,\ldots,K</math>. The matrix <math>\mathbb{H}_\mathrm{e}(\mathbf{R})</math> is diagonal and the nuclear Hamilton matrix is non-diagonal with the following off-diagonal (''vibronic coupling'') terms,
| |
| :<math>
| |
| \big(\mathbb{H}_\mathrm{n}(\mathbf{R})\big)_{k'k} = \langle\chi_{k'}(\mathbf{r};\mathbf{R}) |
| |
| T_\mathrm{n}|\chi_k(\mathbf{r};\mathbf{R})\rangle_{(\mathbf{r})}.
| |
| </math>
| |
| | |
| The vibronic coupling in this approach is through nuclear kinetic energy terms.
| |
| Solution of these coupled equations gives an approximation for energy and wavefunction that goes beyond the Born–Oppenheimer approximation.
| |
| Unfortunately, the off-diagonal kinetic energy terms are usually difficult to handle. This is why often a [[diabatic]] transformation is applied, which retains part of the nuclear kinetic energy terms on the diagonal, removes the kinetic energy terms from the off-diagonal and creates coupling terms between the adiabatic PESs on the off-diagonal.
| |
| | |
| If we can neglect the off-diagonal elements the equations will uncouple and simplify drastically. In order to show when this neglect is justified, we
| |
| suppress the coordinates in the notation and write, by applying the [[Leibniz rule (generalized product rule)|Leibniz rule]] for differentiation, the matrix elements of <math>T_{\textrm{n}}</math> as
| |
| :<math>
| |
| \mathrm{H_n}(\mathbf{R})_{k'k}\equiv
| |
| \big(\mathbb{H}_\mathrm{n}(\mathbf{R})\big)_{k'k}
| |
| = \delta_{k'k} T_{\textrm{n}}
| |
| + \sum_{A,\alpha}\frac{1}{M_A} \langle\chi_{k'}|\big(P_{A\alpha}\chi_k\big)\rangle_{(\mathbf{r})} P_{A\alpha} + \langle\chi_{k'}|\big(T_\mathrm{n}\chi_k\big)\rangle_{(\mathbf{r})}.
| |
| </math>
| |
| | |
| The diagonal (<math>k'=k</math>) matrix elements <math>\langle\chi_{k}|\big(P_{A\alpha}\chi_k\big)\rangle_{(\mathbf{r})}</math> of the operator <math>P_{A\alpha}\,</math> vanish, because we assume time-reversal invariant so <math>\chi_k</math> can be chosen to be always real. The off-diagonal matrix elements satisfy
| |
| | |
| :<math>
| |
| \langle\chi_{k'}|\big(P_{A\alpha}\chi_k\big)\rangle_{(\mathbf{r})}
| |
| =
| |
| \frac{\langle\chi_{k'}
| |
| |\big[P_{A\alpha}, H_\mathrm{e}\big] |
| |
| \chi_k\rangle_{(\mathbf{r})}} {E_{k}(\mathbf{R})- E_{k'}(\mathbf{R})}.
| |
| </math>
| |
| The matrix element in the numerator is
| |
| :<math>
| |
| \langle\chi_{k'}
| |
| |\big[P_{A\alpha}, H_\mathrm{e}\big] |
| |
| \chi_k\rangle_{(\mathbf{r})} =
| |
| iZ_A\sum_i \;\langle\chi_{k'}|\frac{(\mathbf{r}_{iA})_\alpha}{r_{iA}^3}|\chi_k\rangle_{(\mathbf{r})}
| |
| \;\;\mathrm{with}\;\; \mathbf{r}_{iA} \equiv \mathbf{r}_i - \mathbf{R}_A .
| |
| </math>
| |
| The matrix element of the one-electron operator appearing on the right hand side is finite.
| |
| When the two surfaces come close, <math>{E_{k}(\mathbf{R})\approx E_{k'}(\mathbf{R})}</math>, the nuclear momentum coupling term becomes large and is no longer negligible. This is the case where
| |
| the BO approximation breaks down and a coupled set of nuclear motion equations must be considered, instead of the one equation appearing in the second step of the BO approximation.
| |
| | |
| Conversely, if all surfaces are well separated, all off-diagonal terms can be neglected and hence the whole matrix of <math>P^{A}_\alpha</math> is effectively zero. The third term on the right hand side of the expression for the matrix element of ''T''<sub>n</sub> (the ''Born–Oppenheimer diagonal correction'') can approximately be written as the matrix of <math>P^{A}_\alpha</math> squared and, accordingly, is then negligible also. Only the first (diagonal) kinetic energy term in this equation survives in the case of well-separated surfaces and a diagonal, uncoupled, set of nuclear motion equations results,
| |
| :<math>
| |
| \left[ T_\mathrm{n} +E_k(\mathbf{R})\right] \; \phi_k(\mathbf{R}) =
| |
| E \phi_k(\mathbf{R})
| |
| \quad\mathrm{for}\quad k=1,\ldots, K,
| |
| </math>
| |
| which are the normal second-step of the BO equations discussed above. | |
| | |
| We reiterate that when two or more potential energy surfaces approach each other, or even cross, the Born–Oppenheimer
| |
| approximation breaks down and one must fall back on the coupled equations. Usually one
| |
| invokes then the [[diabatic]] approximation.
| |
| | |
| ==The Born–Oppenheimer approximation with the correct symmetry==
| |
| To include the correct symmetry within the Born–Oppenheimer (BO) approximation,<ref name=BornOppie>{{cite journal|author1=Max Born |author2=J. Robert Oppenheimer|title=Zur Quantentheorie der Molekeln|journal=Annalen der Physik|year=1927|volume=389|issue=20|pages=457–484|doi=10.1002/andp.19273892002|url=http://onlinelibrary.wiley.com/doi/10.1002/andp.19273892002/abstract;jsessionid=12CCB60F149391AAE8E51282EACEBF33.d02t04|accessdate=28 May 2013|trans_title=On the Quantum Theory of Molecules|language=German|bibcode = 1927AnP...389..457B }}</ref><ref>M. Born and K. Huang, Dynamical Theory of Crystal Lattices, 1954 (Oxford University Press, New York), Chapter IV</ref> a molecular system presented in terms of (mass-dependent) nuclear coordinates, <math>\mathbf{q}</math>, and formed by the two lowest BO adiabatic potential energy surfaces (PES), <math>u_1(\mathbf{q})</math> and <math>u_2 (\mathbf{q})</math>, is considered. To insure the validity of the BO approximation the energy of the system, E, is assumed to be low enough so that <math>u_2 (\mathbf{q})</math> becomes a closed PES in the region of interest, with the exception of sporadic infinitesimal sites surrounding degeneracy points formed by <math>{u}_{1} (\mathbf{q})</math> and <math>{u}_{2} (\mathbf{q})</math> (designated as (1,2) degeneracy points).
| |
| | |
| The starting point is the nuclear adiabatic BO (matrix) equation written in the form:<ref>M. Baer, Beyond Born–Oppenheimer: Electronic non-Adiabatic Coupling Terms and Conical Intersections, 2006 (Wiley and Sons, Inc., Hoboken, N.J.), Chapter 2</ref>
| |
| | |
| : <math>\frac{\hbar^2}{2m} (\nabla + \tau)^2 \Psi + (\mathbf{u} - E)\Psi = 0 </math>
| |
| | |
| where <math>\Psi(\mathbf{q}) </math> is a column vector that contains the unknown nuclear wave functions <math>\psi_{k}(\mathbf{q})</math>, <math>\mathbf{u}(\mathbf{q})</math> is a diagonal matrix which contains the corresponding adiabatic potential energy surfaces <math>u_{k}(\mathbf{q})</math>, m is the reduced mass of the nuclei, E is the total energy of the system, <math>\nabla</math> is the grad operator with respect to the nuclear coordinates <math>\mathbf{q}</math> and <math>\mathbf{\tau} (\mathbf{q})</math> is a matrix which contains the vectorial Non-Adiabatic Coupling Terms (NACT):
| |
| | |
| <math>\mathbf{\tau}_{jk} = \langle\zeta_{j} | \nabla\zeta_{k}\rangle</math>
| |
| | |
| Here <math>|\zeta_{n}\rangle ; n = j,k</math> are eigenfunctions of the [[electronic Hamiltonian]] assumed to form a complete [[Hilbert space]] in the given region in [[configuration space]].
| |
| | |
| To study the scattering process taking place on the two lowest surfaces one extracts, from the above BO equation, the two corresponding equations: | |
| | |
| : <math>-\frac{\hbar^2}{2m} \nabla^{2}\psi_{1} + (\tilde{u}_{1}- E)\psi_{1} - \frac{\hbar^2}{2m} [2\mathbf{\tau}_{12}\nabla + \nabla \mathbf{\tau}_{12}]\psi_{2} = 0 </math>
| |
| | |
| : <math>-\frac{\hbar^2}{2m} \nabla^{2}\psi_{2} + (\tilde{u}_{2}- E)\psi_{2} + \frac{\hbar^2}{2m} [2\mathbf{\tau}_{12}\nabla + \nabla \mathbf{\tau}_{12}]\psi_{1} = 0 </math>
| |
| | |
| where <math>\tilde{u}_{k}(\mathbf{q}) = u_{k}(\mathbf{q}) + (\hbar^{2}/2m)\tau_{12}^{2} ; k = 1,2,</math> and <math>\mathbf\tau_{12} (=\mathbf\tau_{12}(\mathbf{q}))</math> is the (vectorial) NACT responsible for the coupling between <math>u_{1}(\mathbf{q}) </math> and <math>u_{2}(\mathbf{q}) </math>.
| |
| | |
| Next a new function is introduced:<ref>M. Baer and R. Englman, Chem. Phys. Lett. 265, 105 (1997)</ref>
| |
| | |
| : <math> \chi = \psi_1 + i\psi_2 \, </math>
| |
| | |
| and the corresponding rearrangements are made:
| |
| | |
| (i) Multiplying the second equation by <math> i </math> and combining it with the first equation yields the (complex) equation:
| |
| | |
| : <math>-\frac{\hbar^2}{2m} \nabla^{2}\chi + (\tilde{u}_{1}- E)\chi + i\frac{\hbar^2}{2m}[2 \mathbf{\tau}_{12}\nabla + \nabla\mathbf{\tau}_{12}]\chi + i ({u}_{1} - {u}_{2} )\psi_{2} = 0 </math>
| |
| | |
| (ii) The last term in this equation can be deleted for the following reasons: At those points where <math>{u}_{2} (\mathbf{q})</math> is classically closed, <math> \psi_{2}(\mathbf{q}) </math> ~ 0 by definition, and at those points where <math>{u}_{2} (\mathbf{q})</math> becomes classically allowed (which happens at the vicinity of the (1,2) degeneracy points) this implies that: <math>{u}_{1} (\mathbf{q}) </math> ~ <math>{u}_{2} (\mathbf{q}) </math> or <math>{u}_{1} (\mathbf{q}) </math> - <math>{u}_{2} (\mathbf{q}) </math> ~ 0. Consequently the last term is, indeed, negligibly small at every point in the region of interest and the equation simplifies to become:
| |
| | |
| : <math>-\frac{\hbar^2}{2m} \nabla^{2}\chi + (\tilde{u}_{1}- E)\chi + i\frac{\hbar^2}{2m}[2 \mathbf{\tau}_{12}\nabla + \nabla\mathbf{\tau}_{12}]\chi = 0 </math>
| |
| | |
| In order for this equation to yield a solution with the correct symmetry, it is suggested to apply a perturbation approach based on an elastic potential, <math> {u}_{0} (\mathbf{q})</math>, which coincides with <math> {u}_{1} (\mathbf{q})</math> at the asymptotic region.
| |
| | |
| The equation with an elastic potential can be solved, in a straightforward manner, by substitution. Thus, if <math>\chi_{0}</math> is the solution of this equation, it is presented as:
| |
| | |
| : <math>\chi_{0} (\mathbf{q} | \Gamma) = \xi_{0}(\mathbf{q}) exp [-i \int_{\Gamma} d\mathbf{q} \mathbf{'} \cdot \mathbf{\tau} (\mathbf{q} \mathbf{'}|\Gamma)]</math>
| |
| | |
| where <math>\Gamma</math> is an arbitrary contour and the exponential function contains the relevant symmetry as created while moving along <math>\Gamma</math>.
| |
| | |
| The function <math>\xi_{0}(\mathbf{q})</math> can be shown to be a solution of the (unperturbed/elastic) equation:
| |
| | |
| : <math>-\frac{\hbar^2}{2m} \nabla^{2}\xi_{0} + (u_{0}- E) \xi_{0} = 0 </math>
| |
| | |
| Having <math>\chi_{0}(\mathbf{q} | \Gamma)</math>, the full solution of the above decoupled equation takes the form:
| |
| | |
| : <math>\chi(\mathbf{q} | \Gamma) = \chi_{0}(\mathbf{q} | \Gamma) + \eta(\mathbf{q} | \Gamma)</math>
| |
| | |
| where <math>\eta(\mathbf{q} | \Gamma)</math> satisfies the resulting inhomogeneous equation:
| |
| | |
| : <math>-\frac{\hbar^2}{2m} \nabla^{2}\eta + (\tilde{u}_{1}- E)\eta + i\frac{\hbar^2}{2m}[2 \mathbf{\tau}_{12}\nabla + \nabla\mathbf{\tau}_{12}]\eta = ({u}_{1} - {u}_{0})\chi_{0}</math>
| |
| | |
| In this equation the inhomogeneity ensures the symmetry for the perturbed part of the solution along any contour and therefore for the solution in the required region in configuration space.
| |
| | |
| The relevance of the present approach was demonstrated while studying a two-arrangement-channel model (containing one inelastic channel and one reactive channel) for which the two adiabatic states were coupled via a [[Jahn-Teller]] conical intersection.<ref>(a) R. Baer, D.M. Charutz, R. Kosloff and M. Baer, J. Chem. Phys. 111, 9141 (1996); (b) S. Adhikari and G.D. Billing, J. Chem. Phys. 111, 40 (1999)</ref><ref>D.M. Charutz, R. Baer and M. Baer, Chem. Phys. Lett. 265, 629 (1996)</ref> A nice fit between the symmetry-preserved, single-state, treatment and the corresponding two-state treatment was obtained. This applies in particular to the reactive state-to-state probabilities (see Table III in Ref. 5a and Table III in Ref. 5b) for which the ordinary BO approximation led to erroneous results, whereas the symmetry-preserving BO approximation produced the accurate results as they followed from solving the two coupled equations.
| |
| | |
| ==See also==
| |
| * [[Vibronic coupling]]
| |
| * [[Born–Huang approximation]]
| |
| * [[Adiabatic process (quantum mechanics)]]
| |
| * [[Avoided crossing]]
| |
| * [[Diabatic]]
| |
| * [[Franck–Condon principle]]
| |
| * [[Adiabatic ionization]]
| |
| | |
| ==References ==
| |
| <references/>
| |
| | |
| ===Historical note===
| |
| | |
| Born and Oppenheimer wrote the paper: {{cite journal|doi=10.1002/andp.19273892002|title=Zur Quantentheorie der Molekeln|trans_title=On the Quantum Theory of Molecules|language=German| journal=Annalen der Physik|volume=389|issue= 20|pages= 457–484 |year=1927| first1=Max|last1= Born|first2=J. Robert |last2=Oppenheimer|bibcode = 1927AnP...389..457B }} This paper describes the separation of electronic motion, nuclear vibrations, and molecular rotation. However, in the original paper the presentation of the BO approximation is well hidden in Taylor expansions (in terms of internal and external nuclear coordinates) of (i) electronic wave functions, (ii) potential energy surfaces and (iii) nuclear kinetic energy terms. ''Internal coordinates'' are the relative positions of the nuclei in the molecular equilibrium and their displacements (vibrations) from equilibrium. ''External coordinates'' are the position of the center of mass and the orientation of the molecule. The Taylor expansions complicate the theory and make the derivations
| |
| very hard to follow, in contrast to the explanation and derivation above. Moreover, knowing that the proper separation of vibrations and rotations was not achieved in this paper, but only eight years later [{{ cite journal|doi=10.1103/PhysRev.46.383|first=C. |last=Eckart|title=The Kinetic Energy of Polyatomic Molecules|journal= Physical Review|volume= 46| pages=383–387|year=1935|bibcode = 1934PhRv...46..383E }}] (see [[Eckart conditions]]), there is little need to understand the work by Born and Oppenheimer, however famous it may be. Although the article still collects many citations each year, it is safe to say that it is not read anymore (except perhaps by historians of science).
| |
| | |
| ==External links==
| |
| Resources related to the Born–Oppenheimer approximation:
| |
| * [http://gallica.bnf.fr/ark:/12148/bpt6k15386r The original article] (in German)
| |
| * [http://www.ulb.ac.be/cpm/people/scientists/bsutclif/bornopn_corr.pdf Translation by S. M. Blinder]
| |
| * [http://www.tcm.phy.cam.ac.uk/%7Epdh1001/thesis/node13.html The Born–Oppenheimer approximation], a section from Peter Haynes' doctoral thesis
| |
| | |
| {{DEFAULTSORT:Born-Oppenheimer approximation}}
| |
| [[Category:Quantum chemistry]]
| |
| [[Category:Approximations]]
| |
| | |
| [[it:Approssimazione di Born - Oppenheimer]]
| |
| [[ja:ボルン-オッペンハイマー近似]]
| |
| [[zh:玻恩-奥本海默近似]]
| |