Phase transition theory of pulse formation in passively mode-locked lasers with dispersion and Kerr nonlinearity
|
|
- Warren Summers
- 5 years ago
- Views:
Transcription
1 Optics Counications 223 (2003) Phase transition theory of pulse foration in passively ode-locked lasers with dispersion and Kerr nonlinearity Ariel Gordon, Baruch Fischer * Departent of Electrical Engineering, Technion-Israel Institute of Technology, Haifa 32000, Israel Received 19February 2003; accepted 30 May 2003 Abstract We extend our statistical echanics approach (Phys. Rev. Lett. 89(2002) ) to passively ode-locked lasers to also include dispersion and the Kerr nonlinearity and in particular the soliton laser case. The theory predicts for these general cases a discontinuous first-order phase transition at the onset of pulses that explains the power threshold behavior of ode-locked lasers. Ó 2003 Elsevier Science B.V. All rights reserved. PACS: Re; Fc; Mi; Wd 1. Introduction In the course of the recent three decades, uch progress has been ade in the study of passive ode locking (PML) of lasers. The analytical part of the studies has led to a rather consensual theory [2,3]. There are any nuances of this theory, but in ost of the conteporary analytical descriptions of PML the core reains the sae: an equation of otion (the aster equation), reiniscent to the coplex Landau Ginzburg equation, for the evolution of optical pulses, odeling the action of gain dispersion, group velocity dispersion (GVD), a usually fast saturable absorber, the Kerr nonlinearity and * Corresponding author. Tel.: ; fax: E-ail address: fischer@ee.technion.ac.il (B. Fischer). slow saturable gain. This faily of theories has been used for extensive study of PML: the pulse shape [4], stability issues [5], the threshold behavior [6], noise characteristics [7] of the pulses and any ore. A recent review is given in [3]. Most of the existing theories of the PML threshold behavior [6,8 11] consider noise or soe sort of incoherence to be responsible for it. When the effect of noise on PML is studied through the Langevin aster equation, noise is often assued to be sall, which facilitates linearization of the Langevin aster equation [7]. This way, the lowest oents describing the statistics of the output signal of PML lasers can be calculated. As we find below, crucial features of the syste can be issed in this procedure when the noise is large. In a recent work [1] we presented a statistical echanics approach for a any interacting ode laser syste. We have shown that such systes /03/$ - see front atter Ó 2003 Elsevier Science B.V. All rights reserved. doi: /s (03)
2 152 A. Gordon, B. Fischer / Optics Counications 223 (2003) exhibit phase transition of the first kind, a result of the interplay between entropy and (a generalization of) energy in a syste of any interacting odes, where the role of teperature is taken by noise. Contrary to sall noise in the laser, which causes noise in its output, increasing the noise (or equivalently decreasing the signal) to a certain critical value causes an abrupt breakdown of the pulses and a transition to a noisy CW operation. The pulse foration in lasers is shown to be siilar to elting freezing behavior in solid liquid and paraferroagnet transitions. We have recently experientally deonstrated the predicted first-order phase transition in a passively ode-locked fiber laser [12], showing the therodynaic elting freezing path. In this paper we extend the study of the threshold behavior of PML [1] to include GVD, the Kerr nonlinearity and a ore coon description of the optical saturable gain. We extend the Hailtonian foralis [1] to include gain saturation. We discuss the purely dissipative otion, the otion induced by the gain and losses (linear and nonlinear) alone, and present a dynaical picture of the balance between forces and noise, reflecting the balance between interaction and entropy. We then write the Fokker Planck equation associated with the Langevin equation of otion and look for steady-state solutions which are functions of the Hailtonians alone. We find such solutions when the soliton condition is satisfied [3], which extends our theory fro [1] to this iportant case. Beyond the soliton condition we investigate the Langevin equations of otion nuerically, for general values of dispersion and the Kerr nonlinearity. We find that the equations consistently exhibit the first-order phase transition behavior. We suggest this echanis as an explanation to the passive ode-locking threshold behavior. 2. Hailtonian forulation of the aster equation The coon equation of otion [3] of PML (the aster equation ), which is usually written in the tie doain, can be written in the spectral doain as follows: _a ¼ðg l þ i/ Þa þðic k c s Þ X a j a k a l þ C : j kþl ¼0 ð1þ The boundaries of the cavity adit only discrete axial odes into the laser, the slowly varying aplitudes of which are denoted by a. _a stands for da =ds, where s is a diensionless variable easuring tie in units of the cavity roundtrip tie. g denotes the net gain at the wavelength corresponding to and is signal dependent. / represents the phase acquired by the th ode after a roundtrip, which in particular describes dispersion. l represents the wavelength-dependent linear loss. c s and c k represent the nonlinear self-aplitude (saturable absorber) and self-phase odulation (Kerr nonlinearity) coefficients, and C is a white Gaussian Langevin additive noise ter which in particular describes quantu noise [7]. We assue the CÕs tobe Gaussian, white and uncorrelated C ðt 1 ÞC n ðt 2Þ ¼ 2T dn dðt 1 t 2 Þ; ð2þ hc ðt 1 ÞC n ðt 2 Þi ¼ 0 with hi denoting enseble averaging. T is the spectral power of the noise, which functions as teperature. It is well known that lasers and in particular PML lasers are stable due to gain saturation [5]. When the response of the laser aplifier and the evolution of the wavefor are uch slower than the cavity roundtrip tie, the gain at a given oent is deterined by the whole wavefor all over the cavity. The ost coon odel for gain saturation in that case is g 0 g ¼ ; ð3þ 1 þ P=P sat where P ¼ X a a ð4þ and P sat is the saturation power of the aplifier. Eq. (3) describes a wavelength-independent aplifier: the filtering action of the aplifier is ebedded in l and its signal dependence has been neglected. This coon approxiation is often satisfying [5], but wavelength-dependent saturation behavior is easily incorporated in our foralis by defining
3 A. Gordon, B. Fischer / Optics Counications 223 (2003) g 0 g ¼ ; ð5þ 1 þ G=G sat where G ¼ X g 0 a a : ð6þ Eq. (5) reflects the fact that in a hoogeneously broadened aplifier the gain saturation is deterined by the overlap between the optical spectru and the gain profile [13]. G sat is the analog of the saturation power for that case, and can be related to its fundaental properties [13]. Defining H R ¼ X / a a c X k a j a k 2 a la ð7þ j kþl ¼0 and H I ¼ X l a a c X k a j a k 2 a la j kþl ¼0 G sat lnðg sat þ GÞ; ð8þ it is straight forward to verify that (1) can be rewritten as _a ¼ oðh R ih I Þ þ C ioa : ð9þ Changing variables to the real and iaginary parts of a (a ¼ a R þ iai ) we have _a R ¼ oh R oa I _a I ¼ oh R oa R oh I oa R oh I oa I þ C R ; þ C I ; ð10þ C R and CI are the real and iaginary parts of C. Observing Eq. (10), we see clearly that H R generates a Hailtonian otion and H I generates a purely dissipative otion (a gradient flow ). Eqs. (10) are Langevin equations describing a rando process, which can be equivalently described by the associated Fokker Planck equation [14] _q ¼ X þ X o oa R þ T X oh R o oa I oa I oh R o 2 oa R2 oh I q oa R oa R þ o2 oa I2 oh I q oa I! q; ð11þ where q is the statistical distribution of the odes a and T is half the spectral power of the Langevin force. As shown in [1], when H R ¼ 0 it is straight forward to see that the distribution qða 1 ;...; a N Þ/e H I=T ð12þ is the exact solution of Eq. (11) with _q ¼ 0 (the steady-state solution). The solution (12) has far reaching consequences, since it is identical by structure to the statistical distribution of a canonical enseble. We elaborate this case in the next section (Section 3). The ore general cases that include the Kerr nonlinearity and dispersion (H R 6¼ 0) are discussed afterwards (in Sections 4 and 5). 3. Purely dissipative otion (H R =0) For the case of zero GVD and no Kerr nonlinearity (H R ¼ 0) the dynaics of the syste is gradient flow (_a R ¼ oh I=oa R, _ai ¼ oh I=oa I ). This is useful for additional intuitive understanding of the dynaics induced by the saturable absorber, gain saturation and filtering. Let us exaine the evolution of H I without noise and for H R ¼ 0. We have then " 2 þ dh I ds ¼ X oh I a R which eans that H I oh # 2 I a I ; ð13þ keeps decreasing until a stationary point of H I is reached. If the stationary point is not a iniu, then a sall perturbation will push H I down the hill again. If it is a local iniu, then the point is locally stable: after a sall perturbation the syste will return to the iniu, but a large one will throw the syste away without necessarily going back. Only a global iniu of H I will be a truly stable solution to Eq. (1). It follows therefore that a truly stable solution to H I will exist if and only if H I is bounded fro below. This provides a convenient tool for stability analysis, and a atheatical ebodient to the coonly used intuitive arguents on what lasers like and dislike : lasers tend to iniize H I.
4 154 A. Gordon, B. Fischer / Optics Counications 223 (2003) An iediate result of this stability analysis technique is that the aster equation adits no globally stable solutions: it is easy to see fro Eq. (8) that H I is not bounded fro below. The reason for this is that in the siplest odel of absorption saturation we used, which takes into account only its the lowest nonlinear ter, increasing the optical power enough will eventually turn the absorber into an aplifier, which is physically incorrect. This proble can be fixed by adding to the aster equation a stabilizing ter which will turn the local iniu of the odelocked configuration into a global one. An exaple for such a ter is the next one in the expansion of the saturable absorber. When noise is present, along with sliding down, the syste will be randoly thrown fro state to state. This rando otion is likely to draw the syste away fro rare states (like the iniu of H I ), or in other words, noise tends to increase the entropy of the syste. As entioned above, just like in statistical echanics, this copetition between iniizing H I and axiizing entropy leads to the Gibbs distribution, with the noise power taking the role of teperature. Contrary to linearization techniques, which can provide inforation on the few lowest oents only, Eq. (12) is an exact solution for the entire distribution qða 1 ;...; a N Þ/e H I=T of the wavefors (or odes) in the laser. Needless to say equipped with an exact solution for the distribution one can thoroughly study the noise properties of PML lasers. More interesting is that the threshold behavior of PML is ebedded in Eq. (12): order (ode locking) eerges in the laser upon cooling. This can be shown analytically too, through the ean field approxiation [1]. Eq. (12) establishes an intiate relation between PML lasers and equilibriu statistical echanics and phase transition theory, and indeed it can be shown that the occurrence of ode locking upon reaching threshold is a phase transition of the first kind. The results provided by statistical echanics, both through nuerical siulation (a Monte-Carlo siulation) and the ean field approxiation are shown in Fig. 1. The connection to phase transition theory and statistical echanics is not accidental. In fact, the situation in the PML laser is very siilar to the one in traditional systes exhibiting order disorder transitions like elting solids, ferroagnets and any ore. Just like the latter, the laser is a syste of any interacting degrees of freedo the odes, with an ordering interaction the nonlinear saturable absorber, and with teperature (which in our case is replaced by the noise variance) encouraging gain of entropy which opposes the ordering interaction. The theory of phase transitions is a ature theory with any well-established powerful results, which can be now directly applied to PML. One of the interesting consequences of our analysis is that whenever a saturable absorber is present (for zero GVD and no Kerr nonlinearity) CW is not a stable solution of Eq. (1): it is easy to show that a CW solution is a axiu or a saddle point of H I. As we have seen, what stabilizes a CW solution is noise alone. This reinds other situations in nature, for exaple the liquid phase, which reains liquid due to entropy considerations, although its energy would be lower had it solidified, and nevertheless it reains liquid. 4. Analytical solution for a case with H R 6¼ 0 including the soliton laser case Let us now proceed to H R 6¼ 0, that is, a laser with nonzero dispersion and Kerr coefficient. Inspired by equilibriu statistical echanics, where a very fundaental stateent is that the statistical distribution is a function of the Hailtonian only, and even ore by the solution (12), we can try to seek for a solution to the steady-state distribution which is again a function of the Hailtonians only. We therefore substitute the solution qða 1 ;...; a N Þ/fðH R ða 1 ;...; a N Þ; H I ða 1 ;...; a N ÞÞ ð14þ into Eq. (11) and equating to zero leads via a lengthy but trivial calculation to the result that Eq. (12) is the only possible solution of the for (14), and it solves Eq. (11) as long as
5 A. Gordon, B. Fischer / Optics Counications 223 (2003) Fig. 1. The order paraeter Q, inversely proportional to the pulse duty cycle, as function of the generalized teperature T for several choices of the different paraeters. Graphs 1 4, which alost coincide, are results of nuerical siulations: A Monte-Carlo siulation of Eq. (12) and a direct siulation of the Langevin equation (10), with (graphs 3 and 4) or without (graphs 1 and 2) the Kerr nonlinearity. The thick lines are their analytical ean field based approxiations [12]. Graphs 5 10 siulations of Eq. (10) with GVD, and graphs 11 and 12 are the sae siulations with both GVD and the Kerr effect. Qualitatively the behavior of graphs 5 12 is very siilar to that of graphs 1 4, all exhibiting an abrupt transition to a ode-locked phase. In the siulations we used 100 discrete odes (a 1 ; a 2...; a 100 ), with equal spectral spacing. The stabilization we used is a constraint of constant power P ¼ P ja k j 2. T is in units of c s P 2. Graphs 3 and 4 are with c k =c s ¼ 2 and 3. In graphs 5 12 we took a parabolic dispersion relation for GVD / ¼ c d ðn 50:5Þ 2 with c d ¼ 0:002; 0:003; 0:004; 0:006; 0:009; 0:02 for graphs 5; 6; 7; 8; 9; 10; respectively. In graphs 11 and 12 c k ¼ 3c s and c d is )0.005 and 0.005, respectively. X oh R oa I oh I oa R oh R oa R oh I ¼ 0; oa I ð15þ that is, H I is a constant of otion under the evolution generated by H R alone. Putting the coon parabolic profiles for dispersion and gain dispersion and restricting ourselves to Eq. (3) for odeling gain saturation, it is straight forward to show that Eq. (15) holds in particular for solitons [3]. The right balance between dispersion and the Kerr nonlinearity is kept. The conclusion of the above calculation is twofold. First, we have seen that Eq. (12) and hence all the subsequent analysis [1] holds as it is in a soliton laser. Second, we have seen that the steady-state distribution is not a function of the Hailtonians alone when dispersion and Kerr are nonzero and the soliton condition does not hold. It is therefore an explicit function of the aõs, which is difficult to find and presuably also to analyze. We have therefore turned to nuerical studies. 5. Nuerical study of the general H R 6¼ 0 case The ain question we were addressing in the nuerical study is whether the sharp breakdown of pulses upon increasing noise occurs with the ost general choice of the paraeters in Eq. (1). We have therefore conducted tens of nuerical experients, where we nuerically solved Eq. (1) and calculated the tie average of the noralized order paraeter P j kþl ¼0 ~Q ¼ a ja k a la 2 ð16þ P ja j 2
6 156 A. Gordon, B. Fischer / Optics Counications 223 (2003) as function of the noise level T. This paraeter is inversely proportional to the pulsewidth and directly proportional to the nuber of correlated odes: In the ode-locked phase it is of order of the total nuber of odes in the band, whereas in the CW phase it is of order one. All the siulations presented in this paper were perfored with 100 odes and with the constraint of constant optical power P 0 as explained above. All the nuerical experients exhibited the sae sharp breakdown of the pulses, however the critical T depends on the paraeters. It sees that pulses with higher chirp paraeter are less stable against noise, which is siilar to the behavior of stability against sall perturbations [3]. Additionally the general effect of the unbalanced GVD and Kerr nonlinearity is to reduce the effective nuber of locked odes. Soe results are presented in Fig. 1. With nonzero GVD and Kerr nonlinearity, the behavior of a PML laser with noise is qualitatively the sae as PML in the purely dissipative case (H R ¼ 0), and in particular the phase transition occurs, at least for a wide variety of paraeters. This is not very surprising, since a phase transition of the first kind is a draatic singularity, and it is not expected to be reoved by slight odifications of the dynaics. 6. Conclusion We have extended our statistical echanics approach for any interacting ode systes to general PML lasers that include dispersion and Kerr nonlinearity. Such systes were shown to exhibit a first-order transition. For a specific balance between dispersion and Kerr nonlinearity, that includes the iportant case of soliton lasers, the proof was given by an exact analysis. For the broader case we used a nuerical study. We also ention again our recent [12] experiental deonstration of the phase transition with a variety of therodynaic properties of pulse foration in a fiber PML laser. We have shown the elting freezing transition as the teperature noise was varied. The present paper can provide a further understanding of such systes when solitons are incorporated, as often happens in fiber lasers. We end with an experiental reark. We have analyzed the behavior of PML lasers as the noise power is altered. In an experient, the odelocking threshold is usually observed as the power in the cavity is raised. The interaction energy responsible for the ode locking is quadratic with the power in the cavity (see the last ter in Eq. (8)), while the noise power is close to be linear with it. This explains the existence of a threshold power. Acknowledgeents We thank Shuel Fishan for fruitful discussions. We acknowledge the support of the Division of Research Funds of the Israeli Ministry of Science, and the Fund for Prootion of Research at the Technion. References [1] A. Gordon, B. Fischer, Phys. Rev. Lett. 89(2002) [2] H.A. Haus, J. Appl. Phys. 46 (1975) [3] H.A. Haus, IEEE J. Sel. Top. Quant. 6 (2000) [4] H.A. Haus, C.V. Shank, E.P. Ippen, Opt. Coun. 15 (1) (1975) 29. [5] C.J. Chen, P.K. Wai, C.R. Menyuk, Opt. Lett. 19(3) (1994) 198. [6] C.J. Chen, P.K. Wai, C.R. Menyuk, Opt. Lett. 20 (4) (1995) 350. [7] H.A. Haus, A. Mecozzi, IEEE J. Quantu Electron. 29(3) (1993) 983. [8] F. Krausz, T. Brabec, Ch. Spilann, Opt. Lett. 16 (1991) 235. [9] J. Herrann, Opt. Coun. 98 (1993) 111. [10] E.P. Ippen, L.Y. Liu, H.A. Haus, Opt. Lett. 15 (1990) 737. [11] H.A. Haus, E.P. Ippen, Opt. Lett. 16 (1991) [12] A. Gordon, B. Vodonos, V. Sulakovski, B. Fischer, unpublished. [13] M. Horowitz, C.R. Menyuk, S. Keren, IEEE Photon Technol. Lett. 11 (10) (1999) [14] H. Risken, The Fokker Planck Equation, second ed., Springer, Berlin, 1989.
2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all
Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either
More informationNumerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term
Nuerical Studies of a Nonlinear Heat Equation with Square Root Reaction Ter Ron Bucire, 1 Karl McMurtry, 1 Ronald E. Micens 2 1 Matheatics Departent, Occidental College, Los Angeles, California 90041 2
More informationClassical systems in equilibrium
35 Classical systes in equilibriu Ideal gas Distinguishable particles Here we assue that every particle can be labeled by an index i... and distinguished fro any other particle by its label if not by any
More informationANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER
IEPC 003-0034 ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER A. Bober, M. Guelan Asher Space Research Institute, Technion-Israel Institute of Technology, 3000 Haifa, Israel
More informationma x = -bv x + F rod.
Notes on Dynaical Systes Dynaics is the study of change. The priary ingredients of a dynaical syste are its state and its rule of change (also soeties called the dynaic). Dynaical systes can be continuous
More informationReading from Young & Freedman: For this topic, read the introduction to chapter 25 and sections 25.1 to 25.3 & 25.6.
PHY10 Electricity Topic 6 (Lectures 9 & 10) Electric Current and Resistance n this topic, we will cover: 1) Current in a conductor ) Resistivity 3) Resistance 4) Oh s Law 5) The Drude Model of conduction
More informationδ 12. We find a highly accurate analytic description of the functions δ 11 ( δ 0, n)
Coplete-return spectru for a generalied Rosen-Zener two-state ter-crossing odel T.A. Shahverdyan, D.S. Mogilevtsev, V.M. Red kov, and A.M Ishkhanyan 3 Moscow Institute of Physics and Technology, 47 Dolgoprudni,
More informationLecture #8-3 Oscillations, Simple Harmonic Motion
Lecture #8-3 Oscillations Siple Haronic Motion So far we have considered two basic types of otion: translation and rotation. But these are not the only two types of otion we can observe in every day life.
More informationNUMERICAL MODELLING OF THE TYRE/ROAD CONTACT
NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT PACS REFERENCE: 43.5.LJ Krister Larsson Departent of Applied Acoustics Chalers University of Technology SE-412 96 Sweden Tel: +46 ()31 772 22 Fax: +46 ()31
More informationPh 20.3 Numerical Solution of Ordinary Differential Equations
Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing
More informationExtension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels
Extension of CSRSM for the Paraetric Study of the Face Stability of Pressurized Tunnels Guilhe Mollon 1, Daniel Dias 2, and Abdul-Haid Soubra 3, M.ASCE 1 LGCIE, INSA Lyon, Université de Lyon, Doaine scientifique
More informationGeneral Properties of Radiation Detectors Supplements
Phys. 649: Nuclear Techniques Physics Departent Yarouk University Chapter 4: General Properties of Radiation Detectors Suppleents Dr. Nidal M. Ershaidat Overview Phys. 649: Nuclear Techniques Physics Departent
More informationA simple phenomenologic model for particle transport in spaceperiodic potentials in underdamped systems
A siple phenoenologic odel for particle transport in spaceperiodic potentials in underdaped systes IG MARCHENKO 1,(a,b), II MARCHENKO 3, A ZHIGLO 1 1 NSC Kharov Institute of Physics and Technology, Aadeichesaya
More informationSOLUTIONS. PROBLEM 1. The Hamiltonian of the particle in the gravitational field can be written as, x 0, + U(x), U(x) =
SOLUTIONS PROBLEM 1. The Hailtonian of the particle in the gravitational field can be written as { Ĥ = ˆp2, x 0, + U(x), U(x) = (1) 2 gx, x > 0. The siplest estiate coes fro the uncertainty relation. If
More informationCommun Nonlinear Sci Numer Simulat
Coun Nonlinear Sci Nuer Siulat 17 (2012) 1980 1985 Contents lists available at SciVerse ScienceDirect Coun Nonlinear Sci Nuer Siulat journal hoepage: www.elsevier.co/locate/cnsns Using delayed daping to
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationCh 12: Variations on Backpropagation
Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith
More informationSupporting information for Self-assembly of multicomponent structures in and out of equilibrium
Supporting inforation for Self-assebly of ulticoponent structures in and out of equilibriu Stephen Whitela 1, Rebecca Schulan 2, Lester Hedges 1 1 Molecular Foundry, Lawrence Berkeley National Laboratory,
More informationChapter 1: Basics of Vibrations for Simple Mechanical Systems
Chapter 1: Basics of Vibrations for Siple Mechanical Systes Introduction: The fundaentals of Sound and Vibrations are part of the broader field of echanics, with strong connections to classical echanics,
More informationSPECTRUM sensing is a core concept of cognitive radio
World Acadey of Science, Engineering and Technology International Journal of Electronics and Counication Engineering Vol:6, o:2, 202 Efficient Detection Using Sequential Probability Ratio Test in Mobile
More informationOptical Properties of Plasmas of High-Z Elements
Forschungszentru Karlsruhe Techni und Uwelt Wissenschaftlishe Berichte FZK Optical Properties of Plasas of High-Z Eleents V.Tolach 1, G.Miloshevsy 1, H.Würz Project Kernfusion 1 Heat and Mass Transfer
More informationAnalysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning
Analysis of Ipulsive Natural Phenoena through Finite Difference Methods A MATLAB Coputational Project-Based Learning Nicholas Kuia, Christopher Chariah, Mechatronics Engineering, Vaughn College of Aeronautics
More informationDETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION
DETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION Masaki WAKUI 1 and Jun IYAMA and Tsuyoshi KOYAMA 3 ABSTRACT This paper shows a criteria to detect
More informationSEISMIC FRAGILITY ANALYSIS
9 th ASCE Specialty Conference on Probabilistic Mechanics and Structural Reliability PMC24 SEISMIC FRAGILITY ANALYSIS C. Kafali, Student M. ASCE Cornell University, Ithaca, NY 483 ck22@cornell.edu M. Grigoriu,
More informationModulation of Harmonic Emission Spectra from Intense Laser-Plasma Interactions
Modulation of Haronic Eission Spectra fro Intense Laser-Plasa Interactions T.J.M. Boyd and R. Ondarza-Rovira 2 Centre for Physics, University of Essex, Wivenhoe Park, Colchester CO4 3SQ, U.K. 2 ININ, A.P.
More informationForce and dynamics with a spring, analytic approach
Force and dynaics with a spring, analytic approach It ay strie you as strange that the first force we will discuss will be that of a spring. It is not one of the four Universal forces and we don t use
More informationChapter 2: Introduction to Damping in Free and Forced Vibrations
Chapter 2: Introduction to Daping in Free and Forced Vibrations This chapter ainly deals with the effect of daping in two conditions like free and forced excitation of echanical systes. Daping plays an
More informationFigure 1: Equivalent electric (RC) circuit of a neurons membrane
Exercise: Leaky integrate and fire odel of neural spike generation This exercise investigates a siplified odel of how neurons spike in response to current inputs, one of the ost fundaental properties of
More informationSupporting Information
Supporting Inforation Nash et al. 10.1073/pnas.1507413112 Equation of Motion If a gyroscope is spinning with a large constant angular frequency, ω, around its principal axis, ^l, then its dynaics are deterined
More informationScattering and bound states
Chapter Scattering and bound states In this chapter we give a review of quantu-echanical scattering theory. We focus on the relation between the scattering aplitude of a potential and its bound states
More informationIN modern society that various systems have become more
Developent of Reliability Function in -Coponent Standby Redundant Syste with Priority Based on Maxiu Entropy Principle Ryosuke Hirata, Ikuo Arizono, Ryosuke Toohiro, Satoshi Oigawa, and Yasuhiko Takeoto
More informationESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics
ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS A Thesis Presented to The Faculty of the Departent of Matheatics San Jose State University In Partial Fulfillent of the Requireents
More informationDynamic analysis of frames with viscoelastic dampers: a comparison of damper models
Structural Engineering and Mechanics, Vol. 41, No. 1 (2012) 113-137 113 Dynaic analysis of fraes with viscoelastic dapers: a coparison of daper odels R. Lewandowski*, A. Bartkowiak a and H. Maciejewski
More informationNow multiply the left-hand-side by ω and the right-hand side by dδ/dt (recall ω= dδ/dt) to get:
Equal Area Criterion.0 Developent of equal area criterion As in previous notes, all powers are in per-unit. I want to show you the equal area criterion a little differently than the book does it. Let s
More informationDensity and structure of undercooled liquid titanium
Article Condensed Matter Physics March 2012 Vol.57 No.7: 719 723 doi: 10.1007/s11434-011-4945-6 Density and structure of undercooled liquid titaniu WANG HaiPeng, YANG ShangJing & WEI BingBo * Departent
More informationAll Excuses must be taken to 233 Loomis before 4:15, Monday, April 30.
Miscellaneous Notes he end is near don t get behind. All Excuses ust be taken to 233 Loois before 4:15, Monday, April 30. he PHYS 213 final exa ties are * 8-10 AM, Monday, May 7 * 8-10 AM, uesday, May
More informationKernel Methods and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley ENSIAG 2 / osig 1 Second Seester 2012/2013 Lesson 20 2 ay 2013 Kernel ethods and Support Vector achines Contents Kernel Functions...2 Quadratic
More informationarxiv: v2 [physics.soc-ph] 4 Oct 2018
Analytical and nuerical study of the non-linear noisy voter odel on coplex networks A. F. Peralta, 1, a) A. Carro,, 3 M. San Miguel, 1 and R. Toral 1 1) IFISC (Instituto de Física Interdisciplinar y Sisteas
More informationSolving initial value problems by residual power series method
Theoretical Matheatics & Applications, vol.3, no.1, 13, 199-1 ISSN: 179-9687 (print), 179-979 (online) Scienpress Ltd, 13 Solving initial value probles by residual power series ethod Mohaed H. Al-Sadi
More informationIII.H Zeroth Order Hydrodynamics
III.H Zeroth Order Hydrodynaics As a first approxiation, we shall assue that in local equilibriu, the density f 1 at each point in space can be represented as in eq.iii.56, i.e. f 0 1 p, q, t = n q, t
More informationlecture 37: Linear Multistep Methods: Absolute Stability, Part I lecture 38: Linear Multistep Methods: Absolute Stability, Part II
lecture 37: Linear Multistep Methods: Absolute Stability, Part I lecture 3: Linear Multistep Methods: Absolute Stability, Part II 5.7 Linear ultistep ethods: absolute stability At this point, it ay well
More informationHybrid System Identification: An SDP Approach
49th IEEE Conference on Decision and Control Deceber 15-17, 2010 Hilton Atlanta Hotel, Atlanta, GA, USA Hybrid Syste Identification: An SDP Approach C Feng, C M Lagoa, N Ozay and M Sznaier Abstract The
More informationXiaoming Mao. Department of Physics and Astronomy, University of Pennsylvania. Collaborators: Tom Lubensky, Ning Xu, Anton Souslov, Andrea Liu
Xiaoing Mao Departent of Physics and Astronoy, University of Pennsylvania Collaborators: To Lubensky, Ning Xu, Anton Souslov, Andrea Liu Feb., 009 What is isostaticity? Isostatic systes are at the onset
More informationSynchronization in large directed networks of coupled phase oscillators
CHAOS 16, 015107 2005 Synchronization in large directed networks of coupled phase oscillators Juan G. Restrepo a Institute for Research in Electronics and Applied Physics, University of Maryland, College
More informationKinetic Theory of Gases: Elementary Ideas
Kinetic Theory of Gases: Eleentary Ideas 9th February 011 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion of
More informationSymbolic Analysis as Universal Tool for Deriving Properties of Non-linear Algorithms Case study of EM Algorithm
Acta Polytechnica Hungarica Vol., No., 04 Sybolic Analysis as Universal Tool for Deriving Properties of Non-linear Algoriths Case study of EM Algorith Vladiir Mladenović, Miroslav Lutovac, Dana Porrat
More informationHysteresis model for magnetic materials using the Jiles-Atherton model
Hysteresis odel for agnetic aterials using the Jiles-Atherton odel Predrag Petrovic Technical faculty Svetog Save 65 32 Cacak, pegi@ei.yu Nebojsa itrovic Technical faculty Svetog Save 65 32 Cacak, itar@tfc.tfc.kg.ac.yu
More informationOn the characterization of non-linear diffusion equations. An application in soil mechanics
On the characterization of non-linear diffusion equations. An application in soil echanics GARCÍA-ROS, G., ALHAMA, I., CÁNOVAS, M *. Civil Engineering Departent Universidad Politécnica de Cartagena Paseo
More informationKinetic Theory of Gases: Elementary Ideas
Kinetic Theory of Gases: Eleentary Ideas 17th February 2010 1 Kinetic Theory: A Discussion Based on a Siplified iew of the Motion of Gases 1.1 Pressure: Consul Engel and Reid Ch. 33.1) for a discussion
More informationDynamical thermalization of disordered nonlinear lattices
PHYSICAL REVIEW E 8, 56 9 Dynaical theralization of disordered nonlinear lattices Mario Mulansky, Karsten Ahnert, Arkady Pikovsky,, and Dia L. Shepelyansky, Departent of Physics and Astronoy, Potsda University,
More informationNon-Parametric Non-Line-of-Sight Identification 1
Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,
More informationDispersion. February 12, 2014
Dispersion February 1, 014 In aterials, the dielectric constant and pereability are actually frequency dependent. This does not affect our results for single frequency odes, but when we have a superposition
More informationEfficient dynamic events discrimination technique for fiber distributed Brillouin sensors
Efficient dynaic events discriination technique for fiber distributed rillouin sensors Carlos A. Galindez,* Francisco J. Madruga, and Jose M. Lopez-Higuera Photonics Engineering Group, Universidad de Cantabria,Edif.
More informationMSEC MODELING OF DEGRADATION PROCESSES TO OBTAIN AN OPTIMAL SOLUTION FOR MAINTENANCE AND PERFORMANCE
Proceeding of the ASME 9 International Manufacturing Science and Engineering Conference MSEC9 October 4-7, 9, West Lafayette, Indiana, USA MSEC9-8466 MODELING OF DEGRADATION PROCESSES TO OBTAIN AN OPTIMAL
More informationCHAPTER 19: Single-Loop IMC Control
When I coplete this chapter, I want to be able to do the following. Recognize that other feedback algoriths are possible Understand the IMC structure and how it provides the essential control features
More informationNonmonotonic Networks. a. IRST, I Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I Povo (Trento) Italy
Storage Capacity and Dynaics of Nononotonic Networks Bruno Crespi a and Ignazio Lazzizzera b a. IRST, I-38050 Povo (Trento) Italy, b. Univ. of Trento, Physics Dept., I-38050 Povo (Trento) Italy INFN Gruppo
More informationElectrical Engineering 3BB3: Cellular Bioelectricity (2008) Solutions to Midterm Quiz #1
Electrical Engineering 3BB3: Cellular Bioelectricity (2008) Solutions to Midter Quiz #1 1. In typical excitable cells there will be a net influx of K + through potassiu ion channels if: a. V Vrest >, b.
More informationDepartment of Physics Preliminary Exam January 3 6, 2006
Departent of Physics Preliinary Exa January 3 6, 2006 Day 1: Classical Mechanics Tuesday, January 3, 2006 9:00 a.. 12:00 p.. Instructions: 1. Write the answer to each question on a separate sheet of paper.
More informationIntelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes
More informationAn Improved Particle Filter with Applications in Ballistic Target Tracking
Sensors & ransducers Vol. 72 Issue 6 June 204 pp. 96-20 Sensors & ransducers 204 by IFSA Publishing S. L. http://www.sensorsportal.co An Iproved Particle Filter with Applications in Ballistic arget racing
More informationThe accelerated expansion of the universe is explained by quantum field theory.
The accelerated expansion of the universe is explained by quantu field theory. Abstract. Forulas describing interactions, in fact, use the liiting speed of inforation transfer, and not the speed of light.
More informationElectrical Engineering 3BB3: Cellular Bioelectricity (2013) Solutions to Midterm Quiz #1
Electrical Engineering 3BB3: Cellular Bioelectricity (2013) Solutions to Midter Quiz #1 1. In a typical excitable cell, the ion species with the ost positive Nernst equilibriu potential is: a. potassiu,
More informationSF Chemical Kinetics.
SF Cheical Kinetics. Lecture 5. Microscopic theory of cheical reaction inetics. Microscopic theories of cheical reaction inetics. basic ai is to calculate the rate constant for a cheical reaction fro first
More informationChapter 1 Introduction and Kinetics of Particles
Chapter 1 Introduction and Kinetics of Particles 1.1 Introduction There are two ain approaches in siulating the transport equations (heat, ass, and oentu), continuu and discrete. In continuu approach,
More informationA remark on a success rate model for DPA and CPA
A reark on a success rate odel for DPA and CPA A. Wieers, BSI Version 0.5 andreas.wieers@bsi.bund.de Septeber 5, 2018 Abstract The success rate is the ost coon evaluation etric for easuring the perforance
More information2.141 Modeling and Simulation of Dynamic Systems Assignment #2
2.141 Modeling and Siulation of Dynaic Systes Assignent #2 Out: Wednesday Septeber 20, 2006 Due: Wednesday October 4, 2006 Proble 1 The sketch shows a highly siplified diagra of a dry-dock used in ship
More informationGeneralized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials
Copyright 7 Tech Science Press CMC vol.53 no.4 pp.37-34 7 Generalized Rayleigh Wave Dispersion in a Covered Half-space Made of Viscoelastic Materials S.D. Akbarov and M. Negin 3 Abstract: Dispersion of
More informationA DESIGN GUIDE OF DOUBLE-LAYER CELLULAR CLADDINGS FOR BLAST ALLEVIATION
International Journal of Aerospace and Lightweight Structures Vol. 3, No. 1 (2013) 109 133 c Research Publishing Services DOI: 10.3850/S201042862013000550 A DESIGN GUIDE OF DOUBLE-LAYER CELLULAR CLADDINGS
More informationExperimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis
City University of New York (CUNY) CUNY Acadeic Works International Conference on Hydroinforatics 8-1-2014 Experiental Design For Model Discriination And Precise Paraeter Estiation In WDS Analysis Giovanna
More informationNonuniqueness of canonical ensemble theory. arising from microcanonical basis
onuniueness of canonical enseble theory arising fro icrocanonical basis arxiv:uant-ph/99097 v2 25 Oct 2000 Suiyoshi Abe and A. K. Rajagopal 2 College of Science and Technology, ihon University, Funabashi,
More informationSlanted coupling of one-dimensional arrays of small tunnel junctions
JOURNAL OF APPLIED PHYSICS VOLUME 84, NUMBER 1 15 DECEMBER 1998 Slanted coupling of one-diensional arrays of sall tunnel junctions G. Y. Hu Departent of Physics and Astronoy, Louisiana State University,
More informationRole of rf electric and magnetic fields in heating of micro-protrusions in accelerating structures. Gregory S. Nusinovich and Thomas M. Antonsen, Jr.
Role of rf electric and agnetic fields in heating of icro-protrusions in accelerating structures. Gregory S. Nusinovich and Thoas M. Antonsen, Jr. Abstract It is known that high-gradient operation in etallic
More informationCOULD A VARIABLE MASS OSCILLATOR EXHIBIT THE LATERAL INSTABILITY?
Kragujevac J. Sci. 3 (8) 3-44. UDC 53.35 3 COULD A VARIABLE MASS OSCILLATOR EXHIBIT THE LATERAL INSTABILITY? Nebojša Danilović, Milan Kovačević and Vukota Babović Institute of Physics, Faculty of Science,
More informationPlasma Modulation of Harmonic Emission Spectra from Laser-Dense Plasma Interactions
Plasa Modulation of Haronic Eission Spectra fro Laser-Dense Plasa Interactions R. Ondarza-Rovira and T.J.M. Boyd 2 ININ, A.P. 8-027, México 80, Distrito Federal, Mexico; ondarza@nuclear.inin.x 2 Centre
More informationI. Understand get a conceptual grasp of the problem
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Departent o Physics Physics 81T Fall Ter 4 Class Proble 1: Solution Proble 1 A car is driving at a constant but unknown velocity,, on a straightaway A otorcycle is
More informationThe Transactional Nature of Quantum Information
The Transactional Nature of Quantu Inforation Subhash Kak Departent of Coputer Science Oklahoa State University Stillwater, OK 7478 ABSTRACT Inforation, in its counications sense, is a transactional property.
More information13 Harmonic oscillator revisited: Dirac s approach and introduction to Second Quantization
3 Haronic oscillator revisited: Dirac s approach and introduction to Second Quantization. Dirac cae up with a ore elegant way to solve the haronic oscillator proble. We will now study this approach. The
More informationDesign and Experimental Research of Atomizer Based on Micro Abrasive Ultrasonic Polishing Bang-fu WANG, Yin ZHEN, Juan SONG and A-chun ZHU
217 3rd International Conference on Applied Mechanics and Mechanical Autoation (AMMA 217) ISBN: 978-1-6595-479- Design and Experiental Research of Atoizer Based on Micro Abrasive Ultrasonic Polishing Bang-fu
More informationA Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair
Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving
More informationUsing a De-Convolution Window for Operating Modal Analysis
Using a De-Convolution Window for Operating Modal Analysis Brian Schwarz Vibrant Technology, Inc. Scotts Valley, CA Mark Richardson Vibrant Technology, Inc. Scotts Valley, CA Abstract Operating Modal Analysis
More informationREDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION
ISSN 139 14X INFORMATION TECHNOLOGY AND CONTROL, 008, Vol.37, No.3 REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION Riantas Barauskas, Vidantas Riavičius Departent of Syste Analysis, Kaunas
More informationSupport Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization
Recent Researches in Coputer Science Support Vector Machine Classification of Uncertain and Ibalanced data using Robust Optiization RAGHAV PAT, THEODORE B. TRAFALIS, KASH BARKER School of Industrial Engineering
More informationOptimizing energy potentials for success in protein tertiary structure prediction Ting-Lan Chiu 1 and Richard A Goldstein 1,2
Research Paper 223 Optiizing energy potentials for success in protein tertiary structure prediction Ting-Lan Chiu 1 and Richard A Goldstein 1,2 Background: Success in solving the protein structure prediction
More informationIn this chapter, we consider several graph-theoretic and probabilistic models
THREE ONE GRAPH-THEORETIC AND STATISTICAL MODELS 3.1 INTRODUCTION In this chapter, we consider several graph-theoretic and probabilistic odels for a social network, which we do under different assuptions
More informationEasy Evaluation Method of Self-Compactability of Self-Compacting Concrete
Easy Evaluation Method of Self-Copactability of Self-Copacting Concrete Masanori Maruoka 1 Hiroi Fujiwara 2 Erika Ogura 3 Nobu Watanabe 4 T 11 ABSTRACT The use of self-copacting concrete (SCC) in construction
More informationMeasuring orbital angular momentum superpositions of light by mode transformation
CHAPTER 7 Measuring orbital angular oentu superpositions of light by ode transforation In chapter 6 we reported on a ethod for easuring orbital angular oentu (OAM) states of light based on the transforation
More informationNumerically repeated support splitting and merging phenomena in a porous media equation with strong absorption. Kenji Tomoeda
Journal of Math-for-Industry, Vol. 3 (C-), pp. Nuerically repeated support splitting and erging phenoena in a porous edia equation with strong absorption To the eory of y friend Professor Nakaki. Kenji
More informationRegular article Maximum radius of convergence perturbation theory: test calculations on Be, Ne, H 2 and HF
Theor Che Acc (2003) 110: 185 189 DOI 10.1007/s00214-003-0473-z Regular article Maxiu radius of convergence perturbation theory: test calculations on Be, Ne, H 2 and HF Kotaro Yokoyaa 1, Haruyuki Nakano
More informationThe Solution of One-Phase Inverse Stefan Problem. by Homotopy Analysis Method
Applied Matheatical Sciences, Vol. 8, 214, no. 53, 2635-2644 HIKARI Ltd, www.-hikari.co http://dx.doi.org/1.12988/as.214.43152 The Solution of One-Phase Inverse Stefan Proble by Hootopy Analysis Method
More informationPattern Recognition and Machine Learning. Artificial Neural networks
Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2016 Lessons 7 14 Dec 2016 Outline Artificial Neural networks Notation...2 1. Introduction...3... 3 The Artificial
More informationData-Driven Imaging in Anisotropic Media
18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands
More informationThe source of THz radiation based on dielectric waveguide excited by sequence of electron bunches
Journal of Physics: Conference Series PAPER OPEN ACCESS The source of THz radiation based on dielectric waveguide excited by sequence of electron bunches To cite this article: A M Altark and A D Kanareykin
More informationChapter 4: Hypothesis of Diffusion-Limited Growth
Suary This section derives a useful equation to predict quantu dot size evolution under typical organoetallic synthesis conditions that are used to achieve narrow size distributions. Assuing diffusion-controlled
More informationOBJECTIVES INTRODUCTION
M7 Chapter 3 Section 1 OBJECTIVES Suarize data using easures of central tendency, such as the ean, edian, ode, and idrange. Describe data using the easures of variation, such as the range, variance, and
More information13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices
CS71 Randoness & Coputation Spring 018 Instructor: Alistair Sinclair Lecture 13: February 7 Disclaier: These notes have not been subjected to the usual scrutiny accorded to foral publications. They ay
More informationChapter 12. Quantum gases Microcanonical ensemble
Chapter 2 Quantu gases In classical statistical echanics, we evaluated therodynaic relations often for an ideal gas, which approxiates a real gas in the highly diluted liit. An iportant difference between
More informationHee = ~ dxdy\jj+ (x) 'IJ+ (y) u (x- y) \jj (y) \jj (x), V, = ~ dx 'IJ+ (x) \jj (x) V (x), Hii = Z 2 ~ dx dy cp+ (x) cp+ (y) u (x- y) cp (y) cp (x),
SOVIET PHYSICS JETP VOLUME 14, NUMBER 4 APRIL, 1962 SHIFT OF ATOMIC ENERGY LEVELS IN A PLASMA L. E. PARGAMANIK Khar'kov State University Subitted to JETP editor February 16, 1961; resubitted June 19, 1961
More informationOn Constant Power Water-filling
On Constant Power Water-filling Wei Yu and John M. Cioffi Electrical Engineering Departent Stanford University, Stanford, CA94305, U.S.A. eails: {weiyu,cioffi}@stanford.edu Abstract This paper derives
More informationOptimal nonlinear Bayesian experimental design: an application to amplitude versus offset experiments
Geophys. J. Int. (23) 155, 411 421 Optial nonlinear Bayesian experiental design: an application to aplitude versus offset experients Jojanneke van den Berg, 1, Andrew Curtis 2,3 and Jeannot Trapert 1 1
More informationŞtefan ŞTEFĂNESCU * is the minimum global value for the function h (x)
7Applying Nelder Mead s Optiization Algorith APPLYING NELDER MEAD S OPTIMIZATION ALGORITHM FOR MULTIPLE GLOBAL MINIMA Abstract Ştefan ŞTEFĂNESCU * The iterative deterinistic optiization ethod could not
More information