Questions tagged [geometry-optimization]

Questions about finding the molecular or material geometry that has the lowest energy.

Filter by
Sorted by
Tagged with
2 votes
1 answer
17 views

Large energy differences across successive geometry optimization steps in SIESTA

I am trying to calculate adsorption energies of acetophenone on a 6x6x3 nickel surface with about 12 Angstrom of vacuum above it using SIESTA 5. I have arranged the molecule so that the ring plane is ...
user1420303's user avatar
2 votes
0 answers
42 views

Wrong lattice parameter of CeO2 after structure relaxation with VASP

I am calculating the lattice parameter of CeO2 with VASP, but the optimized lattice parameter is wrong according to many references (Physical Review B 75, 035109 (2007); The Journal of Chemical ...
FaDA's user avatar
  • 121
1 vote
0 answers
20 views

LBFGS -- Incomplete Convergence due to BadIntpln

I am doing energy minimization of a molecular system in TINKER. While doing so, I am getting the error in the title of this question. What might be my problem?
Himanshu's user avatar
4 votes
1 answer
155 views

Computing proton-proton potentials for infinite periodic materials

For electronic structure calculation, when optimizing the wavefunction for an infinite periodic material, the proton-proton energy is needed to monitor the convergence of the SCF loop. How is the ...
mle's user avatar
  • 733
3 votes
0 answers
27 views

Constrained geometry optmization: Molecule interacts with its surroundings

I have a water cluster(suppose) generated from MD run. Now, I want to optimize the geometry of only one water molecule of interest using QC methods (e.g. B3LYP/6-31+G(d,p)). I intend that the water ...
Pro's user avatar
  • 810
3 votes
0 answers
25 views

Discrepancy between Force Convergence in Geometrical Optimization and SCF Calculation in VASP: Reasons and Troubleshooting Suggestions

During a geometrical optimization of ions and lattice parameters in VASP, the forces acting on atoms were minimized to less than 0.01 eV/Å. However, when performing the subsequent self-consistent ...
Jaafar Mehrez's user avatar
1 vote
0 answers
41 views

Structural phase transition under uniaxial strain

I am trying calculate a strain-induced phase transition in some weakly correlated crystalline material using Quantum Espresso (QE). The compound has two (experimentally observed) phases, one of which ...
manju9's user avatar
  • 31
2 votes
1 answer
52 views

Optimal choice between ISIF = 3 and ISIF = 4 for lattice parameter and atomic position relaxation in VASP calculations

In my VASP calculations, I am interested in simultaneously optimizing both the lattice parameters and atomic positions of my system. I understand that using ISIF = 3 allows for relaxation of both ...
Jaafar Mehrez's user avatar
3 votes
1 answer
113 views

How to pre-optimize metal-organic complexes in RDKit using UFF?

I have run into a particular problem while trying to automate my input file generation for ORCA calculations. Apparently, I either do not fully understand how SMILES strings work, or I am demanding ...
J.Doe's user avatar
  • 443
3 votes
0 answers
56 views

Quantum Espresso optimization calculation not converging

I have been trying to optimize the geometry of a copper based metal halide perovskite. I opened the CIF file using BURAI GUI and generated the VC-RELAX input file. I have increased the steps to be 200 ...
Kamal Nassar's user avatar
8 votes
0 answers
89 views

How does the line search in geometry optimisations work?

I am trying to write a Newton-Raphson geometry optimisation script in python (Current version of my code) using the energy, gradient and Hessian from Gaussian 16. In order to reach a minimum, I would ...
evd's user avatar
  • 81
7 votes
1 answer
82 views

How many non-redundant internal coordinates are needed for optimisation of linear molecule?

I have been working to implement optimisation in internal coordinates (for a transition state search software I am working on). However, the implementation details are almost always glossed over in ...
S R Maiti's user avatar
  • 6,591
3 votes
1 answer
93 views

UseSym Keyword in ORCA fails to find appropriate number of SALCs

I am trying to run some geometry optimizations in ORCA 5.0.4 and I am using the UseSym keyword to ensure that my molecules retain Cs Symmetry during optimization. ...
J.Doe's user avatar
  • 443
5 votes
1 answer
157 views

ORCA crashes because of RI approximation?

I am doing a geometry optimization for a cationic imidazolate-zinc complex. I'm trying to run an optimization using this input file. It uses the RI approximation. ORCA seems to be running fine until ...
J.Doe's user avatar
  • 443
5 votes
1 answer
157 views

Is it okay to use semi-empirical methods to optimize structures if I perform high-level energy calculations?

For non-covalent systems comprised of small-molecule monomers, folks are usually interested in understanding the interaction and binding energies of small (gas-phase) clusters ... things like $\ce{X(...
epalos's user avatar
  • 4,719
2 votes
2 answers
95 views

How to select charge and multiplicity for macro-molecules?

In optimising proteins with DFT/B3LYP, how can we detect the charge and multiplicity of the protein?
PriZarah's user avatar
  • 301
4 votes
2 answers
102 views

Cannot optimize hexa aqua Cu(II) and Fe(III)

I try to optimize very simple systems - hexa aqua Cu(II) and hexa aqua Fe(III) - in Gaussian 16. The inputs are provided in the end of the post. Unfortunately, the geometry optimization never finishes ...
farmaceut's user avatar
  • 536
6 votes
1 answer
119 views

Use of genetic algorithm based methods for geometry optimization

What are the benefits or pitfalls in using a GA based optimizer for geometry optimization instead of a more traditional Hessian based algorithm? It can be argued that performing single point energy ...
Hemanth Haridas's user avatar
6 votes
2 answers
199 views

Number of k-points for unit and super cell

I am just new with solid state calculation and would like to know one thing. Why is plane wave basis set is good for 1s and 2s orbitals, but for higher convergence is poor? Why is it better to use ...
Alexander Teryahskin's user avatar
5 votes
1 answer
201 views

Convergence problem in geometry optimization in Gaussian16

I have a quite big set of molecules to optimize and study and, after some research in the literature, I selected B3LYP/6-311G* level of theory. These molecules are quite big and wobbly, so I expected ...
Laura's user avatar
  • 967
10 votes
1 answer
505 views

How can I specify different theory levels for different atoms in Gaussian?

I have a bunch of geometry optimisations I want to run for a cation Metal-organic complex. This is going to be a repeat of calculations that were done in "Turbomole" previously and my Prof. ...
J.Doe's user avatar
  • 443
4 votes
0 answers
42 views

How to compute only part of molecular Hessian accurately (hybrid Hessian or partial Hessian calculation)?

Hessian is required for minimisation and transition state (TS) optimisation algorithms. While the full Hessian at a certain level of theory is calculated usually for TS optimisation, usually ...
S R Maiti's user avatar
  • 6,591
3 votes
1 answer
364 views

Why the relaxation of my material in Quantum ESPRESSO wasn't achieved?

I did a vc-relax calculation using HSE by Quantum ESPRESSO and after 22h time WALL the job was complete and JOB DONE was mentionned in the end of the output file ...
Camilla's user avatar
  • 1,957
5 votes
1 answer
128 views

DFT calculation for conjugated polymers

I am new to DFT calculation. I want to know the structure of conjugated polymers (plane angle, etc.) and their electronic and optical properties (such as HOMO-LUMO, IR, etc.). I am using 6-311+G(d) ...
Duy-Minh Hoang's user avatar
5 votes
1 answer
195 views

Geometry optimization: nuclear coordinates vs lattice parameters

I would like to know whether, at least in principle, minimization of energy with respect to nuclear coordinates leads to the same result as minimization of cell energy with respect to lattice ...
chemdamned's user avatar
4 votes
1 answer
155 views

Is the Broyden's algorithm sufficient to obtain the crystal structure of a hypothetical alloy, or should something else be used?

The Broyden (Broyden—Fletcher—Goldfarb—Shanno or BFGS) algorithm can be known in the matter modeling to optimize a crystal lattice geometry — to clarify positions of nodes of a crystal lattice in ...
SFriendly's user avatar
  • 877
4 votes
0 answers
272 views

VASP geometry optimization not converging [closed]

I have tried everything I could find online but I am not able to relax this structure in VASP. I started out with the basic INCAR ...
Chan's user avatar
  • 629
7 votes
3 answers
187 views

What would be the preferred format for a geometry database?

Introduction (analogy to basis set database files) Many programs have a database of basis sets stored in a file such as the following (so that users can just specify in their input file something like ...
Nike Dattani's user avatar
  • 33.6k
6 votes
1 answer
433 views

How can I predict the time of a simulation using Quantum ESPRESSO?

I performed a vc-relax calculation with hybrid functional on a 2 atoms material using the command ...
Camilla's user avatar
  • 1,957
3 votes
1 answer
144 views

Why is line search not popular in quasi-Newton geometry optimisers?

One of the most important quasi-Newton optimisation methods is the BFGS (which uses an iterated Hessian update scheme to avoid recalculation of the actual Hessian). Most numerical optimisation ...
S R Maiti's user avatar
  • 6,591
5 votes
1 answer
198 views

Constrained optimisation on a hypersphere

I am currently trying to implement a GS2 (Gonzalez-Schlegel second order) IRC algorithm in a python code. I am following the original paper ref(1). The main problem is in the constrained optimisation ...
S R Maiti's user avatar
  • 6,591
3 votes
1 answer
145 views

How to calculate the lattice parameter of a supercell using Phonopy? (For Quantum ESPRESSO)

I found in the internet a method to generate a supercell from the input file using Phonopy code. And I have a question if someone can clarify my doubts please. The method is as following: We make an ...
Camilla's user avatar
  • 1,957
4 votes
1 answer
416 views

What's the appropriate method to create a supercell for Quantum ESPRESSO?

I created a CIF file using the experimental parameters of my material, then I used it to create my input file for the unit cell using Cif2Cell. After that, I ...
Camilla's user avatar
  • 1,957
5 votes
1 answer
82 views

Hessian used for calculating GDIIS error vectors

GDIIS (Geometry optimisation with Direct Inversion of Iterative Subspace) is a well known method used by many QM codes for geometry optimisation. It is inspired from the almost universal DIIS method ...
S R Maiti's user avatar
  • 6,591
3 votes
1 answer
80 views

How do you project the Hessian along a certain vector?

I have been trying to understand the Gaussian optimiser (Berny optimiser algorithm) from what has been written in their manual. So far, I have understood most of it, but one part has been giving me ...
S R Maiti's user avatar
  • 6,591
4 votes
3 answers
615 views

What is the bond length of O2--?

I need the bond length in Å for $\ce{O_2^{--}}$ from this geometry database page (type "O2--" for the input). I get a bunch of possible values. I don't know much about those yet but after ...
user avatar
2 votes
0 answers
38 views

How to handle explicit solvent molecules when calculating energies of reactant molecules in isolation [closed]

When determining the energy barrier, it is recommended to compute the energies of reactant molecules in isolation to account for the entropy change of forming the reactant complex (hence factoring ...
learner254's user avatar
3 votes
1 answer
778 views

ORCA: Getting the error "mpirun noticed that process rank 0 with PID 383161 on node exited on signal 11 (Segmentation fault)."

I am getting a segmentation fault error when running the following input file below. The calculation runs fine in gas phase and CPCM modes, but fails with SMD. What could be the issue? Some are saying ...
learner254's user avatar
4 votes
1 answer
91 views

Is it ok to start from an identity matrix Hessian guess for optimisation?

I am attempting to implement an algorithm for transition state finding (in Python, but that's not important). The algorithm is called BITSS (Binary Image-pair Transition State Search). It defines a ...
S R Maiti's user avatar
  • 6,591
6 votes
1 answer
289 views

Do I have to use the same functional for dielectric function calculations, as the one used for relaxing the structure?

This is a general question regarding the practices of doing optical calculations and the role of the ionic contribution to the dielectric function (with heavy reference to VASP). Based on VASP's ...
CW Tan's user avatar
  • 1,034
2 votes
0 answers
68 views

Perovskite Optimization does not converge using gaussian [closed]

I am trying to optimize the geometry of (C5H14N2)[CuCl4] using g09. However, the calculation is not converging at any level of theory in 250 cycles. I am new to this field and I think I have made a ...
Kamal Sobhy's user avatar
4 votes
1 answer
102 views

Crash of DF-CCSD or DF-CCSD(T) geometry optimization (MRCC package)

I am interested in geometry optimization at the CCSD(T) or DF-CCSD(T) level, with MRCC. The calculation, however, crashes when executing prop: ...
EvGeniy's user avatar
  • 191
6 votes
0 answers
168 views

Optimizing geometry of materials under pressure using VASP [closed]

I am currently trying to reproduce the results of [1], wherein they have pressurized $\ce{Cs_2SnI_6}$ from $0$ to $20$ GPa. I ran calculations with PBEsol and used the PSTRESS tag in order to optimize ...
Chan's user avatar
  • 629
4 votes
1 answer
28 views

Predicting parameters by taking a weighted average of multiple possible structures with energy rather than the structure with least energy

I am currently looking to predict the structure of a material who we know exists under the conditions I am intrusted in(pressure), but the cif file does not exist. Now I have to first make the cif ...
Chan's user avatar
  • 629
4 votes
1 answer
184 views

How does ORCA treat intermediates in Nudged Elastic Band (NEB) calculations

I'm trying to run a quite complicated NEB optimization with (hopefully) one intermediate structure. But before losing time, I was wondering: how does ORCA treat that intermediate? does ORCA optimize ...
Andrea Pellegrini's user avatar
1 vote
0 answers
99 views

maxstep Gaussian keywords

I'm running a Berny optimization and I need to include the keyword maxstep=1. Observing the output, in the convergence summary table I can see that the ...
Andrea Pellegrini's user avatar
5 votes
1 answer
78 views

Change in space group upon application of stress during geometry relaxation in VASP

I am trying to duplicate the results of this study using VASP. Where they have observed that the structure of $\ce{Cs_2SnI_6}$ switches from Fm-3m to I2/m upon application of pressure more than 3.3 ...
Chan's user avatar
  • 629
3 votes
1 answer
111 views

Failure to converge geometry of ringed structure

Hi I'm trying to optimise a molecule by using the b3lyp/6-31G(d) method , the calculus took more than 2 days and the energy never converged . Here is my input code ...
diamond999's user avatar
3 votes
0 answers
63 views

Restart geometry optimization in Gaussian after hangup error [closed]

I am trying to optimize the geometry of a paraquat molecule interacting with a gold bulk. The optimization is stopping after a few steps. The error shown in the .log...
Anastasia Bitter's user avatar
5 votes
1 answer
563 views

How to restart an interrupted TS optimization job in Gaussian16 without the CalcFC option from the chk file?

I have to do the TS optimization of a very large structure. Initially, I used the typical opt(ts,calcfc,noeigentest) at the route section, but during the ...
Duffoon's user avatar
  • 51