predictive modelling: a view from the atomic level - warwick...pvmid pvep ppa pvmp ap-4m2e ap-2m3et...
TRANSCRIPT
![Page 1: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/1.jpg)
Predictive modelling: a view from the atomic level
P.M. Rodger
Department of Chemistry
& Centre for Scientific Computing
University of Warwick
![Page 2: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/2.jpg)
Acknowledgements
• David Quigley, Igor Khovanov, Matt Bano, Yuriy Bushuev, Aaron Finney, Ritchie-Mae Gamot, Yuanwei Xu, Salvatore Cosseddu (Warwick)
• John Harding, Colin Freeman and Riccardo Innocenti Malini (Sheffield)
• Paul Smeets, Heiner Friedrich, Jozua Laven, Nico Sommerdijk(Eindhoven), Wouter Habraken (Postdam), Fabio Nudelman (Edinburgh)
• Mike Allen (Warwick), Dorothy Duffy (UCL), Julian Gale (Curtin)
• Warwick CSC
• HECToR & ARCHER
• MidPlus
• EPSRC (mIb: Materials Interfacing Biology)
![Page 3: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/3.jpg)
Philosophy
• Complex materials
complete model is not possible
• Mechanistic information
“predictions” not reliable without the correct underlying physics/chemistry
• Dimensionality
identify key degrees of freedom
target analysis
drive or simplify modelling
![Page 4: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/4.jpg)
Unexpected Mechanisms: old examples
•
•
![Page 5: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/5.jpg)
Motivation: Biominerals
• Nature:
Exquisite control of crystal morphology & aggregation found in nature
Leads to very well “designed” functional materials
Ability to simulate implies ability to guide synthesis Henriksen, K., S. L. S. Stipp, et al.
American Mineralogist 89, 1709-1716
(2004)
Emiliania huxleyi
coccoliths
Reddy, M. M. and A. R. Hoch.
Journal of Colloid and Interface
Science 235, 365-370 (2001)
• Laboratory:
![Page 6: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/6.jpg)
Motivation: additives for inhibition
• Often need to suppress crystal growth
Scale, wax, hydrate
• “kinetic” inhibitors
Delay nucleation or slow growth
Active at low concentrations
Require molecular understanding of nucleation and growth
![Page 7: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/7.jpg)
• Inhibition by surface adsorption
Mechanistically inspired prediction?
-22.0
-21.0
-20.0
-19.0
-18.0
-17.0
-16.0
-15.0
-14.0
-13.0
-12.0
-11.0
-10.0
-22.0-21.0-20.0-19.0-18.0-17.0-16.0-15.0-14.0-13.0-12.0-11.0-10.0
sH (1 0 0) / kcal mol-1
T2
(1
1 1
) / k
ca
l m
ol-1
PVP
PVMFPVMEP2 PVBA
PVI
PVMLO
PBA
AP-23DMt
AP-3E
PVCap
AP-3M
AP-3M4Ec
AP-4M3Et
AP-2M3Ec
AP-3M2Ec
PAFPVPy
PAPe
PAPPVMID
PVEP
PPA
PVMP
AP-4M2E
AP-2M3Et
PVPK
PVDMP
AP-24DM
AP-2M4E
AP-23DMc
PMA
PVDMPc PDHU
PBCOD
AP-3M2Et
PBL
AP-4M3Ec
PCAP
validate screen
new modes?
![Page 8: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/8.jpg)
Synthesise and test
3500
4000
4500
5000
5500
6000
6500
7000
7500
8000
Q5
Q4P
1
Q3P
2
Q2.
5P2.
5
Q2P
3
Q1P
4 P5
Contro
l
Concentration / (1000 ppm)
Ind
uctio
n T
ime
/s
THF Hydrate Tests
induction times for 0.5% mixtures
of tba3S (the “quat”, or Q) and
PVP (P). The Control contains no
inhibitor. Bars indicate one
standard deviation
tinduction / s Rnucl /(ml
-1 min
-1)
Rgrowth /(m min
-1)
Control -90
(32)
1439
(1031)
1.56
(2.18)
tba3S
(0.1 %)
1705
(422)
89
(–)
0.61
(0.23)
tba3S
(0.5 %)
2678
(231)
116
(13)
0.86
(0.25)
Ethane Hydrate Tests
Induction times, nucleation
rates (Rnucl) and growth rates
(Rgrowth). Standard deviations
are given in parentheses
![Page 9: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/9.jpg)
2nd generation
• 2nd generation inhibitors perform about 4 times better than current commercial inhibitor
Test KI Gas Test P/T Subcooling Induction
No. psia/C /C Time/ hrs
1 JI-C003:1 NG 1023 / 5.3 10.7 5
2 JI-C003:1 NG 1023 / 5.3 10.7 >67
3 JI-C003:1 NG 1015 / 5 11 >67
4 JI-C003:1 NG 1023 / 5 11 >150
Commercial NG 10 16
6 JI-C002:1 NG 1095 / 4.5 12.3 3
7 JI-C002:1 NG 1095 / 4 12.8 40
8 JI-C002:1 NG 1095 / 5 11.8 >42
Commercial
5 JI-C003:1 Methane 1700 / 4.5 10 >150
Methane 10 5
![Page 10: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/10.jpg)
But!
• Experiments not repeatable! Oil companies can’t repeat each other’s screens
Heriot-Watt: activity changed with test-cell stirrer design
Scale of activity with new compound by test-site
• Heriot-Watt > Toulouse > Halliburton > Heriot-Watt
• Similar story with wax inhibition (model oil using well characterisededible oil mixture)
• Experimentalists working with nanotoxicology want legislative standards based on computational tests because experiments are too irreproducible!
• Simplistic model: dynamic interface; nucleation; …
![Page 11: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/11.jpg)
Nucleation from MD?
• Homogeneous nucleation?
Experimental nucleation rates < O(106) nuclei per cm3 per second
MD simulation 10,000 molecules:
• 1 nucleation event every 32 millenia!
• Heterogeneous nucleation?
(some) success with seeding
(some) success with interfaces
![Page 12: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/12.jpg)
Standard MD: modelling surfaces
• Deposition of wax on Fe2O3
C28 from C7
![Page 13: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/13.jpg)
Standard MD: crystal nucleation
• heterogeneous and/or high supersaturation/subcooling
![Page 14: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/14.jpg)
Polydispersity in inhibitors
• molecular weight
amount of hydrate present
PVP
• molecular shape
growth of hydrate clusters
pDMAEMA
0
0.03
0.06
100 200 300 400 500
cluster size
clu
ste
r f
req
ue
nc
y
1ns
2ns
3ns
4ns
![Page 15: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/15.jpg)
Nucleation more generally
• Project onto “important” manifold
Define key dimensions (order parameters / reaction coordinates / collective variables)
Project analysis onto these dimensions
Bias simulations to explore these dimensions
• constrain to portion of manifold (umbrella sampling)
• trap fluctuations that explore the manifold (adaptive bias force)
• disfavour current portion of manifold (Wang-Landau; metadynamics)
![Page 16: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/16.jpg)
Free energy hypersurfaces
•
peptide folding
Bussi, …, Parrinello
3-basin model
Quigley, PMR
KcsN ion channel
Cosseddu, PMR, Khovanov
![Page 17: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/17.jpg)
Limitations and Pitfalls: dimensionality!
• Equilibrium / comprehensive sampling
Typically 1–2D(≤ 6; can do ~100 for basinescape & rough sampling)
projection mergesbasins
• Adiabatic Surfaces
dynamics defined by metastable substates?
KcsN
Potential of mean force
Kramers potential
![Page 18: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/18.jpg)
Nucleation with metadynamics: ice
• 4D bias:
Q4, Q6, ς and potential energy (simaltaneous)
• NPT simulation
density change is spontaneous
• T = 180 K; ca. 11 “ns” simulation
• Results independent of periodicboundaries
![Page 19: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/19.jpg)
N = 147
N = 138
N = 157
N = 161
An ensemble of critical nuclei
ice Icice Ihboth
![Page 20: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/20.jpg)
Nucleation with metadynamics: CaCO3
• 6 order parameters
Q4 for: Ca-Ca, Ca-C, Ca-O, C-C, C-O
• measures orientations of X about Y
energy associated with CaCO3
• ca. 10 ns of MD
Exploring configurationsfor 75 units (water notshown)
![Page 21: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/21.jpg)
Calcium Carbonate Biominerals
• Multi-stage hierarchical formation
•
Colfen & Mann
Angew. Chem, 2003
![Page 22: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/22.jpg)
Early stages of carbonate biomineralisation
ions & ion pairs hACC ACC crystalline
hACC: hydrated amorphous calcium cabonatecomposition ca. CaCO3.H2O
ACC: (anhydrous) amorphous calcium carbonate
Crystalline: calcite, aragonite, vaterite
ions & ion pairs hACC ACC crystalline
![Page 23: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/23.jpg)
Free energy maps: 75 CaCO3 units• NVT
• favours amorphous
DQ and PMR, J. Chem. Phys. 128 (2008) 221101
![Page 24: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/24.jpg)
Free energy maps: 75 CaCO3 units
• NPT
• Dominated by calcite
very low barrier (10–20 kT)
DQ, CLF, JHH and PMR, J. Chem. Phys. 134 (2011) 044703
![Page 25: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/25.jpg)
300 CaCO3 Units in water
• two stable states
• calcite (dominant)
• amorphous
• energy barrier > 102 kT
• multiple transitions
![Page 26: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/26.jpg)
Biomineralisation: egg shells
![Page 27: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/27.jpg)
Structure of an eggshell
![Page 28: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/28.jpg)
G-type lectin-type proteins
• found in
Rhea
Duck
Chicken
Emu
Ostrich
![Page 29: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/29.jpg)
Ovocleidin-17
Reyes-Grajeda, Moreno, Romero;
J. Biol. Chem., 2004
![Page 30: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/30.jpg)
metaDynamics of OC17 + CaCO3• 192 & 300 unit nanoparticles
22,000 water molecules
• 20 different protein/nanoparticle orientations
• Select “best” 4 for long meta-dynamics
• Potentials: due to Freeman et al.
Pavese CaCO3 ; Amber protein; Tip3p water;
cross terms derived from crystals with scaled charges
![Page 31: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/31.jpg)
Free energy landscapes
• 192 units of CaCO3 in water
Initial
metadynamics
Crystallised
• 192 units of CaCO3 bound to Ovocleidin-17 in water
![Page 32: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/32.jpg)
Larger nonaparticles?
• 300 CaCO3 unit particles don’t stay bound on crystallisation
• Mechanism for proteins
Bind to small nanoparticles
Facilitate transformation to calcite
Desorb as crystal grows
• Gives catalytic cycle for polycrystalline mamillary layer
![Page 33: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/33.jpg)
OC17: catalytic cycle?
Cryst.
grow
desorb
adsorb
![Page 34: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/34.jpg)
Surface binding & structured water
• Binding is seen to planar surfaces
Structured water dominates adsorption energy on large crystalline surfaces
Calcite nanoparticles and ACC do not give structured water layer; leads to “flat” adsorption
![Page 35: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/35.jpg)
Also seen with AuBP1 on gold
(100) 5x1
(100)
(111)
![Page 36: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/36.jpg)
hACC nulcleation
• Problem
intermolecular potentials
• Controversy
![Page 37: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/37.jpg)
Aragonite?
• needs a better potential (Gale & co workers)
Pavese, Catti, Parker, Wall, Phys Chem.
Miner. , 1996
Raiteri, Gale, Quigley, Rodger, J.
Phys. Chem. C, 2010
![Page 38: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/38.jpg)
Structure of hACC
• create with several protocols
melt monohydrocalcite; melt and dehydrate ikaite; anneal random distribution …
• Good agreement with expt structure factor
![Page 39: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/39.jpg)
Water in hACC
• zones of stability for hACC
• carbonate-mediated H-bond networks
0.0 0.2 0.4 0.6 0.8 1.0-10
0
10
20
30
40
50
-10
0
10
20
30
40
50
G
, kJ/m
ol
H Calc.
H Expt.
G estim.
G estim.
MHC & ikaite
G estim.
H
, kJ/m
ol
x2, H2O mole fraction
MHC
Ikaite
ACC-ET
ACC-AQ
pH = 13
pH =12.2
![Page 40: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/40.jpg)
Percolating water clusters
• Universal scaling behaviour
• percolation threshold ~ n = 0.9
• underlying hexagonal lattice
![Page 41: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/41.jpg)
correlated water dynamics
• correlated large hops along chains of water
![Page 42: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/42.jpg)
Classical nucleation
-20
-10
0
10
0 5 10r
delt
a-G
• Activated process:
Favourable “bulk” energy
Unfavourable interfacial energy
Critical “cluster” size
• Classical Nucleation Theory
• Fundamentally Stochastic
Favourable
Unfavourable
3 244
3G r r
![Page 43: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/43.jpg)
hACC nucleation: classical?
DOLLOP, Raiteri et al.,
Nature Comms, 2011
pre-nucleation clusters
Gebauer et al., Science 2008
dense liquid phase
Wallace et al., Science, 2014
![Page 44: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/44.jpg)
Nature of the clusters?
• spontaneous clusters dominate?
• Random Structure searches
random arrangement of ions; minimise; repeat
• Intelligent Water Drop algorithm
nature-inspired global optimsation
erosion of soil to define river valleys
applications to binary LJ mixtures and Janus particles
• (not yet carbonates)
![Page 45: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/45.jpg)
Random Structure Searches
• random arrangement of ions within sphere
• Conjugate Gradient optimisation
• ~10,000 initial structures for each cluster size (1–40 formula units)
• variety of cluster properties
– e.g. energy vs radius of gyration
![Page 46: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/46.jpg)
Energetics of solvated clusters
![Page 47: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/47.jpg)
Free energy Landscape
~2 kT
F / eV
nO–Ow
![Page 48: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/48.jpg)
• “slow” dissolution: convergence from dispersion or initial clusters
• dissolution slows with cluster size (concentration)
• (requires free energy barrier)
Long timescale MD: 20–50 mM
![Page 49: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/49.jpg)
• stable “DOLLOP” behaviour; large clusters
Long timescale MD: ca. 0.5 M
![Page 50: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/50.jpg)
And higher concentrations: 0.5–1.5 M
•
0.57 M
stable clusters
DOLLOP
1.1 M
infinite, liquid like
(phase transition?)
![Page 51: Predictive modelling: a view from the atomic level - Warwick...PVMID PVEP PPA PVMP AP-4M2E AP-2M3Et PVPK PVDMP AP-24DM AP-2M4E AP-23DMc PMA PVDMPc PDHU PBCOD AP-3M2Et PBL AP-4M3Ec](https://reader036.vdocument.in/reader036/viewer/2022071501/61206d357ae0173d87452dcd/html5/thumbnails/51.jpg)
Conclusions
• Source of new mechanistic insight
kinetic inhibition; eggshell proteins
does your predictive model have the right features?
• Can capture some elements of intrinsic complexity
inhibitor polydispersity; nanoparticle structure
need for intelligent sampling
• Progress through data dimensionality reduction
excellent methods for a handful of dimensions
need to identify key dimensions
need to recognise and avoid pitfalls of too few (or wrong) dimensions