Questions tagged [computational-chemistry]

For questions about computational chemistry.

Filter by
Sorted by
Tagged with
5
votes
1answer
8 views

What does “Test job not archived” mean?

I am running frequency computation with geom=connectivity in Gaussian. It all converges now but at the very end of the last convergence it says: ...
10
votes
1answer
102 views

How to generate accurate displacements along normal modes?

How can i displace my structure along a particular normal mode, without coupling to other modes ? Rather than a exact mathematical procedure, I am more interested in a software/script/package that can ...
11
votes
2answers
59 views

What is the most appropriate (and easy-to-use) all-atoms force field for liquid electrolytes?

I need to simulate a salt diluted in a liquid electrolyte. My experience is with LAMMPS, but so far I only had coarse grained toy models, and setting up the force field was straightforward. Now I have ...
6
votes
1answer
37 views

Vibrational frequency analysis for NEB result through Maestro

I have ran a Nudged Elastic Band (NEB) calculation using QuantumEspresso through Maestro to find the potential transition state. Within Maestro there is an option where you can select for NEB type ...
10
votes
1answer
101 views

Computing optimized 3D structures in python that take solvent into account

I have been using RDKit to generate 3D coordinates, cleaning-up the structures using a general-purpose force field as follows: ...
10
votes
1answer
140 views

What novel drugs are developed via CADD?

Virtual screening , protein ligand docking, MD simulation, FEP free energy calculations and the QM/MM approach have all been widely applied in optimization for drug discovery. Which drugs that have ...
12
votes
1answer
58 views

Excited state energy gradient during geometry optimization of a target excited state

During the geometry optimization of the ground state of a molecule, we try to optimize the energy gradient with respect to the nuclear positions such that the system remains in its lowest energy state....
6
votes
0answers
54 views

Two-stage unbinding & binding energies with DFT

I'm trying to calculate the stability of the product of a coordinate bonding process of molecules A and B in the polarizable continuum solvent model (PCM, solvent set to water). I'm getting energies ...
5
votes
0answers
50 views

Understanding emission spectra using TDDFT calculation

I did a TDDFT calculation to get the emission spectra of my molecule. Here are the first five roots for the emission calculation: ...
7
votes
1answer
59 views

Formation energy and Energy above convex hull as indicators of stability

I was trying to understand some indicators about the stability of materials (that basically tell you whether it can be synthesized or not.. right?). Specifically, I'm trying to link the knowledge of ...
7
votes
1answer
33 views

The number of long-range interaction pairs in Tinker

I tested MM calculation of a benzene molecule with Tinker. The following is the xyz file used for the input. ...
13
votes
1answer
130 views

Are nuclear ab-initio methods related to materials ab-initio methods?

I just read this synopsis in Physics where scientists used first-principles nuclear physics calculations to predict the stability of 700 isotopes up to iron. I didn't know that this was possible, ...
5
votes
0answers
67 views

What is the relation between enthalpy, ionisation potential, internal energy and HOMO/LUMO?

Is there any relation between enthalpy, ionisation potential, internal energy and HOMO/LUMO? I am dehydrogenating a molecule and I need to study the enthalpy of dehydrogenation. If I get any relation ...
8
votes
1answer
52 views

Error running psi4 on charged molecules

I am running a simple energy calculation on a conformer of 3-Sulfopropyl methacrylate, which is given my C=C(C)C(=O)OCCCS(=O)(=O)[O-]. This is my workflow: ...
11
votes
1answer
179 views

How to get more significant digits from OpenBabel?

I am generating conformers from SMILES string using OpenBabel. The shell script looks something like: ...
8
votes
1answer
94 views

Using OpenBabel or RDKIT Python API to generate conformers from a SMILES string

I have a bunch of monomers in the form of SMILES string and I want to generate conformers from them. At the moment, I am using RDKit to convert these SMILES strings to .sdf files, then using ...
8
votes
2answers
251 views

How do I figure out how many monomer configurations I need for my simulation?

I am running a quantum mechanical simulation on Psi4 for a certain number of monomers to generate data for a database. Some of these monomers are short, while some are long. At the moment, I am ...
10
votes
1answer
134 views

Calculating/Visualizing electronic density from orbital information written in the Molden format

I'm trying to obtain an electron density from orbital information written in the Molden format by the formula, $$ \rho(\mathbf{r}) = \sum_{i}^{occ} |\phi_{i}(\mathbf{r})|^{2}. \tag{1} $$ Is there any ...
9
votes
3answers
368 views

What is the procedure for constructing an ab initio potential energy surface for CH3Cl + Ar?

I would like to construct a model analytical potential function for collision between CH3Cl and Ar. This would be the sum of CH3Cl intramolecular + CH3Cl---Ar intermolecular potential. First I would ...
8
votes
2answers
175 views

What are Pair Natural Orbitals?

Recently, I came across electron correlation methods that use local pair natural orbitals (LPNOs) or domain-based local pair natural orbitals (DLPNO), such as DLPNO-MP2, DLPNO-CCSD etc. The dynamic ...
9
votes
1answer
82 views

Energy cut-off value for conformer calculations

I want to get familiar with the conformer calculation. I found a paper where they did the conformer calculations already for the molecule I work with. The paper can found here. In the experimental ...
8
votes
1answer
44 views

Simulation of electronic potential from nanoparticles interacting with molecules

I work with SERS, where the Raman signal gets enhanced by the electric field of the nanoparticles. Now I like to simulate this kind of spectrum with Gaussian. I found a paper[1], where it looks very ...
11
votes
1answer
222 views

What are some of the biggest theoretical and computational challenges in Condensed Phase Chemistry?

What are some of the biggest open-questions / challenges in Condensed Phase theory and computation? What efforts are being done to answer them? I realize Condensed Phase chemistry is a huge field, but ...
6
votes
1answer
214 views

How to choose starting magnetization while doing spin polarised calculation in Quantum ESPRESSO?

I am trying to run spin polarised DFT calculation in Quantum ESPRESSO on system Co2MnSi. I kept starting magnetization for all these three species as, ...
10
votes
1answer
112 views

What are the trade-offs setting DIRSCF to either .TRUE. or .FALSE. (the default) in GAMESS (US)?

The first software I used to run calculations was Gaussian09. When submitting a input, I used to see a pattern of resource usage where the allocated CPU quickly maxed out, remaining high most of the ...
8
votes
1answer
125 views

Which basis set should be used for a system that contains both anions and cations?

I am trying to model a system that has both anions and cations interacting with each other. (In this case it is calcium ion interacting with two small carboxylic acid molecules). Now, looking through ...
11
votes
0answers
70 views

Bioisosteric replacement using SMARTS (KNIME and RDKit) [closed]

I am trying to create a KNIME workflow that would accept a list of compounds and carry out bioisosteric replacements (we will use the following example here: carboxylic acid to tetrazole) ...
15
votes
2answers
93 views

What is the difference between orbital localisation and NBO?

I am new to computational chemistry, I have used mainly GAMESS and Orca. In GAMESS, there are orbital localisation schemes like Pipek-Mezey, Ruedenberg-Edminston, etc. which generate orbitals that ...
13
votes
2answers
184 views

How to get potential energy surface from GAMESS?

I'm trying to learn GAMESS and performing a simple PES scan of oxygen molecule. After the completion of the calculations I get a .log file with lots of sections like that: ...
7
votes
1answer
487 views

What is the physical meaning of occupations='smearing' in Quantum ESPRESSO?

I want to understand, what is the intuition behind using the occupations='smearing' or 'tetrahedra' in Quantum ESPRESSO input file? I mean how we can interpret what is going to happen inside of ...
9
votes
1answer
297 views

How to do spin polarization calculations using Quantum ESPRESSO?

I am currently new to DFT calculations in Quantum ESPRESSO. Though, I have read some research papers for determining spin polarization and Half-metallicity in certain Heusler alloys. These research ...
11
votes
1answer
122 views

Predicting magnetic moment of a metal complex computationally

Recently, for entrance exams, I have been dealing with a lot of weird compounds. One of them is $\ce{K3[Mn(CN)6]}$. The objective was to predict the magnetic moment of the above complex. I predict ...
9
votes
2answers
313 views

How to choose pseudopotential for DFT calculations in Quantum ESPRESSO?

I am a beginner and newly started running DFT calculations to find out the electronic band structure of certain materials in Quantum ESPRESSO. But we have to select the pseudopotential for running the ...
10
votes
0answers
93 views

I have access to 1000 cores across 10 nodes, what would be the best way to prallelize an MP2, GW, or RPA calculation for a periodic system in VASP? [closed]

What's the best way for parallelization in VASP for periodic systems? In the documentation, it mostly explains it for the local computer and molecule calculation. Apparently, these explanations don't ...
9
votes
1answer
381 views

High precision helium energy

In the paper "High Precision Theory of Atomic Helium", Drake lists the then best-known ground state energy of the Schrödinger equation for helium to 22 digits as: $$\lambda_0 \approx -2....
9
votes
1answer
42 views

How can I determine the maximum value of a spin orbital and its derivative in any direction?

I'd like to be able to recover the maximum value that any of the basis functions can take, as well as their maximum derivative. How can I do this? For context, I need it because I want to calculate ...
10
votes
2answers
85 views

Analytic Hessians for meta-GGA functionals

In many of the free/open-source QM programs like GAMESS, ORCA, NWChem, the calculation of analytic hessians is not possible when a meta-GGA density functional is used, like M06-2X. My first question ...
8
votes
1answer
114 views

Are there other custom chemistry programming languages besides SVL?

I must confess some ignorance... I just became aware of the Scientific Vector Language (SVL) which has been around since 1994 and was created solely for computational chemistry. It has a large ...
10
votes
1answer
90 views

Correcting thermodynamic data for concentration, what concentration is best to use?

In computational chemistry when using transition state theory we often use standard states for all involved structures. However, sometimes people argue that it should be corrected for concentration (...
21
votes
2answers
799 views

Why do people care about reversibility in molecular dynamics simulations?

When carrying out a molecular dynamics simulation, one has to choose an ensemble. Depending on the problem, one usually works in the microcanonical (NVE), canonical (NVT), or grand canonical ensemble (...
8
votes
0answers
68 views

Equilibration of normal modes in water molecule [closed]

When running a molecular dynamics simulation on a singular water molecule in vacuum to analyze its normal modes of vibration from the NVE (microcanonical) ensemble, the bending mode frequency is ...
8
votes
1answer
288 views

How to perform virtual crystal approximation calculations in VASP?

Can someone share appropriate tags and methods to perform virtual crystal approximation calculations in VASP?
8
votes
1answer
752 views

How to fix Error Link9999 in Gaussian

I have a problem with geometry optimization in Gaussian software with M062X method. The structure was smoothly optimized using M062X/Gen. Nevertheless, I need very accurate results so afterwards the ...
8
votes
1answer
146 views

Semi-canonicalisation vs canonicalisation of the Fock matrix and orbitals

I have seen the terms semi-canonicalized and canonicalized used in relation to the Fock matrix, density matrices, and orbitals; however, I am unsure what these terms actually describe. For instance: ...
10
votes
0answers
68 views

What are the methods of modeling Jahn-Teller effects in a solid-state material? [closed]

The Jahn-Teller effect arises from the symmetry breaking (distortion) of the local environment for the crystal field of 3d transitional metal ions.The distortion result in energy saving of electronic ...
6
votes
1answer
131 views

I get an error 'Conversion from Z-matrix to cartesian coordinates failed' in Gaussian09. Any suggestions?

I'm doing my calculations using G09W. My OS is Windows 7 (64-bit). This error occurs in the frequency calculation. What's even weirder is that it seems to happen right after everything is already ...
6
votes
1answer
200 views

I get this error in Gaussian09 even though I have enough memory. “Error termination in NtrErr: NtrErr called from FIOCnC.” Any suggestions?

Whenever read-write files get to about 20Gb, I get this error. I know that, in G09, one rw file is limited to 16Gb. So I have set the software to write several RWFs. But it still keeps producing this ...
16
votes
3answers
430 views

How well can we model chemical synthesis?

The Materials Genome Initiative held a 2017 workshop, that led to this published report: de Pablo et al., "New frontiers for the materials initiative," npj Computational Materials 5, 41 (...
13
votes
4answers
176 views

What is the role of a chemist in computational chemistry?

As someone that wants to learn about computational chemistry, mainly for its use in energy storage materials (and therefore DFT as the most practical method of it), I am wondering what is essential ...
12
votes
1answer
107 views

Calculate spin-spin coupling based on previous GIAO calculation with Gaussian

I used GAUSSIAN/GAUSSVIEW to calculate the GIAO Magnetic shielding tensor (ppm), but I forgot to calculate the Nuclear spin-spin coupling J (Hz). How can I use the .chk file or .log file to calculate ...