[{"data":1,"prerenderedAt":1445},["ShallowReactive",2],{"mdc-eot3mt-key":3},{"data":4,"body":5},{},{"type":6,"children":7},"root",[8,17,23,27,34,41,54,59,96,106,118,136,139,145,150,163,168,200,205,208,214,219,225,230,248,260,263,269,274,297,309,312,318,323,336,341,353,356,362,367,390,393,399,405,410,419,430,435,445,450,458,463,481,491,494,500,510,533,538,556,559,565,570,578,583,596,601,609,617,628,638,641,647,652,660,668,673,683,688,706,709,715,721,726,744,749,757,760,766,771,779,784,789,809,812,818,829,834,841,859,865,883,886,892,897,902,910,915,933,943,946,952,958,963,971,975,983,988,996,999,1005,1010,1023,1028,1031,1037,1042,1050,1055,1068,1073,1076,1082,1087,1138,1143,1156,1159,1165,1170,1186,1191,1196,1201,1209,1212,1218,1223,1241,1249,1254,1257,1263,1291,1296,1299,1305,1310,1328,1334,1339,1347,1352,1355,1361,1366,1384,1387,1393,1398,1416,1419,1425,1430,1435,1440],{"type":9,"tag":10,"props":11,"children":13},"element","h1",{"id":12},"molecular-dynamics-md-simulations-complete-educational-summary",[14],{"type":15,"value":16},"text","🧬 Molecular Dynamics (MD) Simulations — Complete Educational Summary",{"type":9,"tag":18,"props":19,"children":20},"p",{},[21],{"type":15,"value":22},"📄 Source:",{"type":9,"tag":24,"props":25,"children":26},"hr",{},[],{"type":9,"tag":28,"props":29,"children":31},"h2",{"id":30},"what-is-molecular-dynamics-md",[32],{"type":15,"value":33},"⭐ What is Molecular Dynamics (MD)?",{"type":9,"tag":35,"props":36,"children":38},"h3",{"id":37},"core-idea",[39],{"type":15,"value":40},"💡 Core idea",{"type":9,"tag":18,"props":42,"children":43},{},[44,46,52],{"type":15,"value":45},"Molecular Dynamics is a ",{"type":9,"tag":47,"props":48,"children":49},"strong",{},[50],{"type":15,"value":51},"computer simulation method",{"type":15,"value":53}," used to mimic how atoms move in real life.",{"type":9,"tag":18,"props":55,"children":56},{},[57],{"type":15,"value":58},"The main concept:",{"type":9,"tag":60,"props":61,"children":62},"ul",{},[63,74,84],{"type":9,"tag":64,"props":65,"children":66},"li",{},[67,69],{"type":15,"value":68},"Atoms interact via a ",{"type":9,"tag":47,"props":70,"children":71},{},[72],{"type":15,"value":73},"potential energy function",{"type":9,"tag":64,"props":75,"children":76},{},[77,79],{"type":15,"value":78},"From this energy → we can calculate ",{"type":9,"tag":47,"props":80,"children":81},{},[82],{"type":15,"value":83},"forces",{"type":9,"tag":64,"props":85,"children":86},{},[87,89,94],{"type":15,"value":88},"Using ",{"type":9,"tag":47,"props":90,"children":91},{},[92],{"type":15,"value":93},"Newton’s laws",{"type":15,"value":95},", we compute how atoms move over time",{"type":9,"tag":18,"props":97,"children":98},{},[99,101],{"type":15,"value":100},"👉 In simple terms:\n",{"type":9,"tag":47,"props":102,"children":103},{},[104],{"type":15,"value":105},"MD = solving physics equations to create a “movie” of atomic motion.",{"type":9,"tag":18,"props":107,"children":108},{},[109,111,116],{"type":15,"value":110},"This makes MD a kind of ",{"type":9,"tag":47,"props":112,"children":113},{},[114],{"type":15,"value":115},"computational microscope",{"type":15,"value":117},", because:",{"type":9,"tag":60,"props":119,"children":120},{},[121,126],{"type":9,"tag":64,"props":122,"children":123},{},[124],{"type":15,"value":125},"We cannot experimentally observe the motion of every atom in a protein over time",{"type":9,"tag":64,"props":127,"children":128},{},[129,131],{"type":15,"value":130},"MD allows us to visualize this at ",{"type":9,"tag":47,"props":132,"children":133},{},[134],{"type":15,"value":135},"femtosecond resolution",{"type":9,"tag":24,"props":137,"children":138},{},[],{"type":9,"tag":28,"props":140,"children":142},{"id":141},"️-basic-md-algorithm",[143],{"type":15,"value":144},"⏱️ Basic MD Algorithm",{"type":9,"tag":18,"props":146,"children":147},{},[148],{"type":15,"value":149},"Time is divided into very small steps:",{"type":9,"tag":60,"props":151,"children":152},{},[153],{"type":9,"tag":64,"props":154,"children":155},{},[156,158],{"type":15,"value":157},"Typical time step:\n",{"type":9,"tag":47,"props":159,"children":160},{},[161],{"type":15,"value":162},"1–2 femtoseconds (10⁻¹⁵ s)",{"type":9,"tag":18,"props":164,"children":165},{},[166],{"type":15,"value":167},"At each step:",{"type":9,"tag":169,"props":170,"children":171},"ol",{},[172,182],{"type":9,"tag":64,"props":173,"children":174},{},[175,177],{"type":15,"value":176},"Calculate forces on every atom using a ",{"type":9,"tag":47,"props":178,"children":179},{},[180],{"type":15,"value":181},"force field",{"type":9,"tag":64,"props":183,"children":184},{},[185,187],{"type":15,"value":186},"Update:",{"type":9,"tag":60,"props":188,"children":189},{},[190,195],{"type":9,"tag":64,"props":191,"children":192},{},[193],{"type":15,"value":194},"Position",{"type":9,"tag":64,"props":196,"children":197},{},[198],{"type":15,"value":199},"Velocity",{"type":9,"tag":18,"props":201,"children":202},{},[203],{"type":15,"value":204},"Then repeat millions to trillions of times.",{"type":9,"tag":24,"props":206,"children":207},{},[],{"type":9,"tag":28,"props":209,"children":211},{"id":210},"applications-of-md",[212],{"type":15,"value":213},"🔬 Applications of MD",{"type":9,"tag":18,"props":215,"children":216},{},[217],{"type":15,"value":218},"MD is extremely powerful and widely used.",{"type":9,"tag":35,"props":220,"children":222},{"id":221},"drug-binding-and-allostery",[223],{"type":15,"value":224},"💊 Drug binding and allostery",{"type":9,"tag":18,"props":226,"children":227},{},[228],{"type":15,"value":229},"MD can help determine:",{"type":9,"tag":60,"props":231,"children":232},{},[233,238,243],{"type":9,"tag":64,"props":234,"children":235},{},[236],{"type":15,"value":237},"Where a drug binds",{"type":9,"tag":64,"props":239,"children":240},{},[241],{"type":15,"value":242},"How binding changes protein structure",{"type":9,"tag":64,"props":244,"children":245},{},[246],{"type":15,"value":247},"How binding at one site affects another site (allosteric effects)",{"type":9,"tag":18,"props":249,"children":250},{},[251,253,258],{"type":15,"value":252},"This helps design ",{"type":9,"tag":47,"props":254,"children":255},{},[256],{"type":15,"value":257},"better drugs",{"type":15,"value":259},".",{"type":9,"tag":24,"props":261,"children":262},{},[],{"type":9,"tag":35,"props":264,"children":266},{"id":265},"️-understanding-protein-mechanisms",[267],{"type":15,"value":268},"⚙️ Understanding protein mechanisms",{"type":9,"tag":18,"props":270,"children":271},{},[272],{"type":15,"value":273},"MD can simulate:",{"type":9,"tag":60,"props":275,"children":276},{},[277,287,292],{"type":9,"tag":64,"props":278,"children":279},{},[280,282],{"type":15,"value":281},"Transition between ",{"type":9,"tag":47,"props":283,"children":284},{},[285],{"type":15,"value":286},"active ↔ inactive states",{"type":9,"tag":64,"props":288,"children":289},{},[290],{"type":15,"value":291},"Signal propagation through receptors",{"type":9,"tag":64,"props":293,"children":294},{},[295],{"type":15,"value":296},"Structural coupling between domains",{"type":9,"tag":18,"props":298,"children":299},{},[300,302,307],{"type":15,"value":301},"This helps understand ",{"type":9,"tag":47,"props":303,"children":304},{},[305],{"type":15,"value":306},"how proteins work dynamically",{"type":15,"value":308},", not just statically.",{"type":9,"tag":24,"props":310,"children":311},{},[],{"type":9,"tag":35,"props":313,"children":315},{"id":314},"protein-folding",[316],{"type":15,"value":317},"🧩 Protein folding",{"type":9,"tag":18,"props":319,"children":320},{},[321],{"type":15,"value":322},"MD can study:",{"type":9,"tag":60,"props":324,"children":325},{},[326,331],{"type":9,"tag":64,"props":327,"children":328},{},[329],{"type":15,"value":330},"How unfolded proteins find their native folded structure",{"type":9,"tag":64,"props":332,"children":333},{},[334],{"type":15,"value":335},"Folding pathways and intermediates",{"type":9,"tag":18,"props":337,"children":338},{},[339],{"type":15,"value":340},"However:",{"type":9,"tag":18,"props":342,"children":343},{},[344,346,351],{"type":15,"value":345},"⚠️ MD is ",{"type":9,"tag":47,"props":347,"children":348},{},[349],{"type":15,"value":350},"not usually the best method to predict final folded structure",{"type":15,"value":352}," from scratch.",{"type":9,"tag":24,"props":354,"children":355},{},[],{"type":9,"tag":35,"props":357,"children":359},{"id":358},"md-beyond-proteins",[360],{"type":15,"value":361},"🌍 MD beyond proteins",{"type":9,"tag":18,"props":363,"children":364},{},[365],{"type":15,"value":366},"Also used for:",{"type":9,"tag":60,"props":368,"children":369},{},[370,375,380,385],{"type":9,"tag":64,"props":371,"children":372},{},[373],{"type":15,"value":374},"DNA / RNA",{"type":9,"tag":64,"props":376,"children":377},{},[378],{"type":15,"value":379},"Lipid membranes",{"type":9,"tag":64,"props":381,"children":382},{},[383],{"type":15,"value":384},"Carbohydrates",{"type":9,"tag":64,"props":386,"children":387},{},[388],{"type":15,"value":389},"Materials science systems",{"type":9,"tag":24,"props":391,"children":392},{},[],{"type":9,"tag":28,"props":394,"children":396},{"id":395},"how-forces-are-calculated",[397],{"type":15,"value":398},"⚡ How Forces are Calculated",{"type":9,"tag":35,"props":400,"children":402},{"id":401},"potential-energy-function",[403],{"type":15,"value":404},"Potential energy function",{"type":9,"tag":18,"props":406,"children":407},{},[408],{"type":15,"value":409},"The system has total potential energy:",{"type":9,"tag":18,"props":411,"children":412},{},[413],{"type":9,"tag":414,"props":415,"children":416},"span",{},[417],{"type":15,"value":418},"\nU(x)\n",{"type":9,"tag":18,"props":420,"children":421},{},[422,424,429],{"type":15,"value":423},"which depends on positions of ",{"type":9,"tag":47,"props":425,"children":426},{},[427],{"type":15,"value":428},"all atoms",{"type":15,"value":259},{"type":9,"tag":18,"props":431,"children":432},{},[433],{"type":15,"value":434},"This function is called a:",{"type":9,"tag":18,"props":436,"children":437},{},[438,440],{"type":15,"value":439},"👉 ",{"type":9,"tag":47,"props":441,"children":442},{},[443],{"type":15,"value":444},"Force field",{"type":9,"tag":18,"props":446,"children":447},{},[448],{"type":15,"value":449},"From energy we compute force:",{"type":9,"tag":18,"props":451,"children":452},{},[453],{"type":9,"tag":414,"props":454,"children":455},{},[456],{"type":15,"value":457},"\nF(x) = -\nabla U(x)\n",{"type":9,"tag":18,"props":459,"children":460},{},[461],{"type":15,"value":462},"Meaning:",{"type":9,"tag":60,"props":464,"children":465},{},[466,476],{"type":9,"tag":64,"props":467,"children":468},{},[469,471],{"type":15,"value":470},"Force points toward ",{"type":9,"tag":47,"props":472,"children":473},{},[474],{"type":15,"value":475},"lower energy",{"type":9,"tag":64,"props":477,"children":478},{},[479],{"type":15,"value":480},"Larger energy change → stronger force",{"type":9,"tag":18,"props":482,"children":483},{},[484,486],{"type":15,"value":485},"At ",{"type":9,"tag":47,"props":487,"children":488},{},[489],{"type":15,"value":490},"energy minima → force = 0",{"type":9,"tag":24,"props":492,"children":493},{},[],{"type":9,"tag":35,"props":495,"children":497},{"id":496},"force-vector",[498],{"type":15,"value":499},"Force Vector",{"type":9,"tag":18,"props":501,"children":502},{},[503,505],{"type":15,"value":504},"For a system with ",{"type":9,"tag":47,"props":506,"children":507},{},[508],{"type":15,"value":509},"N atoms",{"type":9,"tag":60,"props":511,"children":512},{},[513,523],{"type":9,"tag":64,"props":514,"children":515},{},[516,518],{"type":15,"value":517},"There are ",{"type":9,"tag":47,"props":519,"children":520},{},[521],{"type":15,"value":522},"3N coordinates",{"type":9,"tag":64,"props":524,"children":525},{},[526,528],{"type":15,"value":527},"Therefore also ",{"type":9,"tag":47,"props":529,"children":530},{},[531],{"type":15,"value":532},"3N force components",{"type":9,"tag":18,"props":534,"children":535},{},[536],{"type":15,"value":537},"Example:",{"type":9,"tag":60,"props":539,"children":540},{},[541,546,551],{"type":9,"tag":64,"props":542,"children":543},{},[544],{"type":15,"value":545},"Atom 1 → force in x, y, z",{"type":9,"tag":64,"props":547,"children":548},{},[549],{"type":15,"value":550},"Atom 2 → force in x, y, z",{"type":9,"tag":64,"props":552,"children":553},{},[554],{"type":15,"value":555},"etc.",{"type":9,"tag":24,"props":557,"children":558},{},[],{"type":9,"tag":28,"props":560,"children":562},{"id":561},"equations-of-motion",[563],{"type":15,"value":564},"🧮 Equations of Motion",{"type":9,"tag":18,"props":566,"children":567},{},[568],{"type":15,"value":569},"Newton’s Second Law:",{"type":9,"tag":18,"props":571,"children":572},{},[573],{"type":9,"tag":414,"props":574,"children":575},{},[576],{"type":15,"value":577},"\nF = ma\n",{"type":9,"tag":18,"props":579,"children":580},{},[581],{"type":15,"value":582},"Also:",{"type":9,"tag":60,"props":584,"children":585},{},[586,591],{"type":9,"tag":64,"props":587,"children":588},{},[589],{"type":15,"value":590},"Velocity = derivative of position",{"type":9,"tag":64,"props":592,"children":593},{},[594],{"type":15,"value":595},"Acceleration = derivative of velocity",{"type":9,"tag":18,"props":597,"children":598},{},[599],{"type":15,"value":600},"Thus:",{"type":9,"tag":18,"props":602,"children":603},{},[604],{"type":9,"tag":414,"props":605,"children":606},{},[607],{"type":15,"value":608},"\n\frac{dx}{dt} = v\n",{"type":9,"tag":18,"props":610,"children":611},{},[612],{"type":9,"tag":414,"props":613,"children":614},{},[615],{"type":15,"value":616},"\n\frac{dv}{dt} = \frac{F(x)}{m}\n",{"type":9,"tag":18,"props":618,"children":619},{},[620,622,627],{"type":15,"value":621},"These form a large system of ",{"type":9,"tag":47,"props":623,"children":624},{},[625],{"type":15,"value":626},"ordinary differential equations",{"type":15,"value":259},{"type":9,"tag":18,"props":629,"children":630},{},[631,633],{"type":15,"value":632},"👉 Analytical solutions are impossible → must solve ",{"type":9,"tag":47,"props":634,"children":635},{},[636],{"type":15,"value":637},"numerically.",{"type":9,"tag":24,"props":639,"children":640},{},[],{"type":9,"tag":35,"props":642,"children":644},{"id":643},"numerical-integration",[645],{"type":15,"value":646},"🔁 Numerical Integration",{"type":9,"tag":18,"props":648,"children":649},{},[650],{"type":15,"value":651},"Simplest update:",{"type":9,"tag":18,"props":653,"children":654},{},[655],{"type":9,"tag":414,"props":656,"children":657},{},[658],{"type":15,"value":659},"\nx_{i+1} = x_i + delta t v_i\n",{"type":9,"tag":18,"props":661,"children":662},{},[663],{"type":9,"tag":414,"props":664,"children":665},{},[666],{"type":15,"value":667},"\nv_{i+1} = v_i + delta t \frac{F(x_i)}{m}\n",{"type":9,"tag":18,"props":669,"children":670},{},[671],{"type":15,"value":672},"Better method used in practice:",{"type":9,"tag":18,"props":674,"children":675},{},[676,678],{"type":15,"value":677},"⭐ ",{"type":9,"tag":47,"props":679,"children":680},{},[681],{"type":15,"value":682},"Leapfrog Verlet integration",{"type":9,"tag":18,"props":684,"children":685},{},[686],{"type":15,"value":687},"Why?",{"type":9,"tag":60,"props":689,"children":690},{},[691,696,701],{"type":9,"tag":64,"props":692,"children":693},{},[694],{"type":15,"value":695},"Time-symmetric",{"type":9,"tag":64,"props":697,"children":698},{},[699],{"type":15,"value":700},"More stable for long simulations",{"type":9,"tag":64,"props":702,"children":703},{},[704],{"type":15,"value":705},"Conserves energy better",{"type":9,"tag":24,"props":707,"children":708},{},[],{"type":9,"tag":28,"props":710,"children":712},{"id":711},"key-properties-of-md",[713],{"type":15,"value":714},"🎯 Key Properties of MD",{"type":9,"tag":35,"props":716,"children":718},{"id":717},"atoms-never-stop-moving",[719],{"type":15,"value":720},"🔄 Atoms never stop moving",{"type":9,"tag":18,"props":722,"children":723},{},[724],{"type":15,"value":725},"Even at equilibrium:",{"type":9,"tag":60,"props":727,"children":728},{},[729,734],{"type":9,"tag":64,"props":730,"children":731},{},[732],{"type":15,"value":733},"Atoms vibrate continuously",{"type":9,"tag":64,"props":735,"children":736},{},[737,739],{"type":15,"value":738},"Simulation samples the ",{"type":9,"tag":47,"props":740,"children":741},{},[742],{"type":15,"value":743},"Boltzmann distribution",{"type":9,"tag":18,"props":745,"children":746},{},[747],{"type":15,"value":748},"Probability of a structure:",{"type":9,"tag":18,"props":750,"children":751},{},[752],{"type":9,"tag":414,"props":753,"children":754},{},[755],{"type":15,"value":756},"\np(x) propto e^{-U(x)/k_BT}\n",{"type":9,"tag":24,"props":758,"children":759},{},[],{"type":9,"tag":35,"props":761,"children":763},{"id":762},"energy-conservation",[764],{"type":15,"value":765},"🔥 Energy conservation",{"type":9,"tag":18,"props":767,"children":768},{},[769],{"type":15,"value":770},"Total energy:",{"type":9,"tag":60,"props":772,"children":773},{},[774],{"type":9,"tag":64,"props":775,"children":776},{},[777],{"type":15,"value":778},"Potential + kinetic",{"type":9,"tag":18,"props":780,"children":781},{},[782],{"type":15,"value":783},"Should stay constant.",{"type":9,"tag":18,"props":785,"children":786},{},[787],{"type":15,"value":788},"But:",{"type":9,"tag":60,"props":790,"children":791},{},[792,797],{"type":9,"tag":64,"props":793,"children":794},{},[795],{"type":15,"value":796},"Numerical rounding errors cause slow drift",{"type":9,"tag":64,"props":798,"children":799},{},[800,802,807],{"type":15,"value":801},"Therefore simulations often use ",{"type":9,"tag":47,"props":803,"children":804},{},[805],{"type":15,"value":806},"thermostats",{"type":15,"value":808}," to control temperature.",{"type":9,"tag":24,"props":810,"children":811},{},[],{"type":9,"tag":35,"props":813,"children":815},{"id":814},"importance-of-solvent",[816],{"type":15,"value":817},"💧 Importance of Solvent",{"type":9,"tag":18,"props":819,"children":820},{},[821,823,828],{"type":15,"value":822},"Ignoring water causes ",{"type":9,"tag":47,"props":824,"children":825},{},[826],{"type":15,"value":827},"major artifacts",{"type":15,"value":259},{"type":9,"tag":18,"props":830,"children":831},{},[832],{"type":15,"value":833},"Two approaches:",{"type":9,"tag":835,"props":836,"children":838},"h4",{"id":837},"explicit-solvent",[839],{"type":15,"value":840},"🧊 Explicit solvent",{"type":9,"tag":60,"props":842,"children":843},{},[844,849,854],{"type":9,"tag":64,"props":845,"children":846},{},[847],{"type":15,"value":848},"Simulate water molecules directly",{"type":9,"tag":64,"props":850,"children":851},{},[852],{"type":15,"value":853},"More accurate",{"type":9,"tag":64,"props":855,"children":856},{},[857],{"type":15,"value":858},"More computationally expensive",{"type":9,"tag":835,"props":860,"children":862},{"id":861},"️-implicit-solvent",[863],{"type":15,"value":864},"🌫️ Implicit solvent",{"type":9,"tag":60,"props":866,"children":867},{},[868,873,878],{"type":9,"tag":64,"props":869,"children":870},{},[871],{"type":15,"value":872},"Mathematical approximation",{"type":9,"tag":64,"props":874,"children":875},{},[876],{"type":15,"value":877},"Faster",{"type":9,"tag":64,"props":879,"children":880},{},[881],{"type":15,"value":882},"Less realistic",{"type":9,"tag":24,"props":884,"children":885},{},[],{"type":9,"tag":35,"props":887,"children":889},{"id":888},"️-periodic-boundary-conditions",[890],{"type":15,"value":891},"♾️ Periodic Boundary Conditions",{"type":9,"tag":18,"props":893,"children":894},{},[895],{"type":15,"value":896},"Real systems contain ~10²³ molecules.",{"type":9,"tag":18,"props":898,"children":899},{},[900],{"type":15,"value":901},"Simulations contain only:",{"type":9,"tag":60,"props":903,"children":904},{},[905],{"type":9,"tag":64,"props":906,"children":907},{},[908],{"type":15,"value":909},"100–1,000,000 molecules",{"type":9,"tag":18,"props":911,"children":912},{},[913],{"type":15,"value":914},"To mimic bulk environment:",{"type":9,"tag":60,"props":916,"children":917},{},[918,928],{"type":9,"tag":64,"props":919,"children":920},{},[921,923],{"type":15,"value":922},"Simulation box is ",{"type":9,"tag":47,"props":924,"children":925},{},[926],{"type":15,"value":927},"replicated infinitely",{"type":9,"tag":64,"props":929,"children":930},{},[931],{"type":15,"value":932},"Particle leaving one side enters from opposite side",{"type":9,"tag":18,"props":934,"children":935},{},[936,938],{"type":15,"value":937},"Like ",{"type":9,"tag":47,"props":939,"children":940},{},[941],{"type":15,"value":942},"Pac-Man world 🙂",{"type":9,"tag":24,"props":944,"children":945},{},[],{"type":9,"tag":28,"props":947,"children":949},{"id":948},"limitations-of-md",[950],{"type":15,"value":951},"🚧 Limitations of MD",{"type":9,"tag":35,"props":953,"children":955},{"id":954},"timescale-problem",[956],{"type":15,"value":957},"⏳ Timescale problem",{"type":9,"tag":18,"props":959,"children":960},{},[961],{"type":15,"value":962},"Time step ≈ 2 fs\nProtein events:",{"type":9,"tag":60,"props":964,"children":965},{},[966],{"type":9,"tag":64,"props":967,"children":968},{},[969],{"type":15,"value":970},"ns → μs → ms → seconds",{"type":9,"tag":18,"props":972,"children":973},{},[974],{"type":15,"value":600},{"type":9,"tag":60,"props":976,"children":977},{},[978],{"type":9,"tag":64,"props":979,"children":980},{},[981],{"type":15,"value":982},"Millions → trillions of steps needed",{"type":9,"tag":18,"props":984,"children":985},{},[986],{"type":15,"value":987},"Even today:",{"type":9,"tag":60,"props":989,"children":990},{},[991],{"type":9,"tag":64,"props":992,"children":993},{},[994],{"type":15,"value":995},"Microsecond simulations are computationally demanding.",{"type":9,"tag":24,"props":997,"children":998},{},[],{"type":9,"tag":35,"props":1000,"children":1002},{"id":1001},"force-field-approximations",[1003],{"type":15,"value":1004},"🎯 Force field approximations",{"type":9,"tag":18,"props":1006,"children":1007},{},[1008],{"type":15,"value":1009},"Force fields:",{"type":9,"tag":60,"props":1011,"children":1012},{},[1013,1018],{"type":9,"tag":64,"props":1014,"children":1015},{},[1016],{"type":15,"value":1017},"Are not perfect",{"type":9,"tag":64,"props":1019,"children":1020},{},[1021],{"type":15,"value":1022},"May give incorrect dynamics or stability",{"type":9,"tag":18,"props":1024,"children":1025},{},[1026],{"type":15,"value":1027},"Experience is required to interpret MD results.",{"type":9,"tag":24,"props":1029,"children":1030},{},[],{"type":9,"tag":35,"props":1032,"children":1034},{"id":1033},"no-bond-breaking-in-classical-md",[1035],{"type":15,"value":1036},"🔗 No bond breaking in classical MD",{"type":9,"tag":18,"props":1038,"children":1039},{},[1040],{"type":15,"value":1041},"Standard MD:",{"type":9,"tag":60,"props":1043,"children":1044},{},[1045],{"type":9,"tag":64,"props":1046,"children":1047},{},[1048],{"type":15,"value":1049},"Bonds remain fixed",{"type":9,"tag":18,"props":1051,"children":1052},{},[1053],{"type":15,"value":1054},"But real systems can have:",{"type":9,"tag":60,"props":1056,"children":1057},{},[1058,1063],{"type":9,"tag":64,"props":1059,"children":1060},{},[1061],{"type":15,"value":1062},"Disulfide exchange",{"type":9,"tag":64,"props":1064,"children":1065},{},[1066],{"type":15,"value":1067},"Protonation changes",{"type":9,"tag":18,"props":1069,"children":1070},{},[1071],{"type":15,"value":1072},"Advanced methods exist (QM/MM etc.).",{"type":9,"tag":24,"props":1074,"children":1075},{},[],{"type":9,"tag":28,"props":1077,"children":1079},{"id":1078},"md-software",[1080],{"type":15,"value":1081},"💻 MD Software",{"type":9,"tag":18,"props":1083,"children":1084},{},[1085],{"type":15,"value":1086},"Popular packages:",{"type":9,"tag":60,"props":1088,"children":1089},{},[1090,1098,1106,1114,1122,1130],{"type":9,"tag":64,"props":1091,"children":1092},{},[1093],{"type":9,"tag":47,"props":1094,"children":1095},{},[1096],{"type":15,"value":1097},"GROMACS",{"type":9,"tag":64,"props":1099,"children":1100},{},[1101],{"type":9,"tag":47,"props":1102,"children":1103},{},[1104],{"type":15,"value":1105},"AMBER",{"type":9,"tag":64,"props":1107,"children":1108},{},[1109],{"type":9,"tag":47,"props":1110,"children":1111},{},[1112],{"type":15,"value":1113},"NAMD",{"type":9,"tag":64,"props":1115,"children":1116},{},[1117],{"type":9,"tag":47,"props":1118,"children":1119},{},[1120],{"type":15,"value":1121},"Desmond",{"type":9,"tag":64,"props":1123,"children":1124},{},[1125],{"type":9,"tag":47,"props":1126,"children":1127},{},[1128],{"type":15,"value":1129},"OpenMM",{"type":9,"tag":64,"props":1131,"children":1132},{},[1133],{"type":9,"tag":47,"props":1134,"children":1135},{},[1136],{"type":15,"value":1137},"CHARMM",{"type":9,"tag":18,"props":1139,"children":1140},{},[1141],{"type":15,"value":1142},"Visualization tools:",{"type":9,"tag":60,"props":1144,"children":1145},{},[1146,1151],{"type":9,"tag":64,"props":1147,"children":1148},{},[1149],{"type":15,"value":1150},"⭐ VMD (very common)",{"type":9,"tag":64,"props":1152,"children":1153},{},[1154],{"type":15,"value":1155},"PyMOL (alternative)",{"type":9,"tag":24,"props":1157,"children":1158},{},[],{"type":9,"tag":28,"props":1160,"children":1162},{"id":1161},"force-fields",[1163],{"type":15,"value":1164},"🧱 Force Fields",{"type":9,"tag":18,"props":1166,"children":1167},{},[1168],{"type":15,"value":1169},"Main families:",{"type":9,"tag":60,"props":1171,"children":1172},{},[1173,1177,1181],{"type":9,"tag":64,"props":1174,"children":1175},{},[1176],{"type":15,"value":1137},{"type":9,"tag":64,"props":1178,"children":1179},{},[1180],{"type":15,"value":1105},{"type":9,"tag":64,"props":1182,"children":1183},{},[1184],{"type":15,"value":1185},"OPLS-AA",{"type":9,"tag":18,"props":1187,"children":1188},{},[1189],{"type":15,"value":1190},"Important:",{"type":9,"tag":18,"props":1192,"children":1193},{},[1194],{"type":15,"value":1195},"⚠️ Force field name ≠ software name\n(e.g. CHARMM force field vs CHARMM program)",{"type":9,"tag":18,"props":1197,"children":1198},{},[1199],{"type":15,"value":1200},"New research:",{"type":9,"tag":60,"props":1202,"children":1203},{},[1204],{"type":9,"tag":64,"props":1205,"children":1206},{},[1207],{"type":15,"value":1208},"Neural-network force fields trained on quantum data.",{"type":9,"tag":24,"props":1210,"children":1211},{},[],{"type":9,"tag":28,"props":1213,"children":1215},{"id":1214},"why-md-is-computationally-expensive",[1216],{"type":15,"value":1217},"🚀 Why MD is Computationally Expensive",{"type":9,"tag":18,"props":1219,"children":1220},{},[1221],{"type":15,"value":1222},"Reasons:",{"type":9,"tag":169,"props":1224,"children":1225},{},[1226,1231,1236],{"type":9,"tag":64,"props":1227,"children":1228},{},[1229],{"type":15,"value":1230},"Huge number of time steps",{"type":9,"tag":64,"props":1232,"children":1233},{},[1234],{"type":15,"value":1235},"Many force calculations each step",{"type":9,"tag":64,"props":1237,"children":1238},{},[1239],{"type":15,"value":1240},"Non-bonded interactions scale as:",{"type":9,"tag":18,"props":1242,"children":1243},{},[1244],{"type":9,"tag":414,"props":1245,"children":1246},{},[1247],{"type":15,"value":1248},"\nN^2\n",{"type":9,"tag":18,"props":1250,"children":1251},{},[1252],{"type":15,"value":1253},"This dominates runtime.",{"type":9,"tag":24,"props":1255,"children":1256},{},[],{"type":9,"tag":35,"props":1258,"children":1260},{"id":1259},"speed-up-strategies",[1261],{"type":15,"value":1262},"Speed-up strategies",{"type":9,"tag":60,"props":1264,"children":1265},{},[1266,1271,1276,1281,1286],{"type":9,"tag":64,"props":1267,"children":1268},{},[1269],{"type":15,"value":1270},"Cutoff distance for van der Waals forces",{"type":9,"tag":64,"props":1272,"children":1273},{},[1274],{"type":15,"value":1275},"Parallel computing",{"type":9,"tag":64,"props":1277,"children":1278},{},[1279],{"type":15,"value":1280},"GPU acceleration",{"type":9,"tag":64,"props":1282,"children":1283},{},[1284],{"type":15,"value":1285},"Constraint algorithms (freeze fast vibrations)",{"type":9,"tag":64,"props":1287,"children":1288},{},[1289],{"type":15,"value":1290},"Enhanced sampling methods",{"type":9,"tag":18,"props":1292,"children":1293},{},[1294],{"type":15,"value":1295},"Specialized hardware chips for MD are also being developed.",{"type":9,"tag":24,"props":1297,"children":1298},{},[],{"type":9,"tag":28,"props":1300,"children":1302},{"id":1301},"monte-carlo-simulation-alternative-method",[1303],{"type":15,"value":1304},"🎲 Monte Carlo Simulation — Alternative Method",{"type":9,"tag":18,"props":1306,"children":1307},{},[1308],{"type":15,"value":1309},"Instead of Newton’s laws:",{"type":9,"tag":60,"props":1311,"children":1312},{},[1313,1323],{"type":9,"tag":64,"props":1314,"children":1315},{},[1316,1318],{"type":15,"value":1317},"Make ",{"type":9,"tag":47,"props":1319,"children":1320},{},[1321],{"type":15,"value":1322},"random structural changes",{"type":9,"tag":64,"props":1324,"children":1325},{},[1326],{"type":15,"value":1327},"Accept or reject based on energy.",{"type":9,"tag":35,"props":1329,"children":1331},{"id":1330},"metropolis-criterion",[1332],{"type":15,"value":1333},"Metropolis criterion",{"type":9,"tag":18,"props":1335,"children":1336},{},[1337],{"type":15,"value":1338},"If energy decreases → accept\nIf increases → accept with probability:",{"type":9,"tag":18,"props":1340,"children":1341},{},[1342],{"type":9,"tag":414,"props":1343,"children":1344},{},[1345],{"type":15,"value":1346},"\ne^{-Delta U / k_BT}\n",{"type":9,"tag":18,"props":1348,"children":1349},{},[1350],{"type":15,"value":1351},"After long time → samples Boltzmann distribution.",{"type":9,"tag":24,"props":1353,"children":1354},{},[],{"type":9,"tag":35,"props":1356,"children":1358},{"id":1357},"️-simulated-annealing",[1359],{"type":15,"value":1360},"❄️ Simulated annealing",{"type":9,"tag":18,"props":1362,"children":1363},{},[1364],{"type":15,"value":1365},"Gradually lowering temperature:",{"type":9,"tag":60,"props":1367,"children":1368},{},[1369,1379],{"type":9,"tag":64,"props":1370,"children":1371},{},[1372,1374],{"type":15,"value":1373},"Helps find ",{"type":9,"tag":47,"props":1375,"children":1376},{},[1377],{"type":15,"value":1378},"low-energy structures",{"type":9,"tag":64,"props":1380,"children":1381},{},[1382],{"type":15,"value":1383},"Used for structure optimization.",{"type":9,"tag":24,"props":1385,"children":1386},{},[],{"type":9,"tag":28,"props":1388,"children":1390},{"id":1389},"machine-learning-direction",[1391],{"type":15,"value":1392},"🤖 Machine Learning Direction",{"type":9,"tag":18,"props":1394,"children":1395},{},[1396],{"type":15,"value":1397},"New methods:",{"type":9,"tag":60,"props":1399,"children":1400},{},[1401,1406,1411],{"type":9,"tag":64,"props":1402,"children":1403},{},[1404],{"type":15,"value":1405},"Learn to sample equilibrium states directly",{"type":9,"tag":64,"props":1407,"children":1408},{},[1409],{"type":15,"value":1410},"Can give huge speedups",{"type":9,"tag":64,"props":1412,"children":1413},{},[1414],{"type":15,"value":1415},"But do not provide real dynamical trajectories (no “movie”).",{"type":9,"tag":24,"props":1417,"children":1418},{},[],{"type":9,"tag":10,"props":1420,"children":1422},{"id":1421},"final-big-picture",[1423],{"type":15,"value":1424},"⭐ Final Big Picture",{"type":9,"tag":18,"props":1426,"children":1427},{},[1428],{"type":15,"value":1429},"MD simulations allow us to:",{"type":9,"tag":18,"props":1431,"children":1432},{},[1433],{"type":15,"value":1434},"✅ Observe atomic motion\n✅ Study protein mechanisms\n✅ Explore folding pathways\n✅ Understand drug interactions\n✅ Complement experiments",{"type":9,"tag":18,"props":1436,"children":1437},{},[1438],{"type":15,"value":1439},"But challenges remain:",{"type":9,"tag":18,"props":1441,"children":1442},{},[1443],{"type":15,"value":1444},"⚠️ Timescales\n⚠️ Force-field accuracy\n⚠️ Computational cost",1775084300891]