[{"data":1,"prerenderedAt":1213},["ShallowReactive",2],{"mdc-z97cg7-key":3},{"data":4,"body":5},{},{"type":6,"children":7},"root",[8,17,24,38,89,101,124,143,147,153,164,169,196,203,208,214,232,238,268,279,282,288,299,305,323,329,347,359,362,368,373,391,396,414,419,425,430,439,452,470,481,484,490,496,524,530,566,570,583,586,592,597,606,611,629,634,640,645,663,675,678,684,689,702,707,725,736,739,745,751,774,780,793,797,810,815,828,831,837,842,855,860,868,873,891,894,900,905,923,928,946,949,955,960,986,997,1002,1020,1031,1034,1040,1045,1063,1068,1081,1093,1096,1102,1107,1125,1130,1138,1143,1156,1159,1165,1170],{"type":9,"tag":10,"props":11,"children":13},"element","h1",{"id":12},"molecular-dynamics-md-simulation-theoretical-summary",[14],{"type":15,"value":16},"text","🧬 Molecular Dynamics (MD) Simulation — Theoretical Summary",{"type":9,"tag":18,"props":19,"children":21},"h2",{"id":20},"️-1-newtons-equations-and-motion-of-atoms",[22],{"type":15,"value":23},"⚙️ 1. Newton’s Equations and Motion of Atoms",{"type":9,"tag":25,"props":26,"children":27},"p",{},[28,30,36],{"type":15,"value":29},"Molecular dynamics simulations are based on ",{"type":9,"tag":31,"props":32,"children":33},"strong",{},[34],{"type":15,"value":35},"classical Newtonian mechanics",{"type":15,"value":37},".",{"type":9,"tag":39,"props":40,"children":41},"ul",{},[42,79],{"type":9,"tag":43,"props":44,"children":45},"li",{},[46,48],{"type":15,"value":47},"Each atom has:",{"type":9,"tag":39,"props":49,"children":50},{},[51,59,69],{"type":9,"tag":43,"props":52,"children":53},{},[54],{"type":9,"tag":31,"props":55,"children":56},{},[57],{"type":15,"value":58},"Position (x)",{"type":9,"tag":43,"props":60,"children":61},{},[62,67],{"type":9,"tag":31,"props":63,"children":64},{},[65],{"type":15,"value":66},"Velocity (v)",{"type":15,"value":68}," → derivative of position",{"type":9,"tag":43,"props":70,"children":71},{},[72,77],{"type":9,"tag":31,"props":73,"children":74},{},[75],{"type":15,"value":76},"Acceleration (a)",{"type":15,"value":78}," → derivative of velocity",{"type":9,"tag":43,"props":80,"children":81},{},[82,84],{"type":15,"value":83},"Motion is described by ",{"type":9,"tag":31,"props":85,"children":86},{},[87],{"type":15,"value":88},"Force = mass × acceleration",{"type":9,"tag":25,"props":90,"children":91},{},[92,94,99],{"type":15,"value":93},"This means MD simulation becomes a ",{"type":9,"tag":31,"props":95,"children":96},{},[97],{"type":15,"value":98},"huge system of ordinary differential equations",{"type":15,"value":100}," because:",{"type":9,"tag":39,"props":102,"children":103},{},[104,114],{"type":9,"tag":43,"props":105,"children":106},{},[107,109],{"type":15,"value":108},"For ",{"type":9,"tag":31,"props":110,"children":111},{},[112],{"type":15,"value":113},"M atoms → 3M position coordinates + 3M velocity coordinates",{"type":9,"tag":43,"props":115,"children":116},{},[117,119],{"type":15,"value":118},"So the system exists in a ",{"type":9,"tag":31,"props":120,"children":121},{},[122],{"type":15,"value":123},"6M-dimensional phase space",{"type":9,"tag":25,"props":125,"children":126},{},[127,129,134,136,141],{"type":15,"value":128},"👉 Because this system is extremely complex, an ",{"type":9,"tag":31,"props":130,"children":131},{},[132],{"type":15,"value":133},"analytical (exact) solution is impossible",{"type":15,"value":135},".\nInstead, MD uses ",{"type":9,"tag":31,"props":137,"children":138},{},[139],{"type":15,"value":140},"numerical integration",{"type":15,"value":142}," — moving atoms step-by-step in time.",{"type":9,"tag":144,"props":145,"children":146},"hr",{},[],{"type":9,"tag":18,"props":148,"children":150},{"id":149},"️-2-time-step-the-most-critical-parameter",[151],{"type":15,"value":152},"⏱️ 2. Time Step — The Most Critical Parameter",{"type":9,"tag":25,"props":154,"children":155},{},[156,158,163],{"type":15,"value":157},"MD simulations advance in ",{"type":9,"tag":31,"props":159,"children":160},{},[161],{"type":15,"value":162},"small discrete time steps (Δt)",{"type":15,"value":37},{"type":9,"tag":25,"props":165,"children":166},{},[167],{"type":15,"value":168},"Typical values:",{"type":9,"tag":39,"props":170,"children":171},{},[172,184],{"type":9,"tag":43,"props":173,"children":174},{},[175,177,182],{"type":15,"value":176},"🟢 ",{"type":9,"tag":31,"props":178,"children":179},{},[180],{"type":15,"value":181},"~2 femtoseconds",{"type":15,"value":183}," → standard protein simulations",{"type":9,"tag":43,"props":185,"children":186},{},[187,189,194],{"type":15,"value":188},"🟡 ",{"type":9,"tag":31,"props":190,"children":191},{},[192],{"type":15,"value":193},"~4 femtoseconds",{"type":15,"value":195}," → very long or membrane simulations",{"type":9,"tag":197,"props":198,"children":200},"h3",{"id":199},"why-is-time-step-important",[201],{"type":15,"value":202},"Why is time step important?",{"type":9,"tag":25,"props":204,"children":205},{},[206],{"type":15,"value":207},"If the time step is:",{"type":9,"tag":197,"props":209,"children":211},{"id":210},"too-small",[212],{"type":15,"value":213},"❌ Too small",{"type":9,"tag":39,"props":215,"children":216},{},[217,227],{"type":9,"tag":43,"props":218,"children":219},{},[220,222],{"type":15,"value":221},"Simulation becomes ",{"type":9,"tag":31,"props":223,"children":224},{},[225],{"type":15,"value":226},"very slow",{"type":9,"tag":43,"props":228,"children":229},{},[230],{"type":15,"value":231},"You may never reach biologically relevant timescales (µs or ms)",{"type":9,"tag":197,"props":233,"children":235},{"id":234},"too-large",[236],{"type":15,"value":237},"❌ Too large",{"type":9,"tag":39,"props":239,"children":240},{},[241,253,263],{"type":9,"tag":43,"props":242,"children":243},{},[244,246,251],{"type":15,"value":245},"You may ",{"type":9,"tag":31,"props":247,"children":248},{},[249],{"type":15,"value":250},"miss important events",{"type":15,"value":252}," like collisions",{"type":9,"tag":43,"props":254,"children":255},{},[256,258],{"type":15,"value":257},"System can become ",{"type":9,"tag":31,"props":259,"children":260},{},[261],{"type":15,"value":262},"numerically unstable",{"type":9,"tag":43,"props":264,"children":265},{},[266],{"type":15,"value":267},"Atoms may “jump” outside the simulation box → simulation “explodes”",{"type":9,"tag":25,"props":269,"children":270},{},[271,273,278],{"type":15,"value":272},"Thus, MD requires a ",{"type":9,"tag":31,"props":274,"children":275},{},[276],{"type":15,"value":277},"trade-off between accuracy and computational speed",{"type":15,"value":37},{"type":9,"tag":144,"props":280,"children":281},{},[],{"type":9,"tag":18,"props":283,"children":285},{"id":284},"_3-leapfrog-integration-algorithm",[286],{"type":15,"value":287},"🐸 3. Leapfrog Integration Algorithm",{"type":9,"tag":25,"props":289,"children":290},{},[291,293,298],{"type":15,"value":292},"To propagate motion efficiently, MD often uses the ",{"type":9,"tag":31,"props":294,"children":295},{},[296],{"type":15,"value":297},"leapfrog algorithm",{"type":15,"value":37},{"type":9,"tag":197,"props":300,"children":302},{"id":301},"concept",[303],{"type":15,"value":304},"Concept:",{"type":9,"tag":39,"props":306,"children":307},{},[308,318],{"type":9,"tag":43,"props":309,"children":310},{},[311,313],{"type":15,"value":312},"Velocities and positions are calculated at ",{"type":9,"tag":31,"props":314,"children":315},{},[316],{"type":15,"value":317},"offset time points",{"type":9,"tag":43,"props":319,"children":320},{},[321],{"type":15,"value":322},"Values “jump over” each other like a frog hopping",{"type":9,"tag":197,"props":324,"children":326},{"id":325},"advantages",[327],{"type":15,"value":328},"Advantages:",{"type":9,"tag":39,"props":330,"children":331},{},[332,337,342],{"type":9,"tag":43,"props":333,"children":334},{},[335],{"type":15,"value":336},"Time-symmetric",{"type":9,"tag":43,"props":338,"children":339},{},[340],{"type":15,"value":341},"Numerically stable",{"type":9,"tag":43,"props":343,"children":344},{},[345],{"type":15,"value":346},"Accurate for long simulations",{"type":9,"tag":25,"props":348,"children":349},{},[350,352,357],{"type":15,"value":351},"This makes it one of the ",{"type":9,"tag":31,"props":353,"children":354},{},[355],{"type":15,"value":356},"most widely used integration methods",{"type":15,"value":358}," in MD.",{"type":9,"tag":144,"props":360,"children":361},{},[],{"type":9,"tag":18,"props":363,"children":365},{"id":364},"️-4-atoms-are-always-moving-even-at-low-temperature",[366],{"type":15,"value":367},"🌡️ 4. Atoms Are Always Moving (Even at Low Temperature)",{"type":9,"tag":25,"props":369,"children":370},{},[371],{"type":15,"value":372},"Important physical concept:",{"type":9,"tag":39,"props":374,"children":375},{},[376,386],{"type":9,"tag":43,"props":377,"children":378},{},[379,381],{"type":15,"value":380},"Atoms are ",{"type":9,"tag":31,"props":382,"children":383},{},[384],{"type":15,"value":385},"never static",{"type":9,"tag":43,"props":387,"children":388},{},[389],{"type":15,"value":390},"Even near 0 K → quantum and thermal motion still exist",{"type":9,"tag":25,"props":392,"children":393},{},[394],{"type":15,"value":395},"Thus:",{"type":9,"tag":39,"props":397,"children":398},{},[399,409],{"type":9,"tag":43,"props":400,"children":401},{},[402,404],{"type":15,"value":403},"Proteins do ",{"type":9,"tag":31,"props":405,"children":406},{},[407],{"type":15,"value":408},"not sit permanently in minimum energy",{"type":9,"tag":43,"props":410,"children":411},{},[412],{"type":15,"value":413},"They constantly fluctuate around it",{"type":9,"tag":25,"props":415,"children":416},{},[417],{"type":15,"value":418},"Over long simulation time, atomic configurations follow the:",{"type":9,"tag":18,"props":420,"children":422},{"id":421},"boltzmann-distribution",[423],{"type":15,"value":424},"📊 Boltzmann Distribution",{"type":9,"tag":25,"props":426,"children":427},{},[428],{"type":15,"value":429},"This distribution describes:",{"type":9,"tag":431,"props":432,"children":433},"blockquote",{},[434],{"type":9,"tag":25,"props":435,"children":436},{},[437],{"type":15,"value":438},"The probability of a conformation as a function of its potential energy.",{"type":9,"tag":39,"props":440,"children":441},{},[442,447],{"type":9,"tag":43,"props":443,"children":444},{},[445],{"type":15,"value":446},"Low-energy conformations → frequent",{"type":9,"tag":43,"props":448,"children":449},{},[450],{"type":15,"value":451},"High-energy conformations → rare",{"type":9,"tag":25,"props":453,"children":454},{},[455,457,462,464,469],{"type":15,"value":456},"If MD samples ",{"type":9,"tag":31,"props":458,"children":459},{},[460],{"type":15,"value":461},"all conformations",{"type":15,"value":463},", we can calculate ",{"type":9,"tag":31,"props":465,"children":466},{},[467],{"type":15,"value":468},"thermodynamic properties identical to experiments",{"type":15,"value":37},{"type":9,"tag":25,"props":471,"children":472},{},[473,475,480],{"type":15,"value":474},"But reaching full sampling requires ",{"type":9,"tag":31,"props":476,"children":477},{},[478],{"type":15,"value":479},"very long simulations",{"type":15,"value":37},{"type":9,"tag":144,"props":482,"children":483},{},[],{"type":9,"tag":18,"props":485,"children":487},{"id":486},"_5-md-vs-monte-carlo-sampling",[488],{"type":15,"value":489},"🎰 5. MD vs Monte Carlo Sampling",{"type":9,"tag":197,"props":491,"children":493},{"id":492},"monte-carlo-mc",[494],{"type":15,"value":495},"Monte Carlo (MC)",{"type":9,"tag":39,"props":497,"children":498},{},[499,509,514],{"type":9,"tag":43,"props":500,"children":501},{},[502,504],{"type":15,"value":503},"Samples conformations ",{"type":9,"tag":31,"props":505,"children":506},{},[507],{"type":15,"value":508},"randomly",{"type":9,"tag":43,"props":510,"children":511},{},[512],{"type":15,"value":513},"Immediately shows distribution",{"type":9,"tag":43,"props":515,"children":516},{},[517,519],{"type":15,"value":518},"Good for ",{"type":9,"tag":31,"props":520,"children":521},{},[522],{"type":15,"value":523},"thermodynamic properties",{"type":9,"tag":197,"props":525,"children":527},{"id":526},"molecular-dynamics-md",[528],{"type":15,"value":529},"Molecular Dynamics (MD)",{"type":9,"tag":39,"props":531,"children":532},{},[533,543],{"type":9,"tag":43,"props":534,"children":535},{},[536,538],{"type":15,"value":537},"Follows ",{"type":9,"tag":31,"props":539,"children":540},{},[541],{"type":15,"value":542},"realistic time evolution",{"type":9,"tag":43,"props":544,"children":545},{},[546,548],{"type":15,"value":547},"Captures:",{"type":9,"tag":39,"props":549,"children":550},{},[551,556,561],{"type":9,"tag":43,"props":552,"children":553},{},[554],{"type":15,"value":555},"Fast motions",{"type":9,"tag":43,"props":557,"children":558},{},[559],{"type":15,"value":560},"Dynamic pathways",{"type":9,"tag":43,"props":562,"children":563},{},[564],{"type":15,"value":565},"Mechanistic transitions",{"type":9,"tag":25,"props":567,"children":568},{},[569],{"type":15,"value":395},{"type":9,"tag":39,"props":571,"children":572},{},[573,578],{"type":9,"tag":43,"props":574,"children":575},{},[576],{"type":15,"value":577},"MC → better statistics",{"type":9,"tag":43,"props":579,"children":580},{},[581],{"type":15,"value":582},"MD → better physical realism and motion information",{"type":9,"tag":144,"props":584,"children":585},{},[],{"type":9,"tag":18,"props":587,"children":589},{"id":588},"_6-energy-conservation-and-thermostats",[590],{"type":15,"value":591},"🔋 6. Energy Conservation and Thermostats",{"type":9,"tag":25,"props":593,"children":594},{},[595],{"type":15,"value":596},"Total system energy:",{"type":9,"tag":25,"props":598,"children":599},{},[600],{"type":9,"tag":601,"props":602,"children":603},"span",{},[604],{"type":15,"value":605},"\nE_ = E_ + E_\n",{"type":9,"tag":25,"props":607,"children":608},{},[609],{"type":15,"value":610},"However:",{"type":9,"tag":39,"props":612,"children":613},{},[614,624],{"type":9,"tag":43,"props":615,"children":616},{},[617,619],{"type":15,"value":618},"Real atomic collisions are ",{"type":9,"tag":31,"props":620,"children":621},{},[622],{"type":15,"value":623},"not perfectly elastic",{"type":9,"tag":43,"props":625,"children":626},{},[627],{"type":15,"value":628},"Energy is gradually lost",{"type":9,"tag":25,"props":630,"children":631},{},[632],{"type":15,"value":633},"Therefore MD uses:",{"type":9,"tag":197,"props":635,"children":637},{"id":636},"heat-bath-thermostat",[638],{"type":15,"value":639},"🛁 Heat bath / thermostat",{"type":9,"tag":25,"props":641,"children":642},{},[643],{"type":15,"value":644},"Functions:",{"type":9,"tag":39,"props":646,"children":647},{},[648,653,658],{"type":9,"tag":43,"props":649,"children":650},{},[651],{"type":15,"value":652},"Adds energy if system cools",{"type":9,"tag":43,"props":654,"children":655},{},[656],{"type":15,"value":657},"Removes energy if system overheats",{"type":9,"tag":43,"props":659,"children":660},{},[661],{"type":15,"value":662},"Maintains constant temperature",{"type":9,"tag":25,"props":664,"children":665},{},[666,668,673],{"type":15,"value":667},"Energy curves therefore ",{"type":9,"tag":31,"props":669,"children":670},{},[671],{"type":15,"value":672},"fluctuate around a mean value",{"type":15,"value":674},", not perfectly flat.",{"type":9,"tag":144,"props":676,"children":677},{},[],{"type":9,"tag":18,"props":679,"children":681},{"id":680},"_7-equilibration-phase",[682],{"type":15,"value":683},"🚀 7. Equilibration Phase",{"type":9,"tag":25,"props":685,"children":686},{},[687],{"type":15,"value":688},"At simulation start:",{"type":9,"tag":39,"props":690,"children":691},{},[692,697],{"type":9,"tag":43,"props":693,"children":694},{},[695],{"type":15,"value":696},"All atoms may have identical velocities",{"type":9,"tag":43,"props":698,"children":699},{},[700],{"type":15,"value":701},"This is physically unrealistic",{"type":9,"tag":25,"props":703,"children":704},{},[705],{"type":15,"value":706},"During equilibration:",{"type":9,"tag":39,"props":708,"children":709},{},[710,715,720],{"type":9,"tag":43,"props":711,"children":712},{},[713],{"type":15,"value":714},"Velocities redistribute according to mass",{"type":9,"tag":43,"props":716,"children":717},{},[718],{"type":15,"value":719},"System relaxes to correct temperature distribution",{"type":9,"tag":43,"props":721,"children":722},{},[723],{"type":15,"value":724},"Artifacts from starting conditions disappear",{"type":9,"tag":25,"props":726,"children":727},{},[728,730,735],{"type":15,"value":729},"Only after equilibration → ",{"type":9,"tag":31,"props":731,"children":732},{},[733],{"type":15,"value":734},"production simulation begins",{"type":15,"value":37},{"type":9,"tag":144,"props":737,"children":738},{},[],{"type":9,"tag":18,"props":740,"children":742},{"id":741},"_8-solvent-representation",[743],{"type":15,"value":744},"💧 8. Solvent Representation",{"type":9,"tag":197,"props":746,"children":748},{"id":747},"explicit-solvent",[749],{"type":15,"value":750},"Explicit Solvent",{"type":9,"tag":39,"props":752,"children":753},{},[754,759,764,769],{"type":9,"tag":43,"props":755,"children":756},{},[757],{"type":15,"value":758},"Real water molecules simulated",{"type":9,"tag":43,"props":760,"children":761},{},[762],{"type":15,"value":763},"Most accurate",{"type":9,"tag":43,"props":765,"children":766},{},[767],{"type":15,"value":768},"Very computationally expensive",{"type":9,"tag":43,"props":770,"children":771},{},[772],{"type":15,"value":773},"Can introduce artifacts (water clustering, ion shells)",{"type":9,"tag":197,"props":775,"children":777},{"id":776},"implicit-solvent",[778],{"type":15,"value":779},"Implicit Solvent",{"type":9,"tag":39,"props":781,"children":782},{},[783,788],{"type":9,"tag":43,"props":784,"children":785},{},[786],{"type":15,"value":787},"No water molecules",{"type":9,"tag":43,"props":789,"children":790},{},[791],{"type":15,"value":792},"Environment treated mathematically like water",{"type":9,"tag":25,"props":794,"children":795},{},[796],{"type":15,"value":328},{"type":9,"tag":39,"props":798,"children":799},{},[800,805],{"type":9,"tag":43,"props":801,"children":802},{},[803],{"type":15,"value":804},"Faster",{"type":9,"tag":43,"props":806,"children":807},{},[808],{"type":15,"value":809},"Fewer atoms",{"type":9,"tag":25,"props":811,"children":812},{},[813],{"type":15,"value":814},"Disadvantages:",{"type":9,"tag":39,"props":816,"children":817},{},[818,823],{"type":9,"tag":43,"props":819,"children":820},{},[821],{"type":15,"value":822},"Loss of specific hydrogen-bond interactions",{"type":9,"tag":43,"props":824,"children":825},{},[826],{"type":15,"value":827},"Reduced realism near protein surface",{"type":9,"tag":144,"props":829,"children":830},{},[],{"type":9,"tag":18,"props":832,"children":834},{"id":833},"_9-periodic-boundary-conditions-pbc",[835],{"type":15,"value":836},"📦 9. Periodic Boundary Conditions (PBC)",{"type":9,"tag":25,"props":838,"children":839},{},[840],{"type":15,"value":841},"To avoid atoms leaving the system:",{"type":9,"tag":39,"props":843,"children":844},{},[845],{"type":9,"tag":43,"props":846,"children":847},{},[848,850],{"type":15,"value":849},"Simulation box is surrounded by ",{"type":9,"tag":31,"props":851,"children":852},{},[853],{"type":15,"value":854},"copies of itself",{"type":9,"tag":25,"props":856,"children":857},{},[858],{"type":15,"value":859},"When atom exits one side:",{"type":9,"tag":39,"props":861,"children":862},{},[863],{"type":9,"tag":43,"props":864,"children":865},{},[866],{"type":15,"value":867},"It re-enters from opposite side",{"type":9,"tag":25,"props":869,"children":870},{},[871],{"type":15,"value":872},"Benefits:",{"type":9,"tag":39,"props":874,"children":875},{},[876,881,886],{"type":9,"tag":43,"props":877,"children":878},{},[879],{"type":15,"value":880},"Constant particle number",{"type":9,"tag":43,"props":882,"children":883},{},[884],{"type":15,"value":885},"No artificial wall effects",{"type":9,"tag":43,"props":887,"children":888},{},[889],{"type":15,"value":890},"Mimics infinite bulk environment",{"type":9,"tag":144,"props":892,"children":893},{},[],{"type":9,"tag":18,"props":895,"children":897},{"id":896},"_10-computational-limits-of-md",[898],{"type":15,"value":899},"⚡ 10. Computational Limits of MD",{"type":9,"tag":25,"props":901,"children":902},{},[903],{"type":15,"value":904},"Challenges:",{"type":9,"tag":39,"props":906,"children":907},{},[908,913,918],{"type":9,"tag":43,"props":909,"children":910},{},[911],{"type":15,"value":912},"Very small time steps required",{"type":9,"tag":43,"props":914,"children":915},{},[916],{"type":15,"value":917},"Structural changes may occur on ms scale",{"type":9,"tag":43,"props":919,"children":920},{},[921],{"type":15,"value":922},"Huge number of steps needed",{"type":9,"tag":25,"props":924,"children":925},{},[926],{"type":15,"value":927},"Solutions:",{"type":9,"tag":39,"props":929,"children":930},{},[931,936,941],{"type":9,"tag":43,"props":932,"children":933},{},[934],{"type":15,"value":935},"GPU acceleration",{"type":9,"tag":43,"props":937,"children":938},{},[939],{"type":15,"value":940},"Parallel computing",{"type":9,"tag":43,"props":942,"children":943},{},[944],{"type":15,"value":945},"Algorithm improvements",{"type":9,"tag":144,"props":947,"children":948},{},[],{"type":9,"tag":18,"props":950,"children":952},{"id":951},"_11-force-fields-semi-empirical-models",[953],{"type":15,"value":954},"📚 11. Force Fields — Semi-Empirical Models",{"type":9,"tag":25,"props":956,"children":957},{},[958],{"type":15,"value":959},"Force fields are:",{"type":9,"tag":39,"props":961,"children":962},{},[963,968],{"type":9,"tag":43,"props":964,"children":965},{},[966],{"type":15,"value":967},"Parameter tables describing atomic interactions",{"type":9,"tag":43,"props":969,"children":970},{},[971,973],{"type":15,"value":972},"Based on:",{"type":9,"tag":39,"props":974,"children":975},{},[976,981],{"type":9,"tag":43,"props":977,"children":978},{},[979],{"type":15,"value":980},"Physics equations",{"type":9,"tag":43,"props":982,"children":983},{},[984],{"type":15,"value":985},"Experimental data",{"type":9,"tag":25,"props":987,"children":988},{},[989,991,996],{"type":15,"value":990},"Thus MD is ",{"type":9,"tag":31,"props":992,"children":993},{},[994],{"type":15,"value":995},"semi-empirical",{"type":15,"value":37},{"type":9,"tag":25,"props":998,"children":999},{},[1000],{"type":15,"value":1001},"⚠️ Must use scientific judgement:",{"type":9,"tag":39,"props":1003,"children":1004},{},[1005,1015],{"type":9,"tag":43,"props":1006,"children":1007},{},[1008,1010],{"type":15,"value":1009},"Simulations may produce ",{"type":9,"tag":31,"props":1011,"children":1012},{},[1013],{"type":15,"value":1014},"physically unrealistic minima",{"type":9,"tag":43,"props":1016,"children":1017},{},[1018],{"type":15,"value":1019},"Example: alkane chain passing through benzene ring (artifact)",{"type":9,"tag":25,"props":1021,"children":1022},{},[1023,1025],{"type":15,"value":1024},"Always check:\n👉 ",{"type":9,"tag":1026,"props":1027,"children":1028},"em",{},[1029],{"type":15,"value":1030},"Does the simulation result make chemical sense?",{"type":9,"tag":144,"props":1032,"children":1033},{},[],{"type":9,"tag":18,"props":1035,"children":1037},{"id":1036},"_12-what-md-cannot-do",[1038],{"type":15,"value":1039},"❌ 12. What MD Cannot Do",{"type":9,"tag":25,"props":1041,"children":1042},{},[1043],{"type":15,"value":1044},"Classical MD limitations:",{"type":9,"tag":39,"props":1046,"children":1047},{},[1048,1053,1058],{"type":9,"tag":43,"props":1049,"children":1050},{},[1051],{"type":15,"value":1052},"Cannot form or break covalent bonds",{"type":9,"tag":43,"props":1054,"children":1055},{},[1056],{"type":15,"value":1057},"Cannot change protonation states dynamically",{"type":9,"tag":43,"props":1059,"children":1060},{},[1061],{"type":15,"value":1062},"Cannot simulate enzyme reactions directly",{"type":9,"tag":25,"props":1064,"children":1065},{},[1066],{"type":15,"value":1067},"Reason:",{"type":9,"tag":39,"props":1069,"children":1070},{},[1071,1076],{"type":9,"tag":43,"props":1072,"children":1073},{},[1074],{"type":15,"value":1075},"Electronic structure is fixed at simulation start",{"type":9,"tag":43,"props":1077,"children":1078},{},[1079],{"type":15,"value":1080},"Born-Oppenheimer approximation assumed",{"type":9,"tag":25,"props":1082,"children":1083},{},[1084,1086,1091],{"type":15,"value":1085},"To model reactions → need ",{"type":9,"tag":31,"props":1087,"children":1088},{},[1089],{"type":15,"value":1090},"QM/MM methods",{"type":15,"value":1092}," (quantum mechanics + molecular mechanics).",{"type":9,"tag":144,"props":1094,"children":1095},{},[],{"type":9,"tag":18,"props":1097,"children":1099},{"id":1098},"️-13-software-and-acceleration",[1100],{"type":15,"value":1101},"🖥️ 13. Software and Acceleration",{"type":9,"tag":25,"props":1103,"children":1104},{},[1105],{"type":15,"value":1106},"Common MD packages:",{"type":9,"tag":39,"props":1108,"children":1109},{},[1110,1115,1120],{"type":9,"tag":43,"props":1111,"children":1112},{},[1113],{"type":15,"value":1114},"GROMACS",{"type":9,"tag":43,"props":1116,"children":1117},{},[1118],{"type":15,"value":1119},"AMBER",{"type":9,"tag":43,"props":1121,"children":1122},{},[1123],{"type":15,"value":1124},"CHARMM",{"type":9,"tag":25,"props":1126,"children":1127},{},[1128],{"type":15,"value":1129},"Visualization:",{"type":9,"tag":39,"props":1131,"children":1132},{},[1133],{"type":9,"tag":43,"props":1134,"children":1135},{},[1136],{"type":15,"value":1137},"VMD (Visual Molecular Dynamics)",{"type":9,"tag":25,"props":1139,"children":1140},{},[1141],{"type":15,"value":1142},"Speed improvements:",{"type":9,"tag":39,"props":1144,"children":1145},{},[1146,1151],{"type":9,"tag":43,"props":1147,"children":1148},{},[1149],{"type":15,"value":1150},"Parallelization",{"type":9,"tag":43,"props":1152,"children":1153},{},[1154],{"type":15,"value":1155},"GPU (CUDA cores especially important)",{"type":9,"tag":144,"props":1157,"children":1158},{},[],{"type":9,"tag":10,"props":1160,"children":1162},{"id":1161},"key-conceptual-takeaway",[1163],{"type":15,"value":1164},"⭐ Key Conceptual Takeaway",{"type":9,"tag":25,"props":1166,"children":1167},{},[1168],{"type":15,"value":1169},"Molecular dynamics simulations:",{"type":9,"tag":39,"props":1171,"children":1172},{},[1173,1183,1188,1193,1203],{"type":9,"tag":43,"props":1174,"children":1175},{},[1176,1178],{"type":15,"value":1177},"Provide ",{"type":9,"tag":31,"props":1179,"children":1180},{},[1181],{"type":15,"value":1182},"time-resolved molecular motion",{"type":9,"tag":43,"props":1184,"children":1185},{},[1186],{"type":15,"value":1187},"Sample conformational landscapes",{"type":9,"tag":43,"props":1189,"children":1190},{},[1191],{"type":15,"value":1192},"Require careful parameter choices",{"type":9,"tag":43,"props":1194,"children":1195},{},[1196,1198],{"type":15,"value":1197},"Balance ",{"type":9,"tag":31,"props":1199,"children":1200},{},[1201],{"type":15,"value":1202},"accuracy vs computational feasibility",{"type":9,"tag":43,"props":1204,"children":1205},{},[1206,1208],{"type":15,"value":1207},"Must always be interpreted with ",{"type":9,"tag":31,"props":1209,"children":1210},{},[1211],{"type":15,"value":1212},"chemical and physical intuition",1775084312379]