COGNACのコンセプト \(COarse Grained molecular dynamics program developed by NAgoya Cooperation\)

Size: px
Start display at page:

Download "COGNACのコンセプト \(COarse Grained molecular dynamics program developed by NAgoya Cooperation\)"

Transcription

1 COGNAC (COarse-Grained molecular dynamics program by NAgoya Cooperation) ( ),

2 0 sec -3 msec -6 sec -9 nsec -12 psec -15 fsec GOURMET SUSHI PASTA COGNAC MUFFIN fm pm nm m mm m

3 United atom model (CH 2 ) Gay-Berne potential model Bead-spring model

4

5 Molecular dynamics (MD) Ensembles»NVE» NVT,NPH,NPT (loose-coupling / extended Hamiltonian methods) Langevin dynamics Molecular mechanics (MM) Steepest descent / conjugate gradient methods

6 Bonding 2-body(bond):Harmonic,Morse,FENE,Gaussian, Polynomial,Table 3-body(angle):Theta harmonic,cosine harmonic Theta polynomial,table 4-body(torsion):Cosine polynomial,table Non-bonding pair interaction Lennard-Jones,Gay-Berne,LJ-GB, Table Electrostatic Coulomb interaction(ewald,reaction field) Dipole-dipole interaction (Reaction field)

7 : Gay-Berne - Lennard-Jones hybrid potential C CH 2 H 2 C CH 2 H 2 C CH 3 ncb (4-methyl-4 -cyanobiphenyl) Ellipsoid Sphere Smectic phase(non-polar model) Nematic phase(polar model)

8 SILK (1) SILK COGNAC SILK Python GOURMET SILK

9 SILK (2) name="mol" nummol=10 self.engine.createmolecule(name) for i in range(0, 4): self.engine.addatoms(name, "UA", "UA_PE") for i in range(0, 3): self.engine.addbonds(name, i, i+1, "BOND_PE") for i in range(0, 2): self.engine.addangles(name, i, i+1, i+2, "ANGLE_PE") for i in range(0, 1): self.engine.addtorsions(name, i, i+1, i+2, i+3, "TORSION_PE") for i in range(0, 4): self.engine.addinteractionsites(name, [i], "NB_PE", "PAIR") self.engine.setsystem(name, nummol)

10 SILK (3) name="a20b40a20" nummol=50 key="linear" sequence=[("a",20),("b",40),("a",20)] atomtype={"a":"atom1", "B":"atom2"} bondtype={"a_a":"bond1", "A_B":"bond3", "B_B":"bond2"} interactionsitetype={"a":"sitetype1", "B":"siteType2"} self.engine.makebeadspringpolym(name, nummol, key, sequence, atomtype, bondtype, interactionsitetype)

11 Action SILK gift Action GOUMET SILK Selection of diblock

12

13 COGNAC Random: Amorphous like structures Helix: Helical structures at regular lattice points Crystal: Crystal structures defined by crystal data, i.e. unit lattice, symmetric operation and fractional coordinates Semi-crystalline lamella: Semi-crystalline lamella structures consisting of a crystal phase and an amorphous phase Multi phase structure: Micro/macro phase-separeted structures of block copolymer/polymer blend obtained by SUSHI

14 mol/pdb UDF WebLab ViewerLite (TM) mol GOURMET UDF

15 UDF PDB/car/XYZ GOURMET UDF WebLab ViewerLite (TM) car

16 etc. Lees-Edwards MD»

17 Clay(laponite) - Polymer(PEO) composite clay-polymer Clay

18 20nm 20nm

19

20 Density biased Monte Carlo (DBMC) Density biased potential (DBP) SUSHI Staggered reflective boundary condition (SRBC) Lamella builder

21 ABA triblock copolymer ABA triblock copolymer SUSHI Loop/Bridge

22 ABA triblock copolymer 300% Strain BCC sphere phase εσ

23 A/B εσ ε τ elongation δε δε δε δε

24 6nm elongation

25 COGNAC C++ COGNAC UserBond1, UserAngle1 #include "userbond1.h" double UserBond1::calcforce(const Vector3d& dr, Vector3d& ftmp) { double r,delr,ene,tmp; } r=dr.length(); delr=r-r0; tmp=kconst*delr; ftmp=dr*(tmp/r); ene=0.5*tmp*delr; return ene;

26 DPD Dissipative particle dynamics (DPD) dr dt i dv = vi, dt i = f i f i F C ij = ( C D R F + + ) ij Fij Fij i j aij = 0 ( 1 r ) ( < ) ij rˆ ij rij ( r 1) ij 1, F D ij D = w ( )( ) R R r ( ) ij rˆ ij vij rˆ ij, Fij = w rij ijrˆ ij

27 Action Python molecules/atoms/bonds ABA triblock copolymers A

28 COGNAC Python»»»»»»

29 GOURMET Action GOURMET Action

30 :

31 : -

32 UDF HELP

33 COGNAC UDF unit parameters reduced mass in [amu] reduced energy in [kj/mol] reduced length in [nm].

34 COGNAC: χ MUFFIN SUSHI PASTA COGNAC

35 COGNAC : MD/MM

36 COGNAC JCII

OCTAプロジェクト:物質の多階層シミュレーション

OCTAプロジェクト:物質の多階層シミュレーション SS HPC 2003 2003/10/03 OCTA : www.stat.cse.nagoya-u.ac.jp,, 1 SS HPC 2003 2003/10/03 OCTA Open Computational Tool for Advanced material technology 8 2 SS HPC 2003 2003/10/03 Advanced Material Technology

More information

MUFFIN3

MUFFIN3 MUFFIN - MUltiFarious FIeld simulator for Non-equilibrium system - ( ) MUFFIN WG3 - - JCII, - ( ) - ( ) - ( ) - (JSR) - - MUFFIN sec -3 msec -6 sec GOURMET SUSHI MUFFIN -9 nsec PASTA -1 psec -15 fsec COGNAC

More information

スライド タイトルなし

スライド タイトルなし J-OCTA/OCTA と LAMMPS の連携によるメリット http://www.j-octa.com/ 2016 年 2 月 19 日株式会社 JSOL エンジニアリングビジネス事業部 JSOL について 社員数 1300 人 計算科学分野は 150 人 20 以上のシミュレーション, CAE(Computer Aided Engineering) ソフトウェアミクロからマクロまで 幅広いソリューション

More information

4/15 No.

4/15 No. 4/15 No. 1 4/15 No. 4/15 No. 3 Particle of mass m moving in a potential V(r) V(r) m i ψ t = m ψ(r,t)+v(r)ψ(r,t) ψ(r,t) = ϕ(r)e iωt ψ(r,t) Wave function steady state m ϕ(r)+v(r)ϕ(r) = εϕ(r) Eigenvalue problem

More information

卒業研究報告 題 目 Hamiltonian 指導教員 山本哲也教授 報告者 汐月康則 平成 14 年 2 月 5 日 1

卒業研究報告 題 目 Hamiltonian 指導教員 山本哲也教授 報告者 汐月康則 平成 14 年 2 月 5 日 1 卒業研究報告 題 目 Hamiltonian 指導教員 山本哲也教授 報告者 汐月康則 平成 4 年 月 5 日 .....4.....4......6.. 6.. 6....4. 8.5. 9.6....7... 3..... 3.... 3.... 3.3...4 3.4...5 3.5...5 3.5....6 3.5.... 3.5...... 3.5...... 3 3.5.3..4 3.5.4..5

More information

02-量子力学の復習

02-量子力学の復習 4/17 No. 1 4/17 No. 2 4/17 No. 3 Particle of mass m moving in a potential V(r) V(r) m i ψ t = 2 2m 2 ψ(r,t)+v(r)ψ(r,t) ψ(r,t) Wave function ψ(r,t) = ϕ(r)e iωt steady state 2 2m 2 ϕ(r)+v(r)ϕ(r) = εϕ(r)

More information

1 1.1,,,.. (, ),..,. (Fig. 1.1). Macro theory (e.g. Continuum mechanics) Consideration under the simple concept (e.g. ionic radius, bond valence) Stru

1 1.1,,,.. (, ),..,. (Fig. 1.1). Macro theory (e.g. Continuum mechanics) Consideration under the simple concept (e.g. ionic radius, bond valence) Stru 1. 1-1. 1-. 1-3.. MD -1. -. -3. MD 1 1 1.1,,,.. (, ),..,. (Fig. 1.1). Macro theory (e.g. Continuum mechanics) Consideration under the simple concept (e.g. ionic radius, bond valence) Structural relaxation

More information

19 σ = P/A o σ B Maximum tensile strength σ % 0.2% proof stress σ EL Elastic limit Work hardening coefficient failure necking σ PL Proportional

19 σ = P/A o σ B Maximum tensile strength σ % 0.2% proof stress σ EL Elastic limit Work hardening coefficient failure necking σ PL Proportional 19 σ = P/A o σ B Maximum tensile strength σ 0. 0.% 0.% proof stress σ EL Elastic limit Work hardening coefficient failure necking σ PL Proportional limit ε p = 0.% ε e = σ 0. /E plastic strain ε = ε e

More information

A 99% MS-Free Presentation

A 99% MS-Free Presentation A 99% MS-Free Presentation 2 Galactic Dynamics (Binney & Tremaine 1987, 2008) Dynamics of Galaxies (Bertin 2000) Dynamical Evolution of Globular Clusters (Spitzer 1987) The Gravitational Million-Body Problem

More information

80 4 r ˆρ i (r, t) δ(r x i (t)) (4.1) x i (t) ρ i ˆρ i t = 0 i r 0 t(> 0) j r 0 + r < δ(r 0 x i (0))δ(r 0 + r x j (t)) > (4.2) r r 0 G i j (r, t) dr 0

80 4 r ˆρ i (r, t) δ(r x i (t)) (4.1) x i (t) ρ i ˆρ i t = 0 i r 0 t(> 0) j r 0 + r < δ(r 0 x i (0))δ(r 0 + r x j (t)) > (4.2) r r 0 G i j (r, t) dr 0 79 4 4.1 4.1.1 x i (t) x j (t) O O r 0 + r r r 0 x i (0) r 0 x i (0) 4.1 L. van. Hove 1954 space-time correlation function V N 4.1 ρ 0 = N/V i t 80 4 r ˆρ i (r, t) δ(r x i (t)) (4.1) x i (t) ρ i ˆρ i t

More information

3 3.1 R r r + R R r Rr [ ] ˆn(r) = ˆn(r + R) (3.1) R R = r ˆn(r) = ˆn(0) r 0 R = r C nn (r, r ) = C nn (r + R, r + R) = C nn (r r, 0) (3.2) ( 2.2 ) C

3 3.1 R r r + R R r Rr [ ] ˆn(r) = ˆn(r + R) (3.1) R R = r ˆn(r) = ˆn(0) r 0 R = r C nn (r, r ) = C nn (r + R, r + R) = C nn (r r, 0) (3.2) ( 2.2 ) C 3 3.1 R r r + R R r Rr [ ] ˆn(r) = ˆn(r + R) (3.1) R R = r ˆn(r) = ˆn(0) r 0 R = r C nn (r, r ) = C nn (r + R, r + R) = C nn (r r, 0) (3.2) ( 2.2 ) C nn (r r ) = C nn (R(r r )) [2 ] 2 g(r, r ) ˆn(r) ˆn(r

More information

2012専門分科会_new_4.pptx

2012専門分科会_new_4.pptx d dt L L = 0 q i q i d dt L L = 0 r i i r i r r + Δr Δr δl = 0 dl dt = d dt i L L q i q i + q i i q i = q d L L i + q i i dt q i i q i = i L L q i L = 0, H = q q i L = E i q i i d dt L q q i i L = L(q

More information

untitled

untitled SPring-8 RFgun JASRI/SPring-8 6..7 Contents.. 3.. 5. 6. 7. 8. . 3 cavity γ E A = er 3 πε γ vb r B = v E c r c A B A ( ) F = e E + v B A A A A B dp e( v B+ E) = = m d dt dt ( γ v) dv e ( ) dt v B E v E

More information

高知工科大学電子 光システム工学科

高知工科大学電子 光システム工学科 卒業研究報告 題 目 量子力学に基づいた水素分子の分子軌道法的取り扱いと Hamiltonian 近似法 指導教員 山本哲也 報告者 山中昭徳 平成 14 年 月 5 日 高知工科大学電子 光システム工学科. 3. 4.1 4. 4.3 4.5 6.6 8.7 10.8 11.9 1.10 1 3. 13 3.113 3. 13 3.3 13 3.4 14 3.5 15 3.6 15 3.7 17

More information

CMP Technical Report No. 4 Department of Computational Nanomaterials Design ISIR, Osaka University 2 2................................. 2.2......................... 2 3 3 3................................

More information

Molecule tomic rbital bridied tomic rbital Valence Shell Electron Pair Repulsion Rule Molecular rbital 2 1+ + 1+ 1+ 1+ 2 9+ + 9+ 9+ 9+ 2 1+ 1+ 1s 1s 2 9+ 9+ 2p 2p 9+ () 2 (2p ) 2 (2p ) 2 (2p ) 1 Energ

More information

D v D F v/d F v D F η v D (3.2) (a) F=0 (b) v=const. D F v Newtonian fluid σ ė σ = ηė (2.2) ė kl σ ij = D ijkl ė kl D ijkl (2.14) ė ij (3.3) µ η visco

D v D F v/d F v D F η v D (3.2) (a) F=0 (b) v=const. D F v Newtonian fluid σ ė σ = ηė (2.2) ė kl σ ij = D ijkl ė kl D ijkl (2.14) ė ij (3.3) µ η visco post glacial rebound 3.1 Viscosity and Newtonian fluid f i = kx i σ ij e kl ideal fluid (1.9) irreversible process e ij u k strain rate tensor (3.1) v i u i / t e ij v F 23 D v D F v/d F v D F η v D (3.2)

More information

株式会社ローソン 第26期中間事業報告書

株式会社ローソン 第26期中間事業報告書 1 2 3 4 5 797 429 316 1,881 7,583 204 649 334 1,950 564 459 73.10 4,836 3,535 79.05 5,113 4,042 79.96 5,683 4,544 82.41 6,252 5,152 83.47 6,649 5,550 85.75 7,016 6,016 88.45 7,378 6,526 88.95 7,583 6,745

More information

Introduction 2 / 43

Introduction 2 / 43 Batalin-Vilkoviski ( ) 2016 2 22 at SFT16 based on arxiv:1511.04187 BV Analysis of Tachyon Fluctuation around Multi-brane Solutions in Cubic String Field Theory 1 / 43 Introduction 2 / 43 in Cubic open

More information

PDF

PDF 1 1 1 1-1 1 1-9 1-3 1-1 13-17 -3 6-4 6 3 3-1 35 3-37 3-3 38 4 4-1 39 4- Fe C TEM 41 4-3 C TEM 44 4-4 Fe TEM 46 4-5 5 4-6 5 5 51 6 5 1 1-1 1991 1,1 multiwall nanotube 1993 singlewall nanotube ( 1,) sp 7.4eV

More information

卒業論文

卒業論文 Cu-Ru ...7 1.1....7 1.2....7 1.3....8...9 2.1....9 2.1.1....9 2.1.2....10 2.2.... 11 2.2.1.... 11 2.2.2. ()...12 2.2.3. ()...13 2.2....15...18 3.1. Cu z ...18 3.1.1....18 3.1.2....19 3.1.3....29 3.2.

More information

磁性物理学 - 遷移金属化合物磁性のスピンゆらぎ理論

磁性物理学 - 遷移金属化合物磁性のスピンゆらぎ理論 email: takahash@sci.u-hyogo.ac.jp April 30, 2009 Outline 1. 2. 3. 4. 5. 6. 2 / 260 Today s Lecture: Itinerant Magnetism 60 / 260 Multiplets of Single Atom System HC HSO : L = i l i, S = i s i, J = L +

More information

iBookBob:Users:bob:Documents:CurrentData:flMŠÍ…e…L…X…g:Statistics.dvi

iBookBob:Users:bob:Documents:CurrentData:flMŠÍ…e…L…X…g:Statistics.dvi 4 4 9............................................... 3.3......................... 4.4................. 5.5............................ 7 9..................... 9.............................3................................4..........................5.............................6...........................

More information

IPSJ SIG Technical Report Vol.2014-ARC-213 No.24 Vol.2014-HPC-147 No /12/10 GPU 1,a) 1,b) 1,c) 1,d) GPU GPU Structure Of Array Array Of

IPSJ SIG Technical Report Vol.2014-ARC-213 No.24 Vol.2014-HPC-147 No /12/10 GPU 1,a) 1,b) 1,c) 1,d) GPU GPU Structure Of Array Array Of GPU 1,a) 1,b) 1,c) 1,d) GPU 1 GPU Structure Of Array Array Of Structure 1. MPS(Moving Particle Semi-Implicit) [1] SPH(Smoothed Particle Hydrodynamics) [] DEM(Distinct Element Method)[] [] 1 Tokyo Institute

More information

SAXS Table 1 DSC POM SAXSSAXS PF BL-10C BL-15A Fig. 2 LC12 DSC SAXS 138 C T iso T iso SAXS q=1.4 nm -1 q=(4π/λ)sin(θ/2), λ:, θ: Fig. 3 LC12 T iso Figu

SAXS Table 1 DSC POM SAXSSAXS PF BL-10C BL-15A Fig. 2 LC12 DSC SAXS 138 C T iso T iso SAXS q=1.4 nm -1 q=(4π/λ)sin(θ/2), λ:, θ: Fig. 3 LC12 T iso Figu 1 1 1 1,2 1,2 1 2 Correlation between Microphase Separation and Liquid Crystallization in Structure Formation of Liquid Crystalline Block Copolymers Shin-ichi TANIGUCHI 1, Hiroki TAKESHITA 1, Masamitsu

More information

A2, Vol. 69, No. 2 Vol. 16, I_237-I_246, Analytical Investigation of Shear Force Distribution of Perfobond Strip with Plural Perforations * ** *

A2, Vol. 69, No. 2 Vol. 16, I_237-I_246, Analytical Investigation of Shear Force Distribution of Perfobond Strip with Plural Perforations * ** * A2, Vol. 69, No. 2 Vol. 16, I_237-I_246, 213. Analytical Investigation of Shear Force Distribution of Perfobond Strip with Plural Perforations * ** *** **** ***** Noriyuki KUBO, Takeshi SAKAI, Shinji OHGUCHI,

More information

Outline I. Introduction: II. Pr 2 Ir 2 O 7 Like-charge attraction III.

Outline I. Introduction: II. Pr 2 Ir 2 O 7 Like-charge attraction III. Masafumi Udagawa Dept. of Physics, Gakushuin University Mar. 8, 16 @ in Gakushuin University Reference M. U., L. D. C. Jaubert, C. Castelnovo and R. Moessner, arxiv:1603.02872 Outline I. Introduction:

More information

ADM-Hamiltonian Cheeger-Gromov 3. Penrose

ADM-Hamiltonian Cheeger-Gromov 3. Penrose ADM-Hamiltonian 1. 2. Cheeger-Gromov 3. Penrose 0. ADM-Hamiltonian (M 4, h) Einstein-Hilbert M 4 R h hdx L h = R h h δl h = 0 (Ric h ) αβ 1 2 R hg αβ = 0 (Σ 3, g ij ) (M 4, h ij ) g ij, k ij Σ π ij = g(k

More information

ssp2_fixed.dvi

ssp2_fixed.dvi 13 12 30 2 1 3 1.1... 3 1.2... 4 1.3 Bravais... 4 1.4 Miller... 4 2 X 5 2.1 Bragg... 5 2.2... 5 2.3... 7 3 Brillouin 13 3.1... 13 3.2 Brillouin... 13 3.3 Brillouin... 14 3.4 Bloch... 16 3.5 Bloch... 17

More information

.I.v e pmd

.I.v e pmd Structural Design for Curved Panels by Laminated Composite Materials (Identification of Lamination Parameters Using Modal Testing Method ) Tetsuya NARISAWA, Shohei IWATA Abstract - Using a modal testing

More information

Yuzo Nakamura, Kagoshima Univ., Dept Mech Engr. perfect crystal imperfect crystal point defect vacancy self-interstitial atom substitutional impurity

Yuzo Nakamura, Kagoshima Univ., Dept Mech Engr. perfect crystal imperfect crystal point defect vacancy self-interstitial atom substitutional impurity perfect crystal imperfect crystal point defect vacancy self-interstitial atom substitutional impurity atom interstitial impurity atom line defect dislocation planar defect surface grain boundary interface

More information

1 2 2 (Dielecrics) Maxwell ( ) D H

1 2 2 (Dielecrics) Maxwell ( ) D H 2003.02.13 1 2 2 (Dielecrics) 4 2.1... 4 2.2... 5 2.3... 6 2.4... 6 3 Maxwell ( ) 9 3.1... 9 3.2 D H... 11 3.3... 13 4 14 4.1... 14 4.2... 14 4.3... 17 4.4... 19 5 22 6 THz 24 6.1... 24 6.2... 25 7 26

More information

Sample function Re random process Flutter, Galloping, etc. ensemble (mean value) N 1 µ = lim xk( t1) N k = 1 N autocorrelation function N 1 R( t1, t1

Sample function Re random process Flutter, Galloping, etc. ensemble (mean value) N 1 µ = lim xk( t1) N k = 1 N autocorrelation function N 1 R( t1, t1 Sample function Re random process Flutter, Galloping, etc. ensemble (mean value) µ = lim xk( k = autocorrelation function R( t, t + τ) = lim ( ) ( + τ) xk t xk t k = V p o o R p o, o V S M R realization

More information

42 1 Fig. 2. Li 2 B 4 O 7 crystals with 3inches and 4inches in diameter. Fig. 4. Transmission curve of Li 2 B 4 O 7 crystal. Fig. 5. Refractive index

42 1 Fig. 2. Li 2 B 4 O 7 crystals with 3inches and 4inches in diameter. Fig. 4. Transmission curve of Li 2 B 4 O 7 crystal. Fig. 5. Refractive index MEMOIRS OF SHONAN INSTITUTE OF TECHNOLOGY Vol. 42, No. 1, 2008 Li 2 B 4 O 7 (LBO) *, ** * ** ** Optical Scatterer and Crystal Growth Technology of LBO Single Crystal For Development with Optical Application

More information

¼§À�ÍýÏÀ – Ê×ÎòÅŻҼ§À�¤È¥¹¥Ô¥ó¤æ¤é¤®

¼§À�ÍýÏÀ – Ê×ÎòÅŻҼ§À�¤È¥¹¥Ô¥ó¤æ¤é¤® email: takahash@sci.u-hyogo.ac.jp Spring semester, 2012 Outline 1. 2 / 26 Introduction : (d ) : 4f 1970 ZrZn 2, MnSi, Ni 3 Al, Sc 3 In Stoner-Wohlfarth Moriya-Kawabata (1973) 3 / 26 Properties of Weak

More information

untitled

untitled 173 1 2 3 1 PEFC DMFC 1 CMOS beyond CMOS More than Moore IV III-V 2 MRAM FeRAM RAM ReRAM 3 4 EL FED) 1 - - DDS 4 1 2 3 4 5 5 6 2005 2006 2007 2008 WG 2010 2030 2020 2004 2011 20072011 ISO/TC229IEC/TC113

More information

global global mass region (matter ) & (I) M3Y semi-microscopic int. Ref.: H. N., P. R. C68, ( 03) N. P. A722, 117c ( 03) Proc. of NENS03 (to be

global global mass region (matter ) & (I) M3Y semi-microscopic int. Ref.: H. N., P. R. C68, ( 03) N. P. A722, 117c ( 03) Proc. of NENS03 (to be Gogny hard core spin-isospin property @ RCNP (Mar. 22 24, 2004) Collaborator: M. Sato (Chiba U, ) ( ) global global mass region (matter ) & (I) M3Y semi-microscopic int. Ref.: H. N., P. R. C68, 014316

More information

研究室ガイダンス(H28)福山研.pdf

研究室ガイダンス(H28)福山研.pdf 1 2 3 4 5 4 He M. Roger et al., JLTP 112, 45 (1998) A.F. Andreev and I.M. Lifshitz, Sov. Phys. JETP 29, 1107 (1969) Born in 2004 (hcp 4 He) E. Kim and M.H.W. Chan, Nature 427, 225 (2004); Science 305,

More information

i

i 009 I 1 8 5 i 0 1 0.1..................................... 1 0.................................................. 1 0.3................................. 0.4........................................... 3

More information

·«¤ê¤³¤ß·²¤È¥ß¥ì¥Ë¥¢¥àÌäÂê

·«¤ê¤³¤ß·²¤È¥ß¥ì¥Ë¥¢¥àÌäÂê .. 1 10-11 Nov., 2016 1 email:keiichi.r.ito@gmail.com, ito@kurims.kyoto-u.ac.jp ( ) 10-11 Nov., 2016 1 / 45 Clay Institute.1 Construction of 4D YM Field Theory (Jaffe, Witten) Jaffe, Balaban (1980).2 Solution

More information

( ) ) AGD 2) 7) 1

( ) ) AGD 2) 7) 1 ( 9 5 6 ) ) AGD ) 7) S. ψ (r, t) ψ(r, t) (r, t) Ĥ ψ(r, t) = e iĥt/ħ ψ(r, )e iĥt/ħ ˆn(r, t) = ψ (r, t)ψ(r, t) () : ψ(r, t)ψ (r, t) ψ (r, t)ψ(r, t) = δ(r r ) () ψ(r, t)ψ(r, t) ψ(r, t)ψ(r, t) = (3) ψ (r,

More information

1 1-49 11 5 7064 ...4 1.1...5 1....7 1.3...8...9.1...10.1.1 Lennad-Jones...10.1....1....13.3...14.4...15.4.1 L-J...15.4. Hamonic...16.4.3 Multiple time step...16.5...17.6...19.6.1...19.7...0.7.1...0.7.

More information

BIT -2-

BIT -2- 2004.3.31 10 11 12-1- BIT -2- -3-256 258 932 524 585 -4- -5- A B A B AB A B A B C AB A B AB AB AB AB -6- -7- A B -8- -9- -10- mm -11- fax -12- -13- -14- -15- s58.10.1 1255 4.2 30.10-16- -17- -18- -19-6.12.10

More information

nm (T = K, p = kP a (1atm( )), 1bar = 10 5 P a = atm) 1 ( ) m / m

nm (T = K, p = kP a (1atm( )), 1bar = 10 5 P a = atm) 1 ( ) m / m .1 1nm (T = 73.15K, p = 101.35kP a (1atm( )), 1bar = 10 5 P a = 0.9863atm) 1 ( ).413968 10 3 m 3 1 37. 1/3 3.34.414 10 3 m 3 6.0 10 3 = 3.7 (109 ) 3 (nm) 3 10 6 = 3.7 10 1 (nm) 3 = (3.34nm) 3 ( P = nrt,

More information

untitled

untitled 1 2 3 4 5 130mm 32mm UV-irradiation UV-cationic cure UV-cationic cure UV-cationic cure Thermal cationic Reaction heat cure Thermal cationic Cation Reaction heat cure Cation (a) UV-curing of

More information

The Phase Behavior of Monooleoylglycerol-Water Systems Mivoshi Oil & Fat Co.. Ltd. Faculty of Science and Technology, Science University of Tokyo Inst

The Phase Behavior of Monooleoylglycerol-Water Systems Mivoshi Oil & Fat Co.. Ltd. Faculty of Science and Technology, Science University of Tokyo Inst The Phase Behavior of Monooleoylglycerol-Water Systems Mivoshi Oil & Fat Co.. Ltd. Faculty of Science and Technology, Science University of Tokyo Institute of Colloid and Interface Science, Science University

More information

IPSJ SIG Technical Report NetMAS NetMAS NetMAS One-dimensional Pedestrian Model for Fast Evacuation Simulator Shunsuke Soeda, 1 Tomohisa Yam

IPSJ SIG Technical Report NetMAS NetMAS NetMAS One-dimensional Pedestrian Model for Fast Evacuation Simulator Shunsuke Soeda, 1 Tomohisa Yam 1 1 1 1 1 NetMAS NetMAS NetMAS One-dimensional Model for Fast Evacuation Simulator Shunsuke Soeda, 1 Tomohisa Yamashita, 1 Masaki Onishi, 1 Ikushi Yoda 1 and Itsuki Noda 1 We propose the one-dimentional

More information

CHARACTERISTICS OF LOVE WAVE GENERATED AROUND A DIPPING BASEMENT By Susumu NAKAMURA, Iwao SUETOMI, Shinichi AKIYAMA and Nozomu YOSHIDA Source mechanis

CHARACTERISTICS OF LOVE WAVE GENERATED AROUND A DIPPING BASEMENT By Susumu NAKAMURA, Iwao SUETOMI, Shinichi AKIYAMA and Nozomu YOSHIDA Source mechanis CHARACTERISTICS OF LOVE WAVE GENERATED AROUND A DIPPING BASEMENT By Susumu NAKAMURA, Iwao SUETOMI, Shinichi AKIYAMA and Nozomu YOSHIDA Source mechanism and characteristics of the horizontally propagating

More information

Nosé Hoover 1.2 ( 1) (a) (b) 1:

Nosé Hoover 1.2 ( 1) (a) (b) 1: 1 watanabe@cc.u-tokyo.ac.jp 1 1.1 Nosé Hoover 1. ( 1) (a) (b) 1: T ( f(p x, p y, p z ) exp p x + p y + p ) z (1) mk B T p x p y p = = z = 1 m m m k BT () k B T = 1.3 0.04 0.03 0.0 0.01 0-5 -4-3 - -1 0

More information

第62巻 第1号 平成24年4月/石こうを用いた木材ペレット

第62巻 第1号 平成24年4月/石こうを用いた木材ペレット Bulletin of Japan Association for Fire Science and Engineering Vol. 62. No. 1 (2012) Development of Two-Dimensional Simple Simulation Model and Evaluation of Discharge Ability for Water Discharge of Firefighting

More information

修士論文 物性研究 電子版 Vol. 5, No. 2, (2016 年 5 月号 ) 27-2 F14A001B

修士論文 物性研究 電子版 Vol. 5, No. 2, (2016 年 5 月号 ) 27-2 F14A001B 27-2 F14A001B (Molecular Dynamics: MD) Verlet (Verlet Neighbor List: VNL ) (Cell Linked List: CLL ) (VNL-CLL ) VNL MD 10 VNL ( ) CLL (2 8 26 ) CLL ( ) VNL O(N 2 ) CLL O(N 2 ) O(N) VNL-CLL VNL-CLL Verlet

More information

LLG-R8.Nisus.pdf

LLG-R8.Nisus.pdf d M d t = γ M H + α M d M d t M γ [ 1/ ( Oe sec) ] α γ γ = gµ B h g g µ B h / π γ g = γ = 1.76 10 [ 7 1/ ( Oe sec) ] α α = λ γ λ λ λ α γ α α H α = γ H ω ω H α α H K K H K / M 1 1 > 0 α 1 M > 0 γ α γ =

More information

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,.

ii 3.,. 4. F. (), ,,. 8.,. 1. (75%) (25%) =7 20, =7 21 (. ). 1.,, (). 3.,. 1. ().,.,.,.,.,. () (12 )., (), 0. 2., 1., 0,. 24(2012) (1 C106) 4 11 (2 C206) 4 12 http://www.math.is.tohoku.ac.jp/~obata,.,,,.. 1. 2. 3. 4. 5. 6. 7.,,. 1., 2007 (). 2. P. G. Hoel, 1995. 3... 1... 2.,,. ii 3.,. 4. F. (),.. 5... 6.. 7.,,. 8.,. 1. (75%)

More information

chap1_MDpotentials.ppt

chap1_MDpotentials.ppt simplest Morse : simplest (1 Well-chosen functional form is more useful than elaborate fitting strategies!! Phys. Rev. 34, 57 (1929 ( 2 E ij = D e 1" exp("#(r ij " r e 2 r ij = r i " r j r ij =r e E ij

More information

9 1. (Ti:Al 2 O 3 ) (DCM) (Cr:Al 2 O 3 ) (Cr:BeAl 2 O 4 ) Ĥ0 ψ n (r) ω n Schrödinger Ĥ 0 ψ n (r) = ω n ψ n (r), (1) ω i ψ (r, t) = [Ĥ0 + Ĥint (

9 1. (Ti:Al 2 O 3 ) (DCM) (Cr:Al 2 O 3 ) (Cr:BeAl 2 O 4 ) Ĥ0 ψ n (r) ω n Schrödinger Ĥ 0 ψ n (r) = ω n ψ n (r), (1) ω i ψ (r, t) = [Ĥ0 + Ĥint ( 9 1. (Ti:Al 2 O 3 ) (DCM) (Cr:Al 2 O 3 ) (Cr:BeAl 2 O 4 ) 2. 2.1 Ĥ ψ n (r) ω n Schrödinger Ĥ ψ n (r) = ω n ψ n (r), (1) ω i ψ (r, t) = [Ĥ + Ĥint (t)] ψ (r, t), (2) Ĥ int (t) = eˆxe cos ωt ˆdE cos ωt, (3)

More information

positron 1930 Dirac 1933 Anderson m 22Na(hl=2.6years), 58Co(hl=71days), 64Cu(hl=12hour) 68Ge(hl=288days) MeV : thermalization m psec 100

positron 1930 Dirac 1933 Anderson m 22Na(hl=2.6years), 58Co(hl=71days), 64Cu(hl=12hour) 68Ge(hl=288days) MeV : thermalization m psec 100 positron 1930 Dirac 1933 Anderson m 22Na(hl=2.6years), 58Co(hl=71days), 64Cu(hl=12hour) 68Ge(hl=288days) 0.5 1.5MeV : thermalization 10 100 m psec 100psec nsec E total = 2mc 2 + E e + + E e Ee+ Ee-c mc

More information

* 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) *

* 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) * * 1 1 (i) (ii) Brückner-Hartree-Fock (iii) (HF, BCS, HFB) (iv) (TDHF,TDHFB) (RPA) (QRPA) (v) (vi) *1 2004 1 1 ( ) ( ) 1.1 140 MeV 1.2 ( ) ( ) 1.3 2.6 10 8 s 7.6 10 17 s? Λ 2.5 10 10 s 6 10 24 s 1.4 ( m

More information

2 (March 13, 2010) N Λ a = i,j=1 x i ( d (a) i,j x j ), Λ h = N i,j=1 x i ( d (h) i,j x j ) B a B h B a = N i,j=1 ν i d (a) i,j, B h = x j N i,j=1 ν i

2 (March 13, 2010) N Λ a = i,j=1 x i ( d (a) i,j x j ), Λ h = N i,j=1 x i ( d (h) i,j x j ) B a B h B a = N i,j=1 ν i d (a) i,j, B h = x j N i,j=1 ν i 1. A. M. Turing [18] 60 Turing A. Gierer H. Meinhardt [1] : (GM) ) a t = D a a xx µa + ρ (c a2 h + ρ 0 (0 < x < l, t > 0) h t = D h h xx νh + c ρ a 2 (0 < x < l, t > 0) a x = h x = 0 (x = 0, l) a = a(x,

More information

V(x) m e V 0 cos x π x π V(x) = x < π, x > π V 0 (i) x = 0 (V(x) V 0 (1 x 2 /2)) n n d 2 f dξ 2ξ d f 2 dξ + 2n f = 0 H n (ξ) (ii) H

V(x) m e V 0 cos x π x π V(x) = x < π, x > π V 0 (i) x = 0 (V(x) V 0 (1 x 2 /2)) n n d 2 f dξ 2ξ d f 2 dξ + 2n f = 0 H n (ξ) (ii) H 199 1 1 199 1 1. Vx) m e V cos x π x π Vx) = x < π, x > π V i) x = Vx) V 1 x /)) n n d f dξ ξ d f dξ + n f = H n ξ) ii) H n ξ) = 1) n expξ ) dn dξ n exp ξ )) H n ξ)h m ξ) exp ξ )dξ = π n n!δ n,m x = Vx)

More information

力学的性質

力学的性質 Materials Science And Engineering, An Introduction: by William D. Callister, Jr., John Wiley & Sons, Inc. Mechanical Metallurgy, G.E.Dieter, McGraw Hill, 1987 Fundamentals of Metal Forming, Robert H. Wagoner,

More information

EGunGPU

EGunGPU Super Computing in Accelerator simulations - Electron Gun simulation using GPGPU - K. Ohmi, KEK-Accel Accelerator Physics seminar 2009.11.19 Super computers in KEK HITACHI SR11000 POWER5 16 24GB 16 134GFlops,

More information

High Melt Strength Polyolefin Elastomers (HMS POEs) and Olefin Block Copolymers (OBCs) Michio ONO* and Takahiko OHMURA(Dow Chemical Japan Ltd, Dow Japan Development Center, 8-1, Ukishimacho, Kawasaki-ku,

More information

CP-PACS CP-PACS CP-PACS : 2048PU+128IOU 614GFLOPS peak 128GByte memory 1058GByte disk 1992 1996 SR2201 : 1996 8 9 CP-PACS Top 500 List ranking No. 1 November 1996 Linpack 368.2Gflops No. 24 Novermber 1999

More information

96 7 1m =2 10 7 N 1A 7.1 7.2 a C (1) I (2) A C I A A a A a A A a C C C 7.2: C A C A = = µ 0 2π (1) A C 7.2 AC C A 3 3 µ0 I 2 = 2πa. (2) A C C 7.2 A A

96 7 1m =2 10 7 N 1A 7.1 7.2 a C (1) I (2) A C I A A a A a A A a C C C 7.2: C A C A = = µ 0 2π (1) A C 7.2 AC C A 3 3 µ0 I 2 = 2πa. (2) A C C 7.2 A A 7 Lorentz 7.1 Ampère I 1 I 2 I 2 I 1 L I 1 I 2 21 12 L r 21 = 12 = µ 0 2π I 1 I 2 r L. (7.1) 7.1 µ 0 =4π 10 7 N A 2 (7.2) magnetic permiability I 1 I 2 I 1 I 2 12 21 12 21 7.1: 1m 95 96 7 1m =2 10 7 N

More information

1 2 3 4 5 6 0.4% 58.4% 41.2% 10 65 69 12.0% 9 60 64 13.4% 11 70 12.6% 8 55 59 8.6% 0.1% 1 20 24 3.1% 7 50 54 9.3% 2 25 29 6.0% 3 30 34 7.6% 6 45 49 9.7% 4 35 39 8.5% 5 40 44 9.1% 11 70 11.2% 10 65 69 11.0%

More information

42 3 u = (37) MeV/c 2 (3.4) [1] u amu m p m n [1] m H [2] m p = (4) MeV/c 2 = (13) u m n = (4) MeV/c 2 =

42 3 u = (37) MeV/c 2 (3.4) [1] u amu m p m n [1] m H [2] m p = (4) MeV/c 2 = (13) u m n = (4) MeV/c 2 = 3 3.1 3.1.1 kg m s J = kg m 2 s 2 MeV MeV [1] 1MeV=1 6 ev = 1.62 176 462 (63) 1 13 J (3.1) [1] 1MeV/c 2 =1.782 661 731 (7) 1 3 kg (3.2) c =1 MeV (atomic mass unit) 12 C u = 1 12 M(12 C) (3.3) 41 42 3 u

More information

9 8 7 (x-1.0)*(x-1.0) *(x-1.0) (a) f(a) (b) f(a) Figure 1: f(a) a =1.0 (1) a 1.0 f(1.0)

9 8 7 (x-1.0)*(x-1.0) *(x-1.0) (a) f(a) (b) f(a) Figure 1: f(a) a =1.0 (1) a 1.0 f(1.0) E-mail: takio-kurita@aist.go.jp 1 ( ) CPU ( ) 2 1. a f(a) =(a 1.0) 2 (1) a ( ) 1(a) f(a) a (1) a f(a) a =2(a 1.0) (2) 2 0 a f(a) a =2(a 1.0) = 0 (3) 1 9 8 7 (x-1.0)*(x-1.0) 6 4 2.0*(x-1.0) 6 2 5 4 0 3-2

More information

2 Chapter 4 (f4a). 2. (f4cone) ( θ) () g M. 2. (f4b) T M L P a θ (f4eki) ρ H A a g. v ( ) 2. H(t) ( )

2 Chapter 4 (f4a). 2. (f4cone) ( θ) () g M. 2. (f4b) T M L P a θ (f4eki) ρ H A a g. v ( ) 2. H(t) ( ) http://astr-www.kj.yamagata-u.ac.jp/~shibata f4a f4b 2 f4cone f4eki f4end 4 f5meanfp f6coin () f6a f7a f7b f7d f8a f8b f9a f9b f9c f9kep f0a f0bt version feqmo fvec4 fvec fvec6 fvec2 fvec3 f3a (-D) f3b

More information

003村江.indd

003村江.indd *1 Study on Room ressure Control at Cleanroom art 1 Experiments on Room ressure Fluctuation with Door Operation and Local Ventilation Operation Yukitada MURAE *1 Tamio IWAMURA *2 Hiroyuki NAGAI *3 Shigeru

More information

1

1 GL (a) (b) Ph l P N P h l l Ph Ph Ph Ph l l l l P Ph l P N h l P l .9 αl B βlt D E. 5.5 L r..8 e g s e,e l l W l s l g W W s g l l W W e s g e s g r e l ( s ) l ( l s ) r e l ( s ) l ( l s ) e R e r

More information

HREM Manual37JFAQ

HREM Manual37JFAQ xhrem TM (WinHREM TM /MacHREM TM ) V3.7 !!!!! xhrem Userʼs Guide 2 !!!!!!!!!! Support/Update Email: support@hremresearch.com WEB: www.hremresearch.com xhrem Userʼs Guide 3 ! xhrem Userʼs Guide 4 " " "

More information

JFE.dvi

JFE.dvi ,, Department of Civil Engineering, Chuo University Kasuga 1-13-27, Bunkyo-ku, Tokyo 112 8551, JAPAN E-mail : atsu1005@kc.chuo-u.ac.jp E-mail : kawa@civil.chuo-u.ac.jp SATO KOGYO CO., LTD. 12-20, Nihonbashi-Honcho

More information

5 H Boltzmann Einstein Brown 5.1 Onsager [ ] Tr Tr Tr = dγ (5.1) A(p, q) Â 0 = Tr Âe βĥ0 Tr e βĥ0 = dγ e βh 0(p,q) A(p, q) dγ e βh 0(p,q) (5.2) e βĥ0

5 H Boltzmann Einstein Brown 5.1 Onsager [ ] Tr Tr Tr = dγ (5.1) A(p, q) Â 0 = Tr Âe βĥ0 Tr e βĥ0 = dγ e βh 0(p,q) A(p, q) dγ e βh 0(p,q) (5.2) e βĥ0 5 H Boltzmann Einstein Brown 5.1 Onsager [ ] Tr Tr Tr = dγ (5.1) A(p, q) Â = Tr Âe βĥ Tr e βĥ = dγ e βh (p,q) A(p, q) dγ e βh (p,q) (5.2) e βĥ A(p, q) p q Â(t) = Tr Â(t)e βĥ Tr e βĥ = dγ() e βĥ(p(),q())

More information

1-1 (1) (3) 16.6 6.6 60.9 82.2 32.1 1980 199 1-2 1-3 (2) (2001) 2001 2-1 P A e U ij A ij = e Uij R + e U ij A + e U ij C U ija = a 0 +a 1 c ija +a 2 T ija +a 3 AT ija + Air i Rail Car

More information

t χ 2 F Q t χ 2 F 1 2 µ, σ 2 N(µ, σ 2 ) f(x µ, σ 2 ) = 1 ( exp (x ) µ)2 2πσ 2 2σ 2 0, N(0, 1) (100 α) z(α) t χ 2 *1 2.1 t (i)x N(µ, σ 2 ) x µ σ N(0, 1

t χ 2 F Q t χ 2 F 1 2 µ, σ 2 N(µ, σ 2 ) f(x µ, σ 2 ) = 1 ( exp (x ) µ)2 2πσ 2 2σ 2 0, N(0, 1) (100 α) z(α) t χ 2 *1 2.1 t (i)x N(µ, σ 2 ) x µ σ N(0, 1 t χ F Q t χ F µ, σ N(µ, σ ) f(x µ, σ ) = ( exp (x ) µ) πσ σ 0, N(0, ) (00 α) z(α) t χ *. t (i)x N(µ, σ ) x µ σ N(0, ) (ii)x,, x N(µ, σ ) x = x+ +x N(µ, σ ) (iii) (i),(ii) z = x µ N(0, ) σ N(0, ) ( 9 97.

More information

1 filename=mathformula tex 1 ax 2 + bx + c = 0, x = b ± b 2 4ac, (1.1) 2a x 1 + x 2 = b a, x 1x 2 = c a, (1.2) ax 2 + 2b x + c = 0, x = b ± b 2

1 filename=mathformula tex 1 ax 2 + bx + c = 0, x = b ± b 2 4ac, (1.1) 2a x 1 + x 2 = b a, x 1x 2 = c a, (1.2) ax 2 + 2b x + c = 0, x = b ± b 2 filename=mathformula58.tex ax + bx + c =, x = b ± b 4ac, (.) a x + x = b a, x x = c a, (.) ax + b x + c =, x = b ± b ac. a (.3). sin(a ± B) = sin A cos B ± cos A sin B, (.) cos(a ± B) = cos A cos B sin

More information

δ ij δ ij ˆx ˆx ŷ ŷ ẑ ẑ 0, ˆx ŷ ŷ ˆx ẑ, ŷ ẑ ẑ ŷ ẑ, ẑ ˆx ˆx ẑ ŷ, a b a x ˆx + a y ŷ + a z ẑ b x ˆx + b

δ ij δ ij ˆx ˆx ŷ ŷ ẑ ẑ 0, ˆx ŷ ŷ ˆx ẑ, ŷ ẑ ẑ ŷ ẑ, ẑ ˆx ˆx ẑ ŷ, a b a x ˆx + a y ŷ + a z ẑ b x ˆx + b 23 2 2.1 n n r x, y, z ˆx ŷ ẑ 1 a a x ˆx + a y ŷ + a z ẑ 2.1.1 3 a iˆx i. 2.1.2 i1 i j k e x e y e z 3 a b a i b i i 1, 2, 3 x y z ˆx i ˆx j δ ij, 2.1.3 n a b a i b i a i b i a x b x + a y b y + a z b

More information

I

I I 6 4 10 1 1 1.1............... 1 1................ 1 1.3.................... 1.4............... 1.4.1.............. 1.4................. 1.4.3........... 3 1.4.4.. 3 1.5.......... 3 1.5.1..............

More information

( ) I( ) TA: ( M2)

( ) I( ) TA: ( M2) ( ) I( ) TA: ( M) 015 7 17 , 7 ( ) I( ).., M. (hatomura@spin.phys.s.u-tokyo.ac.jp).,,.. Keywords: 1. (gas-liquid phase transition). (critical point) 3. (lattice gas model) (Ising model) H = ϕ 0 i,j n i

More information

The Evaluation on Impact Strength of Structural Elements by Means of Drop Weight Test Elastic Response and Elastic Limit by Hiroshi Maenaka, Member Sh

The Evaluation on Impact Strength of Structural Elements by Means of Drop Weight Test Elastic Response and Elastic Limit by Hiroshi Maenaka, Member Sh The Evaluation on Impact Strength of Structural Elements by Means of Drop Weight Test Elastic Response and Elastic Limit by Hiroshi Maenaka, Member Shigeru Kitamura, Member Masaaki Sakuma Genya Aoki, Member

More information

Lennard Jones Steele * K N nm imac 2015 Mac mini OS *1 N V T

Lennard Jones Steele * K N nm imac 2015 Mac mini OS *1 N V T B 2016 4 2 2 1 1.1 Lennard Jones Steele *1 1.2 1.3 1. 2. 3. 4. 77 K N 2 2 2.5 nm imac 2015 Mac mini OS *1 N V T 1.3 3 1.3.1 1. imac 2. physicalchemistry 3. simulation simulation simulation simulation2

More information

IV (2)

IV (2) COMPUTATIONAL FLUID DYNAMICS (CFD) IV (2) The Analysis of Numerical Schemes (2) 11. Iterative methods for algebraic systems Reima Iwatsu, e-mail : iwatsu@cck.dendai.ac.jp Winter Semester 2007, Graduate

More information

koji07-02.dvi

koji07-02.dvi 007 I II III 1,, 3, 4, 5, 6, 7 5 4 1 ε-n 1 ε-n ε-n ε-n. {a } =1 a ε N N a a N= a a

More information