upvote
Depends what level of accuracy you want. I just started in a computational chemistry lab so I'll probably get some details wrong, but for small systems, you can use a method called CCSD(T) for up to ~20 atoms, but it scales O(N^7). I've been mainly using DFT for the systems I've been simulating, which scales O(N^3). I've been running a system with about 50 atoms with a decent basis set (how the orbitals are modelled), and it takes about 30 minutes for each optimization step with 24 cores and 48 GB of RAM.

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.

reply
Cool details, thanks. To help me understand your life, what would be like a one year and a five year research goal for you? I never spent time in lab sciences so it’s kind of a black box for me.
reply
Disclaimer: I'm only a freshman, so there's still a ton I don't know :)

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.

reply
Good luck!

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.

reply
Thank you for the kind words! I've been wanting to do this research precisely because of firsthand experience with how hard chronic illness can be, and I'm hoping to attack it with a systems approach.

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!

reply
do you think quantum computers would help simulating this? I've seen contradictory opinions from the experts - it can in theory but not really in practice (even assuming sufficiently large quantum computers will be built)
reply
I'm only a freshman, so I don't feel very qualified to comment on that :) I hope so though!
reply
yea im doing my masters in dft research so ik abt this. depends what u want 2 simulate! chemists more do molecular dynamics type stuff and will use experimental data for fitting data etc. like uh what surface of a metal water will react with from thermodynamics or something. (that isnt my field lol i just know a lot of catalysis guys.)

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.

reply
To complete accuracy, we cannot yet manage one proton.
reply
If you want to get pedantic we can't simulate anything with complete accuracy in the absence of a theory that encompasses all the known forces. Which we don't have. (Damn you gravity. Can't you just get along with the others)

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).

reply
That's a fundamentally different problem and a terribly unfair comparison.
reply
Am I right with my assumption that by "fundamentally different problem", you mean we lack a good simulation model, but that the number of degrees of freedom would actually be manageable?
reply
To simulate a proton you need to solve a strongly coupled highly relativistic SU(3) gauge theory (naturally non-abelian i.e. the force carrier field itself carries charge and is self-interacting at tree order) problem with constituents that have masses orders of magnitude below the relevant energy scales (i.e. you have many matter AND force particles that can pop in and out of existence and they all strongly interact with one another).

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.

reply