Optimal driving waveform for the overdamped, rocking ratchets

Size: px
Start display at page:

Download "Optimal driving waveform for the overdamped, rocking ratchets"

Transcription

1 Optimal driving waveform for the overdamped, rocking ratchets Maria Laura Olivera Instituto de Ciencias de la Universidad Nacional de General Sarmiento. Buenos Aires Argentina Directors: R. Álvarez-Nodarse & N.R. Quintero IMUS Doc-course May 20, 2010 Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

2 Plan 1 Introduction to the Ratchet effect 2 Motivation 3 Working Model 4 Objetives and methodology of the project 5 References Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

3 Introduction to the Ratchet effect Introduction Ratchet transport phenomena means a net directed motion of particles or solitons induced by zero average forces. i.e. presence of the thermal noise and some prothotipical perturbation can drives the system out of the equilibrium (directed transport). The ratchet effect usually occurs when there is a break of at less one of main symmetries of the system: temporal symmetries: forces) shift-symmetry (caused by time dependent spatial symmetries (caused by ratchet potential ) We focus in the breaking of the temporal symetries Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

4 Introduction to the Ratchet effect Introduction Elements present in the physical system: Potential U(x) Friction force ηẋ Force f(t) if any Noise ξ(t) The quantity of foremost in the context of transport in periodic systems is the particle current v: v = ẋ The ratchet optimization is given by the maximun of average current of the particle. We study ratchet current and its relationship with the kind of driving force f(t): the driving force shape that maximizes ẋ. Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

5 Introduction to the Ratchet effect Introduction The simplest system is a particle or soliton moving in a potential under influence of a friction force and noise. According to Newton s second law: mẍ + U (x) + ηẋ f (t) = ξ(t), where the left-hand represent the deterministic, conservative part of the particles dynamics, while the right-hand side account for the effects of the thermal environment. Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

6 Introduction to the Ratchet effect The system is assumed to be overdamped, that is, the inertial term ẍ is negligible and we arrive to: ηẋ + U (x) = f (t) + ξ(t). Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

7 Motivation A few works have investigated the dependence of the ratchet current on amplitudes and phases of forces. We focus in two of them. 1. The first one: conjectures there exist a particular force waveform that optimally enhances directed transported. proposes a measure of degree of symmetry breaking (DSB). He defines the DSB of the symmetries of the force f(t) by the expresion: D s (f ) f (t + T /2) T (1) f (t) uses the biharmonic force to prove that (1) exhibits tha same quantitaive behaviour as a funtion of the symmetry-breaking parameters Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

8 Motivation consideres tha case of the elliptic force: f ellip (t) = ɛf (t; T, m, θ) ɛsn(ωt + Θ; m)cn(ωt + Θ; m), and put a biharmonic aproximation of it in (1). recovers the same optimal values obtained for the biharmonic force Concludes that the biharmonic force optimally enhances directed transport. Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

9 2. The second one: Motivation Claims that the optimal driving waveform among a wide class of admissible funtions for an overdamped, rocking ratchet is dichotomous, with a form: { f L, 0 t d (t) = N, d t T where N 0 is a solution of ν (N ) = ν(l) ν(n ) L N uses: ( ẋ = ν(m)d + ν M d 1 d ) (1 d) put the dichotomous force in it and maximizes it with respect to M and to N, finding the dichotomic optimal force Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

10 Moreover, we think that there is not a general optimal force. Probably there is an optimal force for every family of forces. Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / The second one: Motivation Then fixes the parameters ɛ 1 =L/3 and ɛ 2 =2L/3 valid for for the optimal biharmonic force in the limit of small force amplitud in the dichotomous force. concludes that the dichotomous responses is better than the biharmonic responses for larges L Question: The dichotomous responses is better than the biharmonic responses for larges L but the values of parameters fixed are valided only in limit of small force amplitud for the biharmonic force. Then we suspect that the main conclusion: the optimal force is dichotomous, is not correct.

11 Working Model Working Model We study the ratchet dynamic of an overdamped particle with position x(t) moving in a potential, suject to a force and noise. where: ηẋ U (x) + f (t) + ξ(t) = 0 ηẋ represent the friction, with η the damping coeficient and ẋ the derivate on time of the position, the velocity U (x) is the derivate on position of the symmetric periodic potential ξ(t) is an independent Gaussian white noise Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

12 Working Model Working Model We study the ratchet dynamic of an overdamped particle with position x(t) moving in a potential, suject to a force and noise. where: ηẋ U (x) + f (t) + ξ(t) = 0 ηẋ represent the friction, with η the damping coeficient and ẋ the derivate on time of the position, the velocity U (x) is the derivate on position of the symmetric periodic potential ξ(t) is an independent Gaussian white noise Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

13 Working Model Working Model We study the ratchet dynamic of an overdamped particle with position x(t) moving in a potential, suject to a force and noise. where: ηẋ U (x) + f (t) + ξ(t) = 0 ηẋ represent the friction, with η the damping coeficient and ẋ the derivate on time of the position, the velocity U (x) is the derivate on position of the symmetric periodic potential ξ(t) is an independent Gaussian white noise Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

14 Working Model Working Model We study the ratchet dynamic of an overdamped particle with position x(t) moving in a potential, suject to a force and noise. where: ηẋ U (x) + f (t) + ξ(t) = 0 ηẋ represent the friction, with η the damping coeficient and ẋ the derivate on time of the position, the velocity U (x) is the derivate on position of the symmetric periodic potential ξ(t) is an independent Gaussian white noise Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

15 Working Model Hypothesis: ξ(t) is independent Gaussian withe noise of zero mean, ξ(t) = 0, with ξ(t)ξ(t ) = Dδ(t t ), where D is noise coeficient and δ is dirac s delta f(t) breaks the time-shift symmetry f (t) = f (t + T /2) Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

16 Working Model We take: f(t) biharmonic force: f (t) = ɛ 1 cos(qωt + φ 1 ) + ɛ 2 cos(pωt + φ 2 ) where T=2π/ ω, φ 1 and φ 2 are the phases, p and q are coprimes with p+q odd and the amplitudes ɛ 1 and ɛ 2 are small U(x) periodic potential: We fix: q=1 p=2 ɛ 1 = ɛf 0 ɛ 2 = (1 ɛ)f 0 U(x) = U 0 (1 cos(x)) Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

17 Objetives and methodology of the project Goals We will: study the average ratchet velocity dependence on the force strength check numerically that φ opt = 0 check if there is average velocity ratchet dependence on the frequency ω Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

18 Objetives and methodology of the project Goals We will: check the theorem [Lade, 2008]: If force(s) exist that optimize the current : ẋ = 1 T T 0 v(f (t))dt, then one such optimal force is dichotomous, that is, for all t, f(t)=m or N, for some M and N. Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

19 Objetives and methodology of the project Numerical work Methodology to solve the problem: Numerical simulations Proofs of the analytic results Numerical simulations: We solve the equation: ẋ(t) = ηf (t) U (x) + ξ(t) using f (t) = ɛ 1 cos(qωt + φ 1 ) + ɛ 2 cos(pωt + φ 2 ) U(x) = U 0 (1 cos(x)) with U 0 = 1, ω = 0.1, φ 1 = 0, φ 2 = π, p=2, q=1, η = 1, ɛ = 1/3,and D= 0.32 we obtain x(t) averaging 1000 realizations. Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

20 Objetives and methodology of the project Numerical work ẋ(t) = f (t) sin(x) + ξ(t) f (t) = 1 3 cos(ωt+) f 0cos(2ωt + φ 2 ) The figure shows the position of the particle on the time x(t): Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

21 Objetives and methodology of the project Numerical work ẋ(t) = f (t) sin(x) + ξ(t) f (t) = 1 3 cos(ωt+) f 0cos(2ωt + φ 2 ) The figure shows the of the particle on L: Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

22 Objetives and methodology of the project Numerical work ẋ(t) = f (t) sin(x) + ξ(t) f (t) = 1 3 cos(ωt+) f 0cos(2ωt + φ 2 ) The figure shows the of the particle on L: Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

23 Objetives and methodology of the project Numerical work We are currely working on the numerical simulation in order to take the ẋ on φ optimo. The next step, after having the numerical results,is to check the theorems. Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

24 References References 1 S.J.Lade,J.Phys.A:Math.Theor. 41, (2008) 2 R.Chacon, J.Phys.A:Math.Theor. 40, F413(2007) 3 N.R.Quintero, J. A. Cuesta and R. Alvarez-Nodarse,Phys.Rev. E 81, (2010) 4 M. San Miguel, R. Toral, Stochastic Effects in Physical Systems. In: Instabilities and Nonequilibrium Structures VI, Kluwer Academic Publishers,ISBN (2000), Part I. Maria Laura Olivera (ICI.UNGS) Optimal driving waveform rocking ratchets May 20, / 21

Effects of thermal noise on transport in soliton ratchet systems

Effects of thermal noise on transport in soliton ratchet systems Effects of thermal noise on transport in soliton ratchet systems With L. Morales, V. Stehr, P. Müller (Bayreuth) A. Sanchez (Madrid), N.R. Quintero (Sevilla) A.R. Bishop (Los Alamos) J. Buceta, K. Lindenberg

More information

A path integral approach to the Langevin equation

A path integral approach to the Langevin equation A path integral approach to the Langevin equation - Ashok Das Reference: A path integral approach to the Langevin equation, A. Das, S. Panda and J. R. L. Santos, arxiv:1411.0256 (to be published in Int.

More information

Chapter Three Theoretical Description Of Stochastic Resonance 24

Chapter Three Theoretical Description Of Stochastic Resonance 24 Table of Contents List of Abbreviations and Symbols 5 Chapter One Introduction 8 1.1 The Phenomenon of the Stochastic Resonance 8 1.2 The Purpose of the Study 10 Chapter Two The Experimental Set-up 12

More information

Linear and Nonlinear Oscillators (Lecture 2)

Linear and Nonlinear Oscillators (Lecture 2) Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical

More information

BSc/MSci EXAMINATION. Vibrations and Waves. Date: 4 th May, Time: 14:30-17:00

BSc/MSci EXAMINATION. Vibrations and Waves. Date: 4 th May, Time: 14:30-17:00 BSc/MSci EXAMINATION PHY-217 Vibrations and Waves Time Allowed: 2 hours 30 minutes Date: 4 th May, 2011 Time: 14:30-17:00 Instructions: Answer ALL questions in section A. Answer ONLY TWO questions from

More information

5 Applying the Fokker-Planck equation

5 Applying the Fokker-Planck equation 5 Applying the Fokker-Planck equation We begin with one-dimensional examples, keeping g = constant. Recall: the FPE for the Langevin equation with η(t 1 )η(t ) = κδ(t 1 t ) is = f(x) + g(x)η(t) t = x [f(x)p

More information

Exercises Lecture 15

Exercises Lecture 15 AM1 Mathematical Analysis 1 Oct. 011 Feb. 01 Date: January 7 Exercises Lecture 15 Harmonic Oscillators In classical mechanics, a harmonic oscillator is a system that, when displaced from its equilibrium

More information

There is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth)

There is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth) 82 Introduction Liapunov Functions Besides the Liapunov spectral theorem, there is another basic method of proving stability that is a generalization of the energy method we have seen in the introductory

More information

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) 28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.

More information

Low-Resistant Band-Passing Noise and Its Dynamical Effects

Low-Resistant Band-Passing Noise and Its Dynamical Effects Commun. Theor. Phys. (Beijing, China) 48 (7) pp. 11 116 c International Academic Publishers Vol. 48, No. 1, July 15, 7 Low-Resistant Band-Passing Noise and Its Dynamical Effects BAI Zhan-Wu Department

More information

Linear Response and Onsager Reciprocal Relations

Linear Response and Onsager Reciprocal Relations Linear Response and Onsager Reciprocal Relations Amir Bar January 1, 013 Based on Kittel, Elementary statistical physics, chapters 33-34; Kubo,Toda and Hashitsume, Statistical Physics II, chapter 1; and

More information

F = ma, F R + F S = mx.

F = ma, F R + F S = mx. Mechanical Vibrations As we mentioned in Section 3.1, linear equations with constant coefficients come up in many applications; in this section, we will specifically study spring and shock absorber systems

More information

Effect of Non Gaussian Noises on the Stochastic Resonance-Like Phenomenon in Gated Traps. Abstract

Effect of Non Gaussian Noises on the Stochastic Resonance-Like Phenomenon in Gated Traps. Abstract Effect of Non Gaussian Noises on the Stochastic Resonance-Like Phenomenon in Gated Traps Jorge A. Revelli 1 Alejandro D. Sánchez 2, and Horacio S. Wio 1 1) Grupo de Física Estadística Centro Atómico Bariloche

More information

Application of Simple Harmonics Modeling a Shock

Application of Simple Harmonics Modeling a Shock Undergraduate Journal of Mathematical Modeling: One + Two Volume 8 2017 Fall 2017 Issue 1 Article 1 Application of Simple Harmonics Modeling a Shock Kai Raymond University of South Florida Advisors: Thomas

More information

MAT187H1F Lec0101 Burbulla

MAT187H1F Lec0101 Burbulla Spring 2017 Second Order Linear Homogeneous Differential Equation DE: A(x) d 2 y dx 2 + B(x)dy dx + C(x)y = 0 This equation is called second order because it includes the second derivative of y; it is

More information

THE WAVE EQUATION. F = T (x, t) j + T (x + x, t) j = T (sin(θ(x, t)) + sin(θ(x + x, t)))

THE WAVE EQUATION. F = T (x, t) j + T (x + x, t) j = T (sin(θ(x, t)) + sin(θ(x + x, t))) THE WAVE EQUATION The aim is to derive a mathematical model that describes small vibrations of a tightly stretched flexible string for the one-dimensional case, or of a tightly stretched membrane for the

More information

Effect of various periodic forces on Duffing oscillator

Effect of various periodic forces on Duffing oscillator PRAMANA c Indian Academy of Sciences Vol. 67, No. 2 journal of August 2006 physics pp. 351 356 Effect of various periodic forces on Duffing oscillator V RAVICHANDRAN 1, V CHINNATHAMBI 1, and S RAJASEKAR

More information

Section 8.5. z(t) = be ix(t). (8.5.1) Figure A pendulum. ż = ibẋe ix (8.5.2) (8.5.3) = ( bẋ 2 cos(x) bẍ sin(x)) + i( bẋ 2 sin(x) + bẍ cos(x)).

Section 8.5. z(t) = be ix(t). (8.5.1) Figure A pendulum. ż = ibẋe ix (8.5.2) (8.5.3) = ( bẋ 2 cos(x) bẍ sin(x)) + i( bẋ 2 sin(x) + bẍ cos(x)). Difference Equations to Differential Equations Section 8.5 Applications: Pendulums Mass-Spring Systems In this section we will investigate two applications of our work in Section 8.4. First, we will consider

More information

Path integrals for classical Markov processes

Path integrals for classical Markov processes Path integrals for classical Markov processes Hugo Touchette National Institute for Theoretical Physics (NITheP) Stellenbosch, South Africa Chris Engelbrecht Summer School on Non-Linear Phenomena in Field

More information

Damped & forced oscillators

Damped & forced oscillators SEISMOLOGY I Laurea Magistralis in Physics of the Earth and of the Environment Damped & forced oscillators Fabio ROMANELLI Dept. Earth Sciences Università degli studi di Trieste romanel@dst.units.it Damped

More information

KEELE UNIVERSITY PHYSICS/ASTROPHYSICS MODULE PHY OSCILLATIONS AND WAVES PRACTICE EXAM

KEELE UNIVERSITY PHYSICS/ASTROPHYSICS MODULE PHY OSCILLATIONS AND WAVES PRACTICE EXAM KEELE UNIVERSITY PHYSICS/ASTROPHYSICS MODULE PHY-10012 OSCILLATIONS AND WAVES PRACTICE EXAM Candidates should attempt ALL of PARTS A and B, and TWO questions from PART C. PARTS A and B should be answered

More information

Damped Oscillation Solution

Damped Oscillation Solution Lecture 19 (Chapter 7): Energy Damping, s 1 OverDamped Oscillation Solution Damped Oscillation Solution The last case has β 2 ω 2 0 > 0. In this case we define another real frequency ω 2 = β 2 ω 2 0. In

More information

How many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation?

How many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation? How many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation? (A) 0 (B) 1 (C) 2 (D) more than 2 (E) it depends or don t know How many of

More information

Oscillations. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring of constant k is

Oscillations. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring of constant k is Dr. Alain Brizard College Physics I (PY 10) Oscillations Textbook Reference: Chapter 14 sections 1-8. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring

More information

S13 PHY321: Final May 1, NOTE: Show all your work. No credit for unsupported answers. Turn the front page only when advised by the instructor!

S13 PHY321: Final May 1, NOTE: Show all your work. No credit for unsupported answers. Turn the front page only when advised by the instructor! Name: Student ID: S13 PHY321: Final May 1, 2013 NOTE: Show all your work. No credit for unsupported answers. Turn the front page only when advised by the instructor! The exam consists of 6 problems (60

More information

The correlation between stochastic resonance and the average phase-synchronization time of a bistable system driven by colour-correlated noises

The correlation between stochastic resonance and the average phase-synchronization time of a bistable system driven by colour-correlated noises Chin. Phys. B Vol. 19, No. 1 (010) 01050 The correlation between stochastic resonance and the average phase-synchronization time of a bistable system driven by colour-correlated noises Dong Xiao-Juan(

More information

Thursday, August 4, 2011

Thursday, August 4, 2011 Chapter 16 Thursday, August 4, 2011 16.1 Springs in Motion: Hooke s Law and the Second-Order ODE We have seen alrealdy that differential equations are powerful tools for understanding mechanics and electro-magnetism.

More information

Basic Theory of Dynamical Systems

Basic Theory of Dynamical Systems 1 Basic Theory of Dynamical Systems Page 1 1.1 Introduction and Basic Examples Dynamical systems is concerned with both quantitative and qualitative properties of evolution equations, which are often ordinary

More information

Solutions 2: Simple Harmonic Oscillator and General Oscillations

Solutions 2: Simple Harmonic Oscillator and General Oscillations Massachusetts Institute of Technology MITES 2017 Physics III Solutions 2: Simple Harmonic Oscillator and General Oscillations Due Wednesday June 21, at 9AM under Rene García s door Preface: This problem

More information

How to Implement Super-Twisting Controller based on Sliding Mode Observer?

How to Implement Super-Twisting Controller based on Sliding Mode Observer? How to Implement Super-Twisting Controller based on Sliding Mode Observer? Asif Chalanga 1 Shyam Kamal 2 Prof.L.Fridman 3 Prof.B.Bandyopadhyay 4 and Prof.J.A.Moreno 5 124 Indian Institute of Technology

More information

A Theory of Spatiotemporal Chaos: What s it mean, and how close are we?

A Theory of Spatiotemporal Chaos: What s it mean, and how close are we? A Theory of Spatiotemporal Chaos: What s it mean, and how close are we? Michael Dennin UC Irvine Department of Physics and Astronomy Funded by: NSF DMR9975497 Sloan Foundation Research Corporation Outline

More information

C.-H. Lamarque. University of Lyon/ENTPE/LGCB & LTDS UMR CNRS 5513

C.-H. Lamarque. University of Lyon/ENTPE/LGCB & LTDS UMR CNRS 5513 Nonlinear Dynamics of Smooth and Non-Smooth Systems with Application to Passive Controls 3rd Sperlonga Summer School on Mechanics and Engineering Sciences on Dynamics, Stability and Control of Flexible

More information

The Stochastic KdV - A Review

The Stochastic KdV - A Review The Stochastic KdV - A Review Russell L. Herman Wadati s Work We begin with a review of the results from Wadati s paper 983 paper [4] but with a sign change in the nonlinear term. The key is an analysis

More information

Faculty of Computers and Information. Basic Science Department

Faculty of Computers and Information. Basic Science Department 18--018 FCI 1 Faculty of Computers and Information Basic Science Department 017-018 Prof. Nabila.M.Hassan 18--018 FCI Aims of Course: The graduates have to know the nature of vibration wave motions with

More information

Math 221 Topics since the second exam

Math 221 Topics since the second exam Laplace Transforms. Math 1 Topics since the second exam There is a whole different set of techniques for solving n-th order linear equations, which are based on the Laplace transform of a function. For

More information

EE538 Final Exam Fall :20 pm -5:20 pm PHYS 223 Dec. 17, Cover Sheet

EE538 Final Exam Fall :20 pm -5:20 pm PHYS 223 Dec. 17, Cover Sheet EE538 Final Exam Fall 005 3:0 pm -5:0 pm PHYS 3 Dec. 17, 005 Cover Sheet Test Duration: 10 minutes. Open Book but Closed Notes. Calculators ARE allowed!! This test contains five problems. Each of the five

More information

Math 216 Second Midterm 20 March, 2017

Math 216 Second Midterm 20 March, 2017 Math 216 Second Midterm 20 March, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material

More information

MATH Test II

MATH Test II MATH 3113-008 Test II Dr. Darren Ong March 25, 2015 2:30pm-3:20pm Answer the questions in the spaces provided on the question sheets. If you run out of room for an answer, continue on the back of the page.

More information

Analytical Mechanics ( AM )

Analytical Mechanics ( AM ) Analytical Mechanics ( AM ) Olaf Scholten KVI, kamer v8; tel nr 6-55; email: scholten@kvinl Web page: http://wwwkvinl/ scholten Book: Classical Dynamics of Particles and Systems, Stephen T Thornton & Jerry

More information

Math Assignment 5

Math Assignment 5 Math 2280 - Assignment 5 Dylan Zwick Fall 2013 Section 3.4-1, 5, 18, 21 Section 3.5-1, 11, 23, 28, 35, 47, 56 Section 3.6-1, 2, 9, 17, 24 1 Section 3.4 - Mechanical Vibrations 3.4.1 - Determine the period

More information

Chapter 14 Periodic Motion

Chapter 14 Periodic Motion Chapter 14 Periodic Motion 1 Describing Oscillation First, we want to describe the kinematical and dynamical quantities associated with Simple Harmonic Motion (SHM), for example, x, v x, a x, and F x.

More information

P321(b), Assignement 1

P321(b), Assignement 1 P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates

More information

BSc/MSci MidTerm Test

BSc/MSci MidTerm Test BSc/MSci MidTerm Test PHY-217 Vibrations and Waves Time Allowed: 40 minutes Date: 18 th Nov, 2011 Time: 9:10-9:50 Instructions: Answer ALL questions in section A. Answer ONLY ONE questions from section

More information

Vibrational resonance

Vibrational resonance Published in J. Phys. A: Math. Gen. 33, L433 L438 (2000). LETTER TO THE EDITOR Vibrational resonance P S Landa andpvemcclintock Department of Physics, Lomonosov Moscow State University, 119899 Moscow,

More information

LINEAR RESPONSE THEORY

LINEAR RESPONSE THEORY MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior

More information

Fourier Series. Green - underdamped

Fourier Series. Green - underdamped Harmonic Oscillator Fourier Series Undamped: Green - underdamped Overdamped: Critical: Underdamped: Driven: Calculus of Variations b f {y, y'; x}dx is stationary when f y d f = 0 dx y' a Note that y is

More information

A TLM model of a heavy gantry crane system

A TLM model of a heavy gantry crane system A TLM model of a heavy gantry crane system William J O'Connor * Abstract The 1-D wave equation model of a vibrating string is first adjusted to allow for variable tension. The response is verified by comparison

More information

Math 240: Spring/Mass Systems II

Math 240: Spring/Mass Systems II Math 240: Spring/Mass Systems II Ryan Blair University of Pennsylvania Monday, March 26, 2012 Ryan Blair (U Penn) Math 240: Spring/Mass Systems II Monday, March 26, 2012 1 / 12 Outline 1 Today s Goals

More information

Two-Body Problem. Central Potential. 1D Motion

Two-Body Problem. Central Potential. 1D Motion Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of

More information

Chapter 3. Periodic functions

Chapter 3. Periodic functions Chapter 3. Periodic functions Why do lights flicker? For that matter, why do they give off light at all? They are fed by an alternating current which turns into heat because of the electrical resistance

More information

New ideas in the non-equilibrium statistical physics and the micro approach to transportation flows

New ideas in the non-equilibrium statistical physics and the micro approach to transportation flows New ideas in the non-equilibrium statistical physics and the micro approach to transportation flows Plenary talk on the conference Stochastic and Analytic Methods in Mathematical Physics, Yerevan, Armenia,

More information

1 Lyapunov theory of stability

1 Lyapunov theory of stability M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability

More information

arxiv: v1 [nlin.cd] 22 Aug 2016

arxiv: v1 [nlin.cd] 22 Aug 2016 The Two Dynamical States of a Driven Underdamped Classical Pendulum: An Analog Simulation Experiment Ivan Skhem Sawkmie and Mangal C. Mahato Department of Physics, North-Eastern Hill University, Shillong-793022,

More information

1 Pushing your Friend on a Swing

1 Pushing your Friend on a Swing Massachusetts Institute of Technology MITES 017 Physics III Lecture 05: Driven Oscillations In these notes, we derive the properties of both an undamped and damped harmonic oscillator under the influence

More information

Correlation times in stochastic equations with delayed feedback and multiplicative noise. Abstract

Correlation times in stochastic equations with delayed feedback and multiplicative noise. Abstract Correlation times in stochastic equations with delayed feedback and multiplicative noise Mathieu Gaudreault 1, Juliana Militão Berbert 1,2, and Jorge Viñals 1 1 Department of Physics, McGill University,

More information

Inertial Game Dynamics

Inertial Game Dynamics ... Inertial Game Dynamics R. Laraki P. Mertikopoulos CNRS LAMSADE laboratory CNRS LIG laboratory ADGO'13 Playa Blanca, October 15, 2013 ... Motivation Main Idea: use second order tools to derive efficient

More information

Physics III: Final Solutions (Individual and Group)

Physics III: Final Solutions (Individual and Group) Massachusetts Institute of Technology MITES 7 Physics III First and Last Name: Physics III: Final Solutions (Individual and Group Instructor: Mobolaji Williams Teaching Assistant: Rene García (Tuesday

More information

Transduction Based on Changes in the Energy Stored in an Electrical Field

Transduction Based on Changes in the Energy Stored in an Electrical Field Lecture 6- Transduction Based on Changes in the Energy Stored in an Electrical Field Actuator Examples Microgrippers Normal force driving In-plane force driving» Comb-drive device F = εav d 1 ε oε F rwv

More information

08. Brownian Motion. University of Rhode Island. Gerhard Müller University of Rhode Island,

08. Brownian Motion. University of Rhode Island. Gerhard Müller University of Rhode Island, University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 1-19-215 8. Brownian Motion Gerhard Müller University of Rhode Island, gmuller@uri.edu Follow this

More information

S12 PHY321: Practice Final

S12 PHY321: Practice Final S12 PHY321: Practice Final Contextual information Damped harmonic oscillator equation: ẍ + 2βẋ + ω0x 2 = 0 ( ) ( General solution: x(t) = e [A βt 1 exp β2 ω0t 2 + A 2 exp )] β 2 ω0t 2 Driven harmonic oscillator

More information

Simulation of DART Buoy Data

Simulation of DART Buoy Data Simulation of DART Buoy Data goal: develop procedure for simulating detided DART buoy data that can be used to test inversion algorithm, train users etc. simulated data has three stochastic components:

More information

Turbulent drag reduction by streamwise traveling waves

Turbulent drag reduction by streamwise traveling waves 51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For

More information

Section 3.7: Mechanical and Electrical Vibrations

Section 3.7: Mechanical and Electrical Vibrations Section 3.7: Mechanical and Electrical Vibrations Second order linear equations with constant coefficients serve as mathematical models for mechanical and electrical oscillations. For example, the motion

More information

= U / x. That is, U d = m x x dt FRIDAY'S QUIZ QUESTION AN EXAMPLE IN THE CALCULUS OF VARIATIONS

= U / x. That is, U d = m x x dt FRIDAY'S QUIZ QUESTION AN EXAMPLE IN THE CALCULUS OF VARIATIONS FRIDAY'S QUIZ QUESTION AN EXAMPLE IN THE CALCULUS OF VARIATIONS GIVEN S = [ ½ m x 2 U(x) ] dt, find the x(t) that makes S[x(t)] stationary. SOLUTION It's an Euler problem, so the solution is given by Euler's

More information

Lagrangian Data Assimilation and Its Application to Geophysical Fluid Flows

Lagrangian Data Assimilation and Its Application to Geophysical Fluid Flows Lagrangian Data Assimilation and Its Application to Geophysical Fluid Flows Laura Slivinski June, 3 Laura Slivinski (Brown University) Lagrangian Data Assimilation June, 3 / 3 Data Assimilation Setup:

More information

Dynamics of two coupled particles: comparison of Lennard-Jones and spring forces

Dynamics of two coupled particles: comparison of Lennard-Jones and spring forces Dynamics of two coupled particles: comparison of Lennard-Jones and spring forces Renata Retkute a and James P. Gleeson a a Department of Applied Mathematics, UCC, Cork, Ireland; ABSTRACT The motion of

More information

Ratchet Baryogenesis during Reheating

Ratchet Baryogenesis during Reheating Ratchet Baryogenesis during Reheating Neil D. Barrie Kavli IPMU (WPI) February 17, 2018 K. Bamba, N. B., A. Sugamoto, T. Takeuchi and K. Yamashita, arxiv:1610.03268, and upcoming work. N. D. Barrie (Kavli

More information

Matched Filtering for Generalized Stationary Processes

Matched Filtering for Generalized Stationary Processes Matched Filtering for Generalized Stationary Processes to be submitted to IEEE Transactions on Information Theory Vadim Olshevsky and Lev Sakhnovich Department of Mathematics University of Connecticut

More information

PHY 3221 Fall Homework Problems. Instructor: Yoonseok Lee. Submit only HW s. EX s are additional problems that I encourage you to work on.

PHY 3221 Fall Homework Problems. Instructor: Yoonseok Lee. Submit only HW s. EX s are additional problems that I encourage you to work on. PHY 3221 Fall 2012 Homework Problems Instructor: Yoonseok Lee Submit only HW s. EX s are additional problems that I encourage you to work on. Week 1: August 22-24, Due August 27 (nothing to submit) EX:

More information

Numerical investigation of noise-induced resonance in a semiconductor laser with optical feedback

Numerical investigation of noise-induced resonance in a semiconductor laser with optical feedback Physica D 168 169 (2002) 171 176 Numerical investigation of noise-induced resonance in a semiconductor laser with optical feedback C. Masoller Instituto de Física, Facultad de Ciencias, Universidad de

More information

Symmetries 2 - Rotations in Space

Symmetries 2 - Rotations in Space Symmetries 2 - Rotations in Space This symmetry is about the isotropy of space, i.e. space is the same in all orientations. Thus, if we continuously rotated an entire system in space, we expect the system

More information

SDE Coefficients. March 4, 2008

SDE Coefficients. March 4, 2008 SDE Coefficients March 4, 2008 The following is a summary of GARD sections 3.3 and 6., mainly as an overview of the two main approaches to creating a SDE model. Stochastic Differential Equations (SDE)

More information

Introduction to Fluctuation Theorems

Introduction to Fluctuation Theorems Hyunggyu Park Introduction to Fluctuation Theorems 1. Nonequilibrium processes 2. Brief History of Fluctuation theorems 3. Jarzynski equality & Crooks FT 4. Experiments 5. Probability theory viewpoint

More information

Experimental Fourier Transforms

Experimental Fourier Transforms Chapter 5 Experimental Fourier Transforms 5.1 Sampling and Aliasing Given x(t), we observe only sampled data x s (t) = x(t)s(t; T s ) (Fig. 5.1), where s is called sampling or comb function and can be

More information

Wave Physics PHYS 2023 Tim Freegarde

Wave Physics PHYS 2023 Tim Freegarde Wave Physics PHYS 2023 Tim Freegarde Wave Physics WAVE EQUATIONS & SINUSOIDAL SOLUTIONS WAVE PROPAGATION BEHAVIOUR AT INTERFACES SUPERPOSITIONS FURTHER TOPICS general wave phenomena wave equations, derivations

More information

The Harmonic Oscillator

The Harmonic Oscillator The Harmonic Oscillator Math 4: Ordinary Differential Equations Chris Meyer May 3, 008 Introduction The harmonic oscillator is a common model used in physics because of the wide range of problems it can

More information

Fourier Transform in Image Processing. CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012)

Fourier Transform in Image Processing. CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012) Fourier Transform in Image Processing CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012) Basis Decomposition Write a function as a weighted sum of basis functions f ( x) wibi(

More information

a. What is the angular frequency ω of the block in terms of k, l, and m?

a. What is the angular frequency ω of the block in terms of k, l, and m? 1 Problem 1: (4 pts.) Two spherical planets, each of mass M and Radius R, start out at rest with a distance from center to center of 4R. What is the speed of one of the planets at the moment that their

More information

Fourier Transform for Continuous Functions

Fourier Transform for Continuous Functions Fourier Transform for Continuous Functions Central goal: representing a signal by a set of orthogonal bases that are corresponding to frequencies or spectrum. Fourier series allows to find the spectrum

More information

Goals: Second-order Linear Equations Linear Independence of Solutions and the Wronskian Homogeneous DEs with Constant Coefficients

Goals: Second-order Linear Equations Linear Independence of Solutions and the Wronskian Homogeneous DEs with Constant Coefficients Week #3 : Higher-Order Homogeneous DEs Goals: Second-order Linear Equations Linear Independence of Solutions and the Wronskian Homogeneous DEs with Constant Coefficients 1 Second-Order Linear Equations

More information

Damped Harmonic Oscillator

Damped Harmonic Oscillator Damped Harmonic Oscillator Wednesday, 23 October 213 A simple harmonic oscillator subject to linear damping may oscillate with exponential decay, or it may decay biexponentially without oscillating, or

More information

Noise, AFMs, and Nanomechanical Biosensors

Noise, AFMs, and Nanomechanical Biosensors Noise, AFMs, and Nanomechanical Biosensors: Lancaster University, November, 2005 1 Noise, AFMs, and Nanomechanical Biosensors with Mark Paul (Virginia Tech), and the Caltech BioNEMS Collaboration Support:

More information

Uncertainty quantification and calibration of computer models. February 5th, 2014

Uncertainty quantification and calibration of computer models. February 5th, 2014 Uncertainty quantification and calibration of computer models February 5th, 2014 Physical model Physical model is defined by a set of differential equations. Now, they are usually described by computer

More information

PHY217: Vibrations and Waves

PHY217: Vibrations and Waves Assessed Problem set 1 Issued: 5 November 01 PHY17: Vibrations and Waves Deadline for submission: 5 pm Thursday 15th November, to the V&W pigeon hole in the Physics reception on the 1st floor of the GO

More information

Runga-Kutta Schemes. Exact evolution over a small time step: Expand both sides in a small time increment: d(δt) F (x(t + δt),t+ δt) Ft + FF x ( t)

Runga-Kutta Schemes. Exact evolution over a small time step: Expand both sides in a small time increment: d(δt) F (x(t + δt),t+ δt) Ft + FF x ( t) Runga-Kutta Schemes Exact evolution over a small time step: x(t + t) =x(t)+ t 0 d(δt) F (x(t + δt),t+ δt) Expand both sides in a small time increment: x(t + t) =x(t)+x (t) t + 1 2 x (t)( t) 2 + 1 6 x (t)+

More information

Newton's Laws You should be able to state these laws using both words and equations.

Newton's Laws You should be able to state these laws using both words and equations. Review before first test Physical Mechanics Fall 000 Newton's Laws You should be able to state these laws using both words and equations. The nd law most important for meteorology. Second law: net force

More information

Mesoscale Simulation Methods. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai

Mesoscale Simulation Methods. Ronojoy Adhikari The Institute of Mathematical Sciences Chennai Mesoscale Simulation Methods Ronojoy Adhikari The Institute of Mathematical Sciences Chennai Outline What is mesoscale? Mesoscale statics and dynamics through coarse-graining. Coarse-grained equations

More information

Lecture 2: Open quantum systems

Lecture 2: Open quantum systems Phys 769 Selected Topics in Condensed Matter Physics Summer 21 Lecture 2: Open quantum systems Lecturer: Anthony J. Leggett TA: Bill Coish 1. No (micro- or macro-) system is ever truly isolated U = S +

More information

Contrasting measures of irreversibility in stochastic and deterministic dynamics

Contrasting measures of irreversibility in stochastic and deterministic dynamics Contrasting measures of irreversibility in stochastic and deterministic dynamics Ian Ford Department of Physics and Astronomy and London Centre for Nanotechnology UCL I J Ford, New J. Phys 17 (2015) 075017

More information

Physical Acoustics. Hearing is the result of a complex interaction of physics, physiology, perception and cognition.

Physical Acoustics. Hearing is the result of a complex interaction of physics, physiology, perception and cognition. Physical Acoustics Hearing, auditory perception, or audition is the ability to perceive sound by detecting vibrations, changes in the pressure of the surrounding medium through time, through an organ such

More information

Radiation Damping. 1 Introduction to the Abraham-Lorentz equation

Radiation Damping. 1 Introduction to the Abraham-Lorentz equation Radiation Damping Lecture 18 1 Introduction to the Abraham-Lorentz equation Classically, a charged particle radiates energy if it is accelerated. We have previously obtained the Larmor expression for the

More information

Physics Homework 1

Physics Homework 1 Physics 201 - Homewor 1 Ben Coy January 24, 2011 13.11 We are given the speed at the equilibrium position as 12ms 1. This is the maximum speed and is related to the amplitude by v max = Aω So we need to

More information

4.9 Free Mechanical Vibrations

4.9 Free Mechanical Vibrations 4.9 Free Mechanical Vibrations Spring-Mass Oscillator When the spring is not stretched and the mass m is at rest, the system is at equilibrium. Forces Acting in the System When the mass m is displaced

More information

5.74 Introductory Quantum Mechanics II

5.74 Introductory Quantum Mechanics II MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 9 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,

More information

PHYSICS 110A : CLASSICAL MECHANICS HW 2 SOLUTIONS. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see

PHYSICS 110A : CLASSICAL MECHANICS HW 2 SOLUTIONS. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see PHYSICS 11A : CLASSICAL MECHANICS HW SOLUTIONS (1) Taylor 5. Here is a sketch of the potential with A = 1, R = 1, and S = 1. From the plot we can see 1.5 1 U(r).5.5 1 4 6 8 1 r Figure 1: Plot for problem

More information

Chaotic motion. Phys 420/580 Lecture 10

Chaotic motion. Phys 420/580 Lecture 10 Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t

More information

Simple Harmonic Motion

Simple Harmonic Motion Simple Harmonic Motion (FIZ 101E - Summer 2018) July 29, 2018 Contents 1 Introduction 2 2 The Spring-Mass System 2 3 The Energy in SHM 5 4 The Simple Pendulum 6 5 The Physical Pendulum 8 6 The Damped Oscillations

More information

L = 1 2 a(q) q2 V (q).

L = 1 2 a(q) q2 V (q). Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion

More information

Chapter 14 Oscillations. Copyright 2009 Pearson Education, Inc.

Chapter 14 Oscillations. Copyright 2009 Pearson Education, Inc. Chapter 14 Oscillations Oscillations of a Spring Simple Harmonic Motion Energy in the Simple Harmonic Oscillator Simple Harmonic Motion Related to Uniform Circular Motion The Simple Pendulum The Physical

More information

Collective and Stochastic Effects in Arrays of Submicron Oscillators

Collective and Stochastic Effects in Arrays of Submicron Oscillators DYNAMICS DAYS: Long Beach, 2005 1 Collective and Stochastic Effects in Arrays of Submicron Oscillators Ron Lifshitz (Tel Aviv), Jeff Rogers (HRL, Malibu), Oleg Kogan (Caltech), Yaron Bromberg (Tel Aviv),

More information