Name: computation of atmospheric concentrations of molecular clusters from ab initio thermochemistry
Uploaded: 2020-04-08T15:00:00
Duration: 12 min 11 s
Description: Watch this Scientific Journal Video about Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry at JoVE.com

이 프로젝트는 보조금 CHE-1229354, CHE-1662030, CHE-1721511, 그리고 CHE-1903871 국립 과학 재단에서 지원되었다 (GCS), 아놀드와 메이블 베크만 재단 베크만 장학생 상 (AGG), 배리 M. 골드 워터 장학금 (AGG). 머큐리 컨소시엄(http://www.mercuryconsortium.org)의 고성능 컴퓨팅 리소스가 사용되었습니다.

1. 분리된 글리신과 물의 최소 에너지 구조 찾기
참고: 여기서 의 목표는 두 가지입니다: (i) 유전 알고리즘 구성 샘플링에 사용하기 위해 단리된 물과 글리신 분자의 최소 에너지 구조를 얻고, (ii) 대기 농도 계산에 사용하기 위해 이들 분자의 기체 상 에너지에 대한 열역학적 보정을 계산하는 것입니다.
<ol><li>로컬 컴퓨터에서 Avogadro의 새 세션을 엽니다.<ol>
<li>빌드 > 삽입 > 펩티드를 클릭하고 인서트 펩타이드 창에서 글리신 단량체를 생성하려면 글리시를 선택합니다.</li>
		<li>확장 > 가우시안을 클릭하고 텍스트 상자의 첫 번째 줄을 편집하여 '# pw91pw91/6-311++G** int(Acc2E=12,UltraFine) scf(수렴=12) 옵트(꽉, maxcyc=300) 프런시크'를 읽습니다. 입력 파일을 생성하고 저장을 glycine.com.</li>
		<li>분자가 글리신이55를하는 것처럼 분자가 상당한 변형 유연성을 가지고 있다면, 글로벌 최소 구조 및 기타 저지대 적합자를 식별하기 위해 형태 분석을 수행하는 것이 중요합니다. OpenBabel54는 다양한 알고리즘과 빠른 힘 필드를 활용하는 강력한 형태 검색 도구를 제공합니다. conformers는 GA 및 후속 계산 중에 이완및 상호 변환할 수 있지만, 때로는 여러 GA 계산을 실행해야 하는 경우가 있습니다.</li>
	</ol>
</li>
	<li>로컬 컴퓨터에서 Avogadro에서 새 세션을 엽니다.<ol>
<li>빌드 > 삽입 > 조각 및 삽입 조각 창에서 "물"을 검색하여 물 좌표를 가져옵니다.</li>
		<li>확장 > 가우시안을 클릭하고 텍스트 상자의 첫 번째 줄을 편집하여 '# pw91pw91/6-311++G** int(Acc2E=12,UltraFine) scf(수렴=12) 옵트(꽉, maxcyc=300) 프런시크'를 읽습니다. 입력 파일을 생성하고 저장을 water.com.</li>
	</ol>
</li>
	<li>두 .com 파일을 원격 클러스터로 전송합니다. 원격 클러스터에 로그인하면 일괄 처리 제출 스크립트에서 Gaussian 09를 호출하여 계산을 시작합니다. 계산이 완료되면 OpenBabel을 호출하여 최소 에너지 구조의 카르테시안좌표(.xyz 파일)를 추출합니다. 글리신의 경우 실행할 명령은 다음과 같은 것입니다. 오바벨 -ig09 글리신.log -oxyz > 글리신.xyz 이 두 .xyz 파일은 다음 단계에서 GA 구성 샘플링에서 사용됩니다.</li>
</ol>2. 글리 (H2O)n = 1-5 클러스터의 유전자 알고리즘 기반 구성 샘플링
참고: 여기서 의 목표는 MOPAC47에서구현된 PM746 모델을 사용하여 저렴한 반 경험적 이론 수준에서 글리(H2 O)n=1-5에 대한 저에너지 구조 세트를 얻는 것입니다.2 그림 2와같이 작업 디렉토리에 정확한 조직과 구조가 있는 것이 필수적입니다. 이는 사용자 지정 셸 및 Python 스크립트가 오류 없이 작동하도록 하기 위한 것입니다.
<ol><li>필요한 모든 스크립트를 원격 클러스터에 복사하고 해당 위치를 추가하여 $PATH<ol>
<li>모든 스크립트 및 템플릿 파일을 폴더(예: 스크립트)에 넣고 원격 클러스터에 복사합니다.</li>
		<li>모든 스크립트가 실행 가능한지 확인</li>
		<li>터미널에서 다음 명령을 입력하여 스크립트 디렉터리 위치를 $PATH 환경 변수에 추가합니다. 스크립트의 기본 위치는 $HOME/JoVE-demo/script로 설정되어 있지만 스크립트가 포함된 디렉터리를 가리키는 $SCRIPTS_HOME이라는 환경 변수를 정의하고 경로에 $SCRIPTS_HOME을 추가할 수 있습니다.<ol>
<li>배쉬 쉘: 내보내기 SCRIPTS_HOME=/경로/에/스크립트 내보내기 PATH=${SCRIPTS_HOME}:${PATH}</li>
			<li>Tcsh/Csh 쉘: 세텐프 SCRIPTS_HOME /경로 / 에 / 스크립트 setenv PATH ${SCRIPTS_HOME}:${PATH}</li>
		</ol>
</li>
	</ol>
</li>
	<li>원격 클러스터에서 GA 계산을 설정하고 실행합니다.<ol>
<li>gly-h2o-n이라는 디렉토리를 만들었는데 여기서 n은 물 분자의 수입니다.</li>
		<li>gly-h2o-n 디렉터리 아래에 GA라는 하위 디렉터리를 만들어 유전자 알고리즘 계산을 실행합니다.</li>
		<li>OGOLEM 입력 파일(예: pm7.ogo), 단량체 카르테시안 좌표(예: 글리신.xyz, water.xyz) 및 PBS 일괄 제출 스크립트(예: run.pbs)를 GA 디렉토리에 복사합니다.</li>
		<li>OGOLEM 입력 파일 및 일괄 처리 제출 파일을 필요한 변경합니다.</li>
		<li>계산을 제출합니다. 계산이 시작되면 OGOLEM은 GA 디렉터리에 OGOLEM 입력 파일(예: pm7)의 접두사로 명명된 새 디렉터리를 만들고 새로 생성된 좌표를 저장합니다.</li>
	</ol>
</li>
	<li>계산이 완료되면 에너지 및 회전 상수를 컴파일하고 해당 정보를 사용하여 고유한 저에너지 구조가 결정됩니다.<ol>
<li>글리-h2o-n/GA/pm7로 디렉토리 변경 및</li>
		<li>에너지를 추출하고 명령을 통해 GA 에 최적화된 클러스터의 회전 상수를 계산합니다. getRotConsts-GA.csh N 0 99 여기서 N은 분자 클러스터의 원자 수와 '0 99'는 GA 풀 크기가 100이고 인덱스는 0에서 99까지 실행된다는 것을 나타냅니다. 이렇게 하면 rotConstsData_C 라는 파일이 생성되며, 여기에는 GA에 최적화된 모든 클러스터 구성, 에너지 및 회전 상수의 정렬된 목록이 포함됩니다.</li>
		<li>명령을 실행합니다. similarityAnalysis.py pm7 rotConstsData_C 여기서 pm7은 파일 이름 지정 레이블로 사용되어 고유한 GA 최적화 클러스터를 찾고 저장합니다. 이렇게 하면 고유한 GA 최적화 구성의 정렬된 목록이 포함된 고유구조스-pm7.data라는 파일이 생성됩니다. 이것은 이론의 PM7 수준에서 최적화된 Gly(H2O)n 클러스터에 대한 고유한 로컬 최소 구조 목록이며, 이러한 구조는 이제 DFT를 사용하여 구체화할 준비가 되었습니다.</li>
	</ol>
</li>
	<li>gly-h2o-n/GA 디렉토리로 이동하여 combine-GA.csh 스크립트를 사용하여 여러 유사한 GA 실행의 결과를 결합합니다. 구문은 다음과 입니다. 결합-GA.csh & 레이블> <와 GA 실행 및 > 이 특별한 경우 명령은 다음과 같은 것입니다. 결합 GA.csh pm7 pm7 글리-h2o-n/GA 디렉토리에서'고유 구조-pm7.data'라는이름의 새로운 고유 구조 목록을 생성합니다.</li>
</ol>3. 작은 기초 세트와 QM 방법을 사용하여 정제
참고: 여기서 목표는 더 나은 양자-기계적 설명을 사용하여 Gly(H2O)n=1-5 클러스터 구조의 구성 샘플링을 구체화하여 더 작지만 더 정확한 Gly(H2 O)n=1-5 클러스터 구조 집합을 얻는 것입니다.2 이 단계의 시작 구조는 2단계의 출력입니다.
<ol><li>작은 기준 세트 DFT 계산을 준비하고 실행합니다.<ol>
<li>글리 h2o-n 디렉토리 아래에 QM이라는 하위 디렉토리를 만듭니다. QM 디렉터리에서 pw91-sb라는다른 하위 디렉터리를 만듭니다.</li>
		<li>글리-h2o-n/GA 디렉토리에서 QM/pw91-sb 디렉토리에 고유한 구조목록(고유구조-pm7.data)을복사합니다.</li>
		<li>그 gly-h2o-n/QM/pw91-sb로디렉토리를 변경합니다.</li>
		<li>명령을 사용하여 작은 기준 세트 DFT 구성 샘플링 스크립트를 실행합니다. run-pw91-sb.csh 고유구조스-pm7.data sb QUEUE 10 여기서 sb는 이 계산 집합의 레이블이며 QUEUE는 컴퓨팅 클러스터에서 선호되는 큐이며 10개는 10개의 계산을 하나의 일괄 처리 작업으로 그룹화해야 한다는 것을 나타냅니다. 이 스크립트는 Gaussian 09에 대한 입력을 자동으로 생성하고 모든 계산을 제출합니다. 드라이 런을 수행하려면 'QUEUE'에 대한 '테스트'를 입력합니다.</li>
	</ol>
</li>
	<li>제출된 계산이 완료되면 결과를 추출하고 분석합니다.<ol>
<li>에너지를 추출하고 명령을 사용하여 소규모 로 최적화된 클러스터의 회전 상수를 계산합니다. getRotConsts-dft-sb.csh pw91 N 여기서 pw91은 PW91 밀도 함수가 사용되었고 N은 클러스터의 원자 수임을 나타냅니다. 그렇게 하면 rotConstsData_C라는 파일이 생성됩니다.</li>
		<li>이제 명령을 통해 고유한 구조를 식별합니다. similarityAnalysis.py sb rotConstsData_C 여기서 sb가 파일 이름 지정 레이블로 사용됩니다. 이제 파일 고유 구조-sb.data에저장된 이론의 PW91/6-31+G* 수준에서 최적화된 고유한 구성 목록이 표시됩니다.</li>
	</ol>
</li>
	<li>gly-h2o-n/QM 디렉토리로 이동하여 결합-QM.csh 스크립트를 사용하여 여러 유사한 QM 실행의 결과를 결합합니다. 구문은 다음과 입니다. QM.csh&label>&&lt;QM calcs>gt; 이 특별한 경우 명령은 다음과 같은 것입니다. 결합-QM.csh sb pw91-sb 글리-h2o-n/QM 디렉토리에서'고유 구조-sb.data'라는이름의 새로운 고유 구조 목록을 생성합니다.</li>
</ol>4. 큰 기초 세트와 QM 방법을 사용하여 추가 정제
참고: 여기서 목표는 더 나은 양자-기계적 설명을 사용하여2Gly(H2 O)n=1-5 클러스터의 구성 샘플링을 더욱 구체화하는 것입니다. 이 단계의 시작 구조는 3단계의 출력입니다.
<ol><li>더 큰 기준 집합을 사용하여 보다 신뢰할 수 있는 계산을 제출합니다.<ol>
<li>QM 디렉토리 아래에 pw91-lb라는 하위 디렉토리를 만듭니다.</li>
		<li>글리-h2o-n/QM 디렉토리에서 글리-h2o-n/QM/pw91-lb 디렉토리에 고유한 구조목록(고유Structures-sb.data)을복사하고 해당 디렉토리로 변경합니다.</li>
		<li>다음과 같은 명령으로 대단위 DFT 구성 샘플링 스크립트를 실행합니다. run-pw91-lb.csh 고유구조스-sb.data lb QUEUE 10 여기서 lb가 이 계산 집합의 레이블인 경우 QUEUE는 컴퓨팅 클러스터에서 선호되는 큐이며 10은 10개의 계산을 하나의 일괄 처리 작업으로 그룹화해야 한다는 것을 나타냅니다. 이 스크립트는 Gaussian 09에 대한 입력을 자동으로 생성하고 모든 계산을 제출합니다. 'QUEUE'에 대한 '테스트'를 입력하여 드라이 런 테스트를 수행합니다.</li>
	</ol>
</li>
	<li>제출된 계산이 완료되면 데이터를 추출하고 분석합니다.<ol>
<li>다음 명령을 통해 대용량 최적화 클러스터의 회전 상수를 계산합니다. getRotConsts-dft-lb.csh pw91 N 여기서 pw91은 PW91 밀도 함수가 사용되었고 N은 클러스터의 원자 수임을 나타냅니다.</li>
		<li>이제 명령을 통해 고유한 구조를 식별합니다. similarityAnalysis.py 파운드 rotConstsData_C 여기서 lb가 파일 이름 지정 레이블로 사용됩니다. 이제 파일 고유 구조-lb.data에저장된 이론의 PW91/6-311++G** 수준에서 최적화된 고유한 구성 목록이 있습니다.</li>
	</ol>
</li>
</ol>5. 최종 에너지 및 열역학 보정 계산
참고: 여기서 목표는 원하는 열화학 적 보정을 계산하기 위해 큰 기초 세트와 초미세 통합 그리드를 사용하여 Gly(H2O)n = 1-5 클러스터의 진동 구조와 에너지를 얻는 것입니다.
<ol><li>이전 단계의 결과로 시작하여 보다 신뢰할 수 있는 계산을 제출합니다.<ol>
<li>QM/pw91-lb 디렉토리 아래에 초미세라는 하위 디렉토리를 만듭니다. 그런 다음 QM/pw91-lb 디렉토리에서 QM/pw91-lb/초미세 디렉토리에 고유한 구조목록(고유Structures-lb.data)을복사하고 해당 디렉토리로 변경합니다.</li>
		<li>명령과 함께 초미세 대용량 DFT 스크립트를 제출합니다. run-pw91-lb-ultrafine.csh 고유구조-lb.data uf QUEUE 10 uf가 이 계산 집합의 레이블인 경우 QUEUE는 컴퓨팅 클러스터에서 선호되는 큐이며 10개는 10개의 계산을 하나의 일괄 처리 작업으로 그룹화해야 한다는 것을 나타냅니다. 이 스크립트는 Gaussian 09에 대한 입력을 자동으로 생성하고 모든 계산을 제출합니다. 'QUEUE'에 대한 '테스트'를 입력하여 드라이 런 테스트를 수행합니다.</li>
	</ol>
</li>
	<li>제출된 계산이 완료되면 데이터를 추출하고 분석합니다.<ol>
<li>에너지를 추출하고 명령을 통해 대용량 최적화 클러스터의 회전 상수를 계산합니다. getRotConsts-dft-lb-ultrafine.csh pw91 N 여기서 pw91은 PW91 밀도 함수가 사용되었고 N은 클러스터의 원자 수임을 나타냅니다.</li>
		<li>이제 명령을 통해 고유한 구조를 식별합니다. similarityAnalysis.py uf rotConstsData_C uf가 파일 이름 지정 레이블로 사용되는 경우 이제 파일 고유 구조-uf.data에저장된 이론의 PW91/6-311++G** 수준에서 최적화된 고유한 구성 목록이 있습니다.</li>
	</ol>
</li>
	<li>열역학 적 교정을 계산하는 데 필요한 정보의 최종 추출을 수행합니다. 이 정보를 사용하여 열역학 교정을 계산합니다.<ol>
<li>최종 전자 에너지, 회전 상수 및 진동 주파수를 추출하고 이를 사용하여 명령을 사용하여 열역학 보정을 계산합니다. 런 열-pw91.csh 고유구조-uf.data</li>
		<li>'gly-h2o-n.xlsx'라는이름의 Excel 스프레드시트의'Raw_Energies'시트에 명령줄 출력을 복사/붙여넣습니다. 단량체 (글리신 및 물)뿐만 아니라 각 수화물의 가장 낮은 에너지 부재 (gly-h2o-n, n = 1,2, ...)에 대해이 작업을 수행해야합니다.</li>
		<li>'gly-h2o-n.xlsx'스프레드시트의 첫 번째 시트에 원시 에너지가 추가되면 후속'Binding_Energies'및'Hydrate_Distribution'시트가 자동으로 업데이트됩니다. 특히,'Hydrate_Distribution'시트는 서로 다른 온도(예: 298.15K), 상대 습도(20%, 50%, 100%)에서 수화물의 평형 농도를 산출합니다. 및 물 ([H2O]) 및 글리신 ([글리신])의 초기 농도. 이러한 계산뒤에 있는 이론은 다음 단계에서 설명합니다.</li>
	</ol>
</li>
</ol>6. 해수면에서 실온에서글리(H2O)n=0-5 클러스터의 대기 농도 계산
참고: 이전 단계에서 생성된 열역학 데이터를 스프레드시트에 먼저 복사하고 순차수화의 깁스 프리 에너지를 계산하여 이 작업을 수행합니다. 그런 다음 깁스 프리 에너지는 각 순차 수화에 대한 평형 상수를 계산하는 데 사용됩니다. 마지막으로, 선형 방정식 세트는 단량체, 온도 및 압력의 주어진 농도에 대한 수화물의 평형 농도를 얻기 위해 해결된다.
<ol><li>아래와 같이 글리신의 순차적 수화를 위한 화학적 평형 시스템을 설정하여 시작하십시오. <img alt="Equation 1" src="/files/ftp_upload/60964/60964equ01.jpg">
</li>
	<li>K n =e-ΔG n/(kBT)를사용하여 평형 상수 K n을 계산하고, 여기서 n은 수화 수준이고, ΔGn은 nth 수화 반응의 깁스 자유 에너지 변화, kB는 Boltzmann의 상수이고, T는 온도이다. e <img alt="Equation 2" src="/files/ftp_upload/60964/60964equ02.jpg">
</li>
	<li>수화 및 무수화 글리신 클러스터의 평형 농도의 합이 단리글리신[Gly]0의 초기 농도와 같다는 가정을 이용하여 질량 보존을 위한0방정식을 설정한다. 평형 상수 식의 일부 대수 재배열을 사용하여 6 개의 동시 방정식의이 시스템을 다시 작성하십시오. <img alt="Equation 3" src="/files/ftp_upload/60964/60964equ03.jpg">
</li>
	<li>상기에 나타난 방정식의 시스템을 해결하여 글리(H2O)n=0-5의 평형2농도를 얻어 서 대기 중글리신의 농도에 대한 실험값56,57,,58을 이용하여, [Gly]0 = 2.9 x 106 cm-3,100% 상대 습도 및 298.15 K59의온도에서 대기 중의 물 농도 , [H2O] = 7.7 x 1017 cm-3.,</li>
</ol>

이 프로토콜에 의해 생성된 데이터의 정확도는 주로 세 가지에 따라 달라집니다: (i) 단계 2에 의해 샘플링 된 구성의 다양성, (ii) 시스템의 전자 구조의 정확성, (iii) 열역학 적 교정의 정확성. 포함된 스크립트를 편집하여 메서드를 수정하여 이러한 각 요소를 해결할 수 있습니다. 첫 번째 요소는 임의로 생성된 구조의 더 큰 초기 풀, GA의 더 많은 반복 및 GA에 관련된 기준의 느슨한 정의를 사용하여 쉽게 극복할 수 있습니다. 또한, 상이한 물리적 설명의 효과를 탐구하기 위해 자기 일관된 전하 밀도-기능적 타이트 바인딩(SCC-DFTB)모델 및 유효 단편 전위(EFP)63 모델과 같은 상이한 반경험적 방법을 사용할 수 있다. 여기서 주요 제한 사항은 공유 결합을 형성하거나 끊을 수있는 방법의 무능력이며, 이는 단량체가 동결된다는 것을 의미합니다. GA 절차는 반 경험적 설명에 따라 이러한 냉동 단량체의 가장 안정적인 상대 적 위치를 찾습니다.
시스템의 전자 구조의 정확도는 다양한 방법으로 향상 될 수있다, 그 계산 비용 각각. 하나는 M06-2X64 및 wB97X-V65,또는 이러한 Møller-Plesset66,,67,,68(MPn) 섭동 이론 및 결합 클러스터69(CC) 방법과 같은 더 나은 밀도 기능을 선택할 수 있다.68 기능 계층구조에서는 PW91과 같은 일반화된 그라데이션 근사치(GGA) 기능에서 wB97X-D 및 M06-2X와 같은 메타-GGA 하이브리드 기능과 같은 범위 분리하이브리드 기능으로 전환시 성능이 향상됩니다.
DFT 방법의 단점은 정확한 값을 향한 체계적인 수렴이 불가능하다는 것입니다. 그러나 DFT 방법은 계산 비용이 저렴하며 다양한 응용 프로그램에 대한 다양한 기능이 있습니다.
MP2 및 CCSD(T)와 같은 파형 함수 방법을 사용하여 계산된 에너지는 증가하는 추기경 수의 상관 관계 일관된 기준 집합과 함께([8-]cc-pV[D,T,Q,...] Z) 전체 기준 설정 한도를 체계적으로 수렴하지만 시스템 크기가 증가함에 따라 각 계산의 계산 비용이 엄중해집니다. 전자 구조의 추가 개선은 명시적으로 상관된 기준세트(70)를 사용하고 전체 기초 세트(CBS)한계로 추정함으로써 달성될 수 있다. 우리의 최근 연구는 밀도 장착 명시적으로 상관 2 차 Møller-Plesset (DF-MP2-F12) 교란 접근 MP2 / CBS계산32에접근 에너지를 산출 제안한다. 다른 전자 구조 방법을 사용하는 현재 프로토콜의 수정에는 두 단계가 포함됩니다: (i) 소프트웨어에 의해 주어진 구문에 따라 템플릿 입력 파일을 준비하고, (ii) run-pw91-sb.csh, run-pw91-lb.csh,및 run-pw91-lb-ultrafine.csh 스크립트를 편집하여 올바른 입력 파일 구문뿐만 아니라 올바른 스크립트를 생성합니다.
마지막으로, 열역학 적 교정의 정확도는 전자 구조 방법뿐만 아니라 전 세계적으로 최소 한도에 대한 PES의 설명에 따라 달라집니다. PES에 대한 정확한 설명은 매우 비용이 많이 드는 작업인 분력 필드72,,73(QFF)과 같은 핵 자유도의 변위와 관련하여 PES의 3차 및 고차 유도체의 계산을 요구합니다. 현재 프로토콜은 진동 주파수에 대한 고조파 발진기 근사치를 사용하므로 PES의 최대 두 번째 유도체만 계산해야 합니다. 이 접근법은 실제 PES및 고조파 PES의 큰 차이로 인해 매우 플로피 분자 및 대칭 이중 웰 전위와 같은 높은 조화성을 가진 시스템에서 문제가됩니다. 더욱이, 계산적으로 까다로운 전자 구조 방법으로부터 고품질 PES를 갖는 비용은 진동 주파수 계산에 대한 비용 의 문제를 복합화할 뿐이다. 이를 극복하기 위한 한 가지 방법은 고품질 의 전자 구조 계산과 낮은 품질의 PES에서 계산된 진동 주파수를 사용하여 비용과 정확도 간의 균형을 맞추는 것입니다. 현재 프로토콜은 이전 단락에 설명된 바와 같이 상이한 PES 설명을 사용하도록 수정될 수 있다; 그러나 스크립트 및 템플릿의 진동 주파수 키워드를 편집하여 조화진동 주파수를 계산할 수도 있습니다.
모든 구성 샘플링 프로토콜에 대한 두 가지 중요한 문제는 잠재적인 에너지 표면을 샘플링하는 초기 방법과 각 클러스터를 식별하는 데 사용되는 기준입니다. 우리는 전작에서 다양한 방법을 광범위하게 사용했습니다. 첫 번째 문제, 잠재적 인 에너지 표면을 샘플링하기위한 초기 방법, 우리는 이러한 요인에 따라 반 경험적 방법으로 GA를 사용하는 선택을했다. 화학적 직관26,무작위 샘플링 및 분자 역학(MD)29,,30을이용한 구성 샘플링은 물 클러스터18에대한 우리의 연구에서 관찰한 바와 같이 10단량 이상의 클러스터에 대해 정기적으로 가시적 글로벌 미니마를 찾지 못한다. 우리는 성공적으로(H2O)1174의복잡한 PES를 연구하기 위해 분지 호핑 (BH)을 사용했지만, BH 알고리즘이 찾지 못한 잠재적 인 저에너지 이섬의 수동 포함이 필요했습니다. 물 클러스터의 글로벌 최소 를 찾는 BH 및 GA의 성능의 비교,(H2O)n=10-20 GA가 지속적으로 BH75보다더 빠른 글로벌 최소를 발견것을 입증했다. OGOLEM 및 CLUSTER에서 구현된 GA는 모든 분자 클러스터에 적용할 수 있으며 클래식포스 필드, 반경험적, 밀도 기능성 및 ab initio 기능을 갖춘 방대한 수의 패키지와 인터페이스할 수 있기 때문에 매우 다재다능합니다. PM7의 선택은 속도와 합리적인 정확도에 의해 구동됩니다. 사실상 다른 모든 반 경험적 방법은 계산 비용이 상당히 높을 것입니다.
두 번째 문제에 관해서는, 우리는 전자 에너지, 다이폴 모멘트, 중첩 RMSD 및 회전 상수에 이르기까지 독특한 구조를 식별하기 위해 다른 기준을 사용하여 탐구했다. 다이폴 모멘트 를 사용하는 것은 두 다이폴 모멘트 구성 요소가 분자의 배향에 의존하고 총 다이폴 모멘트는 구조가 동일하거나 고유하다고 결정하는 임계값을 설정하기 어려웠던 방식으로 기하학적 차이에 매우 민감했기 때문에 어려웠습니다. 전자 에너지와 회전 상수의 조합이 가장 유용합니다.
두 구조물을 고유하게 간주하기 위한 현재 기준은 0.10 kcal mol-1의 에너지 차이 임계값과 1%의 회전 상수 차이를 기반으로 합니다. 따라서 에너지가 0.10 kcal mol-1 (~0.00015 a.u.) 이상 차이가 있는 경우 두 구조가 다른 것으로 간주됩니다. 그리고 그들의 세 회전 상수 (A, B, C)는 1 % 이상 다릅니다. 수년에 걸쳐 상당한 내부 벤치 마크는 이러한 임계 값이 합리적인 선택으로 나타났습니다. 당사의 구성 샘플링 접근법 및 스크리닝 방법론은 암모니아 및아민(32)을함유하는 강하게 결합된 삼차 황산염 수화물뿐만 아니라 물76,,77로 복잡한 다방향족 탄화수소와 같은 매우 약하게 결합된 클러스터에 적용되었다. 고려해야 할 다른 프로토네이션 상태가 있는 클러스터의 경우, 가장 좋은 방법은 다양한 PROTONATION 상태의 단량체로 시작하여 다양한 GA 계산을 실행하는 것입니다. 이렇게 하면 다른 프로토네이션 상태를 가진 구조물이 신중하게 고려됩니다. 그러나 낮은 수준의 DFT 계산을 통해 형상 최적화 과정에서 프로토네이션 상태가 변경될 수 있으며, 따라서 시작 형상에 관계없이 가장 안정적인 프로토네이션 상태를 생성할 수 있습니다.
GA 구성 샘플링 방법은 GA 코드가 GA 실행 과정에서 단량체가 다른 구성을 채택할 수 있도록 하는 일반적인 비매개 변수화 방법과 인터페이스되어 있는 한 플로피 분자에서도 잘 작동해야 합니다. 예를 들어, PM7과 GA를 인터페이싱하면 단량체의 구조가 변경될 수 있지만, 프로토네이션 상태가 변경될 때 와 같이 채권이 파손되면 구조물이 허용할 수 없는 후보로 버려질 수 있습니다.
우리는 고조파 근사치의 단점을 수정하는 다른 방법을 고려했다, 특히 낮은 진동 주파수에서 발생하는 사람들. 준 고조파 근사치를 현재 방법론에 통합하는 것은 어렵지 않습니다. 그러나, 준 고조파 방법에 대한 질문이 여전히있다, 그것은 적용 될 아래의 컷오프 주파수에 관해서 특히. 또한, 기존의 지혜가 RRHO 근사치에 비해 개선되어야한다고 제안하더라도 준 RRHO 근사치의 신뢰성을 검토하는 엄격한 벤치마킹 작업은 없습니다.
이렇게 제시된 프로토콜은 비협조적으로 결합된 가스 상 분자 클러스터의 임의의 시스템에 일반화될 수 있다. 또한 스크립트 및 템플릿을 편집함으로써 임의의 반경험적 방법, 전자 구조 방법 및 소프트웨어, 진동 분석 방법 및 소프트웨어를 사용하는 것이 일반화될 수 있다. 이는 사용자가 Linux 명령줄 인터페이스, Python 스크립팅 및 고성능 컴퓨팅에 익숙하다고 가정합니다. 익숙하지 않은 구문과 Linux 운영 체제의 모양과 스크립팅 경험 부족은 이 프로토콜의 가장 큰 함정이며 새로운 학생들이 가장 어려움을 겪는 곳입니다. 이 프로토콜은 주로 에어로졸 형성에 황산과 암모니아의 효과에 초점을 맞추고, 우리 그룹에서 년 동안 구현의 다양한 성공적으로 사용되어왔다. 이 프로토콜의 추가 개선은 더 많은 전자 구조 소프트웨어, 유전 알고리즘의 대체 구현 및 전자 및 진동 에너지의 빠른 계산을 위한 최신 방법의 사용에 대한 보다 강력한 인터페이스를 포함합니다. 이 프로토콜의 현재 응용 프로그램은 현재 대기에서 에어로졸 형성의 초기 단계와 프리바이오틱 환경에서 더 큰 생물학적 분자의 형성에서 아미노산의 중요성을 탐구하고 있습니다.

기후 변화의 대기 연구에서 가장 불확실한 매개 변수는 구름 입자가 들어오는 태양 복사를 반영하는 정확한 정도입니다. 가스에 부유하는 미립자 물질인 에어로졸은 들어오는 방사선을 산란시키는 구름 응축 핵(CCN)이라고 불리는 구름 입자를 형성하여대기의흡수와 후속 가열을 방지합니다 1. 이 순 냉각 효과에 대한 자세한 이해를 위해서는 에어로졸이 CCN으로 증가하는 것을 이해해야 하며, 이로 인해 작은 분자 클러스터가 에어로졸 입자로 증가하는 것을 이해해야 합니다. 최근 연구는 에어로졸 형성이 직경 3 nm 이하의 분자 클러스터에 의해 시작된다는 것을 제안했다2; 그러나,이 크기 정권은 실험 기술3,,4를사용하여 액세스하기 가 어렵습니다. 따라서 이러한 실험적 한계를 극복하기 위해서는 전산 모델링 접근법이 요구된다.
아래에 설명된 모델링 방법을 사용하여 수화 클러스터의 성장을 분석할 수 있습니다. 우리는 생물학적 전 환경에서 작은 성분에서 큰 생물학적 분자의 형성에 물의 역할에 관심이 있기 때문에, 우리는 글리신으로 우리의 접근 방식을 설명합니다. 이러한 연구 질문을 해결하는 데 필요한 과제와 도구는 대기 에어로졸 및 사전 핵 화 클러스터,,5,6,,7,,8,,9,,10,11,,12,13,,14,,15의연구에 관여하는 것과 매우 유사하다. 여기에서는, 우리는 격리된 글리신 분자에서 시작하는 수화한 글리신 클러스터를 검토하고 최대 5개의 물 분자의 일련의 단계별 추가에 선행됩니다. 최종 목표는 해수면에서 실온에서 대기 중의 글리(H2O)n=0-5 클러스터의 평형 농도와 100%의 상대 습도(RH)를 계산하는 것입니다.
이 sub 나노미터 분자 클러스터의 소수는 그밖 증기 분자를 추가하거나 기존 클러스터에 응고하여 메타안정 임계 클러스터 (직경 1-3 nm)로 증가합니다. 이러한 중요한 클러스터는 구름의 침전 효율과 입사광을 반사하는 능력에 직접적인 영향을 미치는 훨씬 더 큰 구름 응축 핵(CCN)의 형성으로 이어지는 유리한 성장 프로필을 가지고 있습니다. 따라서 분자 클러스터의 열역학과 평형 분포를 잘 이해하면 에어로졸이 지구 기후에 미치는 영향을 보다 정확하게 예측할 수 있습니다.
에어로졸 형성의 설명 모델은 분자 클러스터 형성의 정확한 열역학을 필요로합니다. 분자 클러스터 형성의 정확한 열역학의 계산은 클러스터의 잠재적 에너지 표면(PES)16에서글로벌 및 로컬 최소를 찾는 것을 포함하는 가장 안정적인 구성의 식별을 요구한다. 이러한 과정은 구성 샘플링이라고 하며 분자 역학(MD)17,18,,19,,20,몬테카를로(MC)21,,22,및 유전 알고리즘(GA)23,,24,,25를기초로 하는 것을 포함하는 다양한 기술을 통해 달성될 수 있다.,
이론의 높은 수준에서 대기 수화물의 구조와 열역학을 얻기 위해 수년에 걸쳐 다른 프로토콜이 개발되었습니다. 이러한 프로토콜은 (i) 구성 샘플링 방법, (ii) 구성 샘플링에 사용되는 로우 레벨 방법의 특성, 및 (iii) 후속 단계에서 결과를 구체화하는 데 사용되는 상위 레벨 방법의 계층구조의 선택에서 달랐다.
구성 샘플링 방법은 화학 직관26,무작위 샘플링27,,28,분자 역학 (MD)29,30,분지 호핑 (BH)31,및 유전 알고리즘 (GA)24,,25,,32를포함했다., 이러한 샘플링 방법에 사용되는 가장 일반적인 저수준 방법은 PM6, PM7 및 SCC-DFTB와 같은 힘 필드 또는 반 경험적 모델입니다. 이들은 종종 점점 더 큰 기초 세트와 야곱의 사다리(33)의높은 횡단에서 더 신뢰할 수있는 기능DFT 계산 뒤에 있습니다. 경우에 따라, 이들은 MP2, CCSD (T) 및 비용 효율적인 DLPNO-CCSD (T)34,,35와같은 더 높은 수준의 파형 방법에 선행된다.
Kildgaard 등36은 더 큰 클러스터에 대한 후보를 생성하기 위해 더 작은 수화 또는 비수화 클러스터 주위의 피보나치구체 37의 지점에서 물 분자가 추가되는 체계적인 방법을 개발했습니다. 비물리적 및 중복 후보는 가까운 접촉 임계값과 서로 다른 순응자 간의 근평균 제곱 거리에 따라 제거됩니다. PM6 반 경험적 방법과 DFT 및 웨이브 함수 방법의 계층 구조를 사용하여 후속 최적화는 높은 수준의 이론에서 낮은 에너지 적합성 세트를 얻는 데 사용됩니다.
인공 꿀벌 콜로니(ABC)알고리즘(38)은 최근 장 외가 ABCluster39라는프로그램에서 분자 클러스터를 연구하기 위해 구현한 새로운 구성 샘플링 접근법이다. Kubecka 등40은 구성 샘플링을 위해 ABCluster를 사용한 다음, 타이트결합 GFN-xTB 반경험적방법(41)을사용하여 낮은 수준의 재최적화를 수행하였다. 그들은 DFT 방법을 사용하여 구조와 에너지를 더욱 정제한 다음 DLPNO-CCSD(T)를 사용하여 최종 에너지를 추가했습니다.
방법에 관계없이 구성 샘플링은 PES에서 임의로 또는 임의로 생성된 포인트 분포로 시작합니다. 각 점은 해당 분자 클러스터의 특정 형상에 해당하며 샘플링 방법에 의해 생성됩니다. 그런 다음 PES의 "내리막" 방향을 따라 각 점에 대해 가장 가까운 로컬 최소를 찾습니다. 이렇게 발견된 최소의 집합은 적어도 얼마 동안 안정되어 있는 분자 클러스터의 그 기하학에 해당합니다. 여기서, PES의 형상과 표면의 각 지점에서의 에너지의 평가는 보다 정확한 물리적 설명이 더 많은 계산비용이 드는 에너지 계산을 초래하는 시스템의 물리적 설명에 민감할 것이다. 우리는 특히 OGOLEM25 프로그램에서 구현된 GA 방법을 사용할 것이며, 이는 다양한 글로벌 최적화 및 구성 샘플링 문제42,,43,,44,,45,초기 샘플링 포인트 세트를 생성하기 위해 성공적으로 적용되었다. PES는 MOPAC2016 프로그램 47에 구현된 PM7모델(46)에의해 설명될 것이다.46 이 조합은 MD 및 MC 방법에 비해 더 다양한 포인트를 생성하고 PES에 대한 보다 상세한 설명보다 더 빠르게 로컬 미니마를 발견하기 때문에 사용된다.
GA에 최적화된 로컬 미니마 세트는 일련의 스크리닝 단계의 시작 형상으로 사용되며, 이로 인해 최소 에너지가 낮아지지 않습니다. 프로토콜의 이 부분은 작은 기초 세트로 밀도 기능 이론(DFT)을 사용하여 고유한 GA 에 최적화된 구조 세트를 최적화하는 것으로 시작됩니다. 이 최적화 세트는 일반적으로 GA에 최적화된 반 경험적 구조에 비해 더 자세하게 모델링되는 더 작은 고유한 로컬 최소 구조 세트를 제공합니다. 그런 다음 더 큰 기초 집합을 사용하여 이 작은 구조 집합에서 또 다른 DFT 최적화가 수행됩니다. 다시 말하지만, 이 단계는 일반적으로 작은 기초 DFT 단계에 비해 더 자세하게 모델링되는 더 작은 구조세트를 제공할 것이다. 그런 다음 최종 고유 구조 세트가 더 엄격한 수렴에 최적화되고 고조파 진동 주파수가 계산됩니다. 이 단계 후 우리는 대기 중 클러스터의 평형 농도를 계산하는 데 필요한 모든 것을 갖습니다. 전체 적인 접근 방식은 그림 1에서다이어그램으로 요약됩니다. 우리는 PW9148 일반화 그라데이션 근사치 (GGA) 교환 상관 관계 기능 Pople49 50 기초 세트의 두 가지 변형과 함께 DFT의 구현 (작은 기초 단계에 대한 6-311 +G*와 큰 기초 단계에 대한 6-311 +G**)를 사용합니다. 이러한 교환 상관 기능 및 기초 세트의 이러한 특별한 조합은 대기클러스터(51,,52)에대한 정확한 깁스 프리 에너지 형성의 이전 성공으로 인해 선택되었다.
이 프로토콜은 사용자가 휴대용 배치시스템(53)을 사용하여 고성능 컴퓨팅(HPC) 클러스터에 액세스할 수 있다고 가정하고, MOPAC2016(http://openmopac.net/MOPAC2016.html)47,OGOLEM(https://www.ogolem.org)25,가우시안 09(https://gaussian.com)49,및 OpenBabel54(http://openbabel.org/wiki/Main_Page) 소프트웨어가 특정 설치 지시에 따라 설치된 소프트웨어이다. 이 프로토콜의 각 단계는 사용자의 $PATH 환경 변수에 포함된 디렉터리에 저장해야 하는 사내 셸 및 Python 2.7 스크립트 집합을 사용합니다. 위의 모든 프로그램을 실행하는 데 필요한 모든 환경 모듈 및 실행 권한도 사용자의 세션에 로드되어야 합니다. GA 코드(OGOLEM) 및 반경험적 코드(MOPAC)에 의한 디스크 및 메모리 사용량은 최신 컴퓨터 리소스 표준에 따라 매우 작습니다. OGOLEM/MOPAC의 전체 메모리 및 디스크 사용량은 사용하려는 스레드 수에 따라 달라지며, 그렇다 하더라도 대부분의 HPC 시스템의 기능에 비해 리소스 사용량이 작습니다. QM 메서드의 리소스 요구 는 클러스터의 크기와 사용된 이론 의 수준에 따라 달라집니다. 이 프로토콜을 사용하면 이론 의 수준이 달라져 서 에너지 구조의 최종 집합을 계산할 수 있으며, 일반적으로 빠른 계산으로 인해 결과의 정확도가 더 불확실해질 수 있다는 점입니다.
명확성을 위해 사용자의 로컬 컴퓨터는"로컬 컴퓨터"라고하며 액세스 권한이 있는 HPC 클러스터는"원격 클러스터"라고합니다.

<table border="0" cellpadding="0" cellspacing="0" style="width:1639px;" width="1638">
	
	<tbody>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">Avogadro</td>
			<td style="width:483px;">https://avogadro.cc</td>
			<td style="width:123px;"></td>
			<td style="width:617px;">Open-source molecular visualization program</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">Gaussian [09/16] Software</td>
			<td style="width:483px;">http://www.gaussian.com/</td>
			<td style="width:123px;"></td>
			<td>Commercial ab initio electronic structure program</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">MOPAC 2016</td>
			<td style="width:483px;">http://openmopac.net/MOPAC2016.html</td>
			<td style="width:123px;"></td>
			<td>Open-source semi-empirical program</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">OGOLEM Software</td>
			<td style="width:483px;">https://www.ogolem.org</td>
			<td style="width:123px;"></td>
			<td>Genetic algorithm-based global optimization program</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">OpenBabel</td>
			<td style="width:483px;">http://openbabel.org/wiki/Main_Page</td>
			<td style="width:123px;"></td>
			<td>Open-source cheminformatics library</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">calcRotConsts.py</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Python script to compute rotational constants</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">calcSymmetry.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to calculate symmetry number of a molecule given Cartesian coordinates</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">combine-GA.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to combine energy and rotational constants from different GA directories</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">combine-QM.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to combine energy and rotational constants from different QM directories</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">gaussianE.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to extract Gaussian 09 energies</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">gaussianFreqs.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to extract Gaussian 09 vibrational frequencies</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">getrotconsts</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Executable to calculate rotational constants given a molecule&#39;s Cartesian coordinates</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">getRotConsts-dft-lb.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to compute rotational constants for a batch of large basis DFT optimized structures</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">getRotConsts-dft-lb-ultrafine.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to compute rotational constants for a batch of ultrafine DFT optimized structures</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">getRotConsts-dft-sb.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to compute rotational constants for a batch of small basis DFT optimized structures</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">getRotConsts-GA.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to compute rotational constants for a batch of genetic algorithm optimized structures</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">global-minimum-coords.xyz</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Cartesian coordinates of global minimum structures of gly-(h2o)n, where n=0-5</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">make-thermo-gaussian.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to extract data from Gaussian output files and make input files for the thermo.pl script</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">ogolem-input-file.ogo</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Ogolem sample input file</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">ogolem-submit-script.pbs</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>PBS batch submission file for Ogolem calculations</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">README.docx</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Clarifications to help readers use the scripts effectively</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">runogolem.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to run OGOLEM</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">run-pw91-lb.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to run a batch of large basis DFT optimization calculations</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">run-pw91-lb-ultrafine.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to run a batch of ultrafine DFT optimization calculations</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">run-pw91-sb.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to run a batch of small basis DFT optimization calculations</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">run-thermo-pw91.csh</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Shell script to compute the thermodynamic corrections for a batch of DFT optimized structures</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">similarityAnalysis.py</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Python script to determine unique structures based on rotational constants and energies</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">symmetry</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Executable to calculate molecular symmetry given Cartesian coordinates</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">symmetry.c</td>
			<td style="width:483px;">(C) 1996, 2003 S. Patchkovskii, Serguei.Patchkovskii@sympatico.ca</td>
			<td style="width:123px;"></td>
			<td style="width:617px;">C code to determine the molecular symmstry of a molecule given Cartesian coordinates</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">template-marcy.pbs</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Template for a PBS submit script which uses OGOLEM</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">template-pw91.com</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Template Gaussian 09 input</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">template-pw91-HL.com</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Template Gaussian 09 input for ultrafine DFT optimization</td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;width:416px;">thermo.pl</td>
			<td style="width:483px;">https://www.nist.gov/mml/csd/chemical-informatics-research-group/products-and-services/program-computing-ideal-gas</td>
			<td style="width:123px;"></td>
			<td>Perl open-source script to compute ideal gas thermodynamic corrections</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">gly-h2o-n.xlsx</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Excel spreadsheet for the complete protocol</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">table-1.xlsx</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Excel spreadsheet</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">table-2.xlsx</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Excel spreadsheet</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">table-3.xlsx</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Excel spreadsheet</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">table-4.xlsx</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Excel spreadsheet</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">water.xyz</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Cartesian coordinates of water</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:416px;">glycine.xyz</td>
			<td style="width:483px;">Shields Group, Department of Chemistry, Furman University</td>
			<td style="width:123px;"></td>
			<td>Cartesian coordinates of glycine</td>
		</tr>
	</tbody>
</table>

computation of atmospheric concentrations of molecular clusters from ab initio thermochemistry

대기 에어로졸의 형성 과 성장에 대한 계산 연구는 가스 상 전자 구조 및 진동 주파수 계산에서 얻을 수있는 정확한 깁스 자유 에너지 표면을 필요로한다. 이러한 수량은 형상이 잠재적 에너지 표면의 최소값에 해당하는 대기 클러스터에 유효합니다. 최소 에너지 구조의 깁스 자유 에너지는 온도 및 압력과 같은 다양한 조건하에서 클러스터의 대기 농도를 예측하는 데 사용될 수 있다. 우리는 점점 더 정확한 스크리닝 계산의 시리즈 뒤에 유전 알고리즘 기반 구성 샘플링에 내장 된 계산 저렴한 절차를 제시한다. 이 절차는 반 경험적 모델을 사용하여 대규모 구성 세트의 형상을 생성하고 진화한 다음 일련의 높은 수준의 ab initio 수준의 이론에서 결과고유 구조를 구체화하는 것으로 시작합니다. 마지막으로, 열역학 보정은 최소 에너지 구조의 결과 집합에 대해 계산되고 깁스 의 형성, 평형 상수 및 대기 농도의 자유 에너지를 계산하는 데 사용됩니다. 우리는 주변 조건하에서 수화 글리신 클러스터의 연구에이 절차의 응용 프로그램을 제시한다.

Watch this Scientific Journal Video about Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry at JoVE.com

Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry

약하게 결합된 분자 클러스터의 대기 농도는 유전 알고리즘과 반 경험적 및 ab initio 양자 화학을 활용하는 다단계 구성 샘플링 방법을 통해 발견되는 저에너지 구조의 열화학적 특성으로부터 계산할 수 있습니다.

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 ...

당사의 프로토콜은 약하게 결합된 분자 클러스터를 연구하는 유연하고 계산적으로 실현 가능한 접근 방식을 제공하며 구조, 형성 및 풍부에 대한 통찰력을 얻기 위해 쉽게 적용할 수 있습니다. 이 기술의 주요 장점은 빠른 힘 필드와 반 경험적 방법에서 엄격한 양자 기계 방법에 이르기까지 다양한 수준의 이론에서 분자 클러스터를 치료하는 데 있어 효율성과 유연성입니다. 대기 및 에어로졸 화학은 기후 변화의 더 나은 모델로 이어지는이 접근 방식에서 가장 혜택을 누릴 수 있습니다.그러나 분자 클러스터와 관련된 모든 분야는 이 접근법을 악용할 수 있습니다. 이 기술을 수행한 적이 없는 개인의 경우 가장 어려운 단계는 초기 프로그램 및 스크립트 설치 및 로컬 컴퓨팅 환경에 대한 적응입니다. 계산 화학에 익숙하지 않은 학생들은 명시적 단계별 지침의 시각화를 통해 고성능 컴퓨팅 클러스터를 사용하는 가파른 학습 곡선을 극복할 수 있습니다.유전 알고리즘 구성 샘플링에 사용하기 위한 절연 글리신 분자의 최소 에너지 구조를 얻으려면, 아보가드로에서 새로운 세션을 열고 빌드, 삽입, 펩타이드, 글리신 및 펩타이드 삽입을 클릭하여 시각화 창에서 글리신 단조머를 생성한다. 광고 확장 및 가우시안을 클릭하고 표시된 대로 텍스트 상자의 첫 번째 줄을 편집합니다. 명령 파일을 생성하고 저장한 glycine.com 클릭합니다.격리된 물의 최소 에너지 구조를 얻으려면 Avogadro에서 새 세션을 열고 빌드, 삽입 및 조각을 선택합니다. 필터 텍스트 상자에 물을 입력하고, 물 파일을 선택하고 삽입을 클릭합니다. 광고 확장 및 가우시안을 클릭하고 표시된 대로 텍스트 상자의 첫 번째 줄을 편집합니다.명령 파일을 생성하고 저장한 water.com 클릭합니다. 그런 다음 두 com 파일을 컴퓨팅 클러스터로 전송하고 적절한 제출 스크립트를 사용하여 Gaussian 09 계산을 실행합니다. 계산이 컴퓨팅 클러스터에서 완료되면 열린 바벨을 호출하여 지시에 따라 명령을 입력하는 최소 에너지 구조의 xyz 파일을 생성합니다.유전자 알고리즘 기반 구성 샘플링의 경우 모든 스크립트와 템플릿을 폴더에 추가하고 폴더를 원격 클러스터에 복사합니다. 모든 스크립트가 실행 가능한지 확인하고 지시된 대로 명령을 사용하여 경로 환경 변수에 스크립트 디렉터리 위치를 추가합니다. 저렴한 반경험적 수준에서 글리신과 물에 대한 저에너지 구조를 얻으려면 n이 물 분자의 수인 글리-h2o-n이라는 디렉토리를 만들고 유전자 알고리즘 계산을 실행하기 위해 글리-h2o-n 디렉토리하에 GA라는 하위 디렉토리를 생성합니다.ogolem 입력 파일, 모노머 카르테시안 좌표 및 PBS 일괄 제출 스크립트를 GA 디렉토리에 복사하고 적절하게 수정된 실행을 사용하여 GA 계산을 실행합니다. pbs는 스크립트를 제출합니다. 계산이 완료되면 디렉터리를 글리-h2o-n GA pm7로 변경하고 13개가 클러스터의 원자 수이고 0과 9개가 인덱스가 0~9개인 10개의 구조가 있는 위치를 나타내는 대로 getRotConsts 명령을 실행합니다.이렇게 하면 GA 최적화된 클러스터의 회전 상수를 계산하고 모든 GA 최적화된 클러스터 구성, 에너지 및 회전 상수의 정렬 된 목록이 포함된 rotConstsData_C 라는 파일을 생성합니다. 유사성을 실행분석합니다. 고유한 GA 최적화 된 클러스터를 찾고 저장하는 입력으로 rotConstsData_C 파일이있는 py 스크립트.Pm7은 고유한구조-pm7.data라는 파일을 생성하는 파일 이름 지정 레이블로 사용됩니다. 여기에는 고유한 GA 최적화 된 구성의 정렬 된 목록이 포함되어 있습니다. 글리-h2o-n GA 디렉토리에서 결합-GA를 사용합니다.csh 스크립트는 여러 비교 가능한 GA 실행에 대한 결과를 결합하고 고유구조-pm7이라는 새로운 고유 구조 목록을 생성합니다. 글리-h2o-n GA 디렉토리의 데이터입니다. 작업 디렉터리에는 그림과 같이 정확한 조직과 구조가 있어야 합니다.반실적 방법을 기반으로 한 유전자 알고리즘으로부터 글리신 수클러스터의 구조를 보다 정확한 양자 기계적 방법을 사용하여 하나에, 글리-h2o-n 디렉토리 하에서 QM이라는 서브디렉토리를 생성한다. QM 디렉토리에서 pw91-sb라는 또 다른 하위 디렉터리를 만들고 gly-h2o-n GA 디렉토리에서 QM pw91-sb 디렉터리에 고유한Structures 목록을 복사합니다. 디렉토리를 글리-h2o-n QM pw91-sb로 변경하고 작은 기초 세트 밀도 기능 이론 스크립트 런-pw91-sb를 실행합니다.sb가 이 계산 집합에 대한 레이블인 csh, Q는 컴퓨팅 클러스터에서 기본 큐이며 10은 10개의 계산이 하나의 일괄 처리 작업으로 그룹화될 것임을 나타냅니다. 제출된 계산이 완료되면 getRotConsts-dft-sb를 사용합니다. 에너지를 추출하고 작은 기초 최적화 된 클러스터의 회전 상수를 계산하는 csh 스크립트.여기서, pw91은 사용되는 밀도 기능이고 n은 클러스터내원자의 수이다. 유사한 분석을 사용합니다. py 스크립트는 이전과 같이 고유 구조를 식별하지만 SB를 레이블로 사용합니다.pw91 631에 최적화된 고유한 구성 목록과 G 별 수준의 이론이 고유한Structures-sb에 저장됩니다. 데이터 파일입니다. 글리-h2o-n QM 디렉토리에서 결합된 결합 QM을 사용합니다.여러 비교 QM 실행의 결과를 결합하는 csh 스크립트. 결합 QM. csh pw91-sb 명령은 고유구조-sb라는 새로운 고유 구조 목록을 생성합니다.글리-h2o-n QM 디렉토리의 데이터입니다. 더 나은 양자 기계적 설명을 사용하여 글리신 및 물 클러스터의 구조를 더욱 미세 조정하려면 QM 디렉토리 아래에 pw91-lb라는 하위 디렉토리를 만듭니다. QM pw91-sb 디렉토리에서 QM pw91-lb 디렉토리에 고유한 구조 목록을 복사하고 디렉터리를 QM pw91-lb로 변경합니다.큰 기초 밀도 기능 이론 스크립트 런-pw91-lb를 실행합니다. lb가 이 계산 집합에 대한 레이블인 csh, Q는 컴퓨팅 클러스터에서 기본 큐이며 10은 10개의 계산을 하나의 일괄 처리 작업으로 그룹화해야 한다는 것을 나타냅니다. 제출된 계산이 완료되면 getRotConsts-dft-lb를 사용하십시오.csh 명령은 대규모 최적화 된 클러스터의 회전 상수를 계산합니다. 여기서, pw91은 사용되는 밀도 기능이고 n은 클러스터내원자의 수이다. 유사성을 사용합니다Analysis.pw91 6311 플러스 G 스타 이론 수준에 최적화 된 고유 구성의 목록을 생성하고 고유구조 - 파운드에 저장하기 위해 레이블로 지금과 같은 py 스크립트. 데이터 파일입니다. 원하는 열화학 적 교정을 계산하는 데 필요한 글리신 및 물 클러스터의 진동 구조 및 에너지를 얻으려면 QM pw91-lb 디렉토리에서 QM pw91-lb 초미세 디렉토리에 고유 구조 목록을 복사하고 디렉토리를 QM / pw91-lb 울트라 파인으로 변경하십시오.초미세 밀도 기능 이론 스크립트 런-pw91-lb-울트라파인을 실행합니다. uf가 이 계산 집합에 대한 레이블인 csh, Q는 컴퓨팅 클러스터에서 기본 큐이며 10은 10개의 계산을 하나의 일괄 처리 작업으로 그룹화해야 한다는 것을 나타냅니다. 이 스크립트는 Gaussian 09에 대한 입력을 자동으로 생성하고 모든 계산을 제출합니다.제출된 계산이 완료되면 getRotConsts-dft-lb-ultrafine을 사용하십시오. 초미세 최적화 된 클러스터의 회전 상수를 계산하는 csh 명령. 여기서, pw91은 사용되는 밀도 기능이고 n은 클러스터내원자의 수이다.유사성을 사용합니다Analysis. pw91 6311 플러스 고유구조-uf에서 초미세 수렴 기준에 최적화된 고유한 구성 목록을 생성하고 저장하는 레이블로 uf와 함께 지금과 같은 파이 스크립트. 데이터 파일입니다.그런 다음 실행 열-pw91을 실행합니다. 독특한 스크립트구조-uf. 열역학 보정을 계산하는 입력으로 데이터 파일.명령줄 출력을 글리-h2o-n.xls라는 연결된 스프레드시트에 복사하여 붙여넣습니다. 이 계산및 후속 n의 원시 에너지가 글리-h2o-n의 첫 번째 시트에 2, 3, 4 및 5개의 계산과 같기 때문에. xls spredsheet, 다른 온도에서 수화물의 평형 농도를 산출 수화물 분포 시트, 상대 습도, 물과 글리신의 초기 농도가 업데이트됩니다.여기서, 글리신-물 클러스터의 가장 낮은 전자 에너지 이소마를 관찰할 수 있다. 수소 결합 네트워크는 물 분자의 수가 주로 평면 네트워크에서 n에서 3 차원 케이지와 같은 구조로 이동 증가함에 따라 복잡성이 증가하는 방법을 주목하십시오. 이 표에서, 런 써모-pw91의 출력의 예.csh 스크립트가 표시됩니다. 각 클러스터에 대해 pw91 6311 플러스 G 스타의 에너지는 pw91 6311 플러스 G 스타 등급의 가스 상 전자 에너지와 하나무 단위로 초미세 통합 그리드에 계산된 G 스타 등급의 이론과 두더지 당 킬로 칼로리 단위로 0포인트 진동 에너지에 해당한다. 각 온도에서 깁스없는 에너지 형성 델타 G에서 계산된 엔탈피 형성 델타 H는 두더지 당 킬로칼로리 단위로 주어지며 계산된 엔탈피 형성 S는 두더지 당 칼로리 단위로 제공됩니다.이 표에서는 수분 공급 및 순차 수화의 총 깁스 없는 에너지 변화의 대표적인 계산이 표시됩니다. 이러한 데이터를 이용하여, 수화 글리신의 대기 농도를 계산할 수 있다. 올바른 소프트웨어를 설치하고 포함된 스크립트를 추가하여 자체 컴퓨팅 환경을 반영해야 합니다.스크립트의 위치를 경로에 추가하는 것이 중요합니다. 이 기술은 펩티드 결합 형성을 향한 대기수 클러스터의 촉매 활성을 결정하는 데 사용되어 프리바이오틱 화학 분야에 기여했다.

Education

Research

JoVE Core

JoVE Journal

Chemistry

Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry

Thermochemistry

Chemical Equilibrium

SCIENTISTS IN ACTION

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 ...

연구

화학

A Simple Method for the Size Controlled Synthesis of Stable Oligomeric Clusters of Gold Nanoparticles under Ambient Conditions

주변 조건에 따라 금 나노 입자의 안정 올리고머 클러스터의 크기 제어 합성을위한 간단한 방법

Reducing dilute aqueous HAuCl4 with sodium thiocyanate (NaSCN) under alkaline conditions produces 2 to 3 nm diameter nanoparticles. Stable grape-like ...

Construction and Systematical Symmetric Studies of a Series of Supramolecular Clusters with Binary or Ternary Ammonium Triphenylacetates

건설 및 바이너리 또는 삼항 암모늄 Triphenylacetates와 초분자 클러스터의 시리즈의 체계적 대칭 연구

Functions of clusters in nano or sub-nano scale significantly depend on not only kinds of their components but also arrangements, or symmetry, of their ...

Improved Heterojunction Quality in Cu2O-based Solar Cells Through the Optimization of Atmospheric Pressure Spatial Atomic Layer Deposited Zn1-xMgxO

Improved Heterojunction Quality in Cu2O-based Solar Cells Through the Optimization of Atmospheric Pressure Spatial Atomic Layer Deposited Zn1-xMgxO

구리에서 개선 된 이종 접합 품질<sub&gt; 2</sub&gt; 입금 대기 압력 공간 원자 층의 최적화를 통해 태양 전지를 O 기반<br /&gt; 아연<sub&gt; 1-X</sub&gt; 마그네슘<sub&gt; X</sub&gt; O

Atmospheric pressure spatial atomic layer deposition (AP-SALD) was used to deposit n-type ZnO and Zn1-xMgxO thin films onto p-type thermally oxidized Cu2O ...

Analyzing the Photo-oxidation of 2-propanol at Indoor Air Level Concentrations Using Field Asymmetric Ion Mobility Spectrometry

2-프로 판 올 필드 비대칭 이온 기동성 분 광 분석을 사용 하 여 실내 공기 수준 농도에서 사진 산화 분석

We demonstrate a versatile protocol to be used for determining the effectiveness of photocatalysts in degrading indoor air concentration (ppb) volatile ...

Direct Imaging of Laser-driven Ultrafast Molecular Rotation

레이저 기반 초고속 분자 회전의 직접 이미징

We present a method for visualizing laser-induced, ultrafast molecular rotational wave packet dynamics. We have developed a new 2-dimensional Coulomb ...

Fizzy Extraction of Volatile Organic Compounds Combined with Atmospheric Pressure Chemical Ionization Quadrupole Mass Spectrometry

대기압 화학 이온화와 결합 된 휘발성 유기 화합물의 탄산염 추출 4 극자 질량 분광법

Chemical analysis of volatile and semivolatile compounds dissolved in liquid samples can be challenging. The dissolved components need to be brought to ...

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 ...

Light-driven Molecular Motors on Surfaces for Single Molecular Imaging

빛 주도 분자 모터 표면에 단일 분자 이미징

The design and synthesis of a synthetic system that aims for the direct visualization of a synthetic rotary motor at the single molecule level on surfaces ...

Synthesis of In37P20(O2CR)51 Clusters and Their Conversion to InP Quantum Dots

Synthesis of In37P20O2CR51 Clusters and Their Conversion to InP Quantum Dots

37P 20(o2CR)51 클러스터의 합성 및 InP 양자 점으로의 변환

This text presents a method for the synthesis of In37P20(O2C14H27)51 clusters and their conversion to indium phosphide quantum dots. The In37P20(O2CR)51 ...