Any computational chemists have any advice?

  • Thread starter Thread starter gravenewworld
  • Start date Start date
  • Tags Tags
    Computational
AI Thread Summary
The discussion centers on performing QST2 calculations for transition states using Gaussian software, specifically after obtaining geometry-optimized structures through DFT calculations. The user reports that while the QST2 calculation successfully identifies a transition state with one imaginary frequency, optimizing the structures adjacent to this transition state yields results that do not align with the original input structure in terms of energy. The user seeks alternatives to QST2 for finding transition states and methods to ensure the geometry to the left of the transition state optimizes back to the original energy level. Additionally, the user mentions a journal article discussing the B3LYP/6-31++G** basis set for lithium-containing organic compounds and inquires about its effectiveness and computational demands. Responses indicate that while this basis set may provide improved results, it will significantly increase computational time and resource requirements, especially on older hardware, with estimates suggesting calculations could take much longer than standard DFT methods.
gravenewworld
Messages
1,128
Reaction score
27
I am trying to do a QST2 calculation on two geometry optimized structures that I obtained by doing DFT calcs. The QST2 calc runs fine and finds 1 imaginary frequency resulting in a transition state that appears reasonable. However, when I pick points to the "left" and "right" of the transition state structure and optimize them back down, I find that the 1 optimized strucutre matches exactly to the input strucuture I originally put in, even in total energy. The structure to the "left" of the transition state optimizes back down but results in a structure that is similar to the original input structure but is off by almost 6kcal/mol in energy to my original input structure. Are there any alternatives to the QST2 calc that can be done in order to find transition states between 2 optimized geometries? Are there any tricks I can do in order to get the geometry to the "left" of the transition state to optimize back down to my original input structure (energy wise)? BTW I am currently using gaussian to do the calcs.

-Thanks, GNW.
 
Chemistry news on Phys.org
I also just found a journal article by a guy from Yale who uses the B3LYP/6-31++G** basis set for calculations involving lithium containing organic compounds. Is anyone familiar with using this basis set and what kind of results did you get? Will this take much longer than a DFT calc (or any idea of how much cpu power this will chew up)? The reaction I am trying to optimize is the is a Biphenylene ring, which when introduced to Li2, opens up to give biphenyl. One structure I have tried to optimize is the Biphenylene system with 2 Lithiums exactly over the center the biphenylene, with 1 lithium on top and the other on the bottom.
 
I think that basis set is pretty standard for computational chemistry. It's going to be a lot more intensive than a standard DFT calculation. The results will be better, but it's going to take longer and require more CPU power.
 
Great, my worst suspicions came true. The computer I am using is a 10 year old Mac. The DFT calcs took several days, the 6-31g++** calcs will probably take a week .
 
Yeah, sorry to break the news, but that's probably about right.
 
It seems like a simple enough question: what is the solubility of epsom salt in water at 20°C? A graph or table showing how it varies with temperature would be a bonus. But upon searching the internet I have been unable to determine this with confidence. Wikipedia gives the value of 113g/100ml. But other sources disagree and I can't find a definitive source for the information. I even asked chatgpt but it couldn't be sure either. I thought, naively, that this would be easy to look up without...
I was introduced to the Octet Rule recently and make me wonder, why does 8 valence electrons or a full p orbital always make an element inert? What is so special with a full p orbital? Like take Calcium for an example, its outer orbital is filled but its only the s orbital thats filled so its still reactive not so much as the Alkaline metals but still pretty reactive. Can someone explain it to me? Thanks!!
Back
Top