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

Dit onderzoek werd gefinancierd door de National Science Foundation via award nummer CBET-1438325. Fellowship ondersteuning voor CJB door NASA opleiding Grant NX14AN43H wordt dankbaar erkend. Simulaties werden uitgevoerd op de Palmetto Supercomputer-Cluster, dat wordt bijgehouden door de Cyberinfrastructure Technology Group Clemson Universiteit. Wij danken Dr. Paul J. Meza-Morales voor het testen van het protocol.

1. het genereren van de adsorbate-structuur
<ol><li>Maak een bestand van de VASP POSCAR bestaande uit een supercell met periodieke randvoorwaarden, zoals u zou doen als u aan het uitvoeren waren simulaties van adsorbates op metalen oppervlakken onder vacuüm. De supercell dient een eerste schatting van de structuur van de adsorbate en het metalen oppervlak, alsmede vacuüm van de ruimte boven de adsorbate voor het toevoegen van H2O moleculen. Details over hoe dit te doen vindt u in de vorige werk35,44. Opmerking: Het is belangrijk dat de hoogte van de vacuüm ruimte ten minste 12 Å boven de top van de adsorbate worden.</li>
	<li>Ontspan de structuur en zijn energie met behulp van de code VASP minimaliseren. Details over hoe dit te doen vindt u in de vorige werk35,44. Dit zal het produceren van een bestand met de naam CONTCAR, die zal worden gebruikt in de volgende sectie.</li>
</ol>2. Voeg expliciete H2O moleculen
<ol><li>Toevoegen van N watermoleculen aan de vacuüm ruimte in de CONTCAR gemaakt in stap 1.2 met behulp van de code van de MCPliQ, waar N = ρV, ρ is de dichtheid van water, en V = het volume van de vacuüm ruimte boven de adsorbate. Opmerking: ρ moeten worden genomen als de dichtheid van water zoals bepaald door het model TIP3P/CHARMM water op de temperatuur van de simulatie. V zal worden verfijnd in de volgende stap.<ol>
<li>Geef de volgende informatie in het bestand MCPliQ master_input.txt: het aantal H2O moleculen om toe te voegen (N) door het veranderen van het eerste argument in lijn 28, het pad naar het bestand water.txt door het tweede argument in lijn 28, en het minimum te wijzigen en maximale hoogte van de supercell die kan worden bezet door watermoleculen door de Minimum z-coördinaat in lijn 11 en de Maximum z-coördinaat in lijn 12 te wijzigen.</li>
		<li>Het uitvoeren van de code van de MCPliQ door te typen mcpliq van de opdrachtregelinterface voor het invoegen van watermoleculen in het CONTCAR bestand. De code zal de uitgang van een of meer bestanden met de bestandsextensie. POSCAR. Opmerking: Als meer dan één POSCAR bestanden worden geproduceerd, zij zal de naam POSCAR_n.POSCAR. Selecteer het bestand waar de n grootste is.</li>
	</ol>
</li>
	<li>Genereren van de LAMMPS-invoerbestanden voor een simulatie van het NPV en equilibreer van het volume van de cel met behulp van FFMD in het NPV ensemble in LAMMPS.<ol>
<li>Uitvoeren van de script-lmps_bond_angle.py op de. POSCAR bestand gegenereerd in stap 2.1.2 door te typen lmps_bond_angle.py $filename. POSCAR op de opdrachtregelinterface, waar $filename de naam van is de. POSCAR bestand gegenereerd in stap 2.1.2. Dit script maakt een bestand met de naam $filename. POSCAR.bond_angle_info.txt waarin obligaties en hoeken die zal worden gebruikt in het gegevensbestand van de LAMMPS.</li>
		<li>Open VMD en selecteer bestand > Nieuwe Molecule in het hoofdvenster te openen van het venster van De Browser van het dossier van de Molecule . VASP_POSCAR selecteren in het keuzemenu Bestandstype bepalen . Klik op Bladeren en navigeer naar de $filename. POSCAR bestand. Klik op Load om te openen de $filename. POSCAR bestand.</li>
		<li>Open de Console van de Tk binnen VMD door het selecteren van extensies > Tk Console van het hoofdvenster van VMD.</li>
		<li>Het volgende commando in de console Tk: topo writelammpsdata $WDPATH/data.myadsorbate volledige, waar $WDPATH is de map op de computer waar VMD zal het schrijven van het gegevensbestand van de LAMMPS en data.myadsorbate is de naam van het gegevensbestand van de LAMMPS.</li>
		<li>De sectie van de obligaties en hoeken aan de onderkant van het data.myadsorbate-bestand verwijderen. Vervolgens voegen de obligatie- en hoek lijsten in het bestand $filename. POSCAR.bond_angle_info.text in data.myadsorbate. Opmerking: De index voor de obligatie-type O-H en H-O-H hoek type voor de watermoleculen in de $filename. Beide POSCAR.bond_angle_info.txt bestand zijn ingesteld op 1. Dus, obligatie- en hoek typen voor adsorbates beginnen met tellen bij 2.</li>
		<li>Bewerk het data.myadsorbate bestand door de Lennard-Jones-parameters toevoegen aan de sectie Coeffs van het paar en de Coulomb parameters naar de sectie van de atomen . Lennard-Jones en Coulomb parameters voor de H2O moleculen, atomen van de adsorbate en metalen oppervlak atomen moeten worden toegevoegd. Opmerking: De Lennard-Jones-parameters voor Pt atomen, watermoleculen en atomen van de adsorbate in dit protocol worden verkregen uit de UFF38, TIP3P/CHARMM37en OPLS-AA39 krachtvelden, respectievelijk. Coulomb parameters voor water moleculen en adsorbates atomen zijn verkregen uit de TIP3P/CHARMM37 en OPLS-AA39 krachtvelden, respectievelijk. Coulomb parameters voor Pt atomen worden ingesteld op 0 in dit protocol. U kunt ook kunnen berekende gedeeltelijke kosten worden gebruikt voor de Coulomb parameters voor adsorbate atomen en Pt atomen.</li>
		<li>Kopieer de LAMMPS input bestand input.equil naar de map $WDPATH. Bewerk de groep op regel 34 om aan te geven van de atoom type indexen voor het water zuurstof en water waterstofatomen en de variabele van de groep op regel 35 om aan te geven van de atoom type indexen voor de Pt- en adsorbate-atomen.</li>
		<li>Uitvoeren van de LAMMPS-software door te typen mpiexec - np XX lmp_mpi < input.equil op de opdrachtregelinterface, waar XX is het aantal CPU-cores te gebruiken, en lmp_mpi is de naam van het uitvoerbare bestand van LAMMPS. Hierdoor loopt een minimalisering van de energie om te verfijnen de H2O-configuratie, gevolgd door een simulatie van de FFMD uitgevoerd op constante aantal H2O moleculen (N), volume (V) en temperatuur (T) om de water tot de temperatuur van de simulatie, gevolgd door een simulatie van de FFMD lopen bij constante N, druk (P), en temperatuur (T) om te bepalen van de fysiek juiste hoogte van de doos van de simulatie. Output-bestanden die worden gebruikt in sectie 3 heten data.myadsorbate_npt en log.myadsorbate. Opmerking: De duur van de NPV-simulatie moet worden lang genoeg om een "evenwichtsinstelling" omvatten uitgevoerd, waar het volume van de supercell komt steady-state, en een "productie" uitvoert, die wordt gebruikt om te monster die het ensemble gemiddelden (hier, de hoogte van de supercell). Tijdens de duur van de evenwichtsinstelling, moet het volume van de supercell wanneer uitgezet tegen de tijd niveau af op de waarde van een steady-state. Zodra dit gebeurt, kan de NPV simulatie worden gezegd dat in de productieperiode. Controleer of evenwichtsinstelling van de NPV-simulatie door ervoor te zorgen dat de schommelingen in de hoogte van de supercell (lz) minimaal zijn, of hebben geconvergeerd naar een constante waarde. Als grote schommelingen optreden, wordt opnieuw gegenereerd zodra een H2O configuratie door de tijdstap op lijn 92 in het input.equil bestand en herhalende stap 2.2.8 minderen of opnieuw vanaf stap 2.1.1.</li>
	</ol>
</li>
</ol>3. Pak de juiste hoogte voor de supercell
<ol><li>Uitvoeren van de script-get_npt_lz.py op het log.myadsorbate-bestand door te typen get_npt_lz.py log.myadsorbate bij de opdrachtregelinterface. Dit script uitgangen de gemiddelde supercell hoogte uit het gedeelte van de "productie" uitvoeren van de simulatie van de NPV in het avg_lz.txt bestand. Opmerking: Het get_npt_lz.py script wordt ervan uitgegaan dat LAMMPS de lengte van de cel z-dimensie (lz) naar het log.myadsorbate bestand elke 1000 fs (aanpasbaar op regel 20 van het get_npt_lz.py script) schrijft, de standaardinstelling in het meegeleverde input.equil LAMMPS invoerbestand. Het get_npt_lz.py script detecteert en verwijdert de eerste 2 ns (aanpasbare op regel 19 van het get_npt_lz.py script) waarde van lz waarden in het bestand log.myadsorbate, als ze bestaan uit het gedeelte van de evenwichtsinstelling van de simulatie, terwijl de overige 3 ns omvatten de " productie"gedeelte en zijn aldus door het get_npt_lz.py script gebruikt voor het berekenen van de gemiddelde z-dimensie lengte. Naast het avg_lz.txt bestand uitgangen het get_npt_lz.py script een bestand met de naam npt_data.txt, waarin waarden van lz als een functie van de tijdstap, evenals een bestand genaamd npt_plot.png, die de percelen van dezelfde gegevens. Het perceel kan worden gebruikt om te controleren of de evenwichtsinstelling van de NPV-simulatie.</li>
	<li>Reconstrueren de supercell met behulp van de gemiddelde hoogte bepaald NPV.<ol>
<li>Kopieer het bestand data.myadsorbate_npt naar een nieuwe map, waarnaar hier als $WD2PATH, en data.myadsorbate de naam.</li>
		<li>Het nieuwe data.myadsorbate bestand zo bewerken dat de lz hoogte gelijk is aan de output van de gemiddelde waarde van het get_npt_lz.py script is door het veranderen van de argumenten van het zlo en zhi in het data.myadsorbate-bestand zodanig dat zlo 0.0 is en zhi de lz-waarde van het avg_lz.txt-bestand gemaakt in de s is TEP 3.1.</li>
	</ol>
</li>
</ol>4. Genereer configuraties van H2O moleculen
<ol><li>Kopieer de LAMMPS invoerbestand input.prod naar $WD2PATH. Bewerk de groep op regel 32 om aan te geven van de atoom type indexen voor het water zuurstof en water waterstofatomen en de variabele van de groep op lijn 33 om aan te geven van de atoom type indexen voor de Pt- en adsorbate-atomen.</li>
	<li>Uitvoeren van de LAMMPS-software door te typen mpiexec - np XX lmp_mpi < input.prod in de opdrachtregelinterface, waar XX is het aantal CPU-cores te gebruiken, en lmp_mpi is de naam van het uitvoerbare bestand van LAMMPS. Hierdoor zal het draaien van een constante NVT-simulatie op de H2O moleculen. De belangrijkste output file van deze simulatie is het dump.myadsorbate.lammpstrj-bestand. Opmerking: De duur van de NVT simulatie moet worden lang genoeg om te bestaan uit een evenwichtsinstelling uitgevoerd, waar de energie van het systeem komt om stabiele staat, en een productieperiode, waaruit het ensemble gemiddelden (hier, de ruimtelijke posities van de watermoleculen) zijn bemonsterd. Tijdens de duur van de evenwichtsinstelling, moet de energie van het systeem wanneer uitgezet tegen de tijd niveau af op de waarde van een steady-state. Zodra dit gebeurt, kan de NVT simulatie worden gezegd dat in de productieperiode.</li>
</ol>5. Bepaal de waterstofbrug levensduur voor bemonstering van de juiste tijd
<ol><li>Bewerk het script hb_lifetime_dist.py om op te geven: de tijdstap van het eerste frame van het dump.myadsorbate.lammpstrj-bestand door het veranderen van de variabele actualStart op regel 22, hoe vaak frames zijn weggeschreven naar het LAMMPS traject bestand door het veranderen van de tijdstap variabele op lijn 23, het eerste en laatste timesteps het script moet overwegen (dat wil zeggen, het productie gedeelte van het traject) door het veranderen van de variabelen van het N_first en N_last op lijnen 24 en 25, of opeenvolgende frames als worden beschouwd of frames worden overgeslagen door het veranderen van de nevery variabele op regel 26, en het aantal regels per frame sectie van het traject-bestand door het veranderen van de variabele frameLine op lijn 27. Bovendien, het bewerken van lijnen 31 door 35 om op te geven die atoomtypes binnen het data.myadsorbate-bestand behoren tot de adsorbate en die atoomtypes behoren tot de H2O moleculen. Opmerking: Het hb_lifetime_dist.py script analyseert de H2O configuraties in de productiegang en bepaalt of elke H2O moleculen waterstof gebonden aan de adsorbate zijn. Het dan telt het simulatie tijdstip waarop elke waterstof bond blijft intact en rapporteert deze informatie als een verdeling van de waterstofbrug levens in eenheden van ps. De specifieke versie van het script dat is voorzien van dit protocol wordt ervan uitgegaan dat LAMMPS de configuratie van H2O moleculen naar het dump.myadsorbate.lammpstrj bestand elke 1000 fs schrijft, de standaardinstelling in het meegeleverde input.prod LAMMPS invoerbestand. Het detecteert en verwijdert de eerste 2 ns waarde van configuraties in de de dump.myadsorbate.lammpstrj bestand, als ze bestaan uit het gedeelte van de evenwichtsinstelling van de simulatie en gebruik de resterende 3 ns voor het berekenen van de waterstofbrug levens.</li>
	<li>Uitvoeren van de script-hb_lifetime_dist.py op het dump.myadsorbate.lammpstrj-bestand door te typen hb_lifetime_dist.py in de opdrachtregelinterface. Hierdoor ontstaat een bestand met de naam distribution_HB_lifetime.dat.</li>
	<li>De gegevens in het bestand distribution_HB_lifetime.dat om weer te geven van de verdeling van de waterstofbrug levens die hebben plaatsgevonden tijdens de NVT simulatie uitzetten.</li>
	<li>Bepaal het tijdsinterval te gebruiken voor het tijdsinterval voor de bemonstering op basis van de berekende waterstofbrug levens. De beste keuze is de maximale waterstofbrug levensduur; u kunt ook kan een waarde die het betrouwbaarheidsinterval van 95 vangen zou % worden gebruikt.</li>
</ol>6. voorbeelden van configuraties van vloeibare H2O moleculen
<ol><li>Bepaal het aantal configuraties van de productiegang van het NVT FFMD traject voor verdere berekeningen. Het aantal configuraties moet worden ingeschakeld, zodat de minimale tijd tussen configuraties gelijk is aan of groter is dan het tijdsinterval voor de bemonstering geïdentificeerd in sectie 5.<ol>
<li>Redigeren naar de standaardwaarde voor de variabele num_frames op lijn 21 van het lammps_frames.py-script op te geven van het aantal configuraties om uit te pakken.</li>
		<li>Het script lammps_frames.py uitvoeren op het bestand dump.myadsorbate.lammpstrj door te typen lammps_frames.py in de opdrachtregelinterface. Hierdoor zal de uitgang van een lijst van simulatie tijden overeenkomt met de configuraties die dient te worden geëxtraheerd uit het dump.myadsorbate.lammpstrj-bestand. Deze configuraties kunnen worden gebruikt als structuren in AIMD of QM simulaties beginnen. Opmerking: 1) het lammps_frames.py script automatisch detecteert de LAMMPS log dumpbestanden alsmede het productie gedeelte van het traject binnen het crashdump-bestand en het aantal configuraties binnen het crashdump-bestand is verdeeld in 10 groepen. Anderzijds kan de gebruiker het logboekbestand, het crashdump-bestand en het aantal configuraties van de command-line interface met behulp van de -l -d, en -n opties, respectievelijk opgeven. Om dit te doen, de gebruiker moet typen lammps_frames.py - n XX -l $logfilename -d $dumpfilename op de opdrachtregelinterface, waar XX is het gewenste aantal configuraties, $logfilename is de naam van het logbestand van de LAMMPS en $dumpfilename is de naam van de LAMMPS traject (dump) bestand. De tijden van de simulatie die uitgang verwijzen naar de mediane tijd in elke groep. 2) als de configuraties worden berekend in VASP met de vlag van de LDIPOLE ingeschakeld, worden een kleine laag van vacuüm ruimte toegevoegd aan de bovenkant van de supercell boven het water laag. Dit zal de convergentie van de elektronische structuur in de berekening VASP vergemakkelijken. Het toevoegen van een extra 3 Å van vacuüm van de ruimte boven de H2O moleculen succesvol is geweest in de simulaties die hieronder besproken.</li>
	</ol>
</li>
</ol>

De auteurs vermelden geen belangenconflicten.

De methode zoals voorgesteld werd geselecteerd voor het gemak van implementatie, maar meerdere aanpassingen kunnen worden aangebracht. Voor een, kunnen de krachtvelden gebruikt in de FFMD-simulaties worden gewijzigd. Wijzigen van de force veldparameters en/of potentieel kan worden gedaan door de LAMMPS-invoer- en gegevens-bestanden te bewerken. Evenzo, oplosmiddelen dan H2O kon worden ingezet. Om deze verandering, de gewenste oplosmiddel molecule zou moeten worden ingevoegd vanaf stap 2.1.1 en de invoerbestanden LAMMPS zouden moeten worden bewerkt om de juiste mogelijkheden en parameters te nemen. Invoegen van de nieuwe oplosmiddel molecule zou moeten verstrekken de interne coördinaten van het oplosmiddel molecuul in een txt-bestand analoog aan het water.txt bestand.
Een andere wijziging die zou kunnen worden gemaakt is het wijzigen van het gebied van het oppervlak van de plaat. De resultaten besproken in dit manuscript werkzaam 3 Pt x 3 Pt of 4 Pt x 4 Pt oppervlakte platen, die minder dan 120 Å2oppervlakten. Als het oppervlak van de plaat verhoogt, verhoogt de computationele kosten ook. Computationele kosten heeft de grootste impact op punt 5 van dit protocol. Als de stappen van de gegevensverwerking in sectie 5 computationeel onbetaalbaar geworden, post grote gegevens processing strategieën zoals die wordt besproken in Li et al. 201845 kunnen worden ingezet.
Mogelijke bronnen van onzekerheid voor deze procedure zijn de krachtveld werkzaam, de bemonsteringsmethode en de bemonsteringsfrequentie. De structuur van water wordt bepaald door de kracht-veld dat wordt gebruikt, wat betekent dat de keuze van krachtveld invloed kan zijn op de specifieke configuraties van H2O moleculen. Onze fractie heeft onderzocht hoe de keuzes van krachtveld voor H2O moleculen en atomen van de Pt beïnvloeden de interactie energieën berekend in FFMD en vond dat de keuze van krachtveld minder dan 0.1 eV aan deze interactie-energie bijdraagt. Een andere bron van onzekerheid is de bemonsteringsmethode, die invloed op de specifieke configuraties die worden gebruikt voor het berekenen van een hoeveelheid van belang. Onze fractie heeft ten opzichte van de prestaties van de "tijd" bemonsteringsmethode gepresenteerd in dit protocol met een "energie" bemonsteringsmethode, die is bevooroordeeld naar lagere energie configuraties van H2O moleculen, op de interactie energieën berekend in DFT en vond zowel van deze bemonsteringsmethoden waarden geven statistisch gelijk35,46. De bemonsteringsfrequentie kan ook invloed hebben op de resultaten. We hebben beoordeeld hoe verhoging van het aantal configuraties van 10 aan 30.000 beïnvloedt de energieën van de gemiddelde interactie berekend in FFMD voor 40 verschillende C3HxO3 adsorbates en vond dat de sampling-frequentie minder bijdraagt dan 0.1 eV tot en met de gemiddelde interactie energie44.
De belangrijkste beperking van deze methode is dat de adsorbates zijn benaderd door de structuren onder drukvermindering tijdens de FFMD simulaties. In werkelijkheid, zou de adsorbates conformationele veranderingen (bond rek, hoek bochten, torsional bewegingen, enz.) als gevolg van normale thermische bewegingen, met inbegrip van interacties met oplosmiddel moleculen vertonen. Pogingen om conformationele veranderingen van adsorbates opgenomen in de FFMD-simulaties zou vereisen gedetailleerde ontwikkeling van krachtvelden voor katalytische oppervlakte adsorbates, dat wil zeggen, die bestaan uit van termen die de rek van de bond, hoek bochten en torsional voorwaarden beschrijven, onder anderen. Als een toekomstige richting van dit protocol, we zijn het ontwikkelen van dergelijke krachtvelden voor adsorbates op vaste oppervlakken, die we gebruiken zullen om te bepalen in hoeverre die met behulp van stijve adsorbates beïnvloedt de resultaten.

Modelleren van moleculaire verschijnselen die betrokken zijn bij heterogene katalyse vloeibare voorwaarden is noodzakelijk voor begrip katalytische functie; Dit blijft echter uitdagend omdat het vereist een delicaat evenwicht tussen chemische nauwkeurigheid en computationele kosten. In het algemeen, aangezien katalyse impliceert het breken en de vorming van chemische bindingen, moeten quantum mechanica worden gebruikt ten minste enigszins; echter, lange simulaties zijn uitdagend in de kwantummechanica, aangezien zij aanzienlijke computerbronnen vereisen. Aangezien moleculen in de vloeibare fase onder constante thermische beweging, simulaties omvat tevens configurationele bemonstering, dat wil zeggen, ze moeten meerdere ruimtelijke regelingen van de vloeibare moleculen, als elke andere ruimtelijke rangschikking (dat wil zeggen, elk configuratie) heeft een verschillende energie. Dit betekent dat meerdere configuraties van vloeibare moleculen moeten worden nagebootst voor elke katalytische soort van belang. Deze behoeften – kwantummechanica gebruiken en voor het uitvoeren van meerdere bewerkingen per katalytische soorten – kunnen renderen modelleren in heterogene katalyse onder vloeibare fase computationeel hardnekkig. Het doel van de hier beschreven methode is om een computationeel hanteerbare simulaties van verschijnselen in heterogene katalyse onder vloeibare fase.
We zijn vooral geïnteresseerd in ongelijkmatig gekatalyseerde reacties die worden uitgevoerd onder vloeibaar water. Watermoleculen hebben aanzienlijke invloed op de katalytische verschijnselen, zoals de interactie met katalytische soorten (bijvoorbeeld via dispersie krachten en waterstof binding)1,2,3,4,5 ,6,7,8,9,10,11,12,13,14 , 15 , 16 , 17 , 18 , 19 , 20 , 21 , 22 , 23, deelnemen aan de katalytische reacties1,7,8,9,15,21,22,24 ,25,26,27, en het beïnvloeden van reactie trajecten en/of katalytische tarieven1,11,12,15, 18,23,25,27,28,29,30,31. Modellering van deze verschijnselen is uitgevoerd met behulp van QM en/of ab initio moleculaire dynamica (AIMD)1,2,6,7,14,22 ,25,27,28,32,33,34, dwingen veld moleculaire dynamica (FFMD)35 , en kwantummechanica/moleculaire mechanica (QM/MM)10. In AIMD en FFMD, worden de atomen in het systeem verplaatst uit hoofde van Newtons bewegingsvergelijkingen volgens de krachten die op hen. In AIMD, worden de energie van het systeem en de krachten berekend met de kwantummechanica, overwegende dat in de FFMD, de energie van het systeem en de krachten zijn berekend met behulp van kracht velden, die algebraïsche uitdrukkingen zijn parameters die zijn gebaseerd op experimentele of QM gegevens. Het gedeelte van het systeem waar de band te breken en de vorming van plaatsvindt wordt berekend met QM in QM/MM, en de rest van het systeem wordt berekend met MM, waarbij krachtvelden in dienst. Omdat ze rechtstreeks in QM dienst, AIMD en QM/MM zijn beter geschikt voor het vastleggen van het breken van de bond en vormen die optreedt in de waterfase heterogene katalyse; FFMD is echter aanzienlijk meer rekenkundig hanteerbare en dus beter geschikt voor het genereren van de configuraties van vloeistof H2O moleculen. De methode die in dit protocol gepresenteerd saldi chemische nauwkeurigheid en computationele kosten door gebruik te maken van een combinatie van QM en FFMD.
Deze methode gebruikt met name FFMD simulaties voor het genereren van configuraties van vloeistof H2O en QM voor het berekenen van de energie van het systeem. FFMD wordt uitgevoerd met behulp van LAMMPS. 36 de krachtvelden gebruikt in FFMD in dezen dienst Lennard-Jones + Coulomb (LJ + C) mogelijkheden, waar de LJ parameters zijn overgenomen uit de TIP3P/CHARMM model37 voor H2O, de universele krachtveld38 (UFF) voor Pt, en de OPLS-AA krachtveld39 voor katalytische soorten, en de Coulomb parameters zijn overgenomen uit de TIP3P/CHARMM37 model voor H2O en de OPLS-AA krachtveld39 voor katalytische soorten. De Coulomb parameters voor Pt atomen zijn ingesteld op 0. QM berekeningen worden uitgevoerd met behulp van de VASP code40,41,42, dat een dichtheid functionele theorie (DFT) code is. Water molecuul invoegingen zijn uitgevoerd met een code die is ontwikkeld in-house genaamd Monte Carlo Plug-in voor Quantum methoden (MCPliQ). Bestandsconversies van VASP aan LAMMPS in dit protocol worden uitgevoerd met de visuele Molecular Dynamics (VMD) software43.
Het protocol is bedoeld voor het genereren van configuraties van vloeibaar watermoleculen rond katalytische soorten op vlakke overgangsmetalen oppervlakken bij lage dekking. Dekking is aangeduid met θ en gedefinieerd als het aantal adsorbates per oppervlakte metaal atoom (dat wil zeggen, het aantal oppervlakte adsorbates genormaliseerd door het aantal metalen atomen in de bovenste laag van de metalen slab in de katalysator-model). In dit manuscript, lage dekking wordt gedefinieerd als θ ≤ 1/9 enkelgelaagde (ML), waar 1 ML een katalytische soorten per oppervlakte metaal atoom betekent. De modellen van de katalysator moeten worden geplaatst in de vakken van de periodieke simulatie. De simulatie vakken hoeft niet te worden van de kubussen. Dit manuscript toont de toepassing van het protocol voor het genereren van configuraties van vloeibare H2O die kunnen worden gebruikt voor het berekenen van de hoeveelheid belangstelling voor heterogene katalyse van de waterige fase.
Dit protocol vereist dat de gebruiker toegang tot geïnstalleerde en werkende versie van de software VASP, MCPliQ, LAMMPS en VMD heeft. Meer informatie over VASP (https://www.vasp.at/), LAMMPS (https://Lammps.sandia.gov/) en VMD (https://www.ks.uiuc.edu/Research/vmd/) zijn beschikbaar op hun respectieve websites. De software van de MCPliQ wordt beschreven op https://github.com/getman-research-group/JoVE_article, samen met alle invoerbestanden en Python scripts vermeld in dit protocol. Dit protocol wordt ervan uitgegaan dat de uitvoerbare bestanden en scripts vermeld binnen zal worden uitgevoerd op de computer van een hoogwaardig onderzoek en zijn geïnstalleerd in een map in de gebruikershandleiding $PATH variabele. Als een uitvoerbaar bestand of script is geplaatst op een locatie die niet in de user's $PATH, en vervolgens het pad naar het uitvoerbare bestand opgenomen worden moet om het uit te voeren. Uitvoerbare bestanden en scripts worden uitgevoerd in stappen 2.1.2, 2.2.1 2.2.8, 3.1, 4.2, 5.2 en 6.1.2. Bijvoorbeeld, voor het uitvoeren van de code MCPliQ in stap 2.1.2 vanuit een directory die zich niet in de gebruiker de $PATH, de gebruiker zou typt $PATHTOMCPLIQ/mcpliq op de opdrachtregelinterface in plaats van mcpliq, waar $PATHTOMCPLIQ de locatie is waar de mcpliq uitvoerbaar bestand is opgeslagen (bijvoorbeeld $PATHTOMCPLIQ misschien wel ~ / bin). Voordat u begint van dit protocol, alle uitvoerbare bestanden en scripts moeten worden gegeven uitvoerbare permissies (bijvoorbeeld in Linux, dit kan gebeuren door te typen chmod + x mcpliq bij de opdrachtregelinterface van de map waarin de uitvoerbare mcpliq is opgeslagen). Verder modules die nodig zijn door de software of scripts moeten worden geladen (deze afhankelijkheden worden specifieke aan afzonderlijke installaties van de verschillende software en de computer waarop de simulaties zullen worden uitgevoerd).

<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

Een groot aantal chemische processen ongelijkmatig-gekatalyseerde voordoen in vloeibare, maar simuleren katalysator functie onder dergelijke omstandigheden is een uitdaging als het moet de oplosmiddel moleculen worden opgenomen. Het breken van de bond en de vorming van processen die gemodelleerd in deze systemen vereisen het gebruik van quantum chemische methoden. Aangezien moleculen in de vloeibare fase onder constante thermische beweging, omvat simulaties tevens configurationele bemonstering. Dit betekent dat meerdere configuraties van vloeibare moleculen moeten worden nagebootst voor elke katalytische soort van belang. Het doel van het hier gepresenteerde protocol is te genereren en te proeven van de trajecten van configuraties van vloeibaar watermoleculen rond katalytische soorten op vlakke overgangsmetalen oppervlakken op een manier die chemische nauwkeurigheid met computationele kosten balanceert. Specifiek, worden krachtveld moleculaire dynamica (FFMD) simulaties gebruikt voor het genereren van configuraties van vloeibare moleculen die vervolgens kunnen worden gebruikt in de kwantummechanica gebaseerde methoden zoals dichtheid functionele theorie of ab initio moleculaire dynamiek. Om dit te illustreren, in dit manuscript, wordt het protocol gebruikt voor de katalytische tussenproducten die kunnen worden betrokken in het traject voor de ontleding van glycerol (C3H8O3). De structuren die worden gegenereerd met behulp van FFMD worden gemodelleerd in DFT om te kunnen inschatten van de enthalpie van was de katalytische soorten en om te bepalen hoe de H2O moleculen in katalytische decomposities deelnemen.

Een gebruik van dit protocol is voor het berekenen van de energieën van interactie tussen de vloeibare water en katalytische soorten, dat wil zeggen, ΔEint35:
∆Eint=Ekatalytische soorten + H2O+Eschone katalysator oppervlak-Ekatalytische soorten-Eschone katalysator oppervlakte + H2O
waar Ekatalytische soorten + H2Ode energie van is is een configuratie van H2O moleculen rond een katalytische soorten op een metalen oppervlak, Eschone katalysator oppervlak de energie van het oppervlak schoon katalysator in vacuüm, E Katalytische soorten is de energie van de katalytische soorten op een metalen oppervlak in vacuüm, en Eschoon katalysator oppervlakte + H2O is de energie van de configuratie van H2O over het oppervlak van de katalysator met de katalytische soorten verwijderd. De posities van de H2O moleculen gebruikt om te berekenen Ekatalytische soorten + H2O en Eschoon katalysator oppervlakte + H2O moet identiek zijn. Alle waarden van E worden berekend met behulp van de code VASP. De hoeveelheid ΔEint omvat alle van de fysische en chemische interacties tussen de moleculen in het vloeibare water-structuur en de katalytische soort en geeft een redelijke schatting van de enthalpie van het oplossen van de katalytische soorten, die nodig is voor het berekenen van de vrije energie van was en totaal vrije energie. Tabel 1 geeft de waarden voor ΔEint berekend voor soorten op een Pt(111) katalysator oppervlak met chemische formules gelijk is aan CxHyOz in eenheden van eV (1 eV = 96.485 kJ/mol). De waarden werden berekend op dekkingen ≤1/9 ML.35,46 , de vermelde waarden zijn de gemiddelden overgenomen van 10 configuraties van vloeibare H2O en de onzekerheden zijn gerapporteerd als standaarddeviaties. Alle waarden zijn negatief, met vermelding van gunstige interacties met water.
Een andere toepassing van dit protocol is het genereren van startende structuren voor AIMD. Filmpje 1 is een film van een AIMD-traject dat is gestart vanuit een configuratie gegenereerd door dit protocol. Aan het begin van deze film, wordt een COH adsorbate weergegeven op een Pt(111) ondergrond onder een structuur van de vloeibare H2O. Een H2O molecuul wordt benadrukt, die vormden een waterstofbrug met COH. In de loop van de film, dit H2O molecuul abstracts van de proton uit de COH adsorbate en deposito's van een tweede waterstofatoom op het oppervlak van de Pt(111). De H2O molecuul helpt dus om het katalyseren van de reactie COH * + * → CO * + H *, waar de * s geven katalytische sites. Deze simulatie belicht de belangrijkste sterkte en het hoofddoel van de multiscale bemonsteringsmethode hierin beschreven. Vele configuraties van H2O moleculen worden gegenereerd met FFMD, als gevolg van zijn kracht in computationele werkwillig. Een beperking van de FFMD is echter dat het niet kan vangen bond breken en de vorming van tenzij een reactieve krachtveld ten uitvoer wordt gelegd. AIMD kwantummechanica gebruikt voor het berekenen van de energie en dus bond breken en de vorming van kunt vastleggen. AIMD eist echter ook rekenkundig voor het genereren van alle van de formaties van H2O moleculen nodig zijn om voldoende bemonstering is bereikt. Dus, dit protocol combineert de twee methoden.
De structuur van de vloeibare H2O moleculen die worden gegenereerd door deze procedure zijn afhankelijk van de input-instellingen. Deze ten onrechte instellen kan hebben onbedoelde invloeden op de water-structuren. Bijvoorbeeld, kan wanneer de intermoleculaire afstanden te klein geworden of andere parameters in de moleculaire dynamica invoerbestanden onjuist zijn ingesteld of unphysical waarden, de structuur van water worden onredelijk. Onder deze omstandigheden zal de structuur van water "opblazen" unintendedly tijdens het FFMD traject. Figuur 1 toont een voorbeeld van dit. De momentopname aan de linker kant is de startende structuur voor een run van de FFMD, en de momentopname aan de rechterkant is een momentopname genomen binnen 1 ps de simulatie te starten. Zoals kan worden gezien, hebben de H2O moleculen verplaatst ver weg van het oppervlak. Dit wordt veroorzaakt door onjuiste instellingen in de input-bestanden van de simulatie en is niet een structuur die dreigt te ontstaan in werkelijkheid.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/59284/59284fig1v2.jpg"> Figuur 1: Voorbeeld van een negatief resultaat. De simulatie van de moleculaire dynamica krachtveld blaast"" als gevolg van een unphysical instelling of waarde. Afbeelding van de linkerhand: de startende geometrie van het oppervlak Pt(111), adsorbate en vloeibaar water structuur. Afbeelding van de rechterhand: de geometrie van het oppervlak van de Pt(111), adsorbate en vloeibaar water structuur minder dan 1 ps later. In de rechter afbeelding, hebben de H2O moleculen gescheiden van het oppervlak door unphysically grote krachten. <a href="https://www.jove.com/files/ftp_upload/59284/59284fig1v2large.jpg" target="_blank">Klik hier voor een grotere versie van dit cijfer.</a>
<img alt="Movie 1" src="/files/ftp_upload/59284/59284mov1.jpg"> Filmpje 1: Ab initio moleculaire dynamica (AIMD) simulatie gestart vanuit een configuratie gegenereerd in multiscale bemonstering. Een H2O molecuul dat oorspronkelijk is waterstof gebonden aan een COH adsorbate op een Pt(111) ondergrond abstracts van de proton uit COH en deposito's van een tweede waterstof op het oppervlak van de Pt(111). Deze band te breken en de vorming van de gebeurtenis kan worden vastgelegd door AIMD maar niet met krachtveld moleculaire dynamica (FFMD) tenzij een reactieve krachtveld wordt gebruikt. De initiële configuratie van H2O moleculen gebruikt in deze simulatie AIMD werd gegenereerd met FFMD zoals beschreven in dit manuscript. <a href="https://www.jove.com/files/ftp_upload/59284/movie1.mp4" target="_blank">Gelieve Klik hier om deze video te bekijken. (Klik met de rechtermuisknop om te downloaden.)</a>
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always"><tbody>
<tr>
<td>Katalytische soorten</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>CH2OH</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>Tabel 1: Water-katalytische energie resultaten van de soorten interactie. Interactie energieën in eV berekend voor acht CxHyOz adsorbates op Pt(111). Waarden gemeld worden de gemiddelden genomen over meerdere configuraties van vloeistof H2O. De onzekerheden zijn de standaarddeviaties van de gemiddelden. 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

Het doel van het hier gepresenteerde protocol is te genereren en te proeven van de trajecten van configuraties van vloeibaar watermoleculen rond katalytische soorten op een vlakke overgangsmetalen oppervlak. De bemonsterde configuraties kunnen worden gebruikt als startende structuren in de kwantummechanica gebaseerde methoden.

Multiscale bemonstering van een heterogene katalysator-Interface voor Water/metaal is met behulp van dichtheid functionele theorie en krachtveld moleculaire dynamica

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

De protocolmodellen, thermische en kwantumverschijnselen in vloeibare fase heterogene katalyse. Het is voor zover wij weten de kwantummechanica met volledige bemonstering van een expliciete vloeibare omgeving. Configuraties van vloeibare moleculen gegenereerd uit dit protocol vertegenwoordigen die verwacht onder de werkelijke reactieomstandigheden en kunnen worden gebruikt om moleculair niveau verschijnselen die afhankelijk zijn van de ruimtelijke regelingen van moleculen te verkennen.De configuraties van vloeibare moleculen gegenereerd door dit protocol geven inzicht in de rollen die oplosmiddel speelt op vloeibare fase heterogene katalyse. Als u dit protocol voor de eerste keer probeert, raad ik u aan eerst ervoor te zorgen dat u toegang hebt tot en de VASP-, MCPliQ-, VMD- en LAMMPS-software uitvoeren. Het aantonen van de procedure zal Tianjun Xie zijn, een afgestudeerde student van mijn laboratorium.Nadat u de adsorbatestructuur hebt gegenereerd, genereert u de LAMMPS-invoerbestanden voor een NPT-simulatie en 7000 het celvolume met FFMD. Kopieer het LAMMPS-invoerbestand naar de werkmap. Bewerk de groepsvariabele op lijn 34 om de indexen van het atoomtype voor de waterzuurstof- en waterwaterstofatomen aan te geven.In de groepsvariabele op lijn 35, om de indexen van het atoomtype voor de platina- en adsorbateatomen aan te geven. Bewerk op lijn 17 van het invoerbestand de run stap variabele om de duur van de NPT-simulatie in te stellen die lang genoeg is om een equilibratierun en een productierun te omvatten. Voer de LAMMPS-software uit door de opdracht te typen op de command-line-interface, die de informatie bevat van het aantal CPU-kernen dat moet worden gebruikt en de naam van de LAMMPS die kan worden uitgevoerd.De energieminimisatie verfijnt de watermolecuulconfiguratie en wordt gevolgd door een FFMD-simulatie uitgevoerd bij constant aantal watermoleculen, volume en temperatuur om het water op de simulatietemperatuur te brengen. Een andere FFMD simulatie wordt vervolgens uitgevoerd op een constant aantal watermoleculen, druk en temperatuur om de fysiek juiste hoogte van de simulatiebox te bepalen. De uitvoerbestanden worden later gebruikt.Na de NPT-simulatie, plot de hoogte van de supercel tegen de tijd. Het punt waar het niveaus af tot een steady state waarde is het punt in de NPT simulatie wanneer de productie kan beginnen. Controleer de balans van de NPV-simulatie door ervoor te zorgen dat de schommelingen in de hoogte van de supercel minimaal zijn of zijn geconvergeerd tot een constante waarde.Als er grote schommelingen optreden, opent u de invoer. equil-bestand en verlaag de timestep op lijn 92 om de watermolecuulconfiguratie te regenereren en de LAMMPS-software uit te voeren zoals voorheen. Typ om te beginnen bij de opdrachtregelinterface om het script uit te voeren.Dit script voert de gemiddelde supercelhoogte uit van het productieuitvoergedeelte van de NPT-simulatie naar een TXT-bestand. Het script leest de lengte van de z-dimensie van de cel op 1.000 femtosecondeintervallen, het standaardinterval voor het afdrukken van informatie in LAMMPS. Als een ander afdrukinterval gewenst is, kan dit worden gewijzigd door regel 20 van de get_npt_lz te bewerken.py script en lijn 16 van LAMMPS input. equil-bestand. Het script detecteert en verwijdert de eerste twee nanoseconden ter waarde van lz-waarden, omdat ze het evenwichtsgedeelte van de simulatie omvatten.De duur van de equilibratierun kan worden gewijzigd door regel 19 van het bestand te bewerken. De overige drie nanoseconden omvatten het productiegedeelte en worden dus gebruikt om de gemiddelde lengte van de z-dimensie te berekenen. Bovendien wordt in het script een ander TXT-bestand uitgevoerd dat waarden van lz biedt als functie van de timestep en een PNG-bestand, waarin dezelfde gegevens worden geploimd.Het perceel kan worden gebruikt om de balans van de NPV-simulatie te verifiëren. Als u de supercel wilt reconstrueren met behulp van de gemiddelde hoogte die in NPT is bepaald, kopieert u het eerder gegenereerde gegevensbestand in een nieuwe werkmap en wijzigt u deze als data.myadsorbate. Bewerk vervolgens het nieuwe gegevensbestand om zlo te wijzigen in 0,0 en zhi naar de lz-waarde van de gemiddelde waardeuitvoer in het TXT-bestand.Kopieer het LAMMPS-invoerbestand naar de nieuwe werkmap. Bewerk de groepsvariabele op lijn 32 om de indexen van het atoomtype voor de waterzuurstof- en waterwaterstofatomen en de groepsvariabele op lijn 33 aan te geven om de indexen van het atoomtype voor de platina- en adsorbateatomen aan te geven. Bewerk vervolgens op lijn 16 de runStep-variabele zodat deze lang genoeg is om een equilibratierun en een productierun te omvatten.Typ de opdracht om LAMMPS uit te voeren in de command-line interface om de LAMMPS-software uit te voeren. Dit zal een constante NVT simulatie op de watermoleculen draaien, en het belangrijkste outputbestand wordt gegenereerd. De NVT simulatie bestaat uit een evenwichtsgedeelte en een productiegedeelte.Het productiegedeelte begint wanneer de energie van het systeem uitgezet tegen tijdsniveaus af. Om te beginnen, open de waterstof bond levensduur script. Wijzig de werkelijke beginvariabele op lijn 22 om de tijdsstap van het eerste tijdsframe in te stellen.Wijzig de tijdstapvariabele op lijn 23 om in te stellen hoe vaak frames worden geschreven naar het LAMMPS-trajectbestand. Wijzig N_first- en N-laatste variabelen op lijn 24 en 25 om de eerste en laatste timesteps in te stellen en wijzig de nevery-variabele op lijn 26 om in te stellen of opeenvolgende frames worden beschouwd of overgeslagen. Stel het aantal lijnen per framesectie van het trajectbestand in door de frameregelvariabele op lijn 27 te wijzigen.Bovendien, bewerken lijnen 31 tot en met 35, om op te geven welke atoomtypen binnen de gegevens. myadsorbate bestand behoren tot de adsorbate, en welke atoomtypes behoren tot de watermoleculen. Het script analyseert de waterconfiguraties in de productie en bepaalt of watermoleculen aan de adsorbaat zijn gebonden.Het telt dan de simulatietijd dat elke waterstofbinding intact blijft en rapporteert deze informatie als een verdeling van waterstof-binding levensduur in eenheden van picoseconden. LAMMPS schrijft de configuratie van watermoleculen om de 1.000 femtoseconden in het bestand, wat de standaard is in het meegeleverde LAMMPS-invoerbestand. Het detecteert en verwijdert de eerste twee nanoseconden ter waarde van configuraties in het bestand als ze het evenwicht gedeelte van de simulatie en maakt gebruik van de resterende drie nanoseconden om waterstof-obligatie levensduur te berekenen.Als u het script wilt uitvoeren, typt u de opdrachtregelinterface. Het genereert dan een DAT-bestand. Plot de gegevens in het bestand om de verdeling van waterstof-obligatie levens die zich hebben voorgedaan tijdens de NVT simulatie te bekijken.Gebruik een tijdstoename die groter is dan of gelijk is aan de maximale levensduur van waterstofbinding om de tijdstoename voor het tijdsbemonsteringsinterval te bepalen. Bepaal het aantal configuraties uit de productierun van het NVT FFMD-traject om zodanig te bemonsteren dat de minimale tijd tussen configuraties gelijk is aan of groter is dan het eerder vastgestelde tijdsbemonsteringsinterval. Bewerk in het eerder geschreven frameextractiescript de standaardwaarde voor het aantal frames dat variabel is op lijn 21 om het aantal configuraties op te geven dat moet worden geëxtraheerd.Als u het script wilt uitvoeren, typt u de scriptnaam op de opdrachtregelinterface. Hiermee wordt een lijst uitgevoerd met simulatietijden die overeenkomen met de configuraties die uit het NVT-simulatiebestand moeten worden geëxtraheerd. Deze configuraties kunnen worden gebruikt als startstructuren in AIMD- of QM-simulaties.In deze procedure werd FFMD gebruikt om de initiële configuratie van watermoleculen te genereren. De AIMD simulatie toont aan dat een watermolecuul dat oorspronkelijk waterstof-gebonden aan een suiker alcohol adsorbaat op een platina-111 oppervlak onttrekkingen van de waterstof uit de alcohol adsorbaat en deposito's een tweede waterstof op het platina-111 oppervlak. De structuren van vloeibare watermoleculen zijn afhankelijk van invoerinstellingen.Het verkeerd instellen van deze kan onbedoelde invloeden hebben op de waterstructuren. In dit cijfer is de linkerkant de startstructuur voor een FFMD-run. En de rechterkant is binnen een picoseconde van het starten van de simulatie.De FFMD simulatie ontploft als gevolg van unphysically grote kracht instellingen, waardoor watermoleculen te bewegen ver weg van het oppervlak. De configuraties kunnen worden gebruikt in de kwantummechanica, of QMM simulatie, of ze kunnen worden gebruikt om statistieken met betrekking tot de ruimtelijke posities van de moleculen te analyseren. Deze techniek maakt de weg vrij voor onderzoekers om de rollen te verkennen die vloeibare reactieomgevingen hebben op katalyse door het genereren van werkelijke configuraties van vloeibare moleculen op katalytische interfaces.

Research

JoVE Journal

Chemistry

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

Layer-by-layer synthese en de overdracht van Vrijstaande geconjugeerde Microporous 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-1 als heterogene katalysator voor de synthese van vanilline

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 Synthese van Worm-micellen door het zichtbare licht Mediated Dispersion polymerisatie onder toepassing van 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

Synthese and Characterizations van de prestaties van overgangsmetalen enkel atoom katalysator voor elektrochemische CO2 reductie

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

Voorbereiding van functionele Silica met behulp van een methode 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

Heterogene verwijdering van in water oplosbare Ruthenium olefine Metathese katalysator uit waterige Media Via Host-Guest interactie

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 corrosie bij de metaal verf-interface met behulp van tijd-of-Flight secundaire Ionenmassa spectrometrie

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: een algemene en in grote lijnen toepasbare koper-gebaseerde fotoredox katalysator

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

Vibrationele Spectra van een N719-Chromophore/Titania Interface van Empirisch-Potential Molecular-Dynamics Simulation, oplosbaar door een ionische vloeistof van kamertemperatuur

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

Ontwikkeling van heterogene enantioselectieve katalysatoren met behulp van chirale metaal-organische kaders (Mof's)

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