Name: multiscale sampling of a heterogeneous water/metal catalyst interface using density functional theory and force-field molecular dynamics
Uploaded: 2019-04-12T14:45:00
Duration: 10 min 52 s
Description: Watch this Scientific Journal Video about Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics at JoVE.com

Denne forskningen ble finansiert av National Science Foundation gjennom prisen nummer CBET-1438325. Fellesskap støtte for CJB gjennom NASA trening Grant NX14AN43H er takknemlig anerkjent. Simuleringer ble utført på Palmetto superdatamaskin klyngen, som vedlikeholdes av gruppen Cyberinfrastructure teknologi ved Clemson University. Vi takker Dr. Paul J. Meza-Morales for testing protokollen.

1. Generer adsorbate strukturen
<ol><li>Opprette en VASP POSCAR fil som består av en supercelle med periodiske betingelser, som du ville gjort hvis du var å utføre simuleringer av adsorbates på metallflater under vakuum. Supercelle bør inneholde en første gjette i adsorbate strukturen og metalloverflate samt vakuum plass over adsorbate for å legge til H2O molekyler. Informasjon om hvordan du gjør dette er gitt i tidligere arbeid35,44. Merk: Det er viktig at høyden av vakuum plass være minst 12 Å over toppen av adsorbate.</li>
	<li>Slapp av strukturen og minimere sin energi ved hjelp av VASP-koden. Informasjon om hvordan du gjør dette er gitt i tidligere arbeid35,44. Dette vil gi en fil kalt CONTCAR, som vil bli brukt i neste avsnitt.</li>
</ol>2. Legg eksplisitt H2O molekyler
<ol><li>Legger til N vannmolekyler vakuum plass i CONTCAR opprettet i trinn 1.2 bruker MCPliQ koden, der N = ρV, ρ er tettheten av vann og V er volumet av vakuum plass over adsorbate. Merk: ρ bør tas som tetthet av vann som TIP3P/CHARMM vann modellen ved simulering temperatur. V vil bli forbedret i neste trinn.<ol>
<li>Angi følgende informasjon i filen MCPliQ master_input.txt: antall H2O molekyler til (N) ved å endre det første argumentet til linje 28, banen til filen water.txt ved å endre det andre argumentet til linje 28 og minimum og Maksimal høyde supercelle som kan være opptatt av vannmolekyler ved å endre den Minimum z-koordinat i tråd 11 og den høyeste z-koordinaten på linje 12.</li>
		<li>Utføre MCPliQ koden ved å skrive inn mcpliq fra kommandolinjegrensesnittet å sette vann inn i filen CONTCAR. Koden vil produksjonen med filtypen for én eller flere filer. POSCAR. Merk: Hvis flere POSCAR filer er produsert, de navnet POSCAR_n.POSCAR. Velg filen der n er største.</li>
	</ol>
</li>
	<li>Generere LAMMPS inndatafiler for en NPT simulering og equilibrate celle volumet ved hjelp av FFMD i NPT ensemblet i LAMMPS.<ol>
<li>Kjøre skript lmps_bond_angle.py på den. POSCAR fil generert i trinn 2.1.2 ved å skrive lmps_bond_angle.py $filename. POSCAR på kommandolinjegrensesnitt, der $filename er navnet på den. POSCAR filen genereres i trinn 2.1.2. Dette skriptet oppretter en fil kalt $filename. POSCAR.bond_angle_info.txt som viser obligasjoner og vinkler som skal brukes i datafilen LAMMPS.</li>
		<li>Åpne VMD og velg fil > Ny molekyl i hovedvinduet åpne Molecule Filutforsker -vinduet. Velg VASP_POSCAR fra rullegardinmenyen Avgjøre filtype . Klikk Bla gjennom og navigerer til $filename. POSCAR fil. Velg Load for å åpne $filename. POSCAR fil.</li>
		<li>Åpne Tk konsollen innen VMD ved å velge utvidelser > Tk konsollen VMD HOVEDVINDU.</li>
		<li>Kjør følgende kommando i konsollen Tk: topo writelammpsdata $WDPATH/data.myadsorbate full, der $WDPATH er mappen på datamaskinen der VMD vil skrive LAMMPS datafilen og data.myadsorbate er navnet på datafilen LAMMPS.</li>
		<li>Slette delen obligasjoner og vinkler nederst i filen data.myadsorbate. Deretter legge listene bond og vinkel i filen $filename. POSCAR.bond_angle_info.text i data.myadsorbate. Merk: Indeksene for O-H bånd og H-O-H vinkel typen for molekylene i $filename. POSCAR.bond_angle_info.txt filen er begge satt til 1. Dermed bør bond og vinkel typer for adsorbates begynne å telle på 2.</li>
		<li>Rediger filen data.myadsorbate ved å legge til parameterne Lennard-Jones delen Par Coeffs og Coulomb parameterne til delen atomer . Norge-Jones og Coulomb parametere for H2O molekyler, adsorbate atomer og metall overflaten atomer må legges. Merk: Norge-Jones parameterne for Pt atomene og molekylene adsorbate atomer i denne protokollen er Hentet fra UFF38, TIP3P/CHARMM37og OPLS-AA39 kraftfelt, henholdsvis. Coulomb parametere for vann molekyler og adsorbates atomer hentes fra TIP3P/CHARMM37 og OPLS-AA39 kraftfelt, henholdsvis. Coulomb parametere for Pt atomer er satt til 0 i denne protokollen. Beregnet delvise belastninger kan også brukes for Coulomb parameterne for adsorbate atomer og Pt atomer.</li>
		<li>Kopier LAMMPS inndatafilen input.equil i mappen $WDPATH. Rediger gruppe variabelen på line 34 å angi atom type indeksene for vann oksygen og vann hydrogenatomer og variabelen gruppe på linje 35 å angi atom type indeksene for atomene Pt og adsorbate.</li>
		<li>Utføre LAMMPS programvaren ved å skrive mpiexec - np XX lmp_mpi < input.equil på kommandolinjegrensesnitt, der XX er antall CPU-kjerner bruker, og lmp_mpi er navnet på den kjørbare filen LAMMPS. Dette vil kjøre en energi minimering for å avgrense H2O konfigurasjonen, etterfulgt av en FFMD simulering på konstant antall H2O molekyler (N), volum (V) og temperatur (T) å bringe den vann til simulering temperatur, etterfulgt av en FFMD simulering kjøre konstant N, press (P), og temperatur (T) å bestemme fysisk riktig høyden på boksen simulering. Utdatafiler som skal brukes i del 3 kalles data.myadsorbate_npt og log.myadsorbate. Merk: Varigheten av NPT simuleringen skal lenge nok til å omfatte en "balanse" kjøres, der volumet av supercelle gjelder steady state, og en "produksjon" kjøre, som brukes til å prøve ensemblet gjennomsnitt (her, høyden på supercelle). Under balanse kjøre, skal volumet av supercelle når plottet mot tid flate ut en steady state verdi. Når dette skjer, kan NPT simuleringen sies å være i produksjonstiden. Kontroller balanse av NPT simuleringen ved å sikre at svingninger i høyden supercelle (lz) er minimal eller har kommet til en jevn verdi. Hvis store svingninger oppstår, deretter å generere en H2O konfigurasjon ved å redusere timestep på linjen 92 i input.equil-filen og gjentatt trinn 2.2.8 eller starter på nytt fra trinn 2.1.1.</li>
	</ol>
</li>
</ol>3. ekstra supercelle riktig høyde
<ol><li>Kjøre skript get_npt_lz.py på filen log.myadsorbate ved å skrive get_npt_lz.py log.myadsorbate i kommandolinjegrensesnittet. Dette skriptet utganger gjennomsnittlig supercelle høyden fra "produksjon" kjøre delen av NPT simulering i filen avg_lz.txt. Merk: Skriptet get_npt_lz.py forutsetter at LAMMPS skriver lengden på cellen z-dimensjonen (lz) til filen log.myadsorbate hver 1000 fs (passelig på linje 20 av skriptet get_npt_lz.py), som er standard i angitte input.equil LAMMPS inndatafilen. Skriptet get_npt_lz.py oppdager og forkaster 2 første ns (kan tilpasses på linje 19 av skriptet get_npt_lz.py) verd av lz verdier i filen log.myadsorbate, som de utgjør delen balanse av simulering, mens de resterende 3 ns utgjør den " produksjon"delen og dermed brukes av skriptet get_npt_lz.py til å beregne gjennomsnittlig z-dimensjonen lengden. I tillegg til filen avg_lz.txt utganger skriptet get_npt_lz.py en fil kalt npt_data.txt, som gir verdiene for lz som en funksjon av timestep, i tillegg til en fil som heter npt_plot.png, som plott de samme dataene. Tomten kan brukes til å bekrefte balanse av NPT simuleringen.</li>
	<li>Rekonstruere supercelle bruker den gjennomsnittlige høyden i PT.<ol>
<li>Kopier filen data.myadsorbate_npt til en ny katalog, her som $WD2PATH, og gi det data.myadsorbate.</li>
		<li>Rediger den nye data.myadsorbate-filen slik at lz høyden er lik gjennomsnittsverdien utdata fra skriptet get_npt_lz.py ved å endre zlo og zhi argumentene i filen data.myadsorbate slik at zlo er 0,0 og er lz verdien fra filen avg_lz.txt produsert i s tep 3.1.</li>
	</ol>
</li>
</ol>4. generere konfigurasjoner av H2O molekyler
<ol><li>Kopier LAMMPS inndatafilen input.prod til $WD2PATH. Rediger gruppe variabelen på linje 32 å angi atom type indeksene for vann oksygen og vann hydrogenatomer og variabelen gruppe på linje 33 å angi atom type indeksene for atomene Pt og adsorbate.</li>
	<li>Utføre LAMMPS programvaren ved å skrive mpiexec - np XX lmp_mpi < input.prod i kommandolinjegrensesnitt, der XX er antall CPU-kjerner bruker, og lmp_mpi er navnet på den kjørbare filen LAMMPS. Dette vil kjøre en konstant NVT simulering på H2O molekyler. Nøkkel utdatafilen fra denne simuleringen er filen dump.myadsorbate.lammpstrj. Merk: Varigheten av NVT simuleringen skal lenge nok til å omfatte en balanse kjøres, der energien i systemet gjelder steady state, og en produksjon løpe, som ensemblet gjennomsnitt (her, romlig stillingene molekylene) er samplet. Under balanse kjøre, skal energi av systemet når plottet mot tid flate ut en steady state verdi. Når dette skjer, kan NVT simuleringen sies å være i produksjonstiden.</li>
</ol>5. Bestem hydrogen bond levetiden for riktig tid prøvetaking
<ol><li>Redigere skriptet hb_lifetime_dist.py angi: timestep på den første rammen i filen dump.myadsorbate.lammpstrj ved å endre variabelen actualStart på linje 22, hvor ofte rammer er skrevet til filen LAMMPS banen ved å endre variabelen timestep på linje 23, den første og siste time skriptet bør vurdere (dvs. produksjon-delen av banen) ved å endre N_first og N_last-variabler på 24 og 25, om påfølgende rammer anses eller rammer hoppes over ved å endre til nevery variabel på line 26, og antall linjer per ramme-delen av filen banen ved å endre variabelen frameLine på linje 27. I tillegg redigere linjer 31 gjennom 35 angi atom-typer i filen data.myadsorbate tilhører adsorbate og som atom typer tilhører H2O molekyler. Merk: Hb_lifetime_dist.py skriptet analyserer H2O konfigurasjonene i produksjonsperiode og bestemmer hvis noen H2O molekyler er hydrogen limt til adsorbate. Det deretter teller simulering klokkeslettet hver hydrogen bånd intakte og rapporterer denne informasjonen som en fordeling av hydrogen bond levetid i enheter av ps. Den spesifikke versjonen av skriptet som følger med denne protokollen forutsetter at LAMMPS skriver konfigurasjonen av H2O molekyler til filen dump.myadsorbate.lammpstrj hver 1000 AS, som er standard i angitte input.prod LAMMPS inndatafilen. Den oppdager og sletter 2 første ns verd av konfigurasjoner i den i dump.myadsorbate.lammpstrj fil, som de utgjør delen balanse av simulering og bruker de resterende 3 ns å beregne hydrogen bond levetid.</li>
	<li>Kjøre skript hb_lifetime_dist.py på filen dump.myadsorbate.lammpstrj ved å skrive hb_lifetime_dist.py i kommandolinjegrensesnittet. Dette vil gi en fil kalt distribution_HB_lifetime.dat.</li>
	<li>Tegne inn dataene i filen distribution_HB_lifetime.dat å vise fordelingen av hydrogen bond levetid som oppstod under NVT simuleringen.</li>
	<li>Bestemme tidsintervallet bruker for prøvetaking tidsintervallet basert på de beregnede hydrogen bond levetid. Det beste valget er maksimal hydrogen bond levetid; Alternativt kan en verdi som ville fange 95% konfidensintervall brukes.</li>
</ol>6. eksempel konfigurasjoner av flytende H2O molekyler
<ol><li>Bestemme antall konfigurasjoner fra produksjonsperiode av NVT FFMD banen for ytterligere beregninger. Antall konfigurasjoner må velges slik at Minimumstid mellom konfigurasjoner er lik eller større enn prøvetaking tidsintervallet i seksjon 5.<ol>
<li>Redigere verdien for variabelen num_frames på linje 21 av skriptet lammps_frames.py til å angi antall konfigurasjoner trekke ut.</li>
		<li>Kjøre skript lammps_frames.py på filen dump.myadsorbate.lammpstrj ved å skrive lammps_frames.py i kommandolinjegrensesnittet. Dette vil generere en liste over simulering ganger tilsvarer konfigurasjonene som skal trekkes fra filen dump.myadsorbate.lammpstrj. Disse konfigurasjonene kan brukes som utgangspunkt strukturer i AIMD eller QM simuleringer. Merk: 1) lammps_frames.py skriptet automatisk oppdager LAMMPS loggen og dumpfiler samt produksjon delen av banen innen dumpfilen og deler antall konfigurasjoner i dumpfilen i 10 grupper. Brukeren kan også angi loggfilen, dumpfilen og antall konfigurasjoner fra kommandolinjegrensesnittet bruker -l, -d, og -n alternativer, henholdsvis. Gjør brukeren må skrive lammps_frames.py - n XX -l $logfilename -d $dumpfilename på kommandolinjegrensesnitt, der XX er ønsket antall konfigurasjoner, $logfilename er navnet på loggfilen LAMMPS og $dumpfilename er navnet på den LAMMPS bane (dump) fil. Simulering tidspunktene som utgang se median ganger i hver gruppe. 2) hvis konfigurasjonene beregnes i VASP med flagget LDIPOLE er aktivert, bør et lite lag vakuum plass legges til toppen av supercelle over vann laget. Dette vil lette konvergens av elektroniske strukturen i beregningen av VASP. Legger til en ytterligere 3 Å vakuum plass over H2O molekyler har vært vellykket i simuleringene beskrives nedenfor.</li>
	</ol>
</li>
</ol>

Forfatterne avsløre noen interessekonflikter.

Metoden som presenteres ble valgt for sin brukervennlighet implementering, men flere tilpasninger kan gjøres. For kan kraftfelt brukes i FFMD simuleringene endres. Endre kraftfelt parametere og/eller potentials gjøres ved å redigere LAMMPS inngang og datafiler. Tilsvarende kan løsemidler enn H2O brukes. For å gjøre denne endringen, må ønsket løsemiddel molekylet settes fra trinn 2.1.1, og LAMMPS inndatafiler må redigeres for å innlemme den aktuelle potensialer og parametere. Sette inn nye løsemiddel molekylet vil også kreve leverer interne koordinater i løsemiddel molekylet i en txt-fil som er analoge til filen water.txt.
En annen endring som gjøres er å endre området av overflaten skive. Resultatene i dette manuskriptet ansatt 3 Pt x 3 pkt eller 4 Pt x 4 pkt. overflaten plater, som har flater mindre enn 120 Å2. Som areal skive øker, øker også beregningsorientert bekostning. Beregningsorientert utgifter har størst innvirkning på seksjon 5 i denne protokollen. Hvis databehandling trinnene i del 5 blir beregningsmessig prohibitive, legge store data behandling strategier som de i Li et al. 201845 kan være ansatt.
Mulige kilder til usikkerhet for denne prosedyren inkluderer kraftfelt ansatt, metoden prøvetaking og samplingfrekvens. Vann strukturen bestemmes av feltet kraft som brukes, betyr at valg av kraftfelt kan påvirke bestemte konfigurasjoner av H2O molekyler. Vår gruppe har vurdert hvordan valgene av kraftfelt for H2O molekyler og Pt atomer påvirke samhandling energiene beregnet i FFMD og funnet at valg av kraftfelt bidrar mindre enn 0,1 eV å denne samhandlingen energi. En annen kilde til usikkerhet er sampling metoden, som påvirker de bestemte konfigurasjonene som brukes til å beregne et antall interessepunkter. Vår gruppe har sammenlignet arbeidsinnsatsen av metoden "tid prøvetaking" presenteres i denne protokollen med en "energi prøvetaking" metoden, som er partisk lavere energi konfigurasjoner av H2O molekyler, på samspill energier beregnet i DFT og funnet både disse Prøvetaking metoder verdier gir Statistisk like35,46. Samplingsfrekvens kan også påvirke resultatene. Vi har vurdert hvordan øke antall konfigurasjoner fra 10 til 30 000 påvirker gjennomsnittlig samhandling energiene beregnes i FFMD for 40 forskjellige C3HxO3 adsorbates og funnet at samplingfrekvens bidrar mindre enn 0.1 eV å gjennomsnittlig samhandling energi44.
Den viktigste begrensningen til denne metoden er at adsorbates er rundet av strukturer under vakuum under FFMD simuleringene. I virkeligheten ville adsorbates utstilling conformational endringer (bond strekninger, vinkel svinger, torsjonsmessig bevegelser, etc.) på grunn av normal termisk bevegelser, inkludert interaksjoner med løsemiddel molekyler. Forsøk på å inkludere conformational endringer i adsorbates i FFMD simuleringene ville kreve detaljert utvikling av kraftfelt for katalytisk overflate adsorbates, dvs. som utgjør vilkårene som beskriver bond strekninger, vinkel svinger og vridningsstivhet vilkår, blant andre. Som en fremtid retning av denne protokollen, utvikler vi slike force felt for adsorbates på faste flater, som vi bruker til å bestemme grad som bruker stive adsorbates påvirker resultatene.

Modellering molekylære fenomener involvert i heterogene katalyse under flytende forhold er nødvendig for forståelsen katalytisk funksjon; Dette er imidlertid utfordrende fordi det krever en fin balanse mellom kjemiske nøyaktighet og beregningsorientert utgifter. Generelt, siden katalyse innebærer brudd og forming av kjemisk obligasjoner, må kvantemekanikk brukes i minst grad; Imidlertid er lang simuleringer utfordrende i kvantemekanikk, siden de krever betydelige datamaskinressurser. Siden molekyler i flytende fase er under konstant termisk bevegelse, simuleringer må også inkludere configurational prøvetaking, dvs., må de innlemme flere romlige ordninger av flytende molekyler, som hver annen romlige arrangement (dvs. hver konfigurasjon) har en annen energi. Dette betyr at flere konfigurasjoner av flytende molekyler må simuleres for hver katalytisk arter av interesse. Disse behov-bruke kvantemekanikk og utføre flere beregninger per katalytisk arter-kan gjengi modellering i heterogene katalyse under flytende fase beregningsmessig uløselige. Formålet med metoden beskrevet her er å aktivere beregningsmessig medgjørlig simuleringer av fenomener i heterogene katalyse under flytende fase.
Vi er spesielt interessert i heterogeneously catalyzed reaksjoner som utføres under flytende vann. Vannmolekyler har betydelig innflytelse på katalytisk fenomener, slik som samhandler med katalytisk arter (f.eks via spredning styrker og hydrogenbinding)1,2,3,4,5 ,6,7,8,9,10,11,12,13,14 , 15 , 16 , 17 , 18 , 19 , 20 , 21 , 22 , 23, delta i katalytisk reaksjoner1,7,8,9,15,21,22,24 ,25,26,27, og påvirke reaksjon veier og/eller katalytisk priser1,11,12,15, 18,23,25,27,28,29,30,31. Modellering av disse fenomenene er utført med QM og/eller ab initio molekylære dynamikk (AIMD)1,2,6,7,14,22 ,25,27,28,32,33,34, tvinge feltet molekylære dynamikk (FFMD)35 , og kvantemekanikk/molekylær mekanikk (QM/MM)10. I AIMD og FFMD flyttes atomene i systemet etter Newtons ligninger av bevegelse i henhold til kreftene som virker på dem. I AIMD, er systemet energi og styrker beregnet med kvantemekanikk, mens i FFMD, systemet energi og styrker beregnes ved hjelp av force felt som er algebraiske uttrykk som er parameteriserte basert på eksperimentelle eller QM data. Delen av systemet der bond bryte og danner oppstår beregnes med QM i QM/MM, og resten av systemet beregnes med MM, som sysselsetter kraftfelt. Fordi de bruker direkte QM, AIMD og QM/MM er bedre egnet for fange bond bryte og forming som forekommer i vandige fasen heterogene katalyse; FFMD er imidlertid betydelig mer beregningsmessig medgjørlig og dermed bedre egnet til å generere konfigurasjonene av væske H2O molekyler. Metoden som presenteres i denne protokollen balanserer kjemiske nøyaktighet og beregningsorientert utgifter ved å bruke en kombinasjon av QM og FFMD.
Spesielt bruker denne metoden FFMD simuleringer for å generere konfigurasjoner av væske H2O og QM beregne systemet energier. FFMD er utført med LAMMPS. 36 kraftfelt brukt i FFMD i dette arbeidet bruker Lennard-Jones + Coulomb (LJ + C) potensialene, der parameterne LJ er tatt fra TIP3P/CHARMM modell37 H2O, den universelle kraftfelt38 (UFF) for Pt, og OPLS-AA kraftfelt39 for katalytisk arter, og parameterne Coulomb er tatt fra TIP3P/CHARMM37 modellen for H2O og OPLS-AA kraftfelt39 for katalytisk arter. Coulomb parameterne for Pt atomer er angitt til 0. QM beregninger utføres ved hjelp VASP kode40,41,42, en tetthet funksjonelle teori (DFT) kode. Vann molekylet innsettinger utføres med en kode utviklet in-house kalt Monte Carlo Plug-in for Quantum metoder (MCPliQ). Filkonverteringer fra VASP til LAMMPS i denne protokollen utføres med visuelle Molecular Dynamics (VMD) programvare43.
Protokollen er ment å generere konfigurasjoner av flytende vannmolekyler rundt katalytisk art på flat transisjonsmetall overflater på lav dekning. Dekningen er betegnet θ og definert som antall adsorbates per overflaten metal atom (dvs. antall overflate adsorbates normalisert av antall metall atomer i det øverste laget av metall skive i catalyst-modellen). I dette manuskriptet lav dekning er definert som θ ≤ 1/9 monolayer (ML), der 1 ML betyr en katalytisk Art per overflaten metal atom. Catalyst modellene bør plasseres i periodiske simulering bokser. Boksene simulering trenger ikke å være kuber. Dette manuskriptet demonstrerer bruken av protokollen for å generere konfigurasjoner av flytende H2O som kan brukes til å beregne antall interesse i vandige fasen heterogene katalyse.
Denne protokollen krever at brukeren har tilgang til installert og fungere versjoner av VASP, MCPliQ, LAMMPS og VMD programvaren. Mer informasjon om VASP (https://www.vasp.at/), LAMMPS (https://Lammps.sandia.gov/) og VMD (https://www.ks.uiuc.edu/Research/vmd/) er tilgjengelig på deres respektive nettsteder. MCPliQ-programvare er dokumentert på https://github.com/getman-research-group/JoVE_article, sammen med inndatafiler og Python-skript som er nevnt i denne protokollen. Denne protokollen forutsetter at den kjørbare filer og skript nevnt i kjøres på en datamaskin som høy ytelse forskning og installeres i en mappe i brukerens $PATH variabel. Hvis en kjørbar fil eller skript er plassert på et sted som ikke er i brukerens $PATH, og banen til den kjørbare filen må være inkludert å effektuere den. Kjørbare filer og skript kjøres i trinn 2.1.2, 2.2.1, 2.2.8, 3.1, 4.2, 5.2 og 6.1.2. For eksempel å kjøre MCPliQ koden i trinn 2.1.2 fra en katalog som ikke brukeren 's $PATH, brukeren skriver inn $PATHTOMCPLIQ/mcpliq på kommandolinjegrensesnittet i stedet for mcpliq, der $PATHTOMCPLIQ er plasseringen der mcpliq kjørbare er lagret (f.eks $PATHTOMCPLIQ kan være ~ / bin). Før du starter denne protokollen, alle kjørbare filer og skript bør gis kjørbart tillatelser (f.eks i Linux, dette kan gjøres ved å skrive chmod + x mcpliq i kommandolinjegrensesnittet fra katalogen der mcpliq kjørbare er lagret). Videre noen moduler som kreves av noen programvare eller skript skal lastes (disse avhengighetene blir bestemte personlige installasjoner av ulike programvaren og datamaskinen der simuleringene kjøres).

<table>
	<tbody>
		<tr>
			<td>VASP software</td>
			<td>Computational Materials Physics, Dept. of Physics, University of Vienna</td>
			<td>vasp.5.4.4</td>
			<td>Standard parallel VASP executable in the newest version.</td>
		</tr>
		<tr>
			<td>LAMMPS software</td>
			<td>Sandia National Laboratory</td>
			<td>31Mar17-dp</td>
			<td>Double-precision, parallel LAMMPS executable from 31 March 2017.</td>
		</tr>
		<tr>
			<td>VMD software</td>
			<td>Theoretical and Computational Biophysics Group, University of Illinois at Urbana-Champaign</td>
			<td>1.9.3</td>
			<td>Standard VMD executable in the newest version.</td>
		</tr>
		<tr>
			<td>MCPliQ software</td>
			<td>Getman Research Group, Dept. of Chemical and Biomolecular Engineering, Clemson University</td>
			<td></td>
			<td>Executable and input files for the MCPliQ software availabe from the Getman Research Group GitHub page.</td>
		</tr>
		<tr>
			<td>JoVE article scripts</td>
			<td>Getman Research Group, Dept. of Chemical and Biomolecular Engineering, Clemson University</td>
			<td></td>
			<td>Python scripts for this JoVE manuscript available from the Getman Research Group GitHub page.</td>
		</tr>
		<tr>
			<td>H2O PDB file</td>
			<td>Getman Research Group, Dept. of Chemical and Biomolecular Engineering, Clemson University or RCSB Protein Data Bank</td>
			<td></td>
			<td>PDB file for a water molecule, available from the Getman Research Group GitHub page or at http://www.rcsb.org/ligand/HOH.</td>
		</tr>
	</tbody>
</table>

multiscale sampling of a heterogeneous water/metal catalyst interface using density functional theory and force-field molecular dynamics

Et betydelig antall heterogeneously-katalysert kjemiske prosesser oppstå under flytende forhold, men simulere katalysator funksjonen under slike forhold er utfordrende når det er nødvendig å inkludere løsemiddel molekyler. Bond bryte og danner prosesser modellert i disse systemene nødvendiggjør bruk av quantum kjemiske metoder. Siden molekyler i flytende fase er under konstant termisk bevegelse, må simuleringer også inkludere configurational prøvetaking. Dette betyr at flere konfigurasjoner av flytende molekyler må simuleres for hver katalytisk arter av interesse. Målet med protokollen presenteres her er å generere og prøve baner av konfigurasjoner av flytende vannmolekyler rundt katalytisk art på flat transisjonsmetall overflater på en måte som balanserer kjemiske nøyaktighet med beregningsorientert bekostning. Spesielt brukes kraftfelt molekylære dynamikk (FFMD) simuleringer til å generere konfigurasjoner av flytende molekyler som senere kan brukes i quantum mekanikk-baserte metoder som tetthet funksjonelle teori eller ab initio molekylær dynamikk. For å illustrere dette, i dette manuskriptet, brukes protokollen for katalytisk mellomprodukter som kan være involvert i veien for nedbryting av glyserol (C3H8O3). Strukturene som genereres ved hjelp av FFMD er modellert i DFT for å anslå dannelsesentalpien solvation av katalytisk arter og identifisere hvor H2O molekyler delta i katalytiske decompositions.

En bruk av denne protokollen er å beregne energiene av samspillet mellom flytende vann og katalytisk arter, dvs. ΔEint35:
∆Eint=Ekatalytisk arter + H2O+Eren katalysator overflaten-Ekatalytisk arter-Eren katalysator overflaten + H2O
der Ekatalytisk arter + H2Oenergien i er en konfigurasjon av H2O molekyler rundt katalytisk art på en metalloverflate, Eren katalysator overflaten energien av ren katalysator overflaten i vakuum, E Katalytisk arter er energien av katalytisk arter på en metalloverflate i vakuum, og Erydde katalysator overflaten + H2O er energien av konfigurasjonen av H2O over katalysator overflaten med katalytisk arter fjernet. Plasseringen av H2O molekylene brukes til å beregne Ekatalytisk arter + H2O og Erydde katalysator overflaten + H2O må være identisk. Alle verdier av E beregnes ved hjelp av VASP-koden. Antall ΔEint omfatter alle fysiske og kjemiske samspillet mellom alle molekylene i flytende vann strukturen og den katalytiske arten og gir en fornuftig anslå verdi av entalpi av solvation av den katalytiske arter, som er nødvendig for å beregne sin fri energi solvation og totalt fri energi. Tabell 1 gir verdier for ΔEint beregnet for arter på Pt(111) katalysator overflaten med kjemiske formler lik CxHyOz i enheter av eV (1 eV = 96.485 kJ/mol). Verdiene ble beregnet på dekning ≤1/9 ML.35,46 verdiene rapportert er gjennomsnitt tatt over 10 konfigurasjoner av flytende H2O, og usikkerhet er rapportert som standardavvik. Alle verdiene er negative, indikerer gunstige interaksjon med vann.
Et annet program av denne protokollen er å generere Start strukturer for AIMD. Filmen 1 er en film av en AIMD bane som ble startet fra en konfigurasjon som er generert av denne protokollen. I starten av denne filmen vises en COH adsorbate på en Pt(111) overflate under en struktur av flytende H2O. En H2O molekylet er understreket, som dannet en hydrogen bånd med COH. I løpet av filmen, dette H2O molekylet Sammendrag proton fra COH adsorbate og innskudd en andre hydrogen atom på Pt(111) overflaten. H2O molekylet dermed bidrar til å katalysere reaksjonen COH * + * → CO * + H *, hvor den * s angir katalytisk nettsteder. Denne simuleringen høydepunkter hovedstyrken og Hovedhensikten med metoden multiscale prøvetaking beskrevet her. Mange konfigurasjoner av H2O molekyler genereres med FFMD, på grunn av sin styrke i beregningsorientert tractability. Men er en begrensning av FFMD at den ikke kan fange bond bryte og danner med mindre en reaktiv kraftfelt implementeres. AIMD bruker kvantemekanikk til å beregne energier og dermed kan ta bond bryte og forming. Men er AIMD også beregningsmessig krevende å generere alle konfigurasjonene av H2O molekyler nødvendig for å sikre tilstrekkelig utvalg er oppnådd. Derfor kombinerer denne protokollen de to metodene.
Strukturer av flytende H2O molekyler generert av denne prosedyren er avhengig av innstillinger. Disse innstillingene kan ha utilsiktede påvirkninger på vann strukturer. For eksempel kan når intermolekylære avstandene blitt for liten eller andre parametere i molekylære dynamikk inndatafiler er feil angitt eller ta på unphysical verdier, vann strukturen bli urimelig. Under disse omstendigheter vil strukturen i vann "blåse" unintendedly under FFMD banen. Figur 1 viser et eksempel på dette. Bildet til venstre er start strukturen for en FFMD kjøre bildet på høyre side er et øyeblikksbilde tatt innen 1 ps å starte simuleringen. Som kan sees, har H2O molekyler flyttet langt fra overflaten. Dette skyldes feil innstillinger gjort i simulering inndatafiler er ikke en struktur som er sannsynlig å skje i virkeligheten.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/59284/59284fig1v2.jpg"> Figur 1: Eksempel på et negativt resultat. Kraftfelt molecular dynamics simulering blåser"opp" på grunn av en unphysical innstilling eller verdi. Venstre bilde: Start geometrien av Pt(111) overflate, adsorbate og flytende vann struktur. Høyre bilde: geometrien av Pt(111) overflate, adsorbate og flytende vann strukturen mindre enn 1 ps senere. I høyre bildet, har H2O molekyler skilt fra overflaten på grunn av unphysically store styrker. <a href="https://www.jove.com/files/ftp_upload/59284/59284fig1v2large.jpg" target="_blank">Klikk her for å se en større versjon av dette tallet.</a>
<img alt="Movie 1" src="/files/ftp_upload/59284/59284mov1.jpg"> Film 1: Ab initio molekylære dynamikk (AIMD) simulering startet fra en konfigurasjon som genereres i multiscale prøvetaking. Et H2O molekyl som er opprinnelig hydrogen limt til en COH adsorbate på en Pt(111) overflate Sammendrag proton fra COH og innskudd en andre hydrogen på Pt(111) overflaten. Denne bond bryte og danner hendelsen kan hentes av AIMD men ikke med kraftfelt molekylære dynamikk (FFMD) med mindre en reaktiv kraftfelt brukes. Den opprinnelige konfigurasjonen av H2O molekyler i denne AIMD simuleringen ble generert ved hjelp FFMD som beskrevet i dette manuskriptet. <a href="https://www.jove.com/files/ftp_upload/59284/movie1.mp4" target="_blank">Vennligst klikk her å se denne videoen. (Høyreklikk for å laste ned.)</a>
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always"><tbody>
<tr>
<td>Katalytisk arter</td>
			<td>∆Eint (eV)</td>
		</tr>
<tr>
<td>COH</td>
			<td>-0.70 ± 0,07</td>
		</tr>
<tr>
<td>CO</td>
			<td>-0.03 ± 0,03</td>
		</tr>
<tr>
<td>LM2OH</td>
			<td>-0.64 ± 0,12</td>
		</tr>
<tr>
<td>CHO-CHOH-CH2OH</td>
			<td>-0.93 ± 0.22</td>
		</tr>
<tr>
<td>COH-COH-CH2OH</td>
			<td>-0.87 ± 0,23</td>
		</tr>
<tr>
<td>COH-CHOH-COH</td>
			<td>-1.72 ± 0.26</td>
		</tr>
<tr>
<td>CHOH-COH-CO</td>
			<td>-1.57 ± 0,25</td>
		</tr>
<tr>
<td>CHO-CO-CO</td>
			<td>-0.31 ± 0,19</td>
		</tr>
</tbody></table>Tabell 1: Vann-katalytiske arter samhandling energi resultater. Samhandling energier i eV beregnet for åtte CxHyOz adsorbates på Pt(111). Verdiene rapportert tas gjennomsnitt over flere konfigurasjoner av væske H2O. Usikkerheten er standardavvik av gjennomsnitt. 1 eV = 96.485 kJ/mol.

Watch this Scientific Journal Video about Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics at JoVE.com

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

Målet med protokollen presenteres her er å generere og prøve baner av konfigurasjoner av flytende vannmolekyler rundt katalytisk art på flat transisjonsmetall underlag. Samplet konfigurasjoner kan brukes som utgangspunkt strukturer i quantum mekanikk-baserte metoder.

Multiscale utvalg av en heterogen vann/Metal katalysator grensesnitt ved hjelp av tetthet funksjonelle teori og Force-feltet molekylære dynamikk

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

Protokollmodellene, termiske og kvantefenomener i væskefase heterogen katalyse. Det er den første, så vidt vi vet, å innlemme kvantemekanikk med fullstendig prøvetaking av et eksplisitt flytende miljø. Konfigurasjoner av flytende molekyler generert fra denne protokollen representerer de som forventes under faktiske reaksjonsforhold og kan brukes til å utforske molekylære nivåfenomener som er avhengige av romlige ordninger av molekyler.Konfigurasjonene av flytende molekyler generert av denne protokollen gir innsikt i rollene som løsemiddel spiller på flytende fase heterogen katalyse. Hvis du prøver denne protokollen for første gang, anbefaler jeg at du først sikrer at du har tilgang til og kan kjøre VASP, MCPliQ, VMD og LAMMPS programvare. Å demonstrere prosedyren vil være Tianjun Xie, en utdannet student fra laboratoriet mitt.Etter å ha generert adsorbatstrukturen, generer LAMMPS-inndatafilene for en NPT-simulering og likevekt cellevolumet ved hjelp av FFMD. Kopier LAMMPS-inndatafilen til arbeidsmappen. Rediger gruppevariabelen på linje 34 for å angi atomtypeindeksene for vann oksygen og vann hydrogenatomer.I gruppevariabelen på linje 35, for å angi atomtypeindeksene for platina- og adsorbatatomer. På linje 17 i inndatafilen redigerer du kjøretrinnvariabelen for å angi varigheten av NPT-simuleringen lenge nok til å utgjøre en likevektskjøring og en produksjonskjøring. Utfør LAMMPS-programvaren ved å skrive inn kommandoen på kommandolinjegrensesnittet, som inkluderer informasjon om antall CPU-kjerner som skal brukes, og navnet på lammps kjørbare filen.Energiminimeringen foredler vannmolekylkonfigurasjonen og etterfølges av en FFMD-simulering utført ved konstant antall vannmolekyler, volum og temperatur for å bringe vannet til simuleringstemperaturen. En annen FFMD-simulering kjøres deretter ved konstant antall vannmolekyler, trykk og temperatur for å bestemme den fysisk korrekte høyden på simuleringsboksen. Utdatafilene vil bli brukt senere.Etter NPT-simuleringen plotter du høyden på supercellen mot tiden. Punktet der det utnivåer til en jevn tilstandsverdi er punktet i NPT-simuleringen når produksjonen kan begynne. Kontroller likevekt av NPT-simuleringen ved å sikre at svingningene i høyden på supercellen er minimale eller har konvergert til en jevn verdi.Hvis det oppstår store svingninger, åpner du inngangen. equil fil og redusere timestep på linje 92 å regenerere vannmolekylet konfigurasjon og utføre LAMMPS programvare som før. Du kan begynne ved å skrive inn kommandolinjegrensesnittet for å utføre skriptet.Dette skriptet sender ut den gjennomsnittlige supercellehøyden fra produksjonskjøringsdelen av NPT-simuleringen til en TXT-fil. Skriptet leser lengden på cellens z-dimensjon ved 1000 femtosekunders intervaller, som er standardintervallet for utskrift av informasjon i LAMMPS. Hvis det er ønskelig med et annet utskriftsintervall, kan det endres ved å redigere linje 20 på get_npt_lz.py script og linje 16 av LAMMPS inngang. equil-filen. Skriptet oppdager og forkaster de to første nanosekunder verdt lz-verdier, da de utgjør likevektsdelen av simuleringen.Varigheten av likevektskjøringen kan endres ved å redigere linje 19 i filen. De resterende tre nanosekunder utgjør produksjonsdelen, og brukes dermed til å beregne den gjennomsnittlige z-dimensjonslengden. I tillegg utdata skriptet en annen TXT-fil som gir verdier av lz som en funksjon av timestep samt en PNG-fil, som plotter de samme dataene.Plottet kan brukes til å verifisere likevekt av NPT-simuleringen. Hvis du vil rekonstruere supercellen ved hjelp av gjennomsnittshøyden som er bestemt i NPT, kopierer du den tidligere genererte datafilen til en ny arbeidsmappe og gir den nytt navn som data.myadsorbate. Deretter redigerer du den nye datafilen for å endre zlo til 0,0 og zhi til lz-verdien fra gjennomsnittlig verdiutdata i TXT-filen.Kopier LAMMPS-inndatafilen til den nye arbeidskatalogen. Rediger gruppevariabelen på linje 32 for å angi atomtypeindeksene for vannsygengen- og vannhydrogenatomer og gruppevariabelen på linje 33 for å angi atomtypeindeksene for platina- og adsorbatatomer. Deretter på linje 16 redigerer du runStep-variabelen slik at den er lang nok til å utgjøre en likevektskjøring og en produksjonskjøring.Skriv inn kommandoen for å kjøre LAMMPS i kommandolinjegrensesnittet for å kjøre LAMMPS-programvaren. Dette vil kjøre en konstant NVT-simulering på vannmolekylene, og nøkkelutdatafilen genereres. NVT-simuleringen består av en likevektsdel og en produksjonsdel.Produksjonsdelen starter når energien i systemet plottet mot tidsnivåer av. For å begynne, åpne hydrogenbindingsskriptet. Endre den faktiske startvariabelen på linje 22 for å angi tidssteget for den første tidsrammen.Endre tidsstegsvariabelen på linje 23 for å angi hvor ofte rammer skrives til LAMMPS banefilen. Endre N_first og N-siste variabler på linje 24 og 25 for å angi første og siste gangtrinn, og endre den nevery variabelen på linje 26 for å angi om etterfølgende rammer skal vurderes eller hoppes over. Angi antall linjer per rammedel i banefilen ved å endre rammelinjevariabelen på linje 27.Rediger i tillegg linjene 31 til 35 for å angi hvilke atomtyper i dataene. myadsorbatfilen tilhører adsorbatet, og hvilke atomtyper som tilhører vannmolekylene. Skriptet analyserer vannkonfigurasjonene i produksjonskjøringen og avgjør om noen vannmolekyler er hydrogenbindet til adsorbatet.Deretter teller simuleringstiden at hvert hydrogenbinding forblir intakt og rapporterer denne informasjonen som en fordeling av hydrogenbindingslevetiden i enheter av picosekunder. LAMMPS skriver konfigurasjonen av vannmolekyler til filen hver 1000 femtosekunder, som er standard i den medfølgende LAMMPS-inndatafilen. Den oppdager og forkaster de to første nanosekunder verdt konfigurasjoner i filen som de utgjør likevektdelen av simuleringen og bruker de resterende tre nanosekunder til å beregne hydrogen-bond levetid.Hvis du vil kjøre skriptet, skriver du inn på kommandolinjegrensesnittet. Det genererer deretter en DAT-fil. Plott dataene i filen for å vise fordelingen av hydrogenbindingslevetid som oppstod under NVT-simuleringen.Hvis du vil bestemme tidsintervallet som skal brukes for tidsprøvetakingsintervallet, bruker du en tidsintervall som er større enn eller lik den maksimale hydrogenbindingslevetiden. Bestem antall konfigurasjoner fra produksjonskjøringen for NVT FFMD-banen for å prøve slik at minimumstiden mellom konfigurasjoner er lik eller større enn tidsprøveintervallet som er identifisert tidligere. På det tidligere skrevne rammeuttrekksskriptet redigerer du standardverdien for antall rammer variabel på linje 21 for å angi antall konfigurasjoner som skal pakkes ut.Hvis du vil kjøre skriptet, skriver du inn skriptnavnet i kommandolinjegrensesnittet. Dette vil sende ut en liste over simuleringstider som tilsvarer konfigurasjonene som skal trekkes ut fra NVT-simuleringsfilen. Disse konfigurasjonene kan brukes som startstrukturer i AIMD- eller QM-simuleringer.I denne prosedyren ble FFMD brukt til å generere den første konfigurasjonen av vannmolekyler. AIMD-simuleringen viser at et vannmolekyl som opprinnelig er hydrogenbundet til en sukkeralkohol adsorbat på en platina-111 overflate abstraherer hydrogenet fra alkoholadsorbatet og avsetter et annet hydrogen på platina-111-overflaten. Strukturene til flytende vannmolekyler er avhengige av inngangsinnstillinger.Å sette disse feil kan ha utilsiktede påvirkninger på vannstrukturene. I dette tallet er venstre side startstrukturen for en FFMD-kjøring. Og høyre side er innenfor ett picosekund for å starte simuleringen.FFMD-simuleringen blåser opp på grunn av ukystisk store kraftinnstillinger, noe som fører til at vannmolekyler beveger seg langt bort fra overflaten. Konfigurasjonene kan brukes i kvantemekanikken, eller QMM-simulering, eller de kan brukes til å analysere statistikk knyttet til molekylenes romlige posisjoner. Denne teknikken baner vei for forskere å utforske rollene som flytende reaksjonsmiljøer har på katalyse ved å generere faktiske konfigurasjoner av flytende molekyler ved katalytiske grensesnitt.

Research

JoVE Journal

Chemistry

Layer-by-layer Synthesis and Transfer of Freestanding Conjugated Microporous Polymer Nanomembranes

Lag-på-lag Syntese og overføring av Frittstående Conjugated mikroporøse Polymer Nanomembranes

CMP as large surface area materials have attracted growing interest recently, due to their high variability in the incorporation of functional groups in ...

HKUST-1 as a Heterogeneous Catalyst for the Synthesis of Vanillin

HKUST-en som en heterogen katalysator for syntese av Vanillin

Vanillin (4-hydoxy-3-methoxybenzaldehyde) is the main component of the extract of vanilla bean. The natural vanilla scent is a mixture of approximately ...

Facile Synthesis of Worm-like Micelles by Visible Light Mediated Dispersion Polymerization Using Photoredox Catalyst

Facile Synthesis of Worm-lignende Miceller av synlig lys mediert dispersjonspolymeriseringen Bruke Photoredox Catalyst

Presented herein is a protocol for the facile synthesis of worm-like micelles by visible light mediated dispersion polymerization. This approach begins ...

Synthesis and Performance Characterizations of Transition Metal Single Atom Catalyst for Electrochemical CO2 Reduction

Synthesis and Performance Characterizations of Transition Metal Single Atom Catalyst for Electrochemical CO2 Reduction

Syntese og ytelse karakteristikkene av transisjonsmetall enkelt Atom katalysator elektrokjemiske CO2 reduksjon

This protocol presents both the synthesis method of the Ni single atom catalyst, and the electrochemical testing of its catalytic activity and selectivity ...

Preparation of Functional Silica Using a Bioinspired Method

Utarbeidelse av funksjonelle Silica bruke en Bioinspired

The goal of the protocols described herein is to synthesize bioinspired silica materials, perform enzyme encapsulation therein, and partially or totally ...

Heterogeneous Removal of Water-Soluble Ruthenium Olefin Metathesis Catalyst from Aqueous Media Via Host-Guest Interaction

Heterogen fjerning av vannløselige Ruthenium olefiner Metathesis Catalyst vandige Media Via Host-gjest interaksjon

A highly efficient transition metal catalyst removal method is developed. The water-soluble catalyst contains a newly-designed NHC ligand for the catalyst ...

Imaging Corrosion at the Metal-Paint Interface Using Time-of-Flight Secondary Ion Mass Spectrometry

Imaging korrosjon på metal-Paint Interface bruke time-of-Flight sekundær ion Mass massespektrometri

Corrosion developed at the paint and aluminum (Al) metal-paint interface of an aluminum alloy is analyzed using time-of-flight secondary ion mass ...

[(DPEPhos)(bcp)Cu]PF6: A General and Broadly Applicable Copper-Based Photoredox Catalyst

[DPEPhosbcpCu]PF6: A General and Broadly Applicable Copper-Based Photoredox Catalyst

[(DPEPhos) (BCP) Cu] PF6: en generell og bredt gjeldende kobber-baserte Photoredox Catalyst

Our group recently reported the use of [(DPEPhos)(bcp)Cu]PF6 as a general copper-based photoredox catalyst which proved efficient to promote the ...

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

Vibrasjonsspektra av et N719-Chromophore/Titania-grensesnitt fra Empirisk-potensiell molekylær-dynamikksimulering, løst av en romtemperatur ionisk væske

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

Development of Heterogeneous Enantioselective Catalysts using Chiral Metal-Organic Frameworks (MOFs)

Development of Heterogeneous Enantioselective Catalysts using Chiral Metal-Organic Frameworks MOFs

Utvikling av heterogene Enantioselective katalysatorer bruker chiral metal-organiske rammeverk (MOFs)

Substrate size discrimination by the pore size and homogeneity of the chiral environment at the reaction sites are important issues in the validation of ...