All this computing power. Can we even simulate a water molecule yet from scratch with QM?
DFT works in many cases, but in some cases it doesn't estimate the energy right, due to how it bypasses some correlation calculations. Bonds are extremely sensitive to energy calculations, so you need to get super close to the actual energy in order to get useful results.
Anyways, someone with more experience here could probably add more, but that's what I've picked up so far.
Right now the lab is having me get comfortable using software like Gaussian and ORCA by simulating a bifurcating reaction. This is a reaction that, depending on the catalyst's momentum, will change what site it bonds to (it makes either a 6-membered or 7-membered ring). I'm finding the intermediate states (where the molecule is most stable) and transition states (the tipping point), and then running trajectories to see which output is more likely.
Once I've finished simulating that, I should be comfortable enough with the process to jump on the bigger project, which is machine learning interatomic potential (MLIP) model distillation. There's a lot of exciting work around speeding up DFT methods by using machine learning (note this is not generative AI, it's merely predicting the molecule energy based on atomic positions). So my one year goal is to get on that project and start contributing.
My five year goal is to, well, graduate. But then I'll probably do a PhD in computational chemistry, since I'm really interested in ways to speed up and scale existing methods. My big dream is to simulate large biological systems while still having bond formation and breaking, to automatically elucidate biochemical pathways, but there's still a lot of steps in-between.
I assume you are familiar with:
https://matt.might.net/articles/phd-school-in-pictures/
I hope and pray that your research helps to make the world a better place and that the rest of us can use your knowledge to help to make the world a place which merits your research.
I haven't seen that website before, but it sounds pretty accurate from what I've heard. It's insane how high of a mountain needs to be climbed just to catch up to the state-of-the-art, and how much work is needed to push through to figure out something truly new.
Here's to making the world a better place!
truly ab initio methods involve figuring out electronic properties from scratch like ionization energy or bandstructure. the real issue is that we dont have exact relations for the exchange and correlation terms. we can know the kinetic energy and charge screening, but we dont know how the electrons are interacting with each other. generally the xc term is treated as a function of electron density or its gradient (see: lda, gga, meta-gga) but there are so many different ways to approximate that. different models are good for different applications also, like transition metals vs organics. and then theres the issue of basis sets (most people use gaussian basis sets that have been tuned over many years but theres also plane waves and finite element methods) which can also change results. and even once u have a decent approximation of density you can try perturbative methods (GW family, delta scf i count also) to try and improve the approximation. i am rambling and typing this on my phone. essentially yes, but often calculations are a little inaccurate. but more accuracy has a higher computational cost, which makes it hard to run larger simulations. tradeoffs of engineering. hope this was coherent.
To simulate a water molecule you do so with a weakly coupled SU(1) gauge theory (light does not interact with itself at tree order) problem where the masses of all constituents are orders of magnitude above the relevant energy scales (you can think of it as the electrons and nuclei and particles coming in and out of existence are contained in a renormalization scheme).
We have "good simulation models" of both, but the former is extraordinarily complicated compared to the latter for the reasons stated above.
To a useful level of accuracy we can certainly simulate water. And we can do the same for a single proton for some definitions of useful (but not other definitions).