Name: analyzing melts and fluids from ab initio molecular dynamics simulations with the umd package
Uploaded: 2021-09-17T13:00:00
Duration: 6 min 37 s
Description: Watch this Scientific Journal Video about Analyzing Melts and Fluids from Ab Initio Molecular Dynamics Simulations with the UMD Package at JoVE.com

이 작품은 유럽 연합 호라이즌 2020 연구 및 혁신 프로그램 (RC에 대한 보조금 계약 번호 681818 영향)에 따라 유럽 연구 위원회 (ERC)에 의해 지원되었으며, 딥 탄소 천문대의 극단적 인 물리 및 화학 이사회에 의해, 그리고 우수 자금 조달 계획의 센터를 통해 노르웨이연구위원회에 의해 프로젝트 번호 223272. 우리는 STl2816 시리즈 eDARI 컴퓨팅 보조금을 통해 GENCI 슈퍼 컴퓨터에 대한 액세스를 인정, PRACE RA4947 프로젝트를 통해 아이린 AMD 슈퍼 컴퓨터, UNINETT 시그마2 NN9697K를 통해 프램 슈퍼 컴퓨터. FS는 마리 Skłodowska-Curie 프로젝트 (보조금 계약 ABISSE No.750901)에 의해 지원되었다.

<ol>
	<li>Frenkel, D., Smit, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Understanding Molecular Simulation. From Algorithms to Applications. , (2001).</li><li>Allen, M. P., Tildesley, D. J., Allen, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Computer Simulation of Liquids. , (1989).</li><li>Zepeda-Ruiz, L. A., Stukowski, A., Oppelstrup, T., Bulatov, V. V. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Probing+the+limits+of+metal+plasticity+with+molecular-dynamics+simulations.">Probing the limits of metal plasticity with molecular-dynamics simulations.</a> Nature Publishing Group. 550 (7677), 492-495 (2017).</li><li>Humphrey, W., Dalke, A., Schulten, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=VMD:+Visual+molecular+dynamics.">VMD: Visual molecular dynamics.</a> Journal of Molecular Graphics &amp; Modeling. 14 (1), 33-38 (1996).</li><li>Momma, K., Izumi, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=VESTA3+for+three-dimensional+visualization+of+crystal,+volumetric+and+morphology+data.">VESTA3 for three-dimensional visualization of crystal, volumetric and morphology data.</a> Journal of Applied Crystallography. 44 (6), 1272-1276 (2011).</li><li>Brehm, M., Kirchner, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=TRAVIS+-+A+free+Analyzer+and+Visualizer+for+Monte+Carlo+and+Molecular+Dynamics+Trajectories.">TRAVIS - A free Analyzer and Visualizer for Monte Carlo and Molecular Dynamics Trajectories.</a> Journal of Chemical Information and Modeling. 51 (8), 2007-2023 (2011).</li><li>Stixrude, L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Visualization-based+analysis+of+structural+and+dynamical+properties+of+simulated+hydrous+silicate+melt.">Visualization-based analysis of structural and dynamical properties of simulated hydrous silicate melt.</a> Physics and Chemistry of Minerals. 37 (2), 103-117 (2009).</li><li>Kresse, G., Hafner, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Ab+initio+Molecular-Dynamics+for+Liquid-Metals.">Ab initio Molecular-Dynamics for Liquid-Metals.</a> Physical Review B. 47 (1), 558-561 (1993).</li><li>Gygi, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Architecture+of+Qbox:+A+scalable+first-principles+molecular+dynamics+code.">Architecture of Qbox: A scalable first-principles molecular dynamics code.</a> IBM Journal of Research and Development. 52 (1-2), 137-144 (2008).</li><li>Harvey, J. P., Asimow, P. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Current+limitations+of+molecular+dynamic+simulations+as+probes+of+thermo-physical+behavior+of+silicate+melts.">Current limitations of molecular dynamic simulations as probes of thermo-physical behavior of silicate melts.</a> American Mineralogist. 100 (8-9), 1866-1882 (2015).</li><li>Caracas, R., Hirose, K., Nomura, R., Ballmer, M. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Melt-crystal+density+crossover+in+a+deep+magma+ocean.">Melt-crystal density crossover in a deep magma ocean.</a> Earth and Planetary Science Letters. 516, 202-211 (2019).</li><li>Green, M. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Markoff+Random+Processes+and+the+Statistical+Mechanics+of+Time-Dependent+Phenomena.+II.+Irreversible+Processes+in+Fluids.">Markoff Random Processes and the Statistical Mechanics of Time-Dependent Phenomena. II. Irreversible Processes in Fluids.</a> The Journal of Chemical Physics. 22 (3), 398-413 (1954).</li><li>Kubo, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Statistical-Mechanical+Theory+of+Irreversible+Processes.+I.+General+Theory+and+Simple+Applications+to+Magnetic+and+Conduction+Problems.">Statistical-Mechanical Theory of Irreversible Processes. I. General Theory and Simple Applications to Magnetic and Conduction Problems.</a> Journal of the Physical Society of Japan. 12 (6), 570-586 (1957).</li><li>Lin, S. T., Blanco, M., Goddard, W. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+two-phase+model+for+calculating+thermodynamic+properties+of+liquids+from+molecular+dynamics:+Validation+for+the+phase+diagram+of+Lennard-Jones+fluids.">The two-phase model for calculating thermodynamic properties of liquids from molecular dynamics: Validation for the phase diagram of Lennard-Jones fluids.</a> The Journal of Chemical Physics. 119 (22), 11792-11805 (2003).</li><li>Meyer, E. R., Kress, J. D., Collins, L. A., Ticknor, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Effect+of+correlation+on+viscosity+and+diffusion+in+molecular-dynamics+simulations.">Effect of correlation on viscosity and diffusion in molecular-dynamics simulations.</a> Physical Review E. 90 (4), 1198-1212 (2014).</li><li>Soubiran, F., Militzer, B., Driver, K. P., Zhang, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Properties+of+hydrogen,+helium,+and+silicon+dioxide+mixtures+in+giant+planet+interiors.">Properties of hydrogen, helium, and silicon dioxide mixtures in giant planet interiors.</a> Physics of Plasmas. 24 (4), 041401-041407 (2017).</li><li>Flyvbjerg, H., Petersen, H. G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Error+estimates+on+averages+of+correlated+data.">Error estimates on averages of correlated data.</a> The Journal of Chemical Physics. 91 (1), 461-466 (1989).</li><li>Tuckerman, M. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Statistical mechanics: theory and molecular simulation. , (2010).</li><li>McDonough, W. F., Sun, S. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+composition+of+the+Earth.">The composition of the Earth.</a> Chemical Geology. 120, 223-253 (1995).</li><li>Elkins-Taton, L. T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Magma+oceans+in+the+inner+solar+system.">Magma oceans in the inner solar system.</a> Annual Review of Earth and Planetary Sciences. 40, 113-139 (2012).</li><li>Lock, S. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+origin+of+the+Moon+within+a+terrestrial+synestia.">The origin of the Moon within a terrestrial synestia.</a> J. Geophysical Research: Planets. 123, 910-951 (2018).</li><li>Solomatova, N. V., Caracas, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Pressure-induced+coordination+changes+in+a+pyrolitic+silicate+melt+from+ab+initio+molecular+dynamics+simulations.">Pressure-induced coordination changes in a pyrolitic silicate melt from ab initio molecular dynamics simulations.</a> Journal of Geophysical Research: Solid Earth. 124, 11232-11250 (2019).</li></ol>

저자는 공개 할 것이 없습니다.

UMD 패키지는 ab initio 시뮬레이션을 통해 더 잘 작동하도록 설계되었으며, 스냅숏 수는 일반적으로 단위 셀당 수백 개의 원자가 있는 수만 개의 스냅숏으로 제한됩니다. 또한 후처리 실행이 충분한 활성 메모리 리소스를 가지고 있는 기계에 대해 더 큰 시뮬레이션을 견딜 수 있습니다. 이 코드는 계산할 수 있는 다양한 속성과 오픈 소스 라이선스로 구분합니다.
탯.dat 파일은 시뮬레이션 전체에서 변경되지 않은 파티클 수를 유지하는 앙상블에 적합합니다. UMD 패키지는 시뮬레이션 상자의 모양과 볼륨이 다른 계산에서 비롯된 파일을 읽을 수 있습니다. 이는 NVT 및 NPT와 같은 가장 일반적인 계산을 포함하며, 여기서 입자, N, 온도 T, 부피, V 및/또는 압력 P의 수가 일정하게 유지됩니다.
당분간 쌍 분포 기능뿐만 아니라 분석 스크립트와 같은 원자 간 거리를 추정하는 데 필요한 모든 스크립트는 축 사이의 각도가 90 °인 입방, 테트라고날 및 orthorhombic 셀에 대한 의미직기 단위 셀에 대해서만 작동합니다.
버전 2.0의 주요 개발 라인은 거리에 대한 직교 제한을 제거하고 개별 화학 결합을 분석하고, 원자 간 각도를 분석하고, 두 번째 조정 구를 구현하는 등 스펙션 스크립트에 대한 더 많은 기능을 추가하는 것입니다. 외부 협업의 도움으로 더 큰 시스템에서 더 빠른 분석을 위해 코드를 GPU에 이식하는 작업을 진행하고 있습니다.

유체 및 용융은 자연 환경에서 활성 화학 및 물리적 수송 벡터입니다. 원자 확산의 높은 비율은 화학 교류와 반응을 선호, 다양한 부력과 결합 낮은 점도 큰 질량 전송을 선호하고, 크리스탈 용융 밀도 관계는 행성 체 내부 레이어를 선호합니다. 주기적인 격자의 부재, 용융 상태에 도달하는 데 필요한 전형적인 고온, 담금질에 대한 어려움은 밀도, 확산 및 점도와 같은 일련의 명백한 특성의 실험적 결정을 매우 어렵게 만듭니다. 이러한 어려움으로 인해 대체 계산 메서드가 이 재료 클래스를 조사하기 위한 강력하고 유용한 도구를 만듭니다.
컴퓨팅 능력의 출현과 슈퍼컴퓨터의 가용성으로, 현재 비결정 분무 시스템, 몬테 카를로1 및 분자 역학(MD)1,2의 동적 상태를 연구하기 위해 두 가지 주요 수치 원자 시뮬레이션 기술이 사용됩니다. 몬테 카를로 시뮬레이션에서 구성 공간은 무작위로 샘플링됩니다. 모든 샘플링 관측이 서로 독립적이면 Monte Carlo 메서드는 병렬화로 선형 배율을 표시합니다. 결과의 품질은 난수 생성기의 품질과 샘플링의 대표성에 따라 달라집니다. 몬테 카를로 메서드는 샘플링이 서로 독립적인 경우 병렬화로 선형 배율을 표시합니다. 분자 역학(MD)에서 구성 공간은 시간 의존적 원자 궤적에 의해 샘플링된다. 지정된 구성에서 시작하여 원자 궤적은 뉴턴의 모션 방정식을 통합하여 계산됩니다. 원자간 힘은 모델 간 원자 적 잠재력 (고전 MD)을 사용하거나 첫 번째 원칙 방법 (ab initio 또는 첫 번째 원칙, MD)을 사용하여 계산할 수 있습니다. 결과의 품질은 궤적의 길이와 로컬 미니마에 끌리지 않는 능력에 따라 달라집니다.
분자 역학 시뮬레이션에는 시스템의 동적 거동과 관련된 다양한 정보가 포함되어 있습니다. 내부 에너지, 온도 및 압력과 같은 열역학 평균 특성은 계산하는 데 표준이 됩니다. 시뮬레이션의 출력 파일(들)에서 추출하고 평균화할 수 있는 반면 원자의 이동과 상호 관계와 직접적인 관련이 있는 수량은 원자 위치와 속도를 추출한 후 계산해야 합니다.
따라서 결과를 시각화하는 데 많은 노력이 전념했으며, 다양한 플랫폼에서 오픈 소스가 아닌 [Ovito3, VMD4, Vesta5, Travis6 등]에서 다양한 패키지를 사용할 수 있습니다. 이러한 모든 시각화 도구는 상호 원자 거리를 효율적으로 처리하므로 쌍 분포 함수및 확산 계수의 효율적인 계산을 허용합니다. 대규모 분자 역학 시뮬레이션을 수행하는 다양한 그룹은 시뮬레이션, 때로는 쉐어웨어 또는 커뮤니티에 대한 제한된 액세스의 다른 형태로 발생하는 다양한 다른 특성을 분석하는 독점 소프트웨어를 가지고 있으며 때로는 범위와 특정 패키지로 제한됩니다. 이러한 패키지 3,4,5,6,7 등에서 원자 결합, 기하학적 패턴 및 열역학에 대한 정보를 추출하는 정교한 알고리즘이 개발 및 구현됩니다.
여기서 우리는 분자 역학 시뮬레이션의 출력을 분석하기 위해 파이썬으로 작성된 오픈 소스 패키지인 UMD 패키지를 제안합니다. UMD 패키지를 사용하면 광범위한 구조적, 동적 및 열역학적 특성(그림 1)을 계산할 수 있습니다. 이 패키지는 GitHub 웹사이트(https://github.com/rcaracas/UMD_package)를 통해 제공되며 ERC IMPACT 프로젝트의 전용 페이지(http://moonimpact.eu/umd-package/)를 오픈 액세스 패키지로 사용할 수 있습니다.
보편적이고 쉽게 처리하기 위해, 우리의 접근 방식은 먼저 실제 분자 역학 실행의 출력 파일에서 열역학 상태 및 원자 궤적과 관련된 모든 정보를 추출하는 것입니다. 이 정보는 시뮬레이션이 실행된 원래 MD 패키지와 는 별개의 전용 파일에 저장됩니다. 유니버설 분자 역학을 의미하는 이러한 파일 "탯" 파일의 이름을 지정합니다. 이러한 방식으로, 우리의 UMD 패키지는 쉽게 적응의 최소 노력으로, 모든 소프트웨어와 AB 이니티오 그룹에 의해 사용할 수 있습니다. 현재 패키지를 사용하는 유일한 요구 사항은 아직 존재하지 않는 경우 특정 MD 소프트웨어의 출력에서 적절한 파서를 umd 파일 형식으로 작성하는 것입니다. 당분간, 우리는 VASP8 및 QBox9 패키지에 대한 이러한 파스터를 제공합니다.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/61534/61534fig01.jpg"> 그림 1: UMD 라이브러리의 순서도입니다. 물리적 특성은 파란색으로 되어 있으며 주요 Python 스크립트와 옵션은 빨간색입니다. <a href="https://www.jove.com/files/ftp_upload/61534/61534fig01large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>
탯줄 파일은 ASCII 파일입니다. 일반적인 확장은 "탯.dat"이지만 필수는 아닙니다. 모든 분석 구성 요소는 실제 이름 확장에 관계없이 UMD 형식의 ASCII 파일을 읽을 수 있습니다. 그러나 여러 시뮬레이션에서 빠른 대규모 통계를 수행하도록 설계된 일부 자동 스크립트는 특히 탯.dat 확장이 있는 파일을 찾습니다. 각 물리적 특성은 한 줄로 표현됩니다. 모든 선은 키워드로 시작됩니다. 이러한 방식으로 형식은 매우 적응력이 뛰어나며 버전 전반에 걸쳐 가독성을 유지하면서 새로운 속성을 umd 파일에 추가할 수 있습니다. 토론에 사용되는 4.6 GPa 및 3000 K에서 파이롤라이트 시뮬레이션의 탯줄 파일의 처음 30선은 그림 2에 나와 있다.
<img alt="Figure 2" class="xfigimg" src="/files/ftp_upload/61534/61534fig02.jpg"> 그림 2: 4.6 GPa 및 3000 K에서 액체 파이롤라이트의 시뮬레이션을 설명하는 탯줄 파일의 시작입니다. 헤더다음에각 스냅숏에 대한 설명이 뒤따릅니다. 각 속성은 물리적 속성의 이름, 값 및 단위의 이름을 포함하는 한 줄로 작성되며 모두 공간으로 구분됩니다. <a href="https://www.jove.com/files/ftp_upload/61534/61534fig02large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>
모든 탯줄 파일에는 원자, 전자 및 원자 유형의 수와 같은 시뮬레이션 셀의 함량을 설명하는 헤더와 유형, 화학 기호, 원자 전자 수 및 질량과 같은 각 원자에 대한 세부 정보가 포함되어 있습니다. 빈 선은 헤더의 끝을 표시하고 배드 파일의 주부와 분리합니다.
그런 다음 시뮬레이션의 각 단계가 상세합니다. 첫째, 즉각적인 열역학 파라미터는 각각 다른 선상에 각각 제공되며, 에너지, 응력, 동등한 유압, 밀도, 부피, 격자 매개변수 등(ii) 값 및 (iii) 단위와 같은 매개 변수의 이름을 지정합니다. 원자를 설명하는 테이블이 다음에 옵니다. 헤더 라인은 카르테시안 위치, 속도, 요금 및 장치와 같은 다양한 측정값을 제공합니다. 그런 다음 각 원자는 한 줄에 상세합니다. 세 개의 x, y, z 축에 해당하는 3개의 그룹으로, 항목은 감소된 위치, 시뮬레이션 셀로 접힌 카르테시안 위치, 카르테시안 위치(원자가 시뮬레이션 중에 여러 단위 셀을 통과할 수 있다는 사실을 적절히 고려함), 원자 속도 및 원자력입니다. 마지막 두 항목은 스칼라입니다: 충전 및 자기 순간.
두 개의 주요 라이브러리는 전체 패키지의 적절한 작동을 보장합니다. umd_process.py 라이브러리는 읽기 및 인쇄와 같은 탯줄 파일을 다룹니다. crystallography.py 라이브러리는 실제 원자 구조와 관련된 모든 정보를 다룹니다. crystallography.py 라이브러리의 기본 철학은 격자를 벡터 공간으로 취급하는 것입니다. 단위 셀 매개 변수와 방향은 기초 벡터를 나타냅니다. "Space"에는 일련의 스칼라 특성(특정 부피, 밀도, 온도 및 특정 수의 원자), 열역학 적 특성(내부 에너지, 압력, 열 용량 등), 일련의 텐서리얼 특성(응력 및 탄력)이 있습니다. 원자는 이 공간을 채웁니다. "격자" 클래스는 특정 볼륨, 밀도, 직접 적도 등에서 상호 격자를 얻는 등과 같은 몇 가지 짧은 계산과 함께이 앙상블을 정의합니다. "원자" 클래스는 원자를 정의합니다. 그들은 일련의 스칼라 속성 (이름, 기호, 질량, 전자 의 수 등)과 일련의 벡터 특성 (격자 클래스에 설명된 벡터 기준에 비해 공간의 위치 또는 보편적 인 카르테시안 좌표, 속도, 힘 등을 특징으로하는 위치)를 특징으로합니다. 이 두 클래스 외에도 crystallography.py 라이브러리에는 원자 거리 또는 셀 곱셈과 같은 다양한 테스트 및 계산을 수행하는 일련의 기능이 포함되어 있습니다. 요소의 주기적인 테이블도 사전으로 포함됩니다.
umd 패키지의 다양한 구성 요소는 여러 출력 파일을 작성합니다. 일반적으로 모든 ASCII 파일이며 모든 항목은 탭으로 구분되며 가능한 한 자명하게 만들어집니다. 예를 들어 물리적 속성과 단위를 항상 명확하게 나타냅니다. 탯.dat 파일은 이 규칙을 완전히 준수합니다.

<table>
	<tbody>
		<tr>
			<td>getopt library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>glob library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>matplotlib library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>numpy library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>os library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>Python software</td>
			<td>The Python Software Foundation</td>
			<td>Version 2 and 3</td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>random library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>re library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>scipy library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>subprocess library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
		<tr>
			<td>sys library</td>
			<td></td>
			<td></td>
			<td>open-source</td>
		</tr>
	</tbody>
</table>

analyzing melts and fluids from ab initio molecular dynamics simulations with the umd package

우리는 유체의 ab initio 분자 역학 시뮬레이션에서 유래 결과를 분석하기 위해 파이썬 기반 오픈 소스 패키지를 개발했습니다. 이 패키지는 규산염 및 산화물 용융, 수성 유체 및 다양한 초임계 유체와 같은 천연 시스템의 응용 분야에 가장 적합합니다. 패키지는 파일 형식을 처리하고 결정예로 처리 두 개의 주요 라이브러리를 포함하는 파이썬 스크립트의 모음입니다. 모든 스크립트는 명령줄에서 실행됩니다. 우리는 원자 궤적과 보편적 인 분자 역학을 위해 서, UMD 파일에 저장되는 시뮬레이션의 관련 열역학 정보를 저장하는 단순화 된 형식을 제안한다. UMD 패키지를 사용하면 일련의 구조, 운송 및 열역학 적 특성을 계산할 수 있습니다. 결합 길이를 정의하고, 원자간 연결 매트릭스를 구축하고, 결국 화학적 표본을 결정합니다. 화학 종의 수명을 결정하는 것은 전체 통계 분석을 실행할 수 있습니다. 그런 다음 전용 스크립트는 원자뿐만 아니라 화학 종에 대한 평균 제곱 변위를 계산합니다. 원자 속도의 구현된 자기 상관 관계 분석은 확산 계수와 진동 스펙트럼을 산출합니다. 응력에 적용된 동일한 분석은 점도를 산출합니다. 이 패키지는 GitHub 웹 사이트를 통해 ERC IMPACT 프로젝트의 전용 페이지를 오픈 액세스 패키지로 사용할 수 있습니다.

Pyrolite는 다중 구성 요소 규산용용(0.5Na2O 2CaO 1.5Al2O3 4FeO 30MgO 24SiO2)으로, 철기반 코어19를 제외한 지구 화학적 평균 또는 지구 전체의 구성을 가장 근사화하는 모델다. 초기 지구는 일련의 대규모 용융 사건에 의해 지배되었다20, 마지막 하나는 프로토 달 디스크에 대한 응축 후, 전체 행성을 침몰 했을 수 있습니다21. Pyrolite는 이러한 행성 규모의 마그마 바다의 구성에 가장 근사한 을 나타냅니다. 따라서 VASP 구현에서 ab initio 분자 역학 시뮬레이션에서 3,000\u20125,000 K 온도 범위및 0\u2012150 GPa 압력 범위에서 용융의 물리적 특성을 광범위하게 연구했습니다. 이러한 열역학 적 조건은 전적으로 지구에서 가장 극단적 인 마그마 해양 조건을 특성화합니다. 우리의 연구는 melts22의 전체 심층 분석을 위한 UMD 패키지를 성공적으로 사용하는 훌륭한 예입니다. 우리는 분포와 평균 결합 길이를 계산하고, 양이온 산소 조정의 변화를 추적하고, 다양한 조성물의 무정형 규산염에 대한 이전의 실험 및 계산 연구와 결과를 비교했습니다. 우리의 심층 분석은 표준 조정 번호를 기본 성분으로 분해하고, 용융시 이국적인 조정 폴리히드라의 존재를 설명하고, 모든 조정 폴리히드라의 수명을 추출하는 데 도움이되었습니다. 또한 궤적의 길이와 모델링된 시스템에 존재하는 원자의 수 의 관점에서 시뮬레이션에서 샘플링의 중요성을 설명했습니다. 사후 처리에 관해서는 UMD 분석은 이러한 요인과 독립적이지만 UMD 패키지에서 제공하는 결과를 해석할 때 고려해야 합니다. 여기서는 UMD 패키지를 사용하여 용융 용융 용융에 대한 응용 프로그램과 함께 용융의 몇 가지 특징을 추출하는 방법에 대한 몇 가지 예를 보여줍니다.
gofrs_umd.py 스크립트에서 얻은 Si-O 쌍 분포 함수는 g(r) 함수의 첫 번째 최소인 첫 번째 조정 구의 반경이 T = 3000 K 및 P = 4.6 GPa에서 약 2.5 개의 앙스트롬에 있음을 보여줍니다. g(r)의 최대 값은 1.635 Å에 있으며, 이는 굽힘 길이에 가장 가장 가까운 근사치입니다. 긴 꼬리는 온도 때문입니다. 이 제한을 Si-O 본드 거리로 사용하여, 분석 결과 최대 몇 피코초 동안 지속될 수 있는 SiO4 유닛이 용융(그림 5)을 지배한다는 것을 알 수 있습니다. Si2O7과 같은 디머와 Si3Ox 유닛과 같은 트리머의 존재에 의해 반사되는 부분 중합을 보여주는 용융의 중요한 부분이 있습니다. 그들의 해당 수명은 피코초의 순서입니다. 고차 폴리머는 모두 수명이 상당히 짧습니다.
<img alt="Figure 5" class="xfigimg" src="/files/ftp_upload/61534/61534fig05.jpg"> 그림 5: 시오 화학 종의 수명. 이 스펙션은 4.6 GPa 및 3000 K에서 다중 성분 용융에서 식별됩니다. 라벨은 SiO3, SiO4 및 SiO5 단조체 및 다양한 SixOy 폴리머를 표시합니다. <a href="https://www.jove.com/files/ftp_upload/61534/61534fig05large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>
위의 -z 및 -v 플래그에 의해 정의된 수직 및 수평 단계의 다른 값은 MSD의 다양한 샘플링을 생성합니다(그림 6). z와 v의 큰 값조차도 경사면을 정의하기에 충분하므로 다른 원자의 확산 계수를 정의할 수 있습니다. z와 v의 큰 값으로 갈 때 후 처리에 대한 시간의 이득은 놀랍습니다. MSD는 시뮬레이션의 품질에 대한 매우 강력한 검증 기준을 제공합니다. MSD의 확산 부분이 충분히 길지 않은 경우 시뮬레이션이 너무 짧고 통계적 의미에서 유체 상태에 도달하지 못한다는 신호입니다. MSD의 확산 부분에 대한 최소 요구 사항은 시스템에 따라 다릅니다. 모든 원자가 유체10으로 간주되기 위해 용융 구조에서 적어도 한 번 은 자신의 부위를 변경하도록 요구할 수 있습니다. 행성 과학의 응용 프로그램과 훌륭한 예는 복잡한 규산염은 액체 라인11에 가깝거나 아래에 있는 고압에서 녹는 복잡한 예입니다. 주요 네트워크 형성 양이온인 Si 원자는 22개 이상의 picoseconds 후에 사이트를 전환합니다. 이 임계값보다 짧은 시뮬레이션은 가능한 구성 공간을 상당히 과소 샘플링할 수 있습니다. 그러나, 조정 된 음의, 즉 O 원자, 중앙 Si 원자보다 빠르게 이동, 그들은 Si의 느린 이동성의 일부를 보상 할 수 있습니다. 따라서 전체 시스템은 Si 변위에서만 가정하는 것보다 구성 공간의 샘플링을 더 잘 다룰 수 있습니다.
<img alt="Figure 6" class="xfigimg" src="/files/ftp_upload/61534/61534fig06.jpg"> 그림 6: 평균 사각형 변위(MSD). MSD는 다중 성분 규산염 용융의 몇 가지 원자 유형에 대해 설명되어 있습니다. 다양한 수평 및 수직 단계, z 및 v가 있는 샘플링은 일관된 결과를 산출합니다. 솔리드 서클: -z 50 -v 50. 오픈 서클: -z 250 -v 500. <a href="https://www.jove.com/files/ftp_upload/61534/61534fig06large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>
마지막으로, 원자 VAC 기능은 용융의 진동 스펙트럼을 산출합니다. 도 7 은 위와 동일한 압력 및 온도 조건에서 스펙트럼을 나타낸다. 우리는 Mg, Si 및 O 원자의 기여뿐만 아니라 총 가치를 나타냅니다. 0 주파수에서 용융의 확산 특성에 해당하는 스펙트럼의 유한 값이 있습니다. 진동 스펙트럼에서 열역학적 특성을 추출하려면 이 가스와 같은 확산 특성을 0에서 제거해야 하지만 더 높은 주파수에서 부패를 적절히 고려해야 합니다.
<img alt="Figure 7" class="xfigimg" src="/files/ftp_upload/61534/61534fig07.jpg"> 그림 7: 용융물의 진동 스펙트럼. 원자 속도-속도 자기 상관 기능의 Fourier 변환의 실제 부분은 진동 스펙트럼을 산출합니다. 여기서 스펙트럼은 다중 구성 요소 규산염 용융을 위해 계산됩니다. 유체는 0 주파수에서 제로 가스와 같은 확산 특성을 가지고 있습니다. <a href="https://www.jove.com/files/ftp_upload/61534/61534fig07large.jpg" target="_blank">이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.</a>

Watch this Scientific Journal Video about Analyzing Melts and Fluids from Ab Initio Molecular Dynamics Simulations with the UMD Package at JoVE.com

Analyzing Melts and Fluids from Ab Initio Molecular Dynamics Simulations with the UMD Package

용융 및 유체는 자연 시스템에서 대량 수송의 유비쿼터스 벡터입니다. 우리는 이러한 시스템의 ab initio 분자 역학 시뮬레이션을 분석하기 위해 오픈 소스 패키지를 개발했습니다. 우리는 구조 (접합, 클러스터화, 화학 적 분석), 운송 (확산, 점도) 및 열역학 적 특성 (진동 스펙트럼)을 계산합니다.

UMD 패키지로 Ab Initio 분자 역학 시뮬레이션에서 용융 및 유체 분석

We have developed a Python-based open-source package to analyze the results stemming from ab initio molecular-dynamics simulations of fluids. The package ...

1차 분자 역학 시뮬레이션을 분석하면 유체의 물리적 및 화학적 특성을 정확하게 예측할 수 있습니다. 이 방법은 물리학 및 화학의 모든 원자 분석에 적용할 수 있습니다. 먼저 패키지에서 하나 이상의 전용 Python 스크립트를 사용하여 각 특정 물리적 특성 집합을 추출합니다.명령줄에서 모든 스크립트를 실행합니다. 그들은 모두 하나의 스크립트에서 다른 스크립트에 가능한 한 일관된 일련의 플래그를 사용합니다. 첫 번째 원칙 코드에서 수행되는 MD 시뮬레이션의 출력을 UMD 파일로 변환한 다음 탯줄 파일을 xyz 파일로 전송하여 VMD 또는 VESTA와 같은 다양한 다른 패키지에서 시각화를 용이하게 합니다.UMD 파일을 umd2poscar를 사용하여 VASP 형 POSCAR 파일로 전환합니다. 미리 정의된 주파수로 시뮬레이션의 스냅샷을 선택합니다. gofrs_umd 실행합니다.원자 유형 A및 B.의 모든 쌍에 대한 쌍 분포 함수를 계산하는 py 스크립트 출력은 확장 gofrs.dat 구분된 하나의 ASCII 파일 탭으로 작성됩니다. 첫 번째 협착 구체의 반으로 평균 원자 간 결합 거리를 추출합니다. 이를 위해 gofrs를 플로팅하여 쌍 분포 함수의 첫 번째 최대값의 위치를 식별합니다.스프레드시트 응용 프로그램에서 파일을 dat하고 각 원자 쌍의 최대마 및 최소마를 검색한 다음 스프레드시트 소프트웨어를 사용하여 첫 번째 조정 구의 반경을 PDF의 첫 번째 최소값으로 식별합니다. 응성 매트릭스를 획득하고 조정 폴리헤드라 또는 중합을 얻기 위해 스펙션 스크립트를 실행합니다. speciation_umd 실행합니다.첫 번째 레벨에서 연결 그래프를 샘플링하여 조정 폴리헤드라를 식별하는 플래그 r0이 있는 py 스크립트입니다. speciation_umd 실행합니다. 모든 깊이 수준에서 연결 그래프를 샘플링하여 중합화를 얻는 플래그 r1을 가진 py 스크립트.구풀에서 발견되는 시뮬레이션에서 발견되는 모든 화학 종의 각 원자 클러스터의 수명을 플롯합니다. dat 파일. 원자의 평균 제곱 변위 또는 MSD를 자체 확산도를 얻기 위한 시간의 함수로서 추출한 다음 일련의 msd_umd 사용하여 MSD를 계산합니다.py 스크립트를 작성하고 각 원자 유형의 평균 MSD를 계산합니다. 각 원자와 화학 종의 MSD를 계산합니다. 스프레드시트 기반 소프트웨어를 사용하여 MSD를 플롯하고 MSD의 경사에서 확산 계수를 계산합니다.vibr_spectrum_umd 실행합니다. py 스크립트는 각 원자 유형에 대한 원자 속도 자동 상관 또는 VAC 기능을 계산하고 빠른 Fourier 변환을 수행합니다. 진동 스펙트럼을 진동 스펙트럼에서 플롯합니다.스프레드시트와 같은 소프트웨어를 사용하여 파일을 dat합니다. 오메가에서 유한값을 식별하는 것은 0과 같으며, 이는 유한 주파수에서 스펙트럼의 다양한 피크에서 유체의 확산 특성에 해당합니다. 평균을 실행합니다.py는 UMD 파일에서 압력, 온도, 밀도 및 내부 에너지에 대한 평균 값및 스프레드를 추출합니다. 마지막으로 전체 평균을 실행합니다. py 스크립트는 평균의 오류를 포함하여 완전한 통계 분석을 수행합니다.Pyrolite는 대량 규산염 지구의 조성물을 가장 근사화하는 다중 구성 요소 규산염 용융 모델입니다. UMD 패키지는 용융 된 파이롤 라이트의 몇 가지 특징을 추출하는 데 사용되었습니다. 실리콘-산소 쌍 분포 함수의 최대는 1.635 angstrom에 놓여 있으며 이는 굴곡 길이에 가장 가장 근사합니다.이 제한을 실리콘-산소 결합 거리로 사용하여, 표본 분석은 최대 몇 pico초 동안 지속될 수 있는 orthosilicate 단위가 용융을 지배한다는 것을 보여줍니다. Si3Ox 유닛과 같은 디실리케이트 및 트리머와 같은 디머의 존재에 의해 반사되는 부분 중합화를 보여주는 용융의 중요한 부분이 있다. 그들의 해당 수명은 피코초의 순서입니다.높은 순서 폴리머는 모두 수명이 상당히 짧습니다. 수직 및 수평 단계의 다른 값은 MSD의 다양한 샘플링을 생성합니다. Z와 V의 큰 값조차도 경사면을 정의하기에 충분하므로 다른 원자의 확산 계수를 정의할 수 있습니다.Z 및 V.Finally의 큰 값에 대해 처리 후 시간이 크게 증가하며, 원자 속도 자가 상관 기능은 용융의 진동 스펙트럼을 생성합니다. 여기에 마그네슘, 실리콘 및 산소 원자의 기여뿐만 아니라 총 값이 표시됩니다. 이 프로토콜을 시도할 때는 항상 수렴을 확인합니다.원자 궤적이 관심있는 현상을 제대로 포착할 수 있을 만큼 충분히 길다는 것을 확인하십시오. 이 기술은 시뮬레이션 결과의 후처리를 다룹니다. 시뮬레이션과 해당 분석은 병렬로 수행해야 합니다.

Research

JoVE Journal

Chemistry

Isolation and Chemical Characterization of Lipid A from Gram-negative Bacteria

그람 음성 박테리아에서 분리 및 지질의 화학 특성

Lipopolysaccharide (LPS) is the major cell surface molecule of gram-negative bacteria, deposited on the outer leaflet of the outer membrane bilayer. LPS ...

연구

화학

Analyzing Protein Dynamics Using Hydrogen Exchange Mass Spectrometry

수소 환율 질량 분석을 사용하여 단백질의 역학 분석

All cellular processes depend on the functionality of proteins. Although the functionality of a given protein is the direct consequence of its unique ...

Visualizing Lignification Dynamics in Plants with Click Chemistry: Dual Labeling is BLISS!

Lignification 역학 클릭 화학 식물에서 시각화: 행복은 듀얼 라벨!

Lignin is one of the most prevalent biopolymers on the planet and a major component of lignocellulosic biomass. This phenolic polymer plays a vital ...

Combustion Chemistry of Fuels: Quantitative Speciation Data Obtained from an Atmospheric High-temperature Flow Reactor with Coupled Molecular-beam Mass Spectrometer

연료의 연소 화학: 양적 Speciation 데이터와 대기 고 열 흐름 반응 기에서 얻은 결합 분자 빔 질량 분석기

This manuscript describes a high-temperature flow reactor experiment coupled to the powerful molecular beam mass spectrometry (MBMS) technique. This ...

Probing the Structure and Dynamics of Interfacial Water with Scanning Tunneling Microscopy and Spectroscopy

구조와 계면 물 검색과 터널링 현미경 분광학의 역학 조사

Water/solid interfaces are ubiquitous and play a key role in many environmental, biophysical, and technological processes. Resolving the internal ...

A Closed-Type Wireless Nanopore Electrode for Analyzing Single Nanoparticles

단일 나노 입자 분석에 대 한 폐쇄-타입 무선 Nanopore 전극

Measuring the intrinsic features of single nanoparticles by nanoelectrochemistry holds deep fundamental importance and has potential impacts in ...

Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics

조밀도 기능적인 이론 및 분자 역학 포스 필드를 사용 하 여 이기종 물/금속 촉매 인터페이스의 multiscale 샘플링

A significant number of heterogeneously-catalyzed chemical processes occur under liquid conditions, but simulating catalyst function under such conditions ...

Vibrational Spectra of a N719-Chromophore/Titania Interface from Empirical-Potential Molecular-Dynamics Simulation, Solvated by a Room Temperature Ionic Liquid

실온 이온 액체에 의해 용해된 경험적 전위 분자 역학 시뮬레이션의 N719-크로모포어/티타니아 인터페이스의 진동 스펙트럼

The accurate molecular simulation prediction of vibrational spectra, and other structural, energetic and spectral characteristics, of photo-active ...

Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry

Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry

ab initio 열화학에서 분자 클러스터의 대기 농도 계산

The computational study of the formation and growth of atmospheric aerosols requires an accurate Gibbs free energy surface, which can be obtained from gas ...

Interactive Molecular Model Assembly with 3D Printing

3D 프린팅을 갖춘 대화형 분자 모델 어셈블리

With the growth in accessibility of 3D printing, there has been a growing application of and interest in additive manufacturing processes in chemical ...