Open research towards the discovery of room-temperature superconductors.
Discover ways to transform this asset
POST /speech/from-post
How this post is connected to other assets
Discover other posts like this one
Temperature ramping AIMD simulation of NaCL (mp-22851), taken from 0 K to 300 K over 10ps.
We had this idea before too, but cool to see Claude agrees. A lot of what we're trying to accomplish with this project requires a room temperature material. As comprehensive as Materials Project may b
Sharing some things I'm learning as I work on temperature ramping simulations. The goal of these simulations is to learn how a material's lattice changes with temperature, as thermal expansion, decomp
2025-01-03
Temperature ramping AIMD simulation of H2O (mp-697111), taken from 0 K to 300 K over 10ps.
https://next-gen.materialsproject.org/materials/mp-697111
Note on .cif visual above: bonding algo for our material visualizations are not so good. Will work on something better, or at least more configurable.
MLIP model: M3GNet
Friction: 0.1 / units.fs
Time step: 1.0 * units.fs
Steps: 15000 -> 15ps
Temperature range: 0 K to 300 K
Temperature ramp duration: 10ps
Supercell: 3x3x3
Total atoms: 324
Changed from previous simulation (NaCl):
Increased friction from 0.01 / units.fs to 0.1 / units.fs
Increased trajectory sampling rate, from 1 frame per 1000 steps to 1 frame per 100 steps
A little more action this time! Very cool to see the lattice pick up more energy and really start falling out of order. Next time, I'll just simulate right below freezing to right above freezing so we can watch the phase transition, if that works. By the end of the simulation here at 300 K, any resemblance of a crystal structure is gone and it's just free floating H20.
Also, I increased the trajectory sampling rate here so now each frame is 100 fs. Theoretically we can get 1 frame per femtosecond since that's our simulation step size, but we'll see how useful that much data actually is.
Watching the hydrogens on the oxygen is interesting too. They don't perfectly keep the 104.5 degree bond angle but instead move around quite a bit.
Simulating ice into water
Steadier this time. It could be due to the increased friction which helped coupling between the environment and the material, but it could also be that this supercell was much larger. Around 300 atoms compared to around 50 for the NaCl simulation.
Trajectories file
Molecular dynamics simulation temperature ramping H2O 3x3x3 supercell from 0 K to 300 K