arxiv: v1 [hep-lat] 19 Nov 2013

Size: px
Start display at page:

Download "arxiv: v1 [hep-lat] 19 Nov 2013"

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 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 information

Lagrangian and Hamiltonian Mechanics

Lagrangian 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 information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey 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 information

6 General properties of an autonomous system of two first order ODE

6 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 information

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

THE 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 information

SYNCHRONOUS SEQUENTIAL CIRCUITS

SYNCHRONOUS 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 information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The 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 information

Chapter 6: Energy-Momentum Tensors

Chapter 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 information

Quantum Mechanics in Three Dimensions

Quantum 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 information

The Exact Form and General Integrating Factors

The 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 information

Table of Common Derivatives By David Abraham

Table 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 information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-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 information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles 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 information

Expected Value of Partial Perfect Information

Expected 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 information

arxiv: v2 [cond-mat.stat-mech] 11 Nov 2016

arxiv: 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 information

PoS(RAD COR 2007)030. Three-jet DIS final states from k -dependent parton showers. F. Hautmann University of Oxford

PoS(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 information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing 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 information

Generalizing Kronecker Graphs in order to Model Searchable Networks

Generalizing 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 information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 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 information

Lower bounds on Locality Sensitive Hashing

Lower 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 information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching 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 information

Generalization of the persistent random walk to dimensions greater than 1

Generalization 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 information

1 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.

1 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 information

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Semiclassical 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

θ 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 information

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

Thermal 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 information

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods

Hyperbolic 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 information

NOTES 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 (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 information

Schrödinger s equation.

Schrö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 information

The Principle of Least Action

The 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

'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 information

Conservation laws a simple application to the telegraph equation

Conservation 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 information

Qubit channels that achieve capacity with two states

Qubit 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 information

05 The Continuum Limit and the Wave Equation

05 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 information

Role of parameters in the stochastic dynamics of a stick-slip oscillator

Role 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 information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 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 information

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010

Solution 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 information

Making 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 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 information

On the number of isolated eigenvalues of a pair of particles in a quantum wire

On 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 information

THE 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 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 information

State observers and recursive filters in classical feedback control theory

State 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 information

Experimental Robustness Study of a Second-Order Sliding Mode Controller

Experimental 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 information

Stable and compact finite difference schemes

Stable 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 information

arxiv: v1 [cond-mat.stat-mech] 9 Jan 2012

arxiv: 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 information

arxiv:hep-th/ v1 3 Feb 1993

arxiv: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 information

Introduction to variational calculus: Lecture notes 1

Introduction 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 information

The Role of Models in Model-Assisted and Model- Dependent Estimation for Domains and Small Areas

The 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 information

1 Heisenberg Representation

1 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 information

Influence of weight initialization on multilayer perceptron performance

Influence 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 information

Delocalization of boundary states in disordered topological insulators

Delocalization 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 information

The AdS/CFT Correspondence PI It from Qubit Summer School: Mukund Rangamani

The 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 information

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

inflow 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 information

How the potentials in different gauges yield the same retarded electric and magnetic fields

How 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 information

Dot trajectories in the superposition of random screens: analysis and synthesis

Dot 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 information

Quantum mechanical approaches to the virial

Quantum 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 information

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations

Characterizing 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 information

Monte Carlo Methods with Reduced Error

Monte 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 information

NON-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 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 information

TEMPORAL AND TIME-FREQUENCY CORRELATION-BASED BLIND SOURCE SEPARATION METHODS. Yannick DEVILLE

TEMPORAL 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 information

Constrained optimal discrimination designs for Fourier regression models

Constrained 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 information

4. Important theorems in quantum mechanics

4. 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 information

Chapter 4. Electrostatics of Macroscopic Media

Chapter 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 information

Nonlinear Dielectric Response of Periodic Composite Materials

Nonlinear 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 information

Introduction to Markov Processes

Introduction 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 information

A Modification of the Jarque-Bera Test. for Normality

A 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 information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-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 information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture 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 information

Least-Squares Regression on Sparse Spaces

Least-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 information

Scalar Field Theory in the AdS/CFT Correspondence Revisited

Scalar 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 information

Sparse Reconstruction of Systems of Ordinary Differential Equations

Sparse 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 information

Linear First-Order Equations

Linear 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 information

Lower Bounds for the Smoothed Number of Pareto optimal Solutions

Lower 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 information

A Review of Multiple Try MCMC algorithms for Signal Processing

A 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 information

Charge { Vortex Duality. in Double-Layered Josephson Junction Arrays

Charge { 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 information

EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION OF UNIVARIATE TAYLOR SERIES

EVALUATING 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 information

Distribution Theory for Discontinuous Test Functions and Differential Operators with Generalized Coefficients

Distribution 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 information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE 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 information

IERCU. 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 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 information

CONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS

CONSERVATION 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 information

Some New Thoughts on the Multipoint Method for Reactor Physics Applications. Sandra Dulla, Piero Ravetto, Paolo Saracco,

Some 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 information

Acute sets in Euclidean spaces

Acute 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 information

This module is part of the. Memobust Handbook. on Methodology of Modern Business Statistics

This 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 information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture 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 information

Witten s Proof of Morse Inequalities

Witten 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 information

Analytic Scaling Formulas for Crossed Laser Acceleration in Vacuum

Analytic 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 information

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

arxiv: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 information

Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs

Perfect 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 information

PHYS 414 Problem Set 2: Turtles all the way down

PHYS 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 information

Tutorial on Maximum Likelyhood Estimation: Parametric Density Estimation

Tutorial 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 information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d 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 information

APPPHYS 217 Thursday 8 April 2010

APPPHYS 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.

. 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 information

arxiv: v1 [hep-ex] 4 Sep 2018 Simone Ragoni, for the ALICE Collaboration

arxiv: 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 information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon 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 information

In the usual geometric derivation of Bragg s Law one assumes that crystalline

In 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 information

Cascaded redundancy reduction

Cascaded 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 information

Balancing Data Security and Blocking Performance with Spectrum Randomization in Optical Networks

Balancing 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 information

MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING

MULTISCALE 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 information

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k

Robust 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 information

A Randomized Approximate Nearest Neighbors Algorithm - a short version

A 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