The Accuracy of the Pseudopotential Approximation within Density-Functional Theory
|
|
- Phebe Hall
- 5 years ago
- Views:
Transcription
1 D. Porezag et a.: The Accuracy of the Pseudopotentia Approximation 219 phys. stat. so. (b) 217, 219 (2000) Subject cassification: Hx; Mb The Accuracy of the Pseudopotentia Approximation within Density-Functiona Theory D. Porezag (a, b), M. R. Pederson (a), and A. Y. Liu (b) (a) Center for Computationa Materias Science, Nava Research Laboratory, Washington, DC 20375, USA (b) Department of Physics, Georgetown University, Washington, DC 20057, USA (Received August 10, 1999) We have investigated the accuracy of pseudopotentia (PSP) density-functiona cacuations with respect to the corresponding a-eectron (AE) resuts for a variety of atoms and sma moecues. It is found that most of the deviations between ab-initio PSP and AE cacuations are due to the inearization of the exchange-correation functiona within the PSP approach. This probem can be eiminated by appying noninear core corrections (NLCC). We find that a correct description of spin-poarized states requires the NLCC, even for first-row atoms. This is essentia for simuations of magnetic systems and reaction processes which invove radicas. The NLCC is aso essentia for a reaistic description of eements with more ong-range core states such as akai atoms. A further improvement of pseudopotentia accuracy may be achieved by expicity incuding semi-core states in the cacuation. 1. Introduction Recent impementations of eectronic structure methodoogies can be divided into two groups. On one side, there are a-eectron (AE) approaches which consider both core and vaence eectrons expicity in the cacuation. Since core eectrons are strongy ocaized, these methods need to use ocaized basis functions such as augmented pane waves (APW) [1], inearized muffin-tin orbitas (LMTO) [2], Sater-type orbitas (STOs), or gaussian orbitas [3 to 5]. Aternativey, one can empoy pseudopotentias (PSPs) which effectivey project out the core states from the probem whie retaining the physica properties of the vaence region. Extensive discussions and reviews of this approach have appeared in Refs. [6 to 10]. Many pseudopotentia appications use a pane-wave basis set [11 to 13] but a variety of oca-orbita based impementations exist as we [14 to 21]. Currenty, a ot of work is aso aimed at mesh-based approaches [22 to 27]. Nevertheess, many PSP appications use pane-waves. The size of the basis needed for a cacuation of this type depends strongy on the shape (hardness) of the atomic pseudopotentias. Severa approaches have been proposed to construct soft pseudopotentias which minimize the numerica costs whie retaining the accuracy of the approach [28 to 31]. In this work, we use a modification of the scheme proposed by Trouier and Martins [31]. It wi be discussed in Section 2. Athough the aspect of PSP softness is vita for practica appications, the more important probem is pseudopotentia transferabiity, i.e. the abiity of the PSP to dea with eectronic configurations that differ from the one used for its construction. By definition, a PSPs that are created from first principes must reproduce the a-eectron
2 220 D. Porezag et a. eigenvaues of the eectronic reference state. In order to avoid an expicit spin dependence of the PSP, the spin-averaged ground state of the corresponding atom is usuay chosen. However, the actua ground states of most atoms are spin-poarized and the PSP must be transferabe in order to describe these states accuratey. This may not be the case for a PSP with a poor transferabiity. Simiar probems arise when atoms are paced in a moecuar environment and donate or accept eectrons. It is therefore important to systematicay study the resuts of PSP cacuations and compare them with the corresponding a-eectron data. This way, systematic PSP deficiencies can be found, understood, and corrected. Our paper is organized as foows. In Section 2, we wi describe the different eves of pseudopotentia approximations that may be used in practica cacuations and how to impement them. Section 3 focuses on discussing the accuracy of these approximations for cacuations on atoms and sma moecues. The resuts are summarized in Section Method If the density-functiona formaism is used in the usua Kohn-Sham approach [32, 33], the tota energy may be written as: where E ˆ T 0 r r Š dr r r V ext r E H r r Š E xc r r Š ; E H r r Š ˆ 1 dr r r V H r ; 2 r r ˆ Pocc jy i r j 2 ; T 0 r r Š ˆ Pocc i V H r ˆ i dr 0 r r 0 jr r 0 j ; 1 i r2 dr Y r 2 Y i r : 2 Note that atomic units for energy (1 a.u. ˆ 27:211 ev) and ength (1 a B = A) have been used in the above equations. V ext r is the externa potentia, for an atomic cacuation V ext r ˆ Z=jrj. T 0 r r Š is the kinetic energy of an interaction-free eectron gas which may be cacuated easiy using the Kohn-Sham orbitas Y i r. E xc r r Š is the exchange-correation energy which is a noninear functiona of the density r r. The exact expression for E xc is unknown but there are reasonabe approximations to it. In this paper, we have used both oca-density and the Perdew-Burke-Ernzerhof generaized-gradient approximations [34]. A generaized tota energy expression for spin-poarized systems may be obtained by aowing different Kohn-Sham orbitas for up and down eectrons and summing over both spin systems. Minimization of the expression (1) eads to the Kohn-Sham equations r2 2 V ext r V H r r Š V xc r r Š Y i r ˆ e i Y i r ; V xc r r Š ˆ de xc r r Š ; V eff r ˆ V ext r V H r r Š V xc r r Š : 3 dr r In an a-eectron cacuation, both core and vaence states are contained in the set of Kohn-Sham orbitas Y i used in the cacuation. Since core orbitas usuay change very itte if the atom is put in a different environment, it is a good approximation to assume
3 The Accuracy of the Pseudopotentia Approximation 221 that they are fixed for a certain atom type and that their effect may be simuated by a modified potentia (pseudopotentia). Experience has shown that a reiabe pseudopotentia shoud fufi the foowing conditions: 1. The (node-ess) owest-energy eigenfunction (pseudo-wavefunction) of the pseudopotentia has the same eigenvaue as the a-eectron vaence wave function. 2. The pseudo-wavefunction is identica to the a-eectron wave function outside a chosen cut-off radius r cut (norm conservation). There are ways to reax the second condition, for a discussion of these methods, see e.g. Ref. [30]. Since the a-eectron vaence functions Y for different anguar momenta (s, p,...) are different, the corresponding pseudopotentias V ps wi be -dependent. Therefore, they are caed semioca. There is a number of schemes eading to pseudopotentias which obey the two conditions stated above. Here, we wi focus on the approach suggested by Trouier and Martins [31]. Within this method, the pseudo-wavefunction Y ps is written as 8 >< Y ps r ˆ exp m 2 P c n r 2n for r r cut ; nˆ0 4 >: Y r for r r cut : The c n are determined from the foowing conditions: d n dr n Y ps r cut ˆ dn dr n Y r cut 8n ˆ 0; m match first m derivatives ; r cut 0 r 2 Y ps r 2 Y r 2 ˆ 0 norm conservation ; d 2 dr 2 Vps=scr 0 ˆ 0 smoothness condition ; 5 where V ps=scr r is the screened pseudopotentia repacing the effective a-eectron potentia in Eq. (3) and reated to Y ps r by V ps=scr r ˆ e r2 Y ps r 2Y ps r (inverted Kohn-Sham equation). In the origina Trouier-Martins approach, ony the first four derivatives of Y ps r were matched at r ˆ r cut. In order to improve the PSP smoothness in this area, we match the first six derivatives. As a resut, the eigenvaue of the a-eectron vaence wave function is usuay reproduced by the pseudo-wavefunction within 1 mev. As an exampe, Fig. 1 shows pseudo and a-eectron vaence wave functions for the suphur atom. As has been noted earier, the screened PSP V ps=scr r in a pseudopotentia cacuation is equivaent to the effective potentia V eff r in an a-eectron cacuation, r2 2 Vps=scr r Y ps r ˆ e Y ps r : 7 However, what is needed for appications of the formaism to arbitrary systems is the unscreened PSP V ps=unscr r which is equivaent to the externa potentia in an a-eectron cacuation. In the traditiona PSP approach, the pseudo-density r ps r takes the 6
4 222 D. Porezag et a. Fig. 1. Pseudo and a-eectron vaence wave functions for suphur pace of the tota density r r in an a-eectron cacuation. Therefore, the PSP Kohn- Sham equation reads r2 2 Vps=unscr r V H r ps r Š V xc r ps r Š Y ps r ˆ e Y ps r : 8 and hence V ps=unscr r ˆ V ps=scr r V H r ps r Š V xc r ps r Š : 9 However, this approach negects that in contrast to the Couomb potentia V H which is inear in the charge density, E xc and V xc are noninear functionas. Consequenty, core and vaence contributions to the exchange-correation energy cannot be separated in a simpe manner. For this reason, Louie et a. [35] introduced a scheme which takes the noninear nature of the E xc functiona into account by expicity adding a core charge density, i.e. repacing r ps r by r ps r r c r. Hence, the unscreening procedure is modified V ps=unscr r ˆ V ps=scr r V H r ps r Š V xc r ps r r c r Š : 10 This correction to the standard PSP approach is commony known as noninear core correction (NLCC). If r c is chosen as r c r ˆ r r r ps r ˆ r c fu r (fu NLCC), the corresponding PSP shoud give the same resuts as an a-eectron cacuation, provided that core reaxation effects are negigibe. However, such an approach is not feasibe for pane-wave cacuations because of the strongy ocaized nature of the true core charge density. For this reason, Louie et a. [35] suggested a partia NLCC characterized by a pseudo core density that is equivaent to the true core charge density for a radius r r core and smooth for smaer radii. Unfortunatey, the origina ansatz for r c r of Ref. [35] does not ead to a smooth transition at r ˆ r core. Therefore, we have deveoped a different procedure to determine this quantity. We set 8 < r c r ˆ exp m 1 P d n r 2n for r r core ; 11 : nˆ0 r c fu r for r r core ;
5 The Accuracy of the Pseudopotentia Approximation 223 and the d n are determined based on the foowing conditions: d n dr n rc r core ˆ dn dr n rc fu r core 8n ˆ 0; m match first m derivatives ; r core! 2 d 2 dr 2 rc r! min smoothness condition ; 12 0 Figure 2 iustrates the PSP unscreening without the NLCC and with a partia NLCC for the suphur atom. Note that the unscreened potentias depend on the magnitude of the NLCC. After the unscreened PSP has been defined, it can be used in eectronic structure cacuations. A resuts presented here have been obtained using the gaussian-orbita based density functiona code NRLMOL deveoped by Pederson, Porezag, and Jackson [4, 5, 20]. This method is based on an accurate numerica integration scheme. Besides the expansion of the wave functions into gaussian orbitas, no additiona approximations are made. The program aows for mixed a-eectron and pseudopotentia descriptions within the same cacuation, i.e. some atoms may be treated within an a-eectron Fig. 2. Upper part: vaence density and core densities with fu and partia NLCC for suphur. Lower part: semioca PSPs without and with partia NLCC for suphur (a.u. atomic units)
6 224 D. Porezag et a. scheme whie others may be described with PSPs. We have empoyed very arge basis sets and checked that a resuts are converged with respect to basis set size. The generaized-gradient (GGA-PBE) and oca-density (LDA) approximations [34] have been used to describe exchange-correation effects. 3. Resuts and Discussion Before presenting our resuts for moecues, it is worthwhie to ook at the cacuated properties of free atoms. In fact, whie the PSP construction is usuay performed based on spin-unpoarized eectronic states, the actua atomic ground state is often characterized by a non-vanishing tota spin. It is therefore important to compare the reative energies of different eectronic configurations and compare them to the corresponding a-eectron atoms. One of the most important figures is the energy difference between the spin-averaged eectronic state used in the PSP construction and the (frequenty spin-poarized) atomic ground state. This quantity wi be referred to as the atomic spin poarization energy (SPE) in the foowing. Tabe 1 summarizes the cut-off and NLCC core radii used to construct the PSPs for the eements discussed here. The spin-unpoarized eectronic ground state of the corresponding atom was used as the PSP construction reference state. 3.1 Atoms Figure 3 shows the cacuated differences between PSP and a-eectron SPEs with and without the NLCC for the first-row atoms B through Ne. These atoms are characterized by a singe, strongy ocaized 1s core she. Therefore, it is usuay assumed that they can be treated without the NLCC. However, since the vaence states in these systems are aso fairy compact, there is sti a significant overap between vaence and core charge densities. As a resut, the PSP based spin poarization energies show significant errors as compared to the corresponding a-eectron vaues if noninear core corrections are ignored. The deviations are argest for N which has a fuy poarized 2p she and a tota SPE of 3.02 ev within LDA and 3.12 ev within GGA-PBE. If the NLCC is not appied, the SPE errors amount to about 10% of the tota SPE. The deviations can be as arge as 0.3 ev (for N and O within GGA-PBE). Errors of this magnitude wi ead to discrepancies of approximatey 0.6 ev if binding energies for homo-nucear dimers are cacuated. For exampe, the AE GGA-PBE dissociation energy for O 2 is 6.2 ev whereas it is determined to be 5.7 ev in the PSP approach without the NLCC. Ta b e 1 Cut-off and core radii [in a B ] used for the pseudopotentia construction. Mn c corresponds to a conservative Mn pseudopotentia which incudes the 3s and 3p semi-core states in the vaence basis H Li B C N O F S Mn Mn c s p d NLCC
7 The Accuracy of the Pseudopotentia Approximation 225 Fig. 3. Errors of the atomic spin-poarization energies with respect to a-eectron vaues as cacuated with fu NLCC and without NLCC Incuding the NLCC greaty improves the performance of the PSP approach, eading to a maximum deviation of ony 0.02 ev (for the SPE of N). Hence, it can be concuded that SPEs of first-row atoms are very we described if a fu NLCC is appied. The question which remains to be answered is how the SPE error depends on the choice of the core radius if ony a partia NLCC is appied. Tabe 2 shows resuts for a GGA-PBE cacuation on oxygen (the behavior is very simiar for the other first-row atoms as we). For core radii of about 0.4a B or smaer, the spin poarization energies determined with a partia NLCC do not differ from the fu NLCC resuts. Ta b e 2 Atomic GGA-PBE spin poarization energies [in ev] for oxygen as cacuated with aeectron and PSP methods using NLCC and different core radii r core (first ine, given in a B ). r core ˆ 1 corresponds to a negect of NLCC AE SPE error
8 226 D. Porezag et a. Ta b e 3 Atomic LDA spin poarization energies [in ev] for Mn as cacuated with a-eectron and PSP methods using different core radii r core (first ine, given in a B ). r core ˆ 1 corresponds to a negect of NLCC AE The deficiencies of the PSP approach without the NLCC become even more apparent for transition meta atoms. These systems have strongy ocaized d eectrons. In fact, the d charge density in the 3d transition metas overaps substantiay with the 3s and 3p core density. In the foowing we consider the Mn atom which has a fuy spinpoarized 3d she. Tabe 3 shows resuts for the LDA spin poarization energies as cacuated within the a-eectron and PSP approaches using different core radii for the NLCC. It is apparent that the PSP method without the NLCC gives competey unreasonabe resuts. Further, the core radius must be chosen fairy sma (0.43a B ) in order to achieve converged resuts. However, if the NLCC is appied correcty, it eads to a very good agreement with the a-eectron cacuation. Another interesting and important property is the s±d transfer energy (i.e. the tota energy differences between the atomic s 2 d n and s 1 d n 1 eectronic configurations). We have cacuated this quantity and give the resuts in Tabe 4. The a-eectron and converged PSP + NLCC s±d transfer energies are found at 1.04 and 1.20 ev, respectivey. Note that the error increases significanty if r core is chosen too arge. Further, athough the NLCC eads to a much improved description of the s±d eectron transfer, there is sti a difference of 0.16 ev between PSP and a-eectron resuts which needs further discussion. The cacuations presented so far were based on PSPs which assigned the 3s and 3p states to the core. However, when we examine the spatia range of these states, we find that they strongy overap with the spatia range of the the 3d eectrons. Transferring an eectron from the rather diffuse 4s state to the ocaized 3d state wi increase the Couomb repusion in this area and hence ead to a reaxation of the 3s and 3p states. Assigning these states to the core makes such a reaxation impossibe. For this reason, we have constructed a `conservative' pseudopotentia for Mn which incudes the 3s and 3p states in the vaence basis. As can be seen from Tabe 4, the conservative approach eads to resuts which are in exceent agreement with the a-eectron data. Therefore, it can be concuded that the PSP approach can describe transition meta atoms with a simiar accuracy as other atom types if both noninear core±vaence interactions and semi-core reaxation effects are incuded in the method. Ta b e 4 Atomic LDA s±d transfer energies [in ev] for Mn as cacuated with a-eectron and PSP methods using different core radii r core (first ine, given in a B ). r core ˆ 1 corresponds to a negect of NLCC. The Ar core ine refers to a cacuation where the 3sp eectrons are treated as core eectrons AE Ar core Ne core
9 The Accuracy of the Pseudopotentia Approximation 227 Fig. 4. LDA cohesive energy for the LiF moecue as a function of internucear distance. The PSP+NLCC resuts for r core ˆ 1.62a B are identica to the fu NLCC 3.2 LiF diatomic moecue The preceding section was aimed at discussing energy differences induced by simpy changing the occupation numbers of the atomic one-eectron orbitas. However, in many practica appications, the most important change in the atomic eectron density occurs when the atom is incorporated into a soid, custer, or moecue. The argest reative changes wi occur for atoms with a very diffuse vaence density if they are paired with atoms characterized by compact ocaized states. Prominent exampes for such systems are akai-haogenide compounds such as LiF. In order to demonstrate the probems which occur for this combination of eements, we focus on the LiF diatomic moecue. Figure 4 shows the cohesive energy of this system as a function of interatomic distance as cacuated within LDA. The PSP cacuation without the NLCC overestimates the binding energy by about 1 ev and underestimates the bond ength by more than 0.5a B. Appying a fu NLCC corrects this deficiency and eads to a good agreement with the a-eectron resut. Partia and fu NLCC resuts agree for r core 1:62a B. The remaining difference between PSP + NLCC and AE data can be attributed to the negect of the 1s core reaxation. Note that in order to rue out additiona sources of error, the F atom has aso been treated in an a-eectron fashion here. 3.3 Structura and vibrationa properties of sma moecues The preceding sections have estabished that it is often necessary to appy the NLCC in order to obtain reiabe resuts. Whie there are aso many cases where the NLCC has negigiby sma effects, we advocate to aways use this correction since it usuay does not increase the numerica effort of a cacuation very much, especiay not in cacuations with gaussian-orbita basis sets. Consequenty, the resuts presented in the foow-
10 228 D. Porezag et a. Ta b e 5 Bond engths [in a B ] and anges for some moecues as cacuated with the a-eectron (AE) and PSP NLCC approaches. The GGA-PBE exchange-correatioin functiona has been used AE PSP AE PSP H 2 r(hh) LiH r(lih) Li 2 r(lili) LiF r(lif) N 2 r(nn) H 2 O r(oh) O 2 r(oo) (HOH) F 2 r(ff) CH 4 r(ch) CO r(co) SO 2 r(so) HF r(hf) (OSO) ing are based on cacuations with pseudopotentias that empoy the NLCC and use the cut-off and core radii summarized in Tabe 1. For a more detaied discussion on the proper choice of the NLCC core radii pease refer to Ref. [36]. Tabe 5 ists the cacuated bond engths and anges for tweve moecues as cacuated with the a-eectron and pseudopotentia approaches and the GGA-PBE exchange-correation functiona. PSP and AE resuts are in exceent agreement with a root-mean-square (RMS) error of ony 0.003a B for the bond engths. The argest errors occur for LiH and LiF which can be attributed to the negect of the Li core reaxation effects. Tabe 6 contains information about cacuated atomization energies. Again, AE and PSP data are in very good agreement, the deviations do not exceed 0.04 ev per atom. The current PSP resuts show much smaer errors than the ones pubished in an earier paper [20] which were based on BHS pseudopotentias without noninear core corrections. This is especiay true for the moecues containing Li where the previous PSP cacuations gave binding energies with errors of severa ev as compared to the corresponding a-eectron vaues. This iustrates that the NLCC is especiay important for akai atoms. Note that a atomization energies are consistenty underestimated. This is most ikey due to the negect of core reaxation effects for the moecuar species. Tabe 7 summarizes the vibrationa properties of the moecues investigated here. Frequencies, IR intensities, and Raman scattering activities have been determined with a method pubished in Ref. [37]. The PSP based vibrationa frequencies are in exceent agreement with the AE cacuations. The two methodoogies usuay differ by ony a few wavenumbers, or at most 1%. However, the behavior of the IR intensities and Ra- Ta b e 6 Atomization energies [in ev] for some moecues as cacuated with the a-eectron (AE) and PSP NLCC approaches. The GGA-PBE exchange-correation functiona has been used. Zero-point motion has been incuded in the vaues presented here AE PSP D AE PSP D H HF Li LiH N LiF O H 2 O F CH CO SO
11 The Accuracy of the Pseudopotentia Approximation 229 Ta b e 7 Harmonic vibrationa frequencies [in cm 1 ], IR intensities [in (D/A) 2 amu 1 ] and Raman activities [in A 4 amu 1 ] for some moecues as cacuated with the a-eectron (AE) and PSP + NLCC approaches. The GGA-PBE exchange-correation functiona has been used AE PSP freq. IR Ram. freq. IR Ram. H 2 S g Li 2 S g N 2 S g O 2 S g F 2 S g CO S HF S LiH S LiF S H 2 O A B CH 4 T E A SO 2 A B man activities (which depends strongy on subte changes of the moecuar charge distribution as a function of the nucear coordinates) is different. Whie many IR and Raman strengths agree very we (within 5%) there are aso cases with significanty arger deviations. For exampe, the spectra intensities of a moecues containing Li show errors of about 10%, indicating that the effect of the Li core eectrons cannot be competey negected. However, even the argest errors do not exceed the typica differences between theory and experiment which were found to be 10 to 30% in Ref. [37]. 4. Summary and Concusions We have investigated the accuracy of the pseudopotentia approximation as compared to the a-eectron approach for a variety of sma moecues using norm-conserving Trouier-Martins PSPs and a gaussian-orbita based numerica scheme. It is found that noninear core corrections are often necessary to obtain reiabe resuts, therefore, we suggest to aways use them in appications of the method to arger systems. We have shown that the PSP approach with NLCC predicts equiibrium structures, energies, and vibrationa frequencies with a very high accuracy. Resuts for IR and Raman intensities show arger errors but are sti sufficienty accurate for most practica appications. Remaining deviations between the PSP and AE methods can be expained by core reaxation effects. Consequenty, the performance of the PSP approximation can be substantiay improved by using conservative pseudopotentias, i.e. by incuding semi-core states expicity in the cacuation. In summary, it can be concuded that the combined PSP + NLCC approach is very we suited for the cass of systems investigated here.
12 230 D. Porezag et a.: The Accuracy of the Pseudopotentia Approximation Acknowedgements DVP greatfuy acknowedges support from the Aexander-von- Humbodt foundation of Germany. This work was supported in part by the ONR Moecuar Design Institute N WX20709 and by Nationa Science Foundation Grant No. DMR References [1] E. Wimmer, H. Krakauer, M. Weinert, and A. J. Freeman, Phys. Rev. B 24, 864 (1981). [2] O. K. Andersen, Phys. Rev. B 12, 3060 (1975). [3] Gaussian 94, M. J. Frisch, G. W. Trucks, H. B. Schege, P. M. W. Gi, B. G. Johnson, M. A. Robb, J. R. Cheeseman, T. Keith, G. A. Petersson, J. A. Montgomery, K. Raghavachari, M. A. A-Laham, V. G. Zakrzewski, J. V. Ortiz, J. B. Foresman, J. Ciosowski, B. B. Stefanov, A. Nanayakkara, M. Chaacombe, C. Y. Peng, P. Y. Ayaa, W. Chen, M. W. Wong, J. L. Andres, E. S. Repoge, R. Gomperts, R. L. Martin, D. J. Fox, J. S. Binkey, D. J. Defrees, J. Baker, J. P. Stewart, M. Head-Gordon, C. Gonzaez, and J. A. Pope, Gaussian Inc., Pittsburgh PA, [4] M. R. Pederson and K. A. Jackson, Phys. Rev. B 41, 7453 (1990). [5] K. A. Jackson and M. R. Pederson, Phys. Rev. B 42, 3276 (1990). [6] M. L. Cohen and V. Heine, Soid State Phys. 24, 37 (1970). [7] W. E. Pickett, Comp. Phys. Rep. 9, 116 (1989). [8] J. R. Cheikowsky and M. L. Cohen, in: Handbook on Semiconductors, Vo. 1, Ed. P. T. Landsberg, Esevier Pub. Co., Amsterdam 1992 (p. 59). [9] M. C. Payne, M. P. Teter, D. C. Aan, T. A. Arias, and J. D. Joannopouos, Rev. Mod. Phys. 64, 1045 (1992). [10] S. Goedecker and K. Maschke, Phys. Rev. B 45, 88 (1992). [11] R. Car and M. Parrineo, Phys. Rev. Lett. 55, 2471 (1985). [12] R. M. Wentzcovich and J. L. Martins, Soid State Commun. 78, 831 (1991). [13] M. Bockstedte, A. Key, J. Neugebauer, and M. Scheffer, Comput. Phys. Commun. 107, 187 (1997). [14] J. Bernhoc and S. T. Panteides, Phys. Rev. B 18, 1780 (1978). [15] D. Vanderbit and S. G. Louie, Phys. Rev. B 30, 6118 (1984). [16] J. R. Cheikowsky and S. G. Louie, Phys. Rev. B 29, 3470 (1984). [17] P. J. Feibeman, Phys. Rev. B 35, 2626 (1987). [18] O. F. Sankey, D. A. Drabod, and G. B. Adams, Bu. Amer. Phys. Soc. 36, 924 (1991). [19] R. Jones and A. Sayyash, J. Phys. C 19, L653 (1986). [20] A. Briey, M. R. Pederson, K. A. Jackson, D. C. Patton, and D. V. Porezag, Phys. Rev. B 58, 1786 (1998). [21] D. Sanchez Porta, P. Ordejon, E. Artacho, and J. M. Soer, Internat. J. Quantum Chem. 65, 453 (1997). [22] F. Gygi, Europhys. Lett. 19, 617 (1992). [23] D. R. Hamann, Phys. Rev. B 51, 7337 (1995). [24] A. Devenyi, K. Cho, T. A. Arias, and J. D. Joannopouos, Phys. Rev. B 49, (1994). [25] N. A. Modine, G. Zumbach, and E. Kaxiras, Phys. Rev. B 55, (1997). [26] J. R. Cheikowsky, N. Trouier, and K. Wu, Phys. Rev. B 16, (1994). [27] E. L. Briggs, D. J. Suivan, and J. Bernhoc, Phys. Rev B 52, R5471 (1995). [28] D. R. Hamann, M. SchuÈter, and C. Chiang, Phys. Rev. Lett. 43, 1494 (1979). [29] A. M. Rappe, K. M. Rabe, E. Kaxiras, and J. D. Joannopouos, Phys. Rev. B 41, 1227 (1990). [30] D. Vanderbit, Phys. Rev. B 41, R7892 (1990). [31] N. Trouier and J. L. Martins, Phys. Rev. B 43, 1993 (1991). [32] P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). [33] W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). [34] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). [35] S. G. Louie, S. Froyen, and M. L. Cohen, Phys. Rev. B 26, 1738 (1982). [36] D. Porezag, M. R. Pederson, and A. Y. Liu, Phys. Rev. B 60, (1999). [37] D. Porezag and M. R. Pederson, Phys. Rev. B 54, 7830 (1996).
Lecture VIII : The pseudopotential
Lecture VIII : The pseudopotentia I. KOHN-SHAM PROBLEM FOR AN ISOLATED ATOM For a one-eectron atom, the Couombic potentia, V ( r) = V (r) = Z/r is sphericay symmetric. The soutions may then be spit into
More informationIntroduction to LMTO method
1 Introduction to MTO method 24 February 2011; V172 P.Ravindran, FME-course on Ab initio Modeing of soar ce Materias 24 February 2011 Introduction to MTO method Ab initio Eectronic Structure Cacuations
More informationIntroduction to DFT and Density Functionals. by Michel Côté Université de Montréal Département de physique
Introduction to DFT and Density Functionas by Miche Côté Université de Montréa Département de physique Eampes Carbazoe moecue Inside of diamant Réf: Jean-François Brière http://www.phys.umontrea.ca/ ~miche_cote/images_scientifiques/
More informationarxiv: v1 [cond-mat.mtrl-sci] 10 Mar 2015
Constructing optima oca pseudopotentias from first principes arxiv:153.2992v1 [cond-mat.mtr-sci] 1 Mar 215 Wenhui Mi, 1,2 Shoutao Zhang, 1,2 Yanming Ma, 1,2, and Maosheng Miao 3,2, 1 State Key Laboratory
More informationElectron-impact ionization of diatomic molecules using a configuration-average distorted-wave method
PHYSICAL REVIEW A 76, 12714 27 Eectron-impact ionization of diatomic moecues using a configuration-average distorted-wave method M. S. Pindzoa and F. Robicheaux Department of Physics, Auburn University,
More informationJost Function for Singular Potentials
Jost Function for Singuar Potentias S. A. Sofianos, S. A. Rakityansky, and S. E. Massen Physics Department, University of South Africa, P.O.Box 392, Pretoria 0003, South Africa (January 2, 999) An exact
More informationSection 6: Magnetostatics
agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More informationIn-plane shear stiffness of bare steel deck through shell finite element models. G. Bian, B.W. Schafer. June 2017
In-pane shear stiffness of bare stee deck through she finite eement modes G. Bian, B.W. Schafer June 7 COLD-FORMED STEEL RESEARCH CONSORTIUM REPORT SERIES CFSRC R-7- SDII Stee Diaphragm Innovation Initiative
More information$, (2.1) n="# #. (2.2)
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationNonlinear Analysis of Spatial Trusses
Noninear Anaysis of Spatia Trusses João Barrigó October 14 Abstract The present work addresses the noninear behavior of space trusses A formuation for geometrica noninear anaysis is presented, which incudes
More informationCharge Density from X-ray Diffraction. Methodology
Charge Density from X-ray Diffraction. Methodoogy Ignasi Mata imata@icmab.es Master on Crystaography and Crystaization, 2012 Outine I. Charge density in crystas II. The mutipoar refinement III. Methodoogy
More informationMore Scattering: the Partial Wave Expansion
More Scattering: the Partia Wave Expansion Michae Fower /7/8 Pane Waves and Partia Waves We are considering the soution to Schrödinger s equation for scattering of an incoming pane wave in the z-direction
More informationLecture 6: Moderately Large Deflection Theory of Beams
Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey
More informationCombining reaction kinetics to the multi-phase Gibbs energy calculation
7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationInternational Journal of Mass Spectrometry
Internationa Journa of Mass Spectrometry 280 (2009) 179 183 Contents ists avaiabe at ScienceDirect Internationa Journa of Mass Spectrometry journa homepage: www.esevier.com/ocate/ijms Stark mixing by ion-rydberg
More informationAn Algorithm for Pruning Redundant Modules in Min-Max Modular Network
An Agorithm for Pruning Redundant Modues in Min-Max Moduar Network Hui-Cheng Lian and Bao-Liang Lu Department of Computer Science and Engineering, Shanghai Jiao Tong University 1954 Hua Shan Rd., Shanghai
More informationSimplified analysis of EXAFS data and determination of bond lengths
Indian Journa of Pure & Appied Physics Vo. 49, January 0, pp. 5-9 Simpified anaysis of EXAFS data and determination of bond engths A Mishra, N Parsai & B D Shrivastava * Schoo of Physics, Devi Ahiya University,
More informationSupporting Information for Suppressing Klein tunneling in graphene using a one-dimensional array of localized scatterers
Supporting Information for Suppressing Kein tunneing in graphene using a one-dimensiona array of ocaized scatterers Jamie D Was, and Danie Hadad Department of Chemistry, University of Miami, Cora Gabes,
More informationNonperturbative Shell Correction to the Bethe Bloch Formula for the Energy Losses of Fast Charged Particles
ISSN 002-3640, JETP Letters, 20, Vo. 94, No., pp. 5. Peiades Pubishing, Inc., 20. Origina Russian Text V.I. Matveev, D.N. Makarov, 20, pubished in Pis ma v Zhurna Eksperimenta noi i Teoreticheskoi Fiziki,
More informationPrediction of high zt in thermoelectric silicon nanowires with axial germanium heterostructures
June 2 EPL, 94 (2) 67 doi:.29/295-575/94/67 www.epjourna.org Prediction of high zt in thermoeectric siicon nanowires with axia germanium heterostructures M. Sheey and A. A. Mostofi (a) The Thomas Young
More informationI. INTRODUCTION. Periodic boundary conditions in ab initio calculations
PHYSCAL REVEW 8 VOLUME 51, NUMBER 7 15 FEBRUARY 1995- Periodic boundary conditions in ab initio cacuations G. Makov and M. C. Payne Cavendish Laboratory, Madingey Road, Cambridge CB3 OHE, United Kingdom
More informationCopyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU
Copyright information to be inserted by the Pubishers Unspitting BGK-type Schemes for the Shaow Water Equations KUN XU Mathematics Department, Hong Kong University of Science and Technoogy, Cear Water
More informationNuclear Size and Density
Nucear Size and Density How does the imited range of the nucear force affect the size and density of the nucei? Assume a Vecro ba mode, each having radius r, voume V = 4/3π r 3. Then the voume of the entire
More informationMidterm 2 Review. Drew Rollins
Midterm 2 Review Drew Roins 1 Centra Potentias and Spherica Coordinates 1.1 separation of variabes Soving centra force probems in physics (physica systems described by two objects with a force between
More informationA Brief Introduction to Markov Chains and Hidden Markov Models
A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,
More informationFRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)
1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using
More informationProblem set 6 The Perron Frobenius theorem.
Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator
More informationhole h vs. e configurations: l l for N > 2 l + 1 J = H as example of localization, delocalization, tunneling ikx k
Infinite 1-D Lattice CTDL, pages 1156-1168 37-1 LAST TIME: ( ) ( ) + N + 1 N hoe h vs. e configurations: for N > + 1 e rij unchanged ζ( NLS) ζ( NLS) [ ζn unchanged ] Hund s 3rd Rue (Lowest L - S term of
More informationSEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l
Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed
More informationA. Distribution of the test statistic
A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Eectrostatic II Notes: Most of the materia presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartoo, Chap... Mathematica Considerations.. The Fourier series and the Fourier
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationA Novel Learning Method for Elman Neural Network Using Local Search
Neura Information Processing Letters and Reviews Vo. 11, No. 8, August 2007 LETTER A Nove Learning Method for Eman Neura Networ Using Loca Search Facuty of Engineering, Toyama University, Gofuu 3190 Toyama
More informationEffects of energy loss on interaction dynamics of energetic electrons with plasmas. C. K. Li and R. D. Petrasso. 1 November 2008
PSFC/JA-8-3 ffects of energy oss on interaction dynamics of energetic ctrons with pasmas C. K. Li and R. D. Petrasso November 8 Pasma Science and Fusion Center Massachusetts Institute of Technoogy Cambridge,
More informationOn the Computation of (2-2) Three- Center Slater-Type Orbital Integrals of 1/r 12 Using Fourier-Transform-Based Analytical Formulas
On the Computation of (2-2) Three- Center Sater-Type Orbita Integras of /r 2 Using Fourier-Transform-Based Anaytica Formuas DANKO ANTOLOVIC, HARRIS J. SILVERSTONE 2 Pervasive Technoogy Labs, Indiana University,
More informationLecture contents. NNSE 618 Lecture #11
Lecture contents Couped osciators Dispersion reationship Acoustica and optica attice vibrations Acoustica and optica phonons Phonon statistics Acoustica phonon scattering NNSE 68 Lecture # Few concepts
More informationModel Calculation of n + 6 Li Reactions Below 20 MeV
Commun. Theor. Phys. (Beijing, China) 36 (2001) pp. 437 442 c Internationa Academic Pubishers Vo. 36, No. 4, October 15, 2001 Mode Cacuation of n + 6 Li Reactions Beow 20 MeV ZHANG Jing-Shang and HAN Yin-Lu
More informationVIII. Addition of Angular Momenta
VIII Addition of Anguar Momenta a Couped and Uncouped Bae When deaing with two different ource of anguar momentum, Ĵ and Ĵ, there are two obviou bae that one might chooe to work in The firt i caed the
More informationOn Complex Preconditioning for the solution of the Kohn-Sham Equation for Molecular Transport
On Compex Preconditioning for the soution of the Kohn-Sham Equation for Moecuar Transport Danie Osei-Kuffuor Lingzhu Kong Yousef Saad James R. Cheikowsky December 23, 2007 Abstract This paper anayzes the
More informationSupplemental Information
Suppementa Information A Singe-eve Tunne Mode to Account for Eectrica Transport through Singe Moecue- Sef-Assembed Monoayer-based Junctions by A. R. Garrigues,. Yuan,. Wang, E. R. Muccioo, D. Thompson,
More informationAb initio full charge-density study of the atomic volume of -phase Fr, Ra, Ac, Th, Pa, U, Np, and Pu
PHYSICA EVIEW B VOUME 55, NUMBE 8 15 FEBUAY 1997-II Ab initio fu charge-density study of the atomic voume of -phase Fr, a, Ac, Th, Pa, U, Np, and Pu. Vitos and J. Koár esearch Institute for Soid State
More informationModal analysis of a multi-blade system undergoing rotational motion
Journa of Mechanica Science and Technoogy 3 (9) 5~58 Journa of Mechanica Science and Technoogy www.springerin.com/content/738-494x DOI.7/s6-9-43-3 Moda anaysis of a muti-bade system undergoing rotationa
More informationTHE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS
ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Sevie, Spain, -6 June 04 THE OUT-OF-PLANE BEHAVIOUR OF SPREAD-TOW FABRICS M. Wysocki a,b*, M. Szpieg a, P. Heström a and F. Ohsson c a Swerea SICOMP
More informationarxiv:quant-ph/ v3 6 Jan 1995
arxiv:quant-ph/9501001v3 6 Jan 1995 Critique of proposed imit to space time measurement, based on Wigner s cocks and mirrors L. Diósi and B. Lukács KFKI Research Institute for Partice and Nucear Physics
More informationarxiv: v2 [physics.atom-ph] 3 Aug 2014
Positron scattering and annihiation on nobe gas atoms D. G. Green, J. A. Ludow, and G. F. Gribakin Department of Appied Mathematics and Theoretica Physics, Queen s University Befast, Befast BT7 1NN, Northern
More informationA SIMPLIFIED DESIGN OF MULTIDIMENSIONAL TRANSFER FUNCTION MODELS
A SIPLIFIED DESIGN OF ULTIDIENSIONAL TRANSFER FUNCTION ODELS Stefan Petrausch, Rudof Rabenstein utimedia Communications and Signa Procesg, University of Erangen-Nuremberg, Cauerstr. 7, 958 Erangen, GERANY
More informationPhysics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions
Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p
More informationPublished in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics
Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of
More informationChemical Kinetics Part 2
Integrated Rate Laws Chemica Kinetics Part 2 The rate aw we have discussed thus far is the differentia rate aw. Let us consider the very simpe reaction: a A à products The differentia rate reates the rate
More informationIn Coulomb gauge, the vector potential is then given by
Physics 505 Fa 007 Homework Assignment #8 Soutions Textbook probems: Ch. 5: 5.13, 5.14, 5.15, 5.16 5.13 A sphere of raius a carries a uniform surface-charge istribution σ. The sphere is rotate about a
More informationProceedings of Meetings on Acoustics
Proceedings of Meetings on Acoustics Voume 9, 23 http://acousticasociety.org/ ICA 23 Montrea Montrea, Canada 2-7 June 23 Architectura Acoustics Session 4pAAa: Room Acoustics Computer Simuation II 4pAAa9.
More informationPartial permutation decoding for MacDonald codes
Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics
More informationPhysics 235 Chapter 8. Chapter 8 Central-Force Motion
Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion
More informationFerromagnetic Coupling of [Ni(dmit) 2 ] - Anions in. (m-fluoroanilinium)(dicyclohexano[18]crown-6)[ni(dmit) 2 ]
Supporting Information Ferromagnetic Coupling of [Ni(dmit) 2 ] - Anions in (m-fluoroanilinium)(dicyclohexano[18]crown-6)[ni(dmit) 2 ] Tomoyuki Akutagawa, *,, Daisuke Sato, Qiong Ye, Shin-ichiro Noro,,
More informationJoel Broida UCSD Fall 2009 Phys 130B QM II. Homework Set 2 DO ALL WORK BY HAND IN OTHER WORDS, DON T USE MATHEMAT- ICA OR ANY CALCULATORS.
Joe Broida UCSD Fa 009 Phys 30B QM II Homework Set DO ALL WORK BY HAND IN OTHER WORDS, DON T USE MATHEMAT- ICA OR ANY CALCULATORS. You may need to use one or more of these: Y 0 0 = 4π Y 0 = 3 4π cos Y
More informationIntegrating Factor Methods as Exponential Integrators
Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been
More informationA MODEL FOR ESTIMATING THE LATERAL OVERLAP PROBABILITY OF AIRCRAFT WITH RNP ALERTING CAPABILITY IN PARALLEL RNAV ROUTES
6 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES A MODEL FOR ESTIMATING THE LATERAL OVERLAP PROBABILITY OF AIRCRAFT WITH RNP ALERTING CAPABILITY IN PARALLEL RNAV ROUTES Sakae NAGAOKA* *Eectronic
More informationAPPENDIX C FLEXING OF LENGTH BARS
Fexing of ength bars 83 APPENDIX C FLEXING OF LENGTH BARS C.1 FLEXING OF A LENGTH BAR DUE TO ITS OWN WEIGHT Any object ying in a horizonta pane wi sag under its own weight uness it is infinitey stiff or
More informationA nodal collocation approximation for the multidimensional P L equations. 3D applications.
XXI Congreso de Ecuaciones Diferenciaes y Apicaciones XI Congreso de Matemática Apicada Ciudad Rea, 1-5 septiembre 9 (pp. 1 8) A noda coocation approximation for the mutidimensiona P L equations. 3D appications.
More informationFormulas for Angular-Momentum Barrier Factors Version II
BNL PREPRINT BNL-QGS-06-101 brfactor1.tex Formuas for Anguar-Momentum Barrier Factors Version II S. U. Chung Physics Department, Brookhaven Nationa Laboratory, Upton, NY 11973 March 19, 2015 abstract A
More informationCryptanalysis of PKP: A New Approach
Cryptanaysis of PKP: A New Approach Éiane Jaumes and Antoine Joux DCSSI 18, rue du Dr. Zamenhoff F-92131 Issy-es-Mx Cedex France eiane.jaumes@wanadoo.fr Antoine.Joux@ens.fr Abstract. Quite recenty, in
More informationBayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with?
Bayesian Learning A powerfu and growing approach in machine earning We use it in our own decision making a the time You hear a which which coud equay be Thanks or Tanks, which woud you go with? Combine
More informationEfficiency of pseudo-spectrum methods for estimation of Cosmic Microwave Background B-mode power spectrum
Efficiency of pseudo-spectrum methods for estimation of Cosmic Microwave Background B-mode power spectrum For these reasons, many poarization experiments targeting B-modes have been buit or proposed, incuding
More informationElements of Kinetic Theory
Eements of Kinetic Theory Diffusion Mean free path rownian motion Diffusion against a density gradient Drift in a fied Einstein equation aance between diffusion and drift Einstein reation Constancy of
More informationAgenda Administrative Matters Atomic Physics Molecules
Fromm Institute for Lifeong Learning University of San Francisco Modern Physics for Frommies IV The Universe - Sma to Large Lecture 3 Agenda Administrative Matters Atomic Physics Moecues Administrative
More informationHYDROGEN ATOM SELECTION RULES TRANSITION RATES
DOING PHYSICS WITH MATLAB QUANTUM PHYSICS Ian Cooper Schoo of Physics, University of Sydney ian.cooper@sydney.edu.au HYDROGEN ATOM SELECTION RULES TRANSITION RATES DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS
More informationChemical Kinetics Part 2. Chapter 16
Chemica Kinetics Part 2 Chapter 16 Integrated Rate Laws The rate aw we have discussed thus far is the differentia rate aw. Let us consider the very simpe reaction: a A à products The differentia rate reates
More informationA Simple and Efficient Algorithm of 3-D Single-Source Localization with Uniform Cross Array Bing Xue 1 2 a) * Guangyou Fang 1 2 b and Yicai Ji 1 2 c)
A Simpe Efficient Agorithm of 3-D Singe-Source Locaization with Uniform Cross Array Bing Xue a * Guangyou Fang b Yicai Ji c Key Laboratory of Eectromagnetic Radiation Sensing Technoogy, Institute of Eectronics,
More informationQuantum Electrodynamical Basis for Wave. Propagation through Photonic Crystal
Adv. Studies Theor. Phys., Vo. 6, 01, no. 3, 19-133 Quantum Eectrodynamica Basis for Wave Propagation through Photonic Crysta 1 N. Chandrasekar and Har Narayan Upadhyay Schoo of Eectrica and Eectronics
More informationLecture 6 Povh Krane Enge Williams Properties of 2-nucleon potential
Lecture 6 Povh Krane Enge Wiiams Properties of -nuceon potentia 16.1 4.4 3.6 9.9 Meson Theory of Nucear potentia 4.5 3.11 9.10 I recommend Eisberg and Resnik notes as distributed Probems, Lecture 6 1 Consider
More informationQuantum Mechanical Models of Vibration and Rotation of Molecules Chapter 18
Quantum Mechanica Modes of Vibration and Rotation of Moecues Chapter 18 Moecuar Energy Transationa Vibrationa Rotationa Eectronic Moecuar Motions Vibrations of Moecues: Mode approximates moecues to atoms
More informationCandidate Number. General Certificate of Education Advanced Level Examination January 2012
entre Number andidate Number Surname Other Names andidate Signature Genera ertificate of Education dvanced Leve Examination January 212 Physics PHY4/1 Unit 4 Fieds and Further Mechanics Section Tuesday
More informationApproximation and Fast Calculation of Non-local Boundary Conditions for the Time-dependent Schrödinger Equation
Approximation and Fast Cacuation of Non-oca Boundary Conditions for the Time-dependent Schrödinger Equation Anton Arnod, Matthias Ehrhardt 2, and Ivan Sofronov 3 Universität Münster, Institut für Numerische
More informationApplied Nuclear Physics (Fall 2006) Lecture 7 (10/2/06) Overview of Cross Section Calculation
22.101 Appied Nucear Physics (Fa 2006) Lecture 7 (10/2/06) Overview of Cross Section Cacuation References P. Roman, Advanced Quantum Theory (Addison-Wesey, Reading, 1965), Chap 3. A. Foderaro, The Eements
More informationXSAT of linear CNF formulas
XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open
More informationStochastic Variational Inference with Gradient Linearization
Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,
More informationhttps://doi.org/ /epjconf/
HOW TO APPLY THE OPTIMAL ESTIMATION METHOD TO YOUR LIDAR MEASUREMENTS FOR IMPROVED RETRIEVALS OF TEMPERATURE AND COMPOSITION R. J. Sica 1,2,*, A. Haefee 2,1, A. Jaai 1, S. Gamage 1 and G. Farhani 1 1 Department
More informationBohr s atomic model. 1 Ze 2 = mv2. n 2 Z
Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated
More informationis the scattering phase shift associated with the th continuum molecular orbital, and M
Moecuar-orbita decomposition of the ionization continuum for a diatomic moecue by ange- and energy-resoved photoeectron spectroscopy. I. Formaism Hongkun Park and Richard. Zare Department of Chemistry,
More informationarxiv: v2 [astro-ph.co] 17 Oct 2012
CMB EB and TB cross-spectrum estimation via pseudospectrum techniques J. Grain Université Paris-Sud 11, Institut d Astrophysique Spatiae, UMR8617, Orsay, France, F-91405 and CNRS, Orsay, France, F-91405
More informationHomework 05 - H Atom and Electron Configuration
HW05 - H Atom and Eectron Configuration This is a preview of the pubished version of the quiz Started: Sep 25 at 6pm Quiz Instructions Homework 05 - H Atom and Eectron Configuration Question 1 Which of
More informationStatistical Inference, Econometric Analysis and Matrix Algebra
Statistica Inference, Econometric Anaysis and Matrix Agebra Bernhard Schipp Water Krämer Editors Statistica Inference, Econometric Anaysis and Matrix Agebra Festschrift in Honour of Götz Trenker Physica-Verag
More informationRELUCTANCE The resistance of a material to the flow of charge (current) is determined for electric circuits by the equation
INTRODUCTION Magnetism pays an integra part in amost every eectrica device used today in industry, research, or the home. Generators, motors, transformers, circuit breakers, teevisions, computers, tape
More informationEffect of Oxygen Injection into Argon Induction Plasmas on Chemically Non-Equilibrium Conditions
Proceedings of 17th Internationa Symposium on Pasma Chemistry, Toronto, Canada, August 7-12, 25 Effect of Oxygen Injection into Argon Induction Pasmas on Chemicay Non-Equiibrium Conditions Nobuhiko Atsuchi
More informationEXPERIMENT 5 MOLAR CONDUCTIVITIES OF AQUEOUS ELECTROLYTES
EXPERIMENT 5 MOLR CONDUCTIVITIES OF QUEOUS ELECTROLYTES Objective: To determine the conductivity of various acid and the dissociation constant, K for acetic acid a Theory. Eectrica conductivity in soutions
More informationHomework 05 - H Atom and Electron Configuration
HW05 - H Atom and Eectron Configura!on! This is a preview of the pubished version of the quiz Started: Sep 18 at 12:47pm Quiz Instruc!ons Homework 05 - H Atom and Eectron Configuration Question 1 Which
More informationEnhanced thermoelectric performance of carbon nanotubes at. elevated temperature
Enhanced thermoeectric performance of carbon nanotubes at eevated temperature P. H. Jiang, H. J. Liu *, D. D. Fan, L. Cheng, J. Wei, J. Zhang, J. H. Liang, J. Shi Key Laboratory of Artificia Micro- and
More informationInterconnect effects on performance of Field Programmable Analog Array
nterconnect effects on performance of Fied Programmabe Anaog Array D. Anderson,. Bir, O. A. Pausinsi 3, M. Spitz, K. Reiss Motoroa, SPS, Phoenix, Arizona, USA, University of Karsruhe, Karsruhe, Germany,
More informationSequential Decoding of Polar Codes with Arbitrary Binary Kernel
Sequentia Decoding of Poar Codes with Arbitrary Binary Kerne Vera Miosavskaya, Peter Trifonov Saint-Petersburg State Poytechnic University Emai: veram,petert}@dcn.icc.spbstu.ru Abstract The probem of efficient
More informationLegendre Polynomials - Lecture 8
Legendre Poynomias - Lecture 8 Introduction In spherica coordinates the separation of variabes for the function of the poar ange resuts in Legendre s equation when the soution is independent of the azimutha
More informationNEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION
NEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION Hsiao-Chang Chen Dept. of Systems Engineering University of Pennsyvania Phiadephia, PA 904-635, U.S.A. Chun-Hung Chen
More informationSTA 216 Project: Spline Approach to Discrete Survival Analysis
: Spine Approach to Discrete Surviva Anaysis November 4, 005 1 Introduction Athough continuous surviva anaysis differs much from the discrete surviva anaysis, there is certain ink between the two modeing
More informationA simple reliability block diagram method for safety integrity verification
Reiabiity Engineering and System Safety 92 (2007) 1267 1273 www.esevier.com/ocate/ress A simpe reiabiity bock diagram method for safety integrity verification Haitao Guo, Xianhui Yang epartment of Automation,
More information17 Lecture 17: Recombination and Dark Matter Production
PYS 652: Astrophysics 88 17 Lecture 17: Recombination and Dark Matter Production New ideas pass through three periods: It can t be done. It probaby can be done, but it s not worth doing. I knew it was
More informationSydU STAT3014 (2015) Second semester Dr. J. Chan 18
STAT3014/3914 Appied Stat.-Samping C-Stratified rand. sampe Stratified Random Samping.1 Introduction Description The popuation of size N is divided into mutuay excusive and exhaustive subpopuations caed
More informationc 2007 Society for Industrial and Applied Mathematics
SIAM REVIEW Vo. 49,No. 1,pp. 111 1 c 7 Society for Industria and Appied Mathematics Domino Waves C. J. Efthimiou M. D. Johnson Abstract. Motivated by a proposa of Daykin [Probem 71-19*, SIAM Rev., 13 (1971),
More informationPHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased
PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization
More informationTurbo Codes. Coding and Communication Laboratory. Dept. of Electrical Engineering, National Chung Hsing University
Turbo Codes Coding and Communication Laboratory Dept. of Eectrica Engineering, Nationa Chung Hsing University Turbo codes 1 Chapter 12: Turbo Codes 1. Introduction 2. Turbo code encoder 3. Design of intereaver
More information