Motions of an Molecular and A-Helical Polypeptide: Harmonic Dynamics

Total Page:16

File Type:pdf, Size:1020Kb

Motions of an Molecular and A-Helical Polypeptide: Harmonic Dynamics Motions of an a-Helical Polypeptide: Comparison of Molecular and Harmonic Dynamics D. PERAHIA,1 R. M. LEVY2 and M. KARPLUS3 ~Laboratoire d’Enzymologie Physico-Chimique et Mol~culaire, Bat 430, Universit~ de Paris-Sud, 91405 Orsay, France; 2Department of Chemistry, Rutgers, The State University of New Jersey, New Brunswick. New Jersey 08903. USA 3Department of Chemistry, Harvard University, Cambridge, Massachusetts 02138, USA, and Laboratoire de Chimie des Interactions Mol~culaires, Coll~ge de France, 75005 Paris, France SYNOPSIS Molecular and harmonic dynamics simulations have been performed for a decaglycine a-helix. The extent of anharmonicity for various observables is studied by a direct comparison of the two types of simulations, at temperatures ranging from 5 to 300 K. The fluctuations of the cartesian, internal and normal mode coordinates, and their time dependence, are analyzed. The heat capacity of the a-helix is evaluated both from the temperature response of the system to an energy perturbation and from the fluctuations in the temperature of the system. It is shown that the anharmonicity depends on the kind of observable. The root mean square atomic fluctuations have significant anharmonic components at temperatures above 100 K. In contrast, the dihedral angle fluctuations are much closer to being harmonic at all the temperatures considered. The analysis of potentials of mean force experienced by individual atoms shows that atomic displacements have approximately Gaussian distributions from 50 to 300 K, with different force constants at each temperature (quasi-harmonic model). At 300 K, the force constants obtained by molecular dynamics are significantly lower than in the harmonic ease. The time dependence of the projection of the molecular dynamics displacements on the normal mode coordinates shows that mode mixing is important above 100 K. The motions of the helix associated with the low-frequency normal modes are described and illustrated. INTRODUCTION diffraction and Mossbauer studies of protein crys- tals.5-v Further, molecular dynamics analysis of A wide range of experimental and theoretical stud- the time scales of the atomic fluctuations have ies have shown the existence of internal motions in indicated that they have a local subpicosecond proteins)’2 Of particular importance have been component superposed on a larger scale more col- molecular dynamics simulations, which have lective motion that develops over 1-10 or more s,9 demonstrated that the root mean square (rms) of picoseconds, Two extreme models for such con- ° atom displacements in proteins are on the order of densed-phase dynamic behavior can be formulated) 0.5 ~,, with a range of 0.1-2/~,~.3.4 a magnitude in In the first, analogous to a liquid, the atomic overall agreement with estimates from x-ray fluctuations can be described as localized diffusion in which the collisions between atoms play the essential role; in the second, analbgous to a solid, the atomic fluctuations involve the coupling of a 1990 John Wiley & Sons, Inc. CCC 0(D6-3525/90/040645-33 $04.00 large number of normal modes. It is likely that the Biopolymers, Vol. 29, 645-677 (1990) interiors of proteins contain both fluid-like and 646 PE~AHIA, LEVY, AND KARPLUS solid-like regionsS’~°’11; e.g., it has been suggested the classical equations of motions for the atoms of that hydrophobic clusters are more fluid-like, which the helix is composed. Up to and including whereas parts with considerable secondary struc- 100 K (5, 50, and 100 K), a kinetic energy twice the ture (a-helices, fl-sheets) are more solid-like.1° Re- final desired temperature was imparted to the sys- sults on the collective character of the motions tem in a single step; in all cases the temperature indicate that the harmonic model may be useful for dropped to its mean value within 50-100 steps. At analysis, even if anharmonic contributions are pre- temperatures above 100 K, the system was heated sent. Further, it has been suggested that the har- 25-50 K every 200 steps until the desired tempera- monic model may be improved by the quasi- ture was reached. The total kinetic energy im- harmonic approximation, in which the force con- parted to the system for the 20 trajectories with a stants are obtained not by differentiating the static given temperature was the same. The values of the potential energy function, but by averaging over average temperature of individual trajectories var- the potential of mean force at a given ied by up to 20%. In all cases, the initial coordi- temperature.12-’4 nates for the system corresponded to the minimum In this paper we present the detailed results of a energy geometry of the decaglycine helix and the study of the full molecular dynamics of an isolated individual trajectories of a set differed only in the a-helix performed at a series of temperatures (5-300 initial velocity components chosen by using a set of K). The results are compared with those for the random numbers corresponding to the appropriate same system in the harmonic approximation.~5-18 Maxwellian distribution. Once the final tempera- The magnitudes of the atomic fluctuations and ture was reached in each case, the analysis of the their time dependence are analyzed. Also, the po- trajectory was begun. tentials of mean force seen by the individual atoms Twenty short trajectories with different initial at the various temperatures are calculated. The conditions were used in an attempt to cover a magnitudes and the time correlation functions of larger portion of phase space than would a single the fluctuations in cartesian and internal coordi- extended trajectory in this relatively harmonic sys- nates are compared with the corresponding har- tem. A set of 20 trajectories was chosen after a monic dynamics results. The contributions of de- series of test calculations performed at 5 K in phasing and population relaxation to the form of which the system was expected to be harmonic. It the correlation functions are examined. The heat was found that this number gave consistent results capacity of the a-helix is evaluated both from the when repeated with a different distribution of ini- temperature response of the system to an energy tial velocities. Sets of 5-10 trajectories showed perturbation and from the fluctuations in the tem- deviation between different samples that were too perature of the system. The motions of the helix large. We found that the differences between the associated with the low-frequency normal modes rms displacements of the C~ atoms of the helix are described. between sets of 5 trajectories at 5 K were about 0.005 /~, which represents almost 10% of the rms value at this temperature. Similar deviations were found with sets of 10 trajectories. The difference METHODOLOGY decreased to 0.001 ~, when two sets of 20 trajecto- ries were compared. In addition, a single longer The system considered is the model a-helix trajectory was calculated at 300 K; it was equili- decaglycine, with the amino and carboxyl ends brated for 9 ps and the analysis period lasted 10 ps. blocked by the CH3CO and NHCH3 groups, re- The numerical integration of the equations of mo- spectively; the decaglycine helix contains eight hy- tions was carried out with the Verlet algorithm~9 drogen bonds. All atoms and degrees of freedom of with a time step of 5.10’~6 s. The program the molecule are treated explicitly, except that the CHARMM2° was employed for all of the computa- CH2 and CH3 groups are represented as extended tions° atoms. The details of the molecular model and the The mean square amplitude of the fluctuations potential functions used in the present simulation of a given atom, (hr2)D, along the direction D, is have been presented previously.~5 A series of tem- evaluated by peratures between 5 and 300 K defined by the mean kinetic energy of the system were studied. ~_, ((r~D- (r)D)2) For each temperature, 20 trajectories, each 1 ps in (hr2)D= ~-1 (1) length, were calculated by solving simultaneously n MOTIONS OF AN a-HELICAL POI,YPEPTIDE 647 in which riD is the projection of the instantaneous was evaluated. The components of this function coordinate of the atom in the Ith trajectory onto along given directions were also considered. In the the D direction and n is the number of trajectories harmonic case, the time correlation function of in the set; the angle brackets represent time aver- positional fluctuations was evaluated by the equa- ages over individual trajectories and (r)D = tion (ET.~(r~D))/n is the mean value of the displace- ment in the D direction for the 20 trajectories. 3 N - 6 The harmonic approximation to the amplitudes k~T ~ aj~ cos (wit) of atomic vibrations are obtained from the eigen- j=~ values and eigenvectors of the decaglycine a-helix15; <~r(O). ~r(t))h~~ = that is (2) Spectral densities corresponding to the variation of a given va~able in the molecular dynamics simu- lation were obtained by taking the Fourier trans- in which k is the Boltzmann constant, T is the e form of its autocorrelation function. The fast absolute temperature, N is the number of atoms, Fourier transform technique was used for these wj is the frequency of the jth normal mode, and calculationsfl2 In the harmonic case, the spectral aid is the component of the jth eigenvector along densities were obtained by considering the contri- D direction for the atom of interest. To obtain the bution of normal modes to the mean square fluc- normal modes, the energy of the a-helix was mini- tuation of the given variable. mized by using the full potential, and the second- derivative matrix of the potential energy with re- spect to the mass-weighted cartesian coordinate displacements relative to the minimal-energy RESULTS AND DISCUSSION structure was determined.
Recommended publications
  • The Halphen Family of Distributions for Modeling Maximum
    Transactions on Ecology and the Environment vol 6, © 1994 WIT Press, www.witpress.com, ISSN 1743-3541 The Halphen family of distributions for modeling maximum annual flood series: New developments L. Perreault, V. Fortin & B. Bobee Chair on Statistical Hydrology, 2800 rue Einstein, Sainte-Foy, Quebec, Canada, G1V 4C7 Abstract Design of hydraulic structures, optimal operation of reservoirs, and flood forecasting involve estimation of flood flows x^,, with given specified return period T. This estimation is generally achieved by fitting a probability distribution to observed data. In this paper, we briefly present some recent theoretical results of a study on Halphen distributions which were developed in the early forties but have been little used due to the computational complexity. This family of three distributions was motivated by Halphen's desire to construct probability models which met specific requirements of fitting hydrological variables. Halphen probability models have great potential applicability in statistical flood analysis but are not yet familiar to hydrologists. 1 Introduction Engineering activities such as the design of hydraulic structures, the management of water resources systems, and the prevention of flood damage, all require an estimation of flood characteristics. The objective of flood frequency analysis FFA is to estimate the flood x^ exceeded in average once every T year, which may be used to design, for example, spillways or as prior information in real- time flood forcasting models. At gauged locations, x? may be estimated by fitting a parametric distribution F(x; 0) to observed maximum annual floods, where x refer to an observation and 0 to a vector of parameters.
    [Show full text]
  • Harmonic Measurement and Reduction in Power Systems. Rutisurhata Kurniawan Hartana Louisiana State University and Agricultural & Mechanical College
    Louisiana State University LSU Digital Commons LSU Historical Dissertations and Theses Graduate School 1989 Harmonic Measurement and Reduction in Power Systems. Rutisurhata Kurniawan Hartana Louisiana State University and Agricultural & Mechanical College Follow this and additional works at: https://digitalcommons.lsu.edu/gradschool_disstheses Recommended Citation Hartana, Rutisurhata Kurniawan, "Harmonic Measurement and Reduction in Power Systems." (1989). LSU Historical Dissertations and Theses. 4849. https://digitalcommons.lsu.edu/gradschool_disstheses/4849 This Dissertation is brought to you for free and open access by the Graduate School at LSU Digital Commons. It has been accepted for inclusion in LSU Historical Dissertations and Theses by an authorized administrator of LSU Digital Commons. For more information, please contact [email protected]. INFORMATION TO USERS The most advanced technology has been used to photograph and reproduce this manuscript from the microfilm master. UMI films the text directly from the original or copy submitted. Thus, some thesis and dissertation copies are in typewriter face, while others may be from any type of computer printer. The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleedthrough, substandard margins, and improper alignment can adversely afreet reproduction. In the unlikely event that the author did not send UMI a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion. Oversize materials (e.g., maps, drawings, charts) are reproduced by sectioning the original, beginning at the upper left-hand corner and continuing from left to right in equal sections with small overlaps.
    [Show full text]
  • Field Guide to Continuous Probability Distributions
    Field Guide to Continuous Probability Distributions Gavin E. Crooks v 1.0.0 2019 G. E. Crooks – Field Guide to Probability Distributions v 1.0.0 Copyright © 2010-2019 Gavin E. Crooks ISBN: 978-1-7339381-0-5 http://threeplusone.com/fieldguide Berkeley Institute for Theoretical Sciences (BITS) typeset on 2019-04-10 with XeTeX version 0.99999 fonts: Trump Mediaeval (text), Euler (math) 271828182845904 2 G. E. Crooks – Field Guide to Probability Distributions Preface: The search for GUD A common problem is that of describing the probability distribution of a single, continuous variable. A few distributions, such as the normal and exponential, were discovered in the 1800’s or earlier. But about a century ago the great statistician, Karl Pearson, realized that the known probabil- ity distributions were not sufficient to handle all of the phenomena then under investigation, and set out to create new distributions with useful properties. During the 20th century this process continued with abandon and a vast menagerie of distinct mathematical forms were discovered and invented, investigated, analyzed, rediscovered and renamed, all for the purpose of de- scribing the probability of some interesting variable. There are hundreds of named distributions and synonyms in current usage. The apparent diver- sity is unending and disorienting. Fortunately, the situation is less confused than it might at first appear. Most common, continuous, univariate, unimodal distributions can be orga- nized into a small number of distinct families, which are all special cases of a single Grand Unified Distribution. This compendium details these hun- dred or so simple distributions, their properties and their interrelations.
    [Show full text]
  • Research of Harmonics in Power System Signal Using Gaussian's
    International Journal of Recent Technology and Engineering (IJRTE) ISSN: 2277-3878, Volume-8 Issue-2S8, August 2019 Research of Harmonics in Power System Signal using Gaussian’s Distribution Overlapping by Receiver Operating Characteristics (Roc) Curve S.Sangeetha, P.Venkatakrishnan, R.Shirisha the signal and image processing, machine learning area in Abstract— Harmonic analysis of the power system signal is a field [1]-[6]. Further, now emerging trends in the ROC proliferating research in the field of electronics technology. analysis in healthcare related biomedical field [7],, Whenever, we analyze odd and ever harmonics are present in the particularly in the radiology frequency analysis [8], where it signal, imperative operation is needed to transform from the time is used to quantify amount of the accuracy of prediction tests. domain to the frequency domain. Hence, all the researchers are utilizing the Fourier Transform technique is very effective for the Mostly, it is a 2D plot of visual interpretation between the analysis of odd and even harmonics in the frequency domain. In performance and classification of the prediction model [9] – the past two decades, Wavelet Transform is a wonderful technique [12]. This may lead the trade-off between the accuracy of to analyze the harmonics both frequency and time domain as well. prediction and the model order. On linear discrimination The analysis of harmonic and its probability distribution are most theory, this type of curves is devoted in communication important for the purpose to predict the harmonic effects in the theory to quantify the amount of discrimination between present situation. We treated all the harmonics and its corresponding frequency distribution are considered as a zero signal and noise [13].
    [Show full text]
  • Supplementary Information for “Unraveling Electronic Absorption
    Supplementary information for \Unraveling Electronic Absorption Spectra using Nuclear Quantum Effects: Photoactive Yellow Protein and Green Fluorescent Protein Chromophores in Water" Tim J. Zuehlsdorff,1, a) Joseph A. Napoli,2 Joel M. Milanese,1 Thomas E. Markland,2, b) and Christine M. Isborn1, c) 1)Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, USA 2)Department of Chemistry, Stanford University, Stanford, California 94305, USA (Dated: 15 June 2018) a)Electronic mail: tzuehlsdorff@ucmerced.edu b)Electronic mail: [email protected] c)Electronic mail: [email protected] 1 CONTENTS I. Double-counting in the E-ZTFC approach 2 A. Updated model systems 3 B. A double-counting free approach: an average vibronic shape function applied to vertical excitation energies obtained from configurations where the chromophore geometry is optimized in a fixed solvent environment 12 C. Application to the GFP chromophore anion 16 D. Concluding remarks 20 II. Simulation cell sizes and QM regions 20 III. Convergence tests: QM region and MM box 23 IV. Convergence tests: AIMD and AI-PIMD sampling 26 V. Computing generalized moments of a distribution 28 VI. Quantifying non-Gaussian features in the ensemble absorption spectra 29 VII. Structural analysis 31 A. Collective motion 31 B. Nuclear quantum effects on specific bond lengths 34 VIII. Analysis of solvent effects 36 IX. Franck-Condon Shape Functions 39 X. Absorption spectra for the Ensemble and E-ZTFC approaches 47 References 47 I. DOUBLE-COUNTING IN THE E-ZTFC APPROACH The E-ZTFC approach includes chromophore nuclear degrees of freedom both from sam- pling the chromophore motion in the ensemble of vertical excitation energies and from the 2 ground state vibrational wave function used to construct the zero temperature Franck- Condon shape function.
    [Show full text]
  • Harmonic Management of Transmission and Distribution Systems Maty Ghezelayagh University of Wollongong
    University of Wollongong Research Online University of Wollongong Thesis Collection University of Wollongong Thesis Collections 1994 Harmonic management of transmission and distribution systems Maty Ghezelayagh University of Wollongong Recommended Citation Ghezelayagh, Maty, Harmonic management of transmission and distribution systems, Doctor of Philosophy thesis, Department of Electrical and Computer Engineering, University of Wollongong, 1994. http://ro.uow.edu.au/theses/1343 Research Online is the open access institutional repository for the University of Wollongong. For further information contact Manager Repository Services: [email protected]. HARMONIC MANAGEMENT OF TRANSMISSION AND DISTRIBUTION SYSTEMS A thesis submitted for the award of the degree DOCTOR OF PHILOSOPHY from UNIVERSITY OF WOLLONGONG BY M. GHEZELAYAGH, B.Sc, M.Sc in Electrical Engineering UK!V5Rr;iV'Gr WOUC,\"3O::G LI- --Y Department of Electrical and Computer Engineering 1994 ABSTRACT HARMONIC MANAGEMENT OF TRANSMISSION AND DISTRIBUTION SYSTEMS This thesis is concerned with some important aspects of harmonic management in a multi-voltage level power system including transmission, sub-transmission and distribution levels. In the work presented both deterministic and statistical approaches for harmonic analysis have been used. Investigations have been made both in Time and Frequency Domains. An already available computer simulation program in Time Domain has been used to derive frequency domain models of distribution feeders with composite loads. A Frequency Domain computer program has been modified to give the statistical data of harmonic voltage distortion levels of a large power system. The thesis examines some of the fundamental assumptions used in harmonic analysis with the aid of a Time Domain simulation program. In harmonic analysis two of the commonly made assumptions are: (i) harmonic sources are represented with ideal ac current sources having equal phase angles, (ii) non-distorting loads are represented with constant impedances.
    [Show full text]
  • A Framework to Analyze the Stochastic Harmonics and Resonance of Wind Energy Grid Interconnection
    Article A Framework to Analyze the Stochastic Harmonics and Resonance of Wind Energy Grid Interconnection Youngho Cho 1, Choongman Lee 1, Kyeon Hur 1,*, Yong Cheol Kang 2, Eduard Muljadi 3, Sang-Ho Park 4, Young-Do Choy 4 and Gi-Gab Yoon 4 1 School of Electrical and Electronic Engineering, Yonsei University, Seoul 03722, Korea; [email protected] (Y.C.); [email protected] (C.L.) 2 Department of Electrical Engineering, Chonbuk National University, Jeonju 54896, Korea; [email protected] 3 National Renewable Energy Laboratory, Golden, CO 80401, USA; [email protected] 4 Korea Electric Power Research Institute, Daejeon 34056, Korea; [email protected] (S.-H.P.); [email protected] (Y.-D.C.); [email protected] (G.-G.Y.) * Correspondence: [email protected]; Tel.: +82-2-2123-5774 Academic Editor: Paras Mandal Received: 29 May 2016; Accepted: 22 August 2016; Published: 31 August 2016 Abstract: This paper addresses a modeling and analysis methodology for investigating the stochastic harmonics and resonance concerns of wind power plants (WPPs). Wideband harmonics from modern wind turbines (WTs) are observed to be stochastic, associated with real power production, and they may adversely interact with the grid impedance and cause unexpected harmonic resonance, if not comprehensively addressed in the planning and commissioning of the WPPs. These issues should become more critical as wind penetration levels increase. We thus propose a planning study framework comprising the following functional steps: First, the best fitted probability density functions (PDFs) of the harmonic components of interest in the frequency domain are determined. In operations planning, maximum likelihood estimations (MLEs) followed by a chi-square test are used once field measurements or manufacturers’ data are available.
    [Show full text]
  • Exact Solution for the Singlet Density Distributions and Second-Order Correlations of Normal-Mode Coordinates for Hard Rods in One Dimension C
    JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 23 15 JUNE 1999 Exact solution for the singlet density distributions and second-order correlations of normal-mode coordinates for hard rods in one dimension C. Daniel Barnes and David A. Kofke Department of Chemical Engineering, State University of New York at Buffalo, Buffalo, New York 14260-4200 ~Received 8 December 1998; accepted 25 March 1999! We examine the distribution of normal-mode coordinates ~defined via the eigenvectors of a chain of harmonic oscillators! for a system of purely repulsive hard rods in one dimension. We obtain an exact solution for the singlet density distribution, and separately for the covariances of the normal-mode coordinates. The hard-rod behavior is examined in terms of its deviation from the corresponding distributions for the system of harmonic oscillators. All off-diagonal covariances are zero in the hard-rod system, and the ~on-diagonal! variances vary with the normal-mode wave number exactly as in the harmonic system. The detailed singlet normal-mode density distributions are very smooth but nonanalytic, and they differ from the ~Gaussian! distributions of the corresponding harmonic system. However, all of the normal-mode coordinate distributions differ in roughly the same way when properly scaled by the distribution variance, and the differences vanish as 1/N in the thermodynamic limit of an infinite number of particles N.©1999 American Institute of Physics. @S0021-9606~99!51523-7# I. INTRODUCTION of this system can be solved analytically by diagonalizing the quadratic form of the internal energy, which effectively con- The system of hard rods in one dimension, also known verts the system into a collection of independent harmonic 1 as the Tonks gas, is of interest as a prototype for studying oscillators.
    [Show full text]
  • Introducing Harmonic Distribution in Wikipedia
    FACULTAT DE CIENCIES` I BIOCIENCIES` DEPARTAMENT DE MATEMATIQUES` Treball de Final de Grau d'Estad´ıstica Aplicada: INTRODUCING HARMONIC DISTRIBUTION IN WIKIPEDIA Alumne: V´ıctor Mart´ınez Ruiz Tutor: Dr. Pere Puig Metadata T´ıtol treball: Introducing Harmonic Distribution in Wikipedia Alumne: V´ıctor Mart´ınez Ruiz Tutor: Pere Puig Abstract Catal`a Volem estudiar la Distribuci´oHarm`onica, aleshores el que farem amb aquest treball ´es implementar aquesta distribuci´oa Wikipedia on descriurem la seva funci´ode densitat, la funci´ode distribuci´o,els moments i els estimadors dels par`ametres entre d'altres propietats. Castellano Queremos estudiar la Distribuci´onHarm´onica, entonces lo que haremos en este trabajo ser´aimplementar ´esta distribuci´onen Wikipedia donde describiremos su funci´onde densidad , la funci´onde distribuci´on, los momentos y los estimadores de los par´ametros entre otras propiedades. English We want to study the Harmonic Distribution. To do this we will implement this distribution in Wikipedia where we will describe its probability density function, its cumulative distribution function, the moments and some parameter estimators among other properties. Paraules claus: Catal`a Funci´ode Densitat de Probabilitat, Funci´ode Disribuci´o,mediana, estimaci´ode par`ametres. Castellano Funci´onde Densidad de Probabilidades, Funci´onde Distribuci´on,mediana, estimaci´onde par´ametros. English Probability Density Function, Cumulative Distribution Function, median, estimation of parameters. 1 Contents 1 Introduction 3 2 History
    [Show full text]
  • Estimation of Harmonics, Interharmonics and Sub-Harmonics in Motor Drive Systems
    ESTIMATION OF HARMONICS, INTERHARMONICS AND SUB-HARMONICS IN MOTOR DRIVE SYSTEMS by Robert John Kerestes B.S. in Electrical Engineering, University of Pittsburgh, 2010 M.S. in Electrical Engineering, University of Pittsburgh, 2011 Submitted to the Graduate Faculty of Swanson School of Engineering in partial fulfillment of the requirements for the degree of Doctor of Philosophy University of Pittsburgh 2014 UNIVERSITY OF PITTSBURGH SWANSON SCHOOL OF ENGINEERING This dissertation was presented by Robert John Kerestes It was defended on November 15th , 2013 and approved by Zhi-Hong Mao, Ph.D., Associate Professor, Electrical and Computer Engineering Thomas McDermott, Ph.D., Assistant Professor, Electrical and Computer Engineering Chengtao Wen, Ph.D., Emerson Process Management Dissertation Director: Gregory F. Reed, Ph.D., Associate Professor, Electrical and Computer Engineering ii Copyright © by Robert J. Kerestes 2014 iii ESTIMATION OF HARMONICS, INTERHARMONICS AND SUB-HARMONICS IN MOTOR DRIVE SYSTEMS Robert Kerestes, Ph.D. University of Pittsburgh, 2014 Harmonic analysis in power systems is a challenge that is always evolving due to a constantly changing power grid. Advances in power electronics have led to the installation of many new non-linear power loads, such as power converters, arc furnaces and motor drive systems, which in turn change the existing infrastructure and inject new types of harmonics. This study analyzes a static Kramer drive system for its harmonic behavior. The system is broken down component- by-component and analyzed. The analysis includes classical harmonics of integer multiple of the fundamental frequency as well as sub-harmonics and interharmonics. This paper shows that there are significant harmonic issues in a static Kramer drive that are likely to be present in similar drive systems.
    [Show full text]
  • Second Harmonic Revisited: an Analytic Quantum Approach
    applied sciences Article Second Harmonic Revisited: An Analytic Quantum Approach Giovanni Chesi 1,*, Matteo M. Wauters 2 , Nadir Fasola 2 , Alessia Allevi 1,3 and Maria Bondani 3 1 Department of Science and High Technology, University of Insubria, Via Valleggio 11, I-22100 Como, Italy; [email protected] 2 SISSA, Via Bonomea 265, I-34136 Trieste, Italy; [email protected] (M.M.W.); [email protected] (N.F.) 3 Institute for Photonics and Nanotechnologies, CNR, Via Valleggio 11, I-22100 Como, Italy; [email protected] * Correspondence: [email protected]; Tel.: +39-031-238-6252 Received: 29 March 2019; Accepted: 17 April 2019; Published: 24 April 2019 Abstract: We address the second-harmonic generation process in a quantum frame. Starting with a perturbative approach, we show that it is possible to achieve a number of analytic results, ranging from the up-conversion probability to the statistical properties of the generated light. In particular, the moments and the correlations of the photon-number distribution of the second-harmonic light generated by any initial state are retrieved. When possible, a comparison with the results achieved with the classical regime is successfully provided. The nonclassicality of some benchmark states is investigated by inspecting the corresponding autocorrelation function. Keywords: second-harmonic generation; quantum description of interaction of light and matter; photon statistics and coherence theory 1. Introduction Following the first observation of second-harmonic generation (SHG) by Franken et al. [1], the process has been the subject of a huge amount of theoretical and experimental investigations. The exact solution of the equation of motion for the classical field amplitudes, provided by Armstrong et al.
    [Show full text]
  • Arxiv:1805.04695V2 [Nucl-Th] 7 Mar 2019 an Analysis More Sophisticated Than a Fourier Analysis
    Non-Bessel-Gaussianity and Flow Harmonic Fine-Splitting Hadi Mehrabpour1, 2, ∗ and Seyed Farid Taghavi2, 3, y 1Department of Physics, Sharif University of Technology, P.O. Box 11155-9161, Tehran, Iran 2School of Particles and Accelerators, Institute for Research in Fundamental Sciences (IPM), P.O. Box 19395-5531, Tehran, Iran 3Physik Department, Technische Universit¨atM¨unchen,Munich, Germany Both collision geometry and event-by-event fluctuations are encoded in the experimentally ob- served flow harmonic distribution p(vn) and 2k-particle cumulants cnf2kg. In the present study, we systematically connect these observables to each other by employing Gram-Charlier A series. We quantify the deviation of p(vn) from Bessel-Gaussianity in terms of flow harmonic fine-splitting. Sub- sequently, we show that the corrected Bessel-Gaussian distribution can fit the simulated data better than the Bessel-Gaussian distribution in the more peripheral collisions. Inspired by Gram-Charlier A series, we introduce a new set of cumulants qnf2kg that are more natural to study distributions near Bessel-Gaussian. These new cumulants are obtained from cnf2kg where the collision geometry effect is extracted from it. By exploiting q2f2kg, we introduce a new set of estimators for averaged ellipticityv ¯2 which are more accurate compared to v2f2kg for k > 1. As another application of q2f2kg, we show we are able to restrict the phase space of v2f4g, v2f6g and v2f8g by demanding the consistency ofv ¯2 and v2f2kg with q2f2kg equation. The allowed phase space is a region such that v2f4g − v2f6g & 0 and 12v2f6g − 11v2f8g − v2f4g & 0, which is compatible with the experimental observations.
    [Show full text]