Supplementary Material Comparative assessment of the ELBA coarse-grained model for water
|
|
- Mitchell Bennett
- 5 years ago
- Views:
Transcription
1 Supplementary Material Comparative assessment of the oarse-grained model for water Mario Orsi Shool of Engineering & Materials Siene, Queen Mary University of London Mile End Road, London E 4NS, United Kingdom m.orsi@qmul.a.uk Shifted-fore interations The total potential energy U ij of an interating pair of water sites i,j is: U ij = U LJ ij +U dip ij, with U LJ ij the Lennard-Jones term and U dip ij the point dipole term. Shifted-fore Lennard-Jones interation The shifted-fore form of the Lennard-Jones potential is []: { [(σ 2 ( σ ) ] [ ( ) 6 2 ( ) ] 6 ( ) 2 ( ) 2 ( ) } 6 σ σ r σ σ Uij LJ = 4ǫ , r) r where r is the interpartile distane, is the utoff radius, and σ and ǫ have the standard meaning [2]. The orresponding fore is given by: { [ ( σ ) 2 ( σ ) ] [ ( ) 6 2 ( ) ] } 6 σ σ f ij = 48ǫ 24ǫ r r r 2 48ǫ 24ǫ r where r is the pair distane vetor, r = r i r j. Shifted-fore point dipole interation Considering two point dipoles µ i and µ j, the shifted-fore pair potential energy is: [ ( ) 3 ( ) ] 4 [ U dip r r ij = 4 +3 r 3(µ i µ j ) 3 ] r 5(µ i r)(µ j r) where r and r are the pair distane vetor and its magnitude, respetively, and is the utoff distane. The standard fator (4πε 0 ) is assumed, and omitted folarity here and in the following formulae. Note that we use the onvention r = r i r j. Fores and torques were derived following the approah in Appendix C.3 of Allen and Tildesley [2]. In partiular, the pair fores are: f ij = f ji = 3 ( ) ] 4 [ r {[ r 5 (µ i µ j ) 3 ] r 2(µ i r)(µ j r) r + The pair torque τ ij is: τ ij = ( r [ 4 r 3 [ ( ) 3 ( r r 4 +3 The pair torque τ ji is: τ ji = ( r [ 4 r 3 ) 4 ] [ (µ j r)µ i +(µ i r)µ j 2 r 2(µ i r)(µ j r)r ) 3 ( ) [ 4 r +3 ](µ r i µ j )+ 3r ( ) 3 ( ) ] 4 r r (µ r j r)(µ i r) ) 3 ( ) [ 4 r +3 ](µ r j µ i )+ 3r ( ) 3 ( ) ] 4 r r (µ r i r)(µ j r) The expressions above, for both the Lenard-Jones and the point dipole potential, have been implemented in the pair lj sf dipole sf.pp module of the program LAMMPS [3]. r 2 ] }
2 Integration timestep and energy onservation Compared to previous publiations [4, 5], in this work we have redued the inertial features of the water model. In partiular, we have here set a realisti mass of 8gmol, and a prinipal moment of inertia of 30gÅ 2 mol. In the original model [4], values of respetively 40gmol and 00gÅ 2 mol were hosen to permit a stable integration of the equations of motion with a 5fs timestep; this hoie was made to maximise the effiieny of serial simulations, as at the time the in-house software used did not allow parallel alulations. Sine the model is now available in the parallel program LAMMPS [3], maximising the timestep is less ritial, and we have deided to redue the timestep size to permit the setting of a realisti (lower) mass. The timestep hoie was validated by running NVE simulations and monitoring the total energy. In partiular, Figure S reports the time evolution of the total energy for a range of timesteps from 5 to 5fs. Total energy / (kal mol - ) t = 5fs t = 2.5fs t = 0fs t = 7.5fs t = 5fs Simulation time / ps Figure S: Total energy onservation for different integration timesteps. Folarity, the different data sets were shifted up or down to avoid overlaps. It an be seen that the total energy is onserved aurately up to t = 0fs. A larger timestep of 2.5fs ould still be aeptable; over 30 ps, while the energy osillates slightly (within 0.%), there is no net drift. For t = 5fs, the total energy is learly not onserved. The flutuations in potential, kineti, and total energy as a funtion of the t are displayed in Figure S2. 00 Potential energy flutuation Kineti energy flutuation Total energy flutuation Energy flutuation E / (kal mol - ) Integration timestep t / fs Figure S2: Average energy flutuation E = (E E ) 2. For timesteps up to 2.5fs, the average flutuation in the total energy is at least 0times smaller than the average flutuations in both the potential and kineti energy, satisfying reommended riteria [6]. 2
3 Overall, ouhoie of t = 0 fs for prodution simulations should guarantee aurate and robust integration of the equations of motion. Computational timings Simulations of a liquid-vapour interfae omprising 8000 water moleules were run for the single-site oarse-grained model, a representative 3-site model (), and a representative 4-site model (). For eah model, the same parameters reported in the main body of the paper were used; in partiular, the timestep was 0fs for and 2fs for the atomisti models, and the utoff radii were.2,.0 and.3 nm for respetively, and. The alulations were run on a GNU/Linux luster of 2.66 GHz Intel Westmere proessors based on 2 sokets and 2 proessors per node, with QLogi TrueSale interonnet. Figure S3 reports a omparison in terms of sampling speed, defined as number of ns sampled / day, where day is intended in terms of alulation (wall-lok) time. 000 Sampling speed [ns/day] Number of proessors Figure S3: Comparison of sampling speed between representative water models. Simulations were run on 2, 24, 48 and 96 proessors. Note the logarithmi sale used to represent the sampling speed. It an be seen that the speed up fators in the simulations are respetively 40 and 00 ompared to and. It is also lear that suh fators are rather insensitive to the number of proessors employed. Mean square displaement and time-dependent diffusion urves The mean square displaement (M SD) was obtained with: MSD = N N [r i (t) r i (0)] 2 () i= with N the number of water moleules, t the time, and r i the 3-dimensional vetor defining the position of partile i. Figure S4 shows the MSD from the liquid bulk simulations of eah model. The self-diffusion oeffiient D was obtained from the bulk systems using the standard expression [7]: D = MSD 6t where MSD is defined in Eq. and t is the simulation time; see Fig. S5. (2) 3
4 5 0 /E MSD / nm Time / ns Figure S4: Mean square displaement as a funtion of simulation time. 4.5 D / (0-9 m 2 /s) /E Time / ns Figure S5: Self-diffusion oeffiient as a funtion of simulation time. 4
5 Radial distribution funtion: main peaks g OO (r ) r /Å g OO (r 2 ) r 2 /Å /E Experiment [8] Table S: Properties of the radial distribution funtion for bulk water at 298K and atm. Height of first peak: g OO (r ). Loation of first peak: r. Height of seond peak: g OO (r 2 ). Loation of seond peak: r 2. Relative unertainties are < 0.%. Mass density profiles ρ / g m /E Figure S6: Mass density profiles. It an be notied that the TIP4P-like models feature small yet nonnegligible density peaks at the liquid-vapour interfaes. In previous work, suh peaks were observed to disappear with inreasing system size [9]. 5
6 ρ / g m /E e-05 e Figure S7: Mass density profiles. To failitate interpretation in the vapour regions, the density is displayed using a logarithmi sale. For the liquid region, see Figure S6. Eletrostati potential and eletri field profiles φ / V /E Figure S8: Eletrostati potential profiles. Foonveniene, a representative mass density profile is also shown with a dotted line (ordinate values not on sale). 6
7 2.5 E z / 0 9 V m /E Figure S9: Eletri field profiles. Foonveniene, a representative mass density profile is also shown with a dotted line (ordinate values not on sale). Referenes [] S.D. Stoddard and J. Ford, Phys. Rev. A 8, 504 (973). [2] M.P. Allen and D.J. Tildesley, Computer Simulation of Liquids, st ed. (Oxford Siene Publiations, Oxford, 987). [3] S. Plimpton, J. Comp. Phys. 7, (995), [4] M. Orsi and J.W. Essex, PLoS ONE 6, e28637 (20). [5] M. Orsi and J.W. Essex, Faraday Disuss. 6, 249 (203). [6] S. Riniker and W.F. van Gunsteren, J. Chem. Phys. 34, 0840 (20). [7] D.C. Rapaport, The Art of Moleular Dynamis Simulation, 2nd ed. (Cambridge University Press, Cambridge, 2004). [8] A. Soper, Chem. Phys. 258, 2 (2000), [9] J. Alejandre and G.A. Chapela, J. Chem. Phys. 32, 0470 (200). 7
arxiv: v1 [physics.gen-ph] 5 Jan 2018
The Real Quaternion Relativity Viktor Ariel arxiv:1801.03393v1 [physis.gen-ph] 5 Jan 2018 In this work, we use real quaternions and the basi onept of the final speed of light in an attempt to enhane the
More informationA model for measurement of the states in a coupled-dot qubit
A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:
More informationA simple expression for radial distribution functions of pure fluids and mixtures
A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.
More informationSupercritical Water Confined In Graphene Nanochannels
Superritial Water Confined In Graphene Nanohannels J. Sala*, E. Guàrdia, and J. Martí Departament de Físia i Enginyeria Nulear, Universitat Politènia de Catalunya, B4-B5 Campus Nord, 08034 Barelona, Catalonia,
More informationEvaluation of effect of blade internal modes on sensitivity of Advanced LIGO
Evaluation of effet of blade internal modes on sensitivity of Advaned LIGO T0074-00-R Norna A Robertson 5 th Otober 00. Introdution The urrent model used to estimate the isolation ahieved by the quadruple
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 4
Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas
More informationMOLECULAR DYNAMICS SIMULATIONS OF HEMATITE NANOPARTICLE DEPOSITION ONTO A SUBSTRATE
MOLECULAR DYNAMICS SIMULATIONS OF HEMATITE NANOPARTICLE DEPOSITION ONTO A SUBSTRATE By DANA R. LOUIE B.S., Mehanial Engineering, Ohio University, 1989 M.S., Aerospae Engineering, University of Florida,
More informationChapter 9. The excitation process
Chapter 9 The exitation proess qualitative explanation of the formation of negative ion states Ne and He in He-Ne ollisions an be given by using a state orrelation diagram. state orrelation diagram is
More informationEffect of magnetization process on levitation force between a superconducting. disk and a permanent magnet
Effet of magnetization proess on levitation fore between a superonduting disk and a permanent magnet L. Liu, Y. Hou, C.Y. He, Z.X. Gao Department of Physis, State Key Laboratory for Artifiial Mirostruture
More informationSupporting Information
Supporting Information Olsman and Goentoro 10.1073/pnas.1601791113 SI Materials Analysis of the Sensitivity and Error Funtions. We now define the sensitivity funtion Sð, «0 Þ, whih summarizes the steepness
More informationMATHEMATICAL AND NUMERICAL BASIS OF BINARY ALLOY SOLIDIFICATION MODELS WITH SUBSTITUTE THERMAL CAPACITY. PART II
Journal of Applied Mathematis and Computational Mehanis 2014, 13(2), 141-147 MATHEMATICA AND NUMERICA BAI OF BINARY AOY OIDIFICATION MODE WITH UBTITUTE THERMA CAPACITY. PART II Ewa Węgrzyn-krzypzak 1,
More informationSimple Considerations on the Cosmological Redshift
Apeiron, Vol. 5, No. 3, July 8 35 Simple Considerations on the Cosmologial Redshift José Franiso Garía Juliá C/ Dr. Maro Mereniano, 65, 5. 465 Valenia (Spain) E-mail: jose.garia@dival.es Generally, the
More informationELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES.
ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES. All systems with interation of some type have normal modes. One may desribe them as solutions in absene of soures; they are exitations of the system
More informationDeveloping Excel Macros for Solving Heat Diffusion Problems
Session 50 Developing Exel Maros for Solving Heat Diffusion Problems N. N. Sarker and M. A. Ketkar Department of Engineering Tehnology Prairie View A&M University Prairie View, TX 77446 Abstrat This paper
More informationGUIDED WAVE ENERGY DISTRIBUTION ANALYSIS IN INHOMOGENEOUS PLATES
GUDED WAVE ENERGY DSTRBUTON ANALYSS N NHOMOGENEOUS PLATES Krishnan Balasubramaniam and Yuyin Ji Department of Aerospae Engineering and Mehanis Mississippi State University, Mississippi State, MS 39762
More informationGyrokinetic calculations of the neoclassical radial electric field in stellarator plasmas
PHYSICS OF PLASMAS VOLUME 8, NUMBER 6 JUNE 2001 Gyrokineti alulations of the neolassial radial eletri field in stellarator plasmas J. L. V. Lewandowski Plasma Physis Laboratory, Prineton University, P.O.
More informationCopyright 2018 Society of Photo-Optical Instrumentation Engineers (SPIE). One print or electronic copy may be made for personal use only.
Copyright 018 Soiety of Photo-Optial Instrumentation Engineers (SPIE) One print or eletroni opy may be made for personal use only Systemati reprodution and distribution, dupliation of any material in this
More informationElectromagnetic radiation of the travelling spin wave propagating in an antiferromagnetic plate. Exact solution.
arxiv:physis/99536v1 [physis.lass-ph] 15 May 1999 Eletromagneti radiation of the travelling spin wave propagating in an antiferromagneti plate. Exat solution. A.A.Zhmudsky November 19, 16 Abstrat The exat
More informationGluing Potential Energy Surfaces with Rare Event Simulations
This is an open aess artile published under an ACS AuthorChoie Liense, whih permits opying and redistribution of the artile or any adaptations for non-ommerial purposes. pubs.as.org/jctc Gluing Potential
More informationModeling of Threading Dislocation Density Reduction in Heteroepitaxial Layers
A. E. Romanov et al.: Threading Disloation Density Redution in Layers (II) 33 phys. stat. sol. (b) 99, 33 (997) Subjet lassifiation: 6.72.C; 68.55.Ln; S5.; S5.2; S7.; S7.2 Modeling of Threading Disloation
More informationAnalysis of discretization in the direct simulation Monte Carlo
PHYSICS OF FLUIDS VOLUME 1, UMBER 1 OCTOBER Analysis of disretization in the diret simulation Monte Carlo iolas G. Hadjionstantinou a) Department of Mehanial Engineering, Massahusetts Institute of Tehnology,
More informationMillennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion
Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six
More informationEFFECTIVE STRESS LAW FOR THE PERMEABILITY OF CLAY-RICH SANDSTONES
SCA22-5 1/6 EFFECTIVE STRESS LAW FOR THE PERMEABILITY OF CLAY-RICH SANDSTONES Widad Al-Wardy and Robert W. Zimmerman Department of Earth Siene and Engineering Imperial College of Siene, Tehnology and Mediine
More informationWhere as discussed previously we interpret solutions to this partial differential equation in the weak sense: b
Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential
More informationNumerical Tests of Nucleation Theories for the Ising Models. Abstract
to be submitted to Physial Review E Numerial Tests of Nuleation Theories for the Ising Models Seunghwa Ryu 1 and Wei Cai 2 1 Department of Physis, Stanford University, Stanford, California 94305 2 Department
More informationChapter 26 Lecture Notes
Chapter 26 Leture Notes Physis 2424 - Strauss Formulas: t = t0 1 v L = L0 1 v m = m0 1 v E = m 0 2 + KE = m 2 KE = m 2 -m 0 2 mv 0 p= mv = 1 v E 2 = p 2 2 + m 2 0 4 v + u u = 2 1 + vu There were two revolutions
More informationthe following action R of T on T n+1 : for each θ T, R θ : T n+1 T n+1 is defined by stated, we assume that all the curves in this paper are defined
How should a snake turn on ie: A ase study of the asymptoti isoholonomi problem Jianghai Hu, Slobodan N. Simić, and Shankar Sastry Department of Eletrial Engineering and Computer Sienes University of California
More informationA NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE. W. J.
A NEW FLEXIBLE BODY DYNAMIC FORMULATION FOR BEAM STRUCTURES UNDERGOING LARGE OVERALL MOTION IIE THREE-DIMENSIONAL CASE W. J. Haering* Senior Projet Engineer General Motors Corporation Warren, Mihigan R.
More informationTime Domain Method of Moments
Time Domain Method of Moments Massahusetts Institute of Tehnology 6.635 leture notes 1 Introdution The Method of Moments (MoM) introdued in the previous leture is widely used for solving integral equations
More information23.1 Tuning controllers, in the large view Quoting from Section 16.7:
Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output
More informationUTC. Engineering 329. Proportional Controller Design. Speed System. John Beverly. Green Team. John Beverly Keith Skiles John Barker.
UTC Engineering 329 Proportional Controller Design for Speed System By John Beverly Green Team John Beverly Keith Skiles John Barker 24 Mar 2006 Introdution This experiment is intended test the variable
More informationRelativity in Classical Physics
Relativity in Classial Physis Main Points Introdution Galilean (Newtonian) Relativity Relativity & Eletromagnetism Mihelson-Morley Experiment Introdution The theory of relativity deals with the study of
More informationMost results in this section are stated without proof.
Leture 8 Level 4 v2 he Expliit formula. Most results in this setion are stated without proof. Reall that we have shown that ζ (s has only one pole, a simple one at s =. It has trivial zeros at the negative
More informationThe Hanging Chain. John McCuan. January 19, 2006
The Hanging Chain John MCuan January 19, 2006 1 Introdution We onsider a hain of length L attahed to two points (a, u a and (b, u b in the plane. It is assumed that the hain hangs in the plane under a
More information11.4 Molecular Orbital Description of the Hydrogen Molecule Electron Configurations of Homonuclear Diatomic Molecules
Chap Moleular Eletroni Struture Table of Contents. The orn-oppenheimer pproximation -. The Hydrogen Moleule Ion.3 Calulation of the Energy of the Hydrogen Moleule Ion.4 Moleular Orbital Desription of the
More informationRESEARCH ON RANDOM FOURIER WAVE-NUMBER SPECTRUM OF FLUCTUATING WIND SPEED
The Seventh Asia-Paifi Conferene on Wind Engineering, November 8-1, 9, Taipei, Taiwan RESEARCH ON RANDOM FORIER WAVE-NMBER SPECTRM OF FLCTATING WIND SPEED Qi Yan 1, Jie Li 1 Ph D. andidate, Department
More informationTemperature-Gradient-Driven Tearing Modes
1 TH/S Temperature-Gradient-Driven Tearing Modes A. Botrugno 1), P. Buratti 1), B. Coppi ) 1) EURATOM-ENEA Fusion Assoiation, Frasati (RM), Italy ) Massahussets Institute of Tehnology, Cambridge (MA),
More informationTutorial 8: Solutions
Tutorial 8: Solutions 1. * (a) Light from the Sun arrives at the Earth, an average of 1.5 10 11 m away, at the rate 1.4 10 3 Watts/m of area perpendiular to the diretion of the light. Assume that sunlight
More informationChapter Outline The Relativity of Time and Time Dilation The Relativistic Addition of Velocities Relativistic Energy and E= mc 2
Chapter 9 Relativeity Chapter Outline 9-1 The Postulate t of Speial Relativity it 9- The Relativity of Time and Time Dilation 9-3 The Relativity of Length and Length Contration 9-4 The Relativisti Addition
More informationOn the Quantum Theory of Radiation.
Physikalishe Zeitshrift, Band 18, Seite 121-128 1917) On the Quantum Theory of Radiation. Albert Einstein The formal similarity between the hromati distribution urve for thermal radiation and the Maxwell
More informationRelativistic Dynamics
Chapter 7 Relativisti Dynamis 7.1 General Priniples of Dynamis 7.2 Relativisti Ation As stated in Setion A.2, all of dynamis is derived from the priniple of least ation. Thus it is our hore to find a suitable
More informationPure Component Phase Diagram. Definitions. Definitions (cont.) Class 17 Non-Ideal Gases
Class 17 Non-Ideal Gases Definitions Critial emperature, ressure Vapor Gas Van der Waals EOS Other Equations of State Compressibility Fator riniple of Corresponding States Kay s Rule Water hase Change
More informationH=250 Oe. Supplementary Figure 1 Magnetic domains: Room temperature 4 x 4 µm 2 MFM phase
1 Supplementary Information Supplementary Figures (b) () 1.6± 1 µm 1 µm 1 µm H290 Oe (d) (e) (f) H410 Oe -1.5± 1 µm 1 µm H250 Oe 1 µm H590 Oe Supplementary Figure 1 Magneti domains: Room temperature 4
More informationThe Working Man s Guide to Obtaining Self Diffusion Coefficients from Molecular Dynamics Simulations
The Working Man s Guide to Obtaining Self Diffusion Coeffiients from Moleular Dynamis Simulations David Keffer Department of Chemial Engineering University of Tennessee, Knoxville Begun: November 11, 2001
More informationCRITICAL EXPONENTS TAKING INTO ACCOUNT DYNAMIC SCALING FOR ADSORPTION ON SMALL-SIZE ONE-DIMENSIONAL CLUSTERS
Russian Physis Journal, Vol. 48, No. 8, 5 CRITICAL EXPONENTS TAKING INTO ACCOUNT DYNAMIC SCALING FOR ADSORPTION ON SMALL-SIZE ONE-DIMENSIONAL CLUSTERS A. N. Taskin, V. N. Udodov, and A. I. Potekaev UDC
More informationEXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION
Journal of Mathematial Sienes: Advanes and Appliations Volume 3, 05, Pages -3 EXACT TRAVELLING WAVE SOLUTIONS FOR THE GENERALIZED KURAMOTO-SIVASHINSKY EQUATION JIAN YANG, XIAOJUAN LU and SHENGQIANG TANG
More informationNon-Markovian study of the relativistic magnetic-dipole spontaneous emission process of hydrogen-like atoms
NSTTUTE OF PHYSCS PUBLSHNG JOURNAL OF PHYSCS B: ATOMC, MOLECULAR AND OPTCAL PHYSCS J. Phys. B: At. Mol. Opt. Phys. 39 ) 7 85 doi:.88/953-75/39/8/ Non-Markovian study of the relativisti magneti-dipole spontaneous
More informationNatural Convection Experiment Measurements from a Vertical Surface
OBJECTIVE Natural Convetion Experiment Measurements from a Vertial Surfae 1. To demonstrate te basi priniples of natural onvetion eat transfer inluding determination of te onvetive eat transfer oeffiient.
More information8.333: Statistical Mechanics I Problem Set # 4 Due: 11/13/13 Non-interacting particles
8.333: Statistial Mehanis I Problem Set # 4 Due: 11/13/13 Non-interating partiles 1. Rotating gas: Consider a gas of N idential atoms onfined to a spherial harmoni trap in three dimensions, i.e. the partiles
More informationMOLECULAR ORBITAL THEORY- PART I
5.6 Physial Chemistry Leture #24-25 MOLECULAR ORBITAL THEORY- PART I At this point, we have nearly ompleted our rash-ourse introdution to quantum mehanis and we re finally ready to deal with moleules.
More informationIN-PLANE VIBRATIONS OF CURVED BEAMS WITH VARIABLE CROSS-SECTIONS CARRYING ADDITIONAL MASS
11 th International Conferene on Vibration Problems Z. Dimitrovová et al. (eds.) Lisbon, Portugal, 9-1 September 013 IN-PLANE VIBRATIONS OF CURVED BEAMS WITH VARIABLE CROSS-SECTIONS CARRYING ADDITIONAL
More informationMassively Parallel Simulations of Diffusion in Dense Polymeric Structures
Massively Parallel Simulations of Diffusion in Dense Polymeri Strutures Jean-Loup Faulon Computational Materials Siene Sandia National Laboratories MS 1111 Albuquerque, New Mexio 87 185-1111 Jean-Lmp Faulon
More informationApplication of the Dyson-type boson mapping for low-lying electron excited states in molecules
Prog. Theor. Exp. Phys. 05, 063I0 ( pages DOI: 0.093/ptep/ptv068 Appliation of the Dyson-type boson mapping for low-lying eletron exited states in moleules adao Ohkido, and Makoto Takahashi Teaher-training
More informationWe consider the nonrelativistic regime so no pair production or annihilation.the hamiltonian for interaction of fields and sources is 1 (p
.. RADIATIVE TRANSITIONS Marh 3, 5 Leture XXIV Quantization of the E-M field. Radiative transitions We onsider the nonrelativisti regime so no pair prodution or annihilation.the hamiltonian for interation
More informationExperiment 03: Work and Energy
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physis Department Physis 8.01 Purpose of the Experiment: Experiment 03: Work and Energy In this experiment you allow a art to roll down an inlined ramp and run into
More informationA generalized equation of state with an application to the Earth s Mantle
Geofísia Internaional 49 (), 77-8 (010) A generalized equation of state with an appliation to the Earth s Mantle J. A. Robles-Gutiérrez 1, J. M. A. Robles-Domínguez 1 and C. Lomnitz 1 Universidad Autónoma
More informationSlenderness Effects for Concrete Columns in Sway Frame - Moment Magnification Method
Slenderness Effets for Conrete Columns in Sway Frame - Moment Magnifiation Method Slender Conrete Column Design in Sway Frame Buildings Evaluate slenderness effet for olumns in a sway frame multistory
More informationThe Corpuscular Structure of Matter, the Interaction of Material Particles, and Quantum Phenomena as a Consequence of Selfvariations.
The Corpusular Struture of Matter, the Interation of Material Partiles, and Quantum Phenomena as a Consequene of Selfvariations. Emmanuil Manousos APM Institute for the Advanement of Physis and Mathematis,
More informationLECTURE 2 Geometrical Properties of Rod Cross Sections (Part 2) 1 Moments of Inertia Transformation with Parallel Transfer of Axes.
V. DEMENKO MECHNCS OF MTERLS 05 LECTURE Geometrial Properties of Rod Cross Setions (Part ) Moments of nertia Transformation with Parallel Transfer of xes. Parallel-xes Theorems S Given: a b = S = 0. z
More informationQUANTUM MECHANICS II PHYS 517. Solutions to Problem Set # 1
QUANTUM MECHANICS II PHYS 57 Solutions to Problem Set #. The hamiltonian for a lassial harmoni osillator an be written in many different forms, suh as use ω = k/m H = p m + kx H = P + Q hω a. Find a anonial
More informationLecture 4. Kinetics and Transport in Sorption (1)
Leture 4. Kinetis and Transport in Sorption () Mass-Transport Proesses in Adsorption Convetion-Dispersion Model Modes of Time-Dependent Sorption Operation Chromatogram Solute Conentration Distriution Separation
More informationFailure Assessment Diagram Analysis of Creep Crack Initiation in 316H Stainless Steel
Failure Assessment Diagram Analysis of Creep Crak Initiation in 316H Stainless Steel C. M. Davies *, N. P. O Dowd, D. W. Dean, K. M. Nikbin, R. A. Ainsworth Department of Mehanial Engineering, Imperial
More information1 sin 2 r = 1 n 2 sin 2 i
Physis 505 Fall 005 Homework Assignment #11 Solutions Textbook problems: Ch. 7: 7.3, 7.5, 7.8, 7.16 7.3 Two plane semi-infinite slabs of the same uniform, isotropi, nonpermeable, lossless dieletri with
More informationarxiv:gr-qc/ v2 6 Feb 2004
Hubble Red Shift and the Anomalous Aeleration of Pioneer 0 and arxiv:gr-q/0402024v2 6 Feb 2004 Kostadin Trenčevski Faulty of Natural Sienes and Mathematis, P.O.Box 62, 000 Skopje, Maedonia Abstrat It this
More informationQCLAS Sensor for Purity Monitoring in Medical Gas Supply Lines
DOI.56/sensoren6/P3. QLAS Sensor for Purity Monitoring in Medial Gas Supply Lines Henrik Zimmermann, Mathias Wiese, Alessandro Ragnoni neoplas ontrol GmbH, Walther-Rathenau-Str. 49a, 7489 Greifswald, Germany
More information( ) ( ) Volumetric Properties of Pure Fluids, part 4. The generic cubic equation of state:
CE304, Spring 2004 Leture 6 Volumetri roperties of ure Fluids, part 4 The generi ubi equation of state: There are many possible equations of state (and many have been proposed) that have the same general
More information( ) which is a direct consequence of the relativistic postulate. Its proof does not involve light signals. [8]
The Speed of Light under the Generalized Transformations, Inertial Transformations, Everyday Clok Synhronization and the Lorentz- Einstein Transformations Bernhard Rothenstein Abstrat. Starting with Edwards
More informationCollinear Equilibrium Points in the Relativistic R3BP when the Bigger Primary is a Triaxial Rigid Body Nakone Bello 1,a and Aminu Abubakar Hussain 2,b
International Frontier Siene Letters Submitted: 6-- ISSN: 9-8, Vol., pp -6 Aepted: -- doi:.8/www.sipress.om/ifsl.. Online: --8 SiPress Ltd., Switzerland Collinear Equilibrium Points in the Relativisti
More informationDiffusion and compensating subsidence as limiting cases of a single flux parameterization. David Randall
! Revised September 014 9:49 AM! 1 Diffusion and ompensating subsidene as limiting ases of a single flux parameterization David Randall The tendeny of the average of a onservative variable h due to a vertial
More informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More informationThe Concept of Mass as Interfering Photons, and the Originating Mechanism of Gravitation D.T. Froedge
The Conept of Mass as Interfering Photons, and the Originating Mehanism of Gravitation D.T. Froedge V04 Formerly Auburn University Phys-dtfroedge@glasgow-ky.om Abstrat For most purposes in physis the onept
More informationComparison of PD and LQR Methods for Spacecraft Attitude Control Using Star Trackers
Comparison of PD and LQ Methods for Spaeraft Attitude Control Using Star raers Sott Beatt, Universit of New Meio, United States, sott.beatt@pangeateh.om ABSAC he wor ontained herein is a omparison of spaeraft
More informationPhysics; Watching the Game From the Outside
Physis; Wathing the Game From the Outside Roald C. Maximo Feb It is a good thing to have two ways of looking at a subjet, and also admit that there are two ways of looking at it. James Clerk Maxwell, on
More informationDiscontinuity of Shannon information entropy for two-electron atoms
Chemial Physis 39 (25) 127 131 www.elsevier.om/loate/hemphys Disontinuity of Shannon information entropy for two-eletron atoms Qiun Shi, Sabre Kais * Department of Chemistry, Purdue University, 1393 Brown
More informationClassical Diamagnetism and the Satellite Paradox
Classial Diamagnetism and the Satellite Paradox 1 Problem Kirk T. MDonald Joseph Henry Laboratories, Prineton University, Prineton, NJ 08544 (November 1, 008) In typial models of lassial diamagnetism (see,
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION The mode-loking transition of random lasers Maro Leonetti, Claudio Conti, 2 and Cefe Lopez Instituto de Cienia de Materiales de Madrid (CSIC) and Unidad Asoiada CSIC-UVigo, Cantoblano
More informationDynamics of the Electromagnetic Fields
Chapter 3 Dynamis of the Eletromagneti Fields 3.1 Maxwell Displaement Current In the early 1860s (during the Amerian ivil war!) eletriity inluding indution was well established experimentally. A big row
More informationBeams on Elastic Foundation
Professor Terje Haukaas University of British Columbia, Vanouver www.inrisk.ub.a Beams on Elasti Foundation Beams on elasti foundation, suh as that in Figure 1, appear in building foundations, floating
More informationAccelerator Physics Particle Acceleration. G. A. Krafft Old Dominion University Jefferson Lab Lecture 4
Aelerator Physis Partile Aeleration G. A. Krafft Old Dominion University Jefferson Lab Leture 4 Graduate Aelerator Physis Fall 15 Clarifiations from Last Time On Crest, RI 1 RI a 1 1 Pg RL Pg L V Pg RL
More informationExperimental Investigation and FE Analysis of Fiber Woven Layered Composites under Dynamic Loading
2th International LS-DYNA Users Conferene Constitutive Modeling(2) xperimental Investigation and F Analysis of Fiber Woven Layered Composites under Dynami Loading Pavel A. Mossakovsky, Fedor K. Antonov,
More informationV. Interacting Particles
V. Interating Partiles V.A The Cumulant Expansion The examples studied in the previous setion involve non-interating partiles. It is preisely the lak of interations that renders these problems exatly solvable.
More informationChapter 3. Volumetric Properties of Pure Fluids
Chapter 3. olumetri roperties of ure Fluids Introdution hermodynami properties (U, H and thus Q, W) are alulated from data data are important for sizing vessels and pipelines Subjets behavior of pure fluids
More informationPHYSICS 432/532: Cosmology Midterm Exam Solution Key (2018) 1. [40 points] Short answer (8 points each)
PHYSICS 432/532: Cosmology Midterm Exam Solution Key (2018) 1. [40 points] Short answer (8 points eah) (a) A galaxy is observed with a redshift of 0.02. How far away is the galaxy, and what is its lookbak
More informationSTUDY OF INTERFACIAL BEHAVIOR OF CNT/POLYMER COMPOSITE BY CFE METHOD
THE 19TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS STUDY OF INTERFACIAL BEHAVIOR OF CNT/POLYMER COMPOSITE BY CFE METHOD Q. S. Yang*, X. Liu, L. D. Su Department of Engineering Mehanis, Beijing University
More informationFINITE WORD LENGTH EFFECTS IN DSP
FINITE WORD LENGTH EFFECTS IN DSP PREPARED BY GUIDED BY Snehal Gor Dr. Srianth T. ABSTRACT We now that omputers store numbers not with infinite preision but rather in some approximation that an be paed
More informationELECTRODYNAMICS: PHYS 30441
. Relativisti Eletromagnetism. Eletromagneti Field Tensor How do E and B fields transform under a LT? They annot be 4-vetors, but what are they? We again re-write the fields in terms of the salar and vetor
More informationFinal Review. A Puzzle... Special Relativity. Direction of the Force. Moving at the Speed of Light
Final Review A Puzzle... Diretion of the Fore A point harge q is loated a fixed height h above an infinite horizontal onduting plane. Another point harge q is loated a height z (with z > h) above the plane.
More informationPHY 108: Optical Physics. Solution to Midterm Test
PHY 108: Optial Physis Solution to Midterm Test TA: Xun Jia 1 May 14, 2008 1 Email: jiaxun@physis.ula.edu Spring 2008 Physis 108 Xun Jia (May 14, 2008) Problem #1 For a two mirror resonant avity, the resonane
More informationMaximum Entropy and Exponential Families
Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It
More informationWhat s New in ChemSep TM 6.8
What s New in ChemSep TM 6.8 January 2011 (Updated Marh 2011) Harry Kooijman and Ross Taylor In this doument we identify and desribe the most important new features in ChemSep. 1. New: GUI an diretly load
More informationLecture 3 - Lorentz Transformations
Leture - Lorentz Transformations A Puzzle... Example A ruler is positioned perpendiular to a wall. A stik of length L flies by at speed v. It travels in front of the ruler, so that it obsures part of the
More informationEffect of Rotation, Magnetic Field and Initial Stresses on Propagation of Plane Waves in Transversely Isotropic Dissipative Half Space
Applied Mathematis 4 7- http://dx.doi.org/.436/am..48a5 Published Online August (http://www.sirp.org/journal/am) Effet of otation Magneti Field and Initial Stresses on Propagation of Plane Waves in Transversely
More informationA Preliminary Explanation for the Pentaquark P found by LHCb
A Preliminary Explanation for the Pentaquark P found by HCb Mario Everaldo de Souza Departamento de Físia, Universidade Federal de Sergipe, Av. Marehal Rondon, s/n, Rosa Elze, 49100-000 São Cristóvão,
More informationSupporting Information for
Eletroni Supplementary Material (ESI) for Nanosale This journal is The Royal Soiety of Chemistry 013 Supporting Information for Exitation polarization modulation in loalization mirosopy allows to resolve
More informationEE 321 Project Spring 2018
EE 21 Projet Spring 2018 This ourse projet is intended to be an individual effort projet. The student is required to omplete the work individually, without help from anyone else. (The student may, however,
More informationUPPER-TRUNCATED POWER LAW DISTRIBUTIONS
Fratals, Vol. 9, No. (00) 09 World Sientifi Publishing Company UPPER-TRUNCATED POWER LAW DISTRIBUTIONS STEPHEN M. BURROUGHS and SARAH F. TEBBENS College of Marine Siene, University of South Florida, St.
More informationExamination of 4He droplets and droplets containing impurities at zero Kelvin using a density functional approach
University of Tennessee, Knoxville Trae: Tennessee Researh and Creative Exhange Masters Theses Graduate Shool 8-2011 Examination of 4He droplets and droplets ontaining impurities at zero Kelvin using a
More informationThe Gravitational Potential for a Moving Observer, Mercury s Perihelion, Photon Deflection and Time Delay of a Solar Grazing Photon
Albuquerque, NM 0 POCEEDINGS of the NPA 457 The Gravitational Potential for a Moving Observer, Merury s Perihelion, Photon Defletion and Time Delay of a Solar Grazing Photon Curtis E. enshaw Tele-Consultants,
More information9 Geophysics and Radio-Astronomy: VLBI VeryLongBaseInterferometry
9 Geophysis and Radio-Astronomy: VLBI VeryLongBaseInterferometry VLBI is an interferometry tehnique used in radio astronomy, in whih two or more signals, oming from the same astronomial objet, are reeived
More informationINVESTIGATION AND IMPROVEMENT OF CRITICALITY CALCULATIONS IN MCNP5 INVOLVING SHANNON ENTROPY CONVERGENCE
INVESTIGATION AND IMPROVEMENT OF CRITICALITY CALCULATIONS IN MCNP5 INVOLVING SHANNON ENTROPY CONVERGENCE A Thesis Presented to The Aademi Faulty by David Jeffrey Koh In Partial Fulfillment of the Requirements
More information