TURBOMOLE Users Forum
TURBOMOLE Modules => Jobex: Structure Optimization and Molecular Dynamics => Topic started by: Krzmbrzl on February 16, 2021, 06:10:37 PM
-
I am trying to optimize a couple of molecular structures using DFT B3LYP/cc-pVDZ in solution (COSMO; epsilon=7.6).
I start by making a single-point calculation via dscf, followed by the geometry optimization via jobex. And finally want to to perform a frequency analysis using aoforce in order to verify that the found structure is indeed a minimum structure.
The problem is that it seems no matter what I do, I always end up with at least one imaginary frequency in the order of i10 cm^{-1}.
Here's what I have tried to get rid of that:
- Use stricter convergence criteria for the SCF energy and for the convergence criteria of the geometry optimization (stp submenu in define)
- Recalculate the hessian at every step in the optimization (did cause the calculation to not converge)
- Use a bigger grid for the DFT calculation (up to grid size of 7)
- Use screwer to distort the molecule along the imagainy mode
The molecule I have tried all these things on is an imino ester with 24 atoms (H, C, N, O).
Before starting to fiddle with all these parameters I (accidentally) performed a separate optimization using the BP86 functional that TurboMole suggests by default. With that I did not have any imaginary modes.
My question basically is: How can I get rid of these imaginary frequencies? Does this indicate that I should use a different functional than B3LYP?
EDIT: I am using Turbomole 7.4.1
-
From my experience, small numbers of imaginary frequencies may exist due to numerical reasons. From you details I think you performed geometry optimization which gives few negative frequencies. I would suggest,
1. Take that geometry, with tight grids do an scf calculation. If imaginary frequencies occur again, do,
2. change the basis, again if you see those,
3. I will see and give further updates. Please reply if the problem is solved or not.