arxiv: v1 [hep-lat] 19 Nov 2013
|
|
- Deborah Quinn
- 6 years ago
- Views:
Transcription
1 HU-EP-13/69 SFB/CPP DESY Applicability of Quasi-Monte Carlo for lattice systems arxiv: v1 [hep-lat] 19 ov 2013, a,b Tobias Hartung, c Karl Jansen, b Hernan Leovey, Anreas Griewank an Michael Müller-Preussker a a Humbolt-University Berlin, Department of Physics Unter en Linen 6, D Berlin, Germany b IC, DESY Zeuthen Platanenallee 6, D Zeuthen, Germany c King s College Lonon, Department of Mathematics Stran, Lonon WC2R 2LS, Unite Kingom Humbolt-University Berlin, Department of Mathematics Unter en Linen 6, D Berlin, Germany Anreas.Ammon@esy.e, tobias.hartung@kcl.ac.uk, Karl.Jansen@esy.e, leovey@math.hu-berlin.e, griewank@math.hu-berlin.e, mmp@physik.hu-berlin.e This project investigates the applicability of quasi-monte Carlo methos to Eucliean lattice systems in orer to improve the asymptotic error scaling of observables for such theories. The error of an observable calculate by averaging over ranom observations generate from orinary Monte Carlo simulations scales like 1/2, where is the number of observations. By means of quasi- Monte Carlo methos it is possible to improve this scaling for certain problems to 1, or even further if the problems are regular enough. We aapte an applie this approach to simple systems like the quantum harmonic an anharmonic oscillator an verifie an improve error scaling of all investigate observables in both cases. 31st International Symposium on Lattice Fiel Theory - LATTICE 2013 July 29 - August 3, 2013 Mainz, Germany Speaker. c Copyright owne by the author(s) uner the terms of the Creative Commons Attribution-onCommercial-ShareAlike Licence.
2 1. Motivation The quasi-monte Carlo (QMC) metho an their ranomizations (RQMC) are highly interesting for the omain of lattice fiel theory. It offers the possibility to improve tremenously the asymptotic error scaling of observables obtaine from Monte Carlo (MC) simulations of lattice path integrals. Substantial reuctions in computing time coul be achieve if the QMC approach coul eventually be applie to lattice-qcd (quantum chromoynamics in its lattice regularize form). A mathematical review of the QMC approach can be foun in [1]. The major part of this contribution is base on our paper [2] (cf. also [3]). The reaer intereste in more etails is referre to this reference at any point of the following iscussion. In orer to better unerstan the point where the QMC approach sets in with its improvement, we want to outline the typical workflow uring the treatment of a general lattice problem with conventional methos. Such a lattice system might be escribe by the partition function Z = Dx e S[x] given the action S. An observable O coul be calculate by O = Z 1 Dx e S[x] O[x]. Dx stans for the path integral measure of all ynamic fiels relevant to the moel uner consieration. This coul be for example the gauge fiel measure for lattice gauge theories or simply a particle path measure for quantum mechanical problems. It is originate in the high imensionality of the lattice path integral that it can naturally only be treate by means of MC methos, provie that analytic solutions are missing. Within the variety of MC approaches the Markov chain-monte Carlo (Mc-MC) approach turns out to be the most efficient, as it allows for importance sampling. During a Mc-MC simulation a number of fiel configurations (x i ) i=1... is generate successively, each of them base on its preecessor an istribute (after the thermalization) accoring to the Boltzmann weight. Then for each sample x i the observable O is measure, leaing to samples O i. Often, these samples are istribute normal, at least to a goo approximation in most cases. Then, the asymptotic error of the mean Ō = 1 i=1 O i scales like 1 2. A fixe-factor increase of the error often arises from correlations between successive observations O i, which originates in the nature of Mc-MC methos. Both features, the crue asymptotic error scaling an the possibly strong auto-correlation, lea to a necessity to generate a large amount of samples to reach a given error level. In many cases it is even impossible to reach the target accuracy ue to the lack of sufficient computing resources. The QMC approach provies the potential to circumvent the aforementione problems, as it exhibits certain favorable properties. Most importantly, it is able to improve the error scaling to 1, given that certain conitions are met (see [2]). It is further encouraging to realize that the QMC technique has alreay been applie successfully in other fiels like financial mathematics [4] for example. But before coming to specific emonstrations of this fascinating approach, we want to take a closer look to a prominent feature of QMC samples. 2. Quasi-Monte Carlo point sets are more uniform Most point sets constructe through QMC techniques fulfill a so-calle low-iscrepancy property (see [2], sections 3 an 4) also referre to uniformity or more uniform (than conventional Monte Carlo point sets). This property can be illustrate through a simple example in two imensions, but whose results can be generalize to arbitrary many imensions. For this experiment the unit square [0,1] [0,1], subivie into 8 8 small squares of equal size, shoul be consiere. ow, 512 points are generate pseuo-ranomly 1 an uniformly in this unit square. Then, for each of the 8 8 = 64 little squares the number of containe points are counte. An example of the outcome of such an experiment is shown in the upper iagram of figure 1. 1 We use the Mersenne Twister pseuo ranom number generator [5]. 2
3 Pseuo ranom 2 point set Histogram of counts count frequency Figure 1: The pseuo-ranom sampling of 512 points in a unit square (upper plot). For each little square the number of containing points is counte an inicate through a colour coe. The meaning of each color coe can be seen from the lower iagram, where a histogram of the counts is shown. The color of each square correspons to the number of points it contains. In the lower part of this figure we have plotte a histogram of the counts. We can clearly see that the istribution of counts is rather broa. This means in practice that many squares contain significantly more or less samples than one woul expect on average, namely 8. If this set of points (x i,y i ) i= woul be use for a Monte Carlo approximation of a two-imensional integral f (x,y)xy i=1 f (x i,y i ) function values in squares with very many or very few points contribute too much or too less respectively to the overall average. This can be seen as a possibly avoiable source of extraneous fluctuations which have nothing to o with the nature of the problem, the integral of f. Hence, it is highly esirable if the filling of squares with samples woul happen more evenly. We will see in the following repetition of the experiment with QMC samples, that this is exactly what can be provie by the QMC metho. We want to mention, that the istribution on the bottom of figure 1 can be escribe theoretically by the Poisson istribution, in the limit of infinitely many little squares keeping the average count per little square fixe (at 8). We repeat now the experiment with exactly the same parameters but instea of a pseuo-ranom number generator we employ the Sobol approach [6], a special QMC metho, for the generation of points. The result is shown in figure 2. As can be seen in the upper plot, the filling of squares in fact happens completely even. Each little square contains exactly 8 points, an this leas to a elta histogram (shown on the lower part of figure 2). If again the points are use in the approximation of a two-imensional integral the function values from each square contributing to the average are represente much better with respect to the area they cover, an hence, smaller stochastic fluctuations are expecte, leaing very likely to smaller errors of this approximation. Through this illustration we might get an unerstaning on how the more evenly istribute samples from QMC methos coul help to ecrease the natural statistical fluctuations of stochastic approximations. 3. Lattice harmonic an anharmonic oscillator We want to briefly introuce the quantum mechanical harmonic an anharmonic oscillator quantize through the lattice path integral, which we will investigate numerically later on. These systems have been investigate in etail alreay in [7] using the Metropolis algorithm, which is consiere as a Markov chain-monte Carlo metho. The unerlying action S = a i=1 ( M0 2 (x i+1 x i ) 2 a 2 + µ2 2 x2 i + λx 4 i ) (3.1) 3
4 with the perioic bounary conition x +1 = x 1 is obtaine from the iscretization of the classical mechanical action of a particle with mass M 0 passing along the path x(t) consiere in Eucliean time on an equiistant finite time lattice with lattice spacing a an lattice points (extent T = a). RQMC point set Histogram of counts count frequency Figure 2: Distribution of 512 Sobol points generate uniformly in the unit square (upper plot) an histogram of counts (lower plot). See also escription of figure 1. The time erivative ẋ(t) is replace by the forwar finite ifference 1 a (x i+1 x i ). λ controls the strength of the anharmonic term x 4 i. Hence, the harmonic oscillator is obtaine for λ = 0 an a spring constant µ 2 > 0. This conition has to be met for a convergent path integral. The anharmonic oscillator can be simulate with λ > 0 an µ 2 R, both being finite. The case µ 2 < 0, to which we restrict ourselves in the following, is particularly interesting, as the potential exhibits two minima in this case (cf. ouble-well potential). The quantization is performe through the partition function Z = e S(x) x 1...x. An observable O of the so quantize system can be expresse as O = Z 1 O(x) e S(x) x 1...x. The primary physical observables are X 2 = 1 i=1x 2 i, X 4 = 1 i=1x 4 i, an Γ(τ) = 1 i=1 x i x i+ τ a. (3.2) The groun state energy E 0 = µ 2 X 2 + 3λX 4 + µ2 16 an the energy gap E = E 1 E 0 between the groun an first excite state can be erive from them. The latter is etermine from a non-linear regression of the moel Γ(τ) C ( e Eτ + e E(T τ)) to the ata for the correlator Γ(τ), efine in (3.2), in a range of times τ where the influence of higher-than-the-first excite states is negligible. 4. Gaussian Sampling As the action of the harmonic oscillator is at most quaratic in the variables x i, the corresponing partition function can be expresse as a multivariate Gaussian integral Z = exp ( 1 2 xt C 1 x ), where the components of C 1 are given by (C 1 ) i, j = 2M 0 a ((1 + a2 µ 2 2M 0 )δ i, j 1 2 (δ i, j+1 + δ i+1, j ) ) (obtaine from comparing: 1 2 xt C 1 x = S(x)). C is calle the covariance matrix. For such a case, the sampling of lattice paths x is particularly straightforwar, an can be base on samples z, which are generate uniformly in the -imensional unit cube. This is particularly useful, as most RQMC methos provie samples in this form. Hence, our algorithm aiming at the generation of properly istribute samples x starts in the first step with 1. the generation of a uniform sample z = (z 1,...,z ) t [0,1]. This is either, as mentione above, a RQMC sample stemming from a scramble (ranomize) Sobol point set, using irection numbers from F. Kuo s page or a sample obtaine from a separate uniform sampling of each imension with a pseuoranom number generator. The latter case will be referre to as (conventional) Monte Carlo (MC) sampling in the following. 2. In the next step, univariate Gaussian samples w = (w 1,...,w ) t are generate by applying the inverse stanar normal istribution function Φ 1 to the z i an multiplying the result with 4
5 the square root of the eigenvalues λ i of C: w i = ( ) λ i Φ 1 z π 1 (i). (4.1) ( ) The eigenvalues are given in a close form as λ i = 2M0 1. a (u cos(2πi/) As inicate through the (inverse of the) permutation π, the orer of imensions in z has to be moifie such that the component z 1 comes upon the largest eigenvalue, z 2 comes upon the secon largest eigenvalue an so on, until the last component z meets the smallest eigenvalue. This can be achieve by etermining a permutation π which brings the eigenvalues in ecreasing orer (λ π(1) λ π(2)... λ π() ) an calculating π 1 as the inverse of this permutation (fulfilling π 1 (π(i)) = i). 3. Finally, the multivariate Gaussian variables x i are generate from the orthonormal transformation x = Gw, where G = R(F) + I(F) is the iscrete Hartley transform an F k,l = 1 e 2πikl/ the iscrete Fourier transform. 5. Results The algorithm we just iscusse was use to generate lattice paths an corresponing samples of the observable X 2 for a harmonic oscillator with the parameters M 0 = 0.5, a = 0.5, µ 2 = 2.0, an = 100 for a fixe number of samples = 2 7,10,13,16,19. For each the experiment is repeate with 300 scramblings (see section 5 in [2]) of the Sobol sequence in the RQMC simulation to allow the approximation of the error an the variance of the error. For the MC simulation 300 ifferent sees have been use to initialize the iniviual runs. A fit of the moel X 2 C α to the etermine errors of X 2 yiels an exponent α = 1.008(15) for the RQMC simulation case an α = 0.49(1) for the MC simulation. A plot illustrating the results is shown in figure 3. The outcome of this investigation basically proofs the full functioning of the (R)QMC metho for a real, even though trivial, physical moel. To stuy the scaling for a less trivial moel, we passe on to the anharmonic oscillator with the parameters λ = 1.0, M 0 = 0.5, µ 2 = 16, a = 1.5/. The experiment was performe for = 100 an 1000 imensions. The sampling of lattice paths is less straightforwar in the present case, but can be realize on the basis of the sampling metho we use before for the harmonic oscillator. This happens with the ai of the reweighting approach (cf. section 7 in [2]), which was first use in [8]. Aitionally to the physical action S, escribing the anharmonic oscillator, an artificial action S is introuce, which is constructe exactly as a harmonic oscillator action with a ifferent set of parameters M 0, a an µ. Then, harmonic oscillator paths (x i ) i=1... are generate corresponing to this unphysical action S. Approximations of observables O of the anharmonic oscillator (escribe by the physical action S) are obtaine from the weighte mean O ( i=1 O(x i)w(x i ) ) / ( i=1 W(x i) ), where the weight function W is given by W(x) = exp( S(x) + S (x)). ow, it remains to fin reasonable parameters with the objective of reucing the fluctuations of the weights W(x i ) as much as possible, leaing most likely to the smallest possible error of the observables. We foun that only the moification of the parameter µ leas alreay to satisfying results, such that M 0 = M 0 an a = a coul be left unchange. A heuristic optimization approach le to a value of µ = We aopte the proceure for the error etermination as well as the regression (fits) for the exponents of the error scaling from the harmonic oscillator experiment. The results are shown in table 1. Plots showing the fits for 1000 imensions are shown in figure 4. From table 1 we can observe that the error scaling of the observables X 2, X 4 an E 0 is significantly improve, although less than in the harmonic oscillator 5
6 Error of <x 2 > for the Harmonic Oscillator Fit QMC, α = Error of <x 2 > 1e 07 1e 05 1e 03 MC QMC 10 i 10 i i = 5 : umber of samples log( <X 2 >) χ 2 = ; of = Figure 3: Error of X 2 for the RQMC an MC simulation of the harmonic oscillator (left). The right plot shows the fit of the asymptotic error scaling for the RQMC simulation. O α logc χ 2 /of X (8) 2.0(1) 7.9 / 6 = 100 X (8) 4.0(1) 13.2 / 6 E (9) 4.0(1) 8.3 / 6 X (14) 2.0(2) 5.0 / 4 = 1000 X (14) 4.0(2) 5.7 / 4 E (13) 4.0(2) 4.0 / 4 Table 1: Results for the error scaling of the observalbes X 2, X 4 an E 0 for the moel of the anharmonic oscillator, simulate through reweighting. Observable errors were fitte to the moel O C α. case. More specifically, we can conclue that, within statistical uncertainties, the error scaling is of O( 3 4 ) in all consiere cases. In a further effort we investigate also the energy gap. In orer to be able to measure this quantity, we increase µ 2 from 16 to 4 (keeping the other parameters fixe), leaing theoretically (in the limit T an a 0) to a change of E from to [9]. For = 100 imensions an sample sizes of = 2 5,8,11,14 we obtain an exponent of α = 0.735(13). These results are very interesting in that the imporve error scaling seems to be rather inepenent of the observable uner consieration. 6. Outlook & conclusions In this contribution we showe a first successful application of one specific RQMC metho to Eucliean lattice moels. We verifie a perfect error scaling of O( 1 ) for the harmonic oscillator an a strongly improve error scaling of O( 3 4 ) for the anharmonic oscillator with a ouble-well potential. The latter investigation also inclue the energy gap, which can be consiere as a rather non-trivial observable, as it is obtaine from the correlator using a non-linear proceure. A better 6
7 Error of <X 2 >, = 1000 Error of <X 4 >, = 1000 Error of E 0, = 1000 <X 2 > 5e 04 2e <X 4 > E Figure 4: Double-log plot of the error scaling of the observables X 2 (left), X 4 (mile) an E 0 (right) with the number of samples in the RQMC approach. The ashe line shows the fit of the moel O C α to the ata. umerical results can be seen in table 1. unerstaning of this in-between behavior of 3 4 is planne for the future, an shoul lea to a better theoretical unerstaning on how the QMC metho behaves when applie to non-trivial problems. Further improvements in the sampling of the anharmonic oscillator are assume when an optimally tune, more generalize covariance matrix is use in the Gaussian sampling step. Furthermore, other promising non-gaussian sampling approaches, like inverse sampling [10], are investigate at the moment, aiming at a better escription of the anharmonic action an probably involving the potential to be applicable to a much broaer class of lattice problems; though, it will be interesting to see in the future how efficient these techniques are in practice. As a next step towars lattice gauge theories we are currently consiering a one-imensional spin like moel, escribe by the action S = ai i=1 1 cos(φ a 2 i+1 φ i ), where I is the moment of inertia, a the lattice spacing, an φ i are angular variables with perioic bounary conitions (φ = φ 0 ). This moel exhibits topological features, visible through the non-vanishing of the wining number a feature not present in the previously consiere oscillator moels but in other lattice gauge theories like QCD. Having manage this moel with generalizable methos it coul be envisage that also generic gauge theories coul be aresse in the future. References [1] F. Kuo, C. Schwab an I. Sloan, AZIAM Journal, 53 (01) (2012). [2] K. Jansen, H. Leovey, A. Ammon, A. Griewank, an M. Müller-Preussker, Comput.Phys.Commun. (2013), arxiv: [hep-lat]. [3] K. Jansen, H. Leovey, A. ube, A. Griewank an M. Müller-Preussker J. Phys.: Conf. Ser. 454 (2013) , arxiv: [hep-lat]. [4] P. Glasserman, Springer-Verlag, ew-york [5] M. Matsumoto an T. ishimura, ACM Trans. Moel. Comput. Simul., 8 no. 1 (1998) [6] I. M. Sobol, U.S.S.R. Comput. Math. an Math. Phys., 7 no. 4 (1967) [7] M. Creutz an B. Freeman, Annals Phys., 132 (1981) 427. [8] A.M. Ferrenberg an R.H. Swensen, Phys. Rev. Lett., 61 no. 23 (1988) [9] R. Blankenbecler an T. A. DeGran, Thomas A. an R. L. Sugar, Phys.Rev., D21 (1980) [10] L. Devroye, Springer-Verlag, ew-york
Applicability of Quasi-Monte Carlo for lattice systems
Applicability of Quasi-Monte Carlo for lattice systems Andreas Ammon 1,2, Tobias Hartung 1,2,Karl Jansen 2, Hernan Leovey 3, Andreas Griewank 3, Michael Müller-Preussker 1 1 Humboldt-University Berlin,
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationSurvey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013
Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y
Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay
More informationExpected Value of Partial Perfect Information
Expecte Value of Partial Perfect Information Mike Giles 1, Takashi Goa 2, Howar Thom 3 Wei Fang 1, Zhenru Wang 1 1 Mathematical Institute, University of Oxfor 2 School of Engineering, University of Tokyo
More informationarxiv: v2 [cond-mat.stat-mech] 11 Nov 2016
Noname manuscript No. (will be inserte by the eitor) Scaling properties of the number of ranom sequential asorption iterations neee to generate saturate ranom packing arxiv:607.06668v2 [con-mat.stat-mech]
More informationPoS(RAD COR 2007)030. Three-jet DIS final states from k -dependent parton showers. F. Hautmann University of Oxford
OUTP-8--P Three-jet DIS final states from k -epenent parton showers University of Oxfor E-mail: hautmann@thphys.ox.ac.uk H. Jung DESY E-mail: jung@mail.esy.e Experimental measurements of angular correlations
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationGeneralizing Kronecker Graphs in order to Model Searchable Networks
Generalizing Kronecker Graphs in orer to Moel Searchable Networks Elizabeth Boine, Babak Hassibi, Aam Wierman California Institute of Technology Pasaena, CA 925 Email: {eaboine, hassibi, aamw}@caltecheu
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationLower bounds on Locality Sensitive Hashing
Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More information1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.
Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency
More informationSemiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom
PHYSICAL REVIEW A 69, 063409 (2004) Semiclassical analysis of long-wavelength multiphoton processes: The Ryberg atom Luz V. Vela-Arevalo* an Ronal F. Fox Center for Nonlinear Sciences an School of Physics,
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationQubit channels that achieve capacity with two states
Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationRole of parameters in the stochastic dynamics of a stick-slip oscillator
Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationSolution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010
NTNU Page of 6 Institutt for fysikk Fakultet for fysikk, informatikk og matematikk This solution consists of 6 pages. Solution to the exam in TFY423 STATISTICAL PHYSICS Wenesay ecember, 2 Problem. Particles
More informationMaking a Wavefunctional representation of physical states congruent with the false vacuum hypothesis of Sidney Coleman
Making a Wavefunctional representation of physical states congruent with the false vacuum hypothesis of Siney Coleman A. W. Beckwith Department of Physics an Texas Center for Superconuctivity an Avance
More informationOn the number of isolated eigenvalues of a pair of particles in a quantum wire
On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in
More informationTHE EFFICIENCIES OF THE SPATIAL MEDIAN AND SPATIAL SIGN COVARIANCE MATRIX FOR ELLIPTICALLY SYMMETRIC DISTRIBUTIONS
THE EFFICIENCIES OF THE SPATIAL MEDIAN AND SPATIAL SIGN COVARIANCE MATRIX FOR ELLIPTICALLY SYMMETRIC DISTRIBUTIONS BY ANDREW F. MAGYAR A issertation submitte to the Grauate School New Brunswick Rutgers,
More informationState observers and recursive filters in classical feedback control theory
State observers an recursive filters in classical feeback control theory State-feeback control example: secon-orer system Consier the riven secon-orer system q q q u x q x q x x x x Here u coul represent
More informationExperimental Robustness Study of a Second-Order Sliding Mode Controller
Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationarxiv: v1 [cond-mat.stat-mech] 9 Jan 2012
arxiv:1201.1836v1 [con-mat.stat-mech] 9 Jan 2012 Externally riven one-imensional Ising moel Amir Aghamohammai a 1, Cina Aghamohammai b 2, & Mohamma Khorrami a 3 a Department of Physics, Alzahra University,
More informationarxiv:hep-th/ v1 3 Feb 1993
NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,
More informationIntroduction to variational calculus: Lecture notes 1
October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational
More informationThe Role of Models in Model-Assisted and Model- Dependent Estimation for Domains and Small Areas
The Role of Moels in Moel-Assiste an Moel- Depenent Estimation for Domains an Small Areas Risto Lehtonen University of Helsini Mio Myrsylä University of Pennsylvania Carl-Eri Särnal University of Montreal
More information1 Heisenberg Representation
1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.
More informationInfluence of weight initialization on multilayer perceptron performance
Influence of weight initialization on multilayer perceptron performance M. Karouia (1,2) T. Denœux (1) R. Lengellé (1) (1) Université e Compiègne U.R.A. CNRS 817 Heuiasyc BP 649 - F-66 Compiègne ceex -
More informationDelocalization of boundary states in disordered topological insulators
Journal of Physics A: Mathematical an Theoretical J. Phys. A: Math. Theor. 48 (05) FT0 (pp) oi:0.088/75-83/48//ft0 Fast Track Communication Delocalization of bounary states in isorere topological insulators
More informationThe AdS/CFT Correspondence PI It from Qubit Summer School: Mukund Rangamani
Lecture The AS/CFT Corresponence PI It from Qubit Summer School: Mukun Rangamani Q. Large N expansion: Consier the following Lagrangian for a zero imensional fiel theory (matrix moel): L = (Tr (ΦΦ) + Tr
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationHow the potentials in different gauges yield the same retarded electric and magnetic fields
How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department
More informationDot trajectories in the superposition of random screens: analysis and synthesis
1472 J. Opt. Soc. Am. A/ Vol. 21, No. 8/ August 2004 Isaac Amiror Dot trajectories in the superposition of ranom screens: analysis an synthesis Isaac Amiror Laboratoire e Systèmes Périphériques, Ecole
More informationQuantum mechanical approaches to the virial
Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from
More informationCharacterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations
Characterizing Real-Value Multivariate Complex Polynomials an Their Symmetric Tensor Representations Bo JIANG Zhening LI Shuzhong ZHANG December 31, 2014 Abstract In this paper we stuy multivariate polynomial
More informationMonte Carlo Methods with Reduced Error
Monte Carlo Methos with Reuce Error As has been shown, the probable error in Monte Carlo algorithms when no information about the smoothness of the function is use is Dξ r N = c N. It is important for
More informationNON-SMOOTH DYNAMICS USING DIFFERENTIAL-ALGEBRAIC EQUATIONS PERSPECTIVE: MODELING AND NUMERICAL SOLUTIONS. A Thesis PRIYANKA GOTIKA
NON-SMOOTH DYNAMICS USING DIFFERENTIAL-ALGEBRAIC EQUATIONS PERSPECTIVE: MODELING AND NUMERICAL SOLUTIONS A Thesis by PRIYANKA GOTIKA Submitte to the Office of Grauate Stuies of Texas A&M University in
More informationTEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE
TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS Yannick DEVILLE Université Paul Sabatier Laboratoire Acoustique, Métrologie, Instrumentation Bât. 3RB2, 8 Route e Narbonne,
More informationConstrained optimal discrimination designs for Fourier regression models
Ann Inst Stat Math (29) 6:43 57 DOI.7/s463-7-33-5 Constraine optimal iscrimination esigns for Fourier regression moels Stefanie Bieermann Holger Dette Philipp Hoffmann Receive: 26 June 26 / Revise: 9 March
More information4. Important theorems in quantum mechanics
TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationNonlinear Dielectric Response of Periodic Composite Materials
onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics
More informationIntroduction to Markov Processes
Introuction to Markov Processes Connexions moule m44014 Zzis law Gustav) Meglicki, Jr Office of the VP for Information Technology Iniana University RCS: Section-2.tex,v 1.24 2012/12/21 18:03:08 gustav
More informationA Modification of the Jarque-Bera Test. for Normality
Int. J. Contemp. Math. Sciences, Vol. 8, 01, no. 17, 84-85 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1988/ijcms.01.9106 A Moification of the Jarque-Bera Test for Normality Moawa El-Fallah Ab El-Salam
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationLecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations
Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationScalar Field Theory in the AdS/CFT Correspondence Revisited
Scalar Fiel Theory in the AS/CFT Corresponence Revisite arxiv:hep-th/9907079v4 Feb 000 Pablo Minces an Victor. Rivelles Universiae e São Paulo, Instituto e Física Caixa Postal 66.38 - CEP 0535-970 - São
More informationSparse Reconstruction of Systems of Ordinary Differential Equations
Sparse Reconstruction of Systems of Orinary Differential Equations Manuel Mai a, Mark D. Shattuck b,c, Corey S. O Hern c,a,,e, a Department of Physics, Yale University, New Haven, Connecticut 06520, USA
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationLower Bounds for the Smoothed Number of Pareto optimal Solutions
Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.
More informationA Review of Multiple Try MCMC algorithms for Signal Processing
A Review of Multiple Try MCMC algorithms for Signal Processing Luca Martino Image Processing Lab., Universitat e València (Spain) Universia Carlos III e Mari, Leganes (Spain) Abstract Many applications
More informationCharge { Vortex Duality. in Double-Layered Josephson Junction Arrays
Charge { Vortex Duality in Double-Layere Josephson Junction Arrays Ya. M. Blanter a;b an Ger Schon c a Institut fur Theorie er Konensierten Materie, Universitat Karlsruhe, 76 Karlsruhe, Germany b Department
More informationEVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION OF UNIVARIATE TAYLOR SERIES
MATHEMATICS OF COMPUTATION Volume 69, Number 231, Pages 1117 1130 S 0025-5718(00)01120-0 Article electronically publishe on February 17, 2000 EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION
More informationDistribution Theory for Discontinuous Test Functions and Differential Operators with Generalized Coefficients
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 1, 9733 1996 ARTICLE NO 56 Distribution Theory for Discontinuous Test Functions an Differential Operators with Generalize Coefficients P Kurasov* Department
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationIERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210
IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.210 Discrete an Continuous Dynamics in Nonlinear Monopolies Akio Matsumoto Chuo University Ferenc
More informationCONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN
More informationSome New Thoughts on the Multipoint Method for Reactor Physics Applications. Sandra Dulla, Piero Ravetto, Paolo Saracco,
Jeju, Korea, April 16-20, 2017, on USB 2017 Some New Thoughts on the Multipoint Metho for Reactor Physics Applications Sanra Dulla, Piero Ravetto, Paolo Saracco, Politecnico i Torino, Dipartimento Energia,
More informationAcute sets in Euclidean spaces
Acute sets in Eucliean spaces Viktor Harangi April, 011 Abstract A finite set H in R is calle an acute set if any angle etermine by three points of H is acute. We examine the maximal carinality α() of
More informationThis module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics
This moule is part of the Memobust Hanbook on Methoology of Moern Business Statistics 26 March 2014 Metho: Balance Sampling for Multi-Way Stratification Contents General section... 3 1. Summary... 3 2.
More informationLecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012
CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration
More informationWitten s Proof of Morse Inequalities
Witten s Proof of Morse Inequalities by Igor Prokhorenkov Let M be a smooth, compact, oriente manifol with imension n. A Morse function is a smooth function f : M R such that all of its critical points
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More informationarxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003
Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep
More informationPerfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs
Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs Ashish Goel Michael Kapralov Sanjeev Khanna Abstract We consier the well-stuie problem of fining a perfect matching in -regular bipartite
More informationPHYS 414 Problem Set 2: Turtles all the way down
PHYS 414 Problem Set 2: Turtles all the way own This problem set explores the common structure of ynamical theories in statistical physics as you pass from one length an time scale to another. Brownian
More informationTutorial on Maximum Likelyhood Estimation: Parametric Density Estimation
Tutorial on Maximum Likelyhoo Estimation: Parametric Density Estimation Suhir B Kylasa 03/13/2014 1 Motivation Suppose one wishes to etermine just how biase an unfair coin is. Call the probability of tossing
More informationd dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1
Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of
More informationAPPPHYS 217 Thursday 8 April 2010
APPPHYS 7 Thursay 8 April A&M example 6: The ouble integrator Consier the motion of a point particle in D with the applie force as a control input This is simply Newton s equation F ma with F u : t q q
More information. Using a multinomial model gives us the following equation for P d. , with respect to same length term sequences.
S 63 Lecture 8 2/2/26 Lecturer Lillian Lee Scribes Peter Babinski, Davi Lin Basic Language Moeling Approach I. Special ase of LM-base Approach a. Recap of Formulas an Terms b. Fixing θ? c. About that Multinomial
More informationarxiv: v1 [hep-ex] 4 Sep 2018 Simone Ragoni, for the ALICE Collaboration
Prouction of pions, kaons an protons in Xe Xe collisions at s =. ev arxiv:09.0v [hep-ex] Sep 0, for the ALICE Collaboration Università i Bologna an INFN (Bologna) E-mail: simone.ragoni@cern.ch In late
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationIn the usual geometric derivation of Bragg s Law one assumes that crystalline
Diffraction Principles In the usual geometric erivation of ragg s Law one assumes that crystalline arrays of atoms iffract X-rays just as the regularly etche lines of a grating iffract light. While this
More informationCascaded redundancy reduction
Network: Comput. Neural Syst. 9 (1998) 73 84. Printe in the UK PII: S0954-898X(98)88342-5 Cascae reunancy reuction Virginia R e Sa an Geoffrey E Hinton Department of Computer Science, University of Toronto,
More informationBalancing Data Security and Blocking Performance with Spectrum Randomization in Optical Networks
Balancing Data Security an Blocking Performance with Spectrum Ranomization in Optical Networks Saneep Kumar Singh, Wolfgang Bziuk, an Amela Jukan Technische Universität Carolo-Wilhelmina zu Braunschweig,
More informationMULTISCALE FRICTION MODELING FOR SHEET METAL FORMING
MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING Authors J. HOL 1, M.V. CID ALFARO 2, M.B. DE ROOIJ 3 AND T. MEINDERS 4 1 Materials innovation institute (M2i) 2 Corus Research Centre 3 University of
More informationRobust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k
A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine
More informationA Randomized Approximate Nearest Neighbors Algorithm - a short version
We present a ranomize algorithm for the approximate nearest neighbor problem in - imensional Eucliean space. Given N points {x } in R, the algorithm attempts to fin k nearest neighbors for each of x, where
More information