Questions tagged [computational-chemistry]
For questions about computational chemistry.
252
questions
4
votes
0
answers
44
views
How do you run Gaussian's cubegen with multiple processes/threads?
The cubegen manual gives the signature of cubegen as:
cubegen nprocs kind fchkfile cubefile npts format cubefile2
Where:
nprocs
Number of shared memory ...
2
votes
1
answer
65
views
Absorption enthaply at different temprature in VASP software via ab initio molecular dynamics
Can we calculate the absorption enthalpy at different temperature in VASP software via ab initio molecular dynamics? If yes, how?
2
votes
0
answers
23
views
How do you determine the Pull Direction for a GROMACS umbrella sampling simulation to generate the most accurate free energy value?
I am trying to calculate the free energy between a protein and aptamer using GROMACS. Literature has indicated that one of the best way to do so is by performing an umbrella sampling MD simulation ...
1
vote
0
answers
19
views
When optimizing AMBER molecular dynamics parameters for protein-ligand interactions, how do I balance accuracy and efficiency?
Are there established practices for fine-tuning AMBER MD parameters in scenarios like mine, involving specific factors like mutant proteins or non-standard ligands? Has anyone successfully tackled ...
6
votes
0
answers
62
views
Updated book/chapter or review article with broad perspective on computational chemistry
I really like Frank Jensen's "Introduction to Computational Chemistry" because it provides a broad perspective, starting from first principles and extending to surveys of current methods for ...
3
votes
0
answers
25
views
Difference between pi-pi* transition and intramolecular charge transfer ICT in organic solar cells?
Intramolecular charge transfer (ICT( occurs from donors to acceptors in molecules. This phenomenon can be observed by visualizing the charge distribution patterns of HOMO and LUMO. While studying D-A ...
7
votes
0
answers
148
views
How do you calculate a full quartic force field for anharmonic frequencies?
I am an experimental chemist interested in anharmonicity—I am particularly interested in how to calculate anharmonic frequencies. The method I have the most experience with is (generalised) second ...
3
votes
0
answers
36
views
Running molecular docking simulation between two proteins
i did run docking simulation between two proteins on HADDOCK server
but now i want to run a simulation of complex (protein+ligand and protein) is it possible with HADDOCK and how can i use it
best ...
5
votes
1
answer
33
views
Database on calculated molecular vibrations including coordination complexes?
Related to the question Are there dashboards displaying experimental properties of molecular materials?
My question is: Are there any good, free, online databases on calculated vibrational modes of ...
4
votes
0
answers
41
views
OOM Error Calculating ESP Charges with Psi4
I'm trying to create a script to generate partial charges for a given .dat file.
As an example of my process, I'm currently using gromacs to run minimization on a box of propane. It outputs a min.gro ...
5
votes
1
answer
48
views
Identifying Key Chemical Series, Pharmacophores, and Scaffolds in a Large SDF Dataset
Hello fellow computational chemists!
I have encountered a challenge that I'd like to discuss and seek advice on. Suppose I have a substantial dataset in SDF format containing thousands of molecules. ...
2
votes
0
answers
43
views
How much should you shift fragment energies to create an energy level diagram?
I am creating an energy level diagram for the orbitals of [XeF₅] - molecule produced by two fragments: [Xe] 4+ and [F₅] 5-. There is a very large energy difference between the two fragments, so I ...
5
votes
0
answers
100
views
How can I properly plot this .cube file in Python?
Motivation for Post: I saw a question about how to plot a .cube file recently on the Mathematica & Wolfram Language Stack Exchange. It was marked "off-...
2
votes
1
answer
57
views
Energy of formation Vs Formation Energy Vs Heat of Formation Vs Energy Above Hull differences
I work in the field of materials science. While exploring the concepts related to material stability, I have come across terms such as "Energy of Formation," "Formation Energy," &...
2
votes
0
answers
67
views
Why the total energy calculated with GGA and metaGGA=SCAN are so different?
Why the total energy calculated with GGA and metaGGA=SCAN are so different?
My INNCAR file is
...
2
votes
0
answers
30
views
Error Simulating bulk diemthyl-ether with angle constrain in gromacs
I am trying to simulate bulk dimethyl ether with total constraints, both bonds and angles.
I am using lincs with constraint = all-angles in *.mdp file in gromacs.
I am attaching the error I ...
5
votes
1
answer
759
views
Why are Coulomb interactions modeled with just r, not r squared?
As far as I understand (as a complete beginner to molecular simulation and modelling), Coulomb's law for two atoms i and j is:
$$F=k\frac{q_iq_j}{(r_{ij})^2}$$
(the force is proportional to the ...
6
votes
2
answers
81
views
Analytical derivative of nuclear-electron attraction integrals over gaussian functions
I'm currently writing a program which evaluates the first derivatives of molecular integrals (over gaussian basis sets) with respect to the cartesian coordinates of the nucleus. I'm somewhat stuck ...
2
votes
0
answers
159
views
How to fit Smooth Overlap of Atomic Positions (SOAP) into Kernel Ridge Regression (KRR)?
I have a dataset of 64 chemical structures and their calculated energies, which I want to predict using machine learning. I need to use SOAP as a descriptor so that each structure would be an ...
3
votes
1
answer
59
views
Basis set of the overlap, kinetic, electronic matrices for RKS, UKS and ROKS code: O2 and NO2 cases
This question is a follow-up of this post.
The basis set to approximate the overlap, kinetic, electron repulsion and proton-electron attraction matrices is to be choosed by the user depending on the ...
4
votes
1
answer
73
views
What is the Atomic Orbital (AO) basis of the (overlap, kinetic, electric) matrices of the RKS, UKS and ROHF?
I am trying to code a simple UKS program. For RKS code, the AO basis functions are set of the spatial part of the wavefunction of every atom involved in the molecule. For ${H_2}$, this is the two ...
2
votes
0
answers
54
views
Help setting the input geometry into standard orientation for a single point g09 calculation
I am trying to run a single point g09 calculation for a pentagonal planar molecule with D5h symmetry using the B3LYP/LANL2DZ level or theory.
I would like the calculation to use symmetry, however the ...
3
votes
1
answer
222
views
Advice on converging a restricted open-shell DFT (rob3lyp/LANL2DZ) single point Gaussian09 calculation for a triplet spin-state molecule?
I am running a single point calculation on a fragment of a molecule using restricted open-shell theory. The fragment is in the triplet spin-state and of pentagonal planar structure (D5h). I ran the ...
3
votes
0
answers
59
views
How to choose the orientation of basis functions in a Gaussian calculation?
I am trying to run a g09 calculation for a D5h symmetric molecule using the HF/LANL2DZ level of theory. The molecule is pentagonal planar in geometry, and I would like to visualize the canonical ...
1
vote
0
answers
16
views
Complex formation site between polymer and external compound
I am studying about conjugated polymers in sensors applications. I want to conduct structural of complex and its properties (such as HOMO-LUMO etc.). Is there any software for docking that external ...
9
votes
1
answer
771
views
Why is chemical bond dissociation difficult to simulate?
I'm a student doing quantum computing research trying to understand better what kinds of chemical systems are difficult to simulate on a classical computer and for which quantum simulation could be ...
1
vote
0
answers
48
views
Resource recommendation for learning matter modeling for a novice [closed]
I am a starting researcher in condensed matter physics. Hence, i have all the necessary background regarding quantum physics, and condensed matter physics. But what i want to do is to learn different ...
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 ...
3
votes
1
answer
106
views
VASP - suggestion on optimization of interfaced slabs
I am currently trying to create a CeO2 (100) / CoO (200) interface for a bulk oxide.
Unlike what I usually see, slabs stacked on the z-axis ("vertically"), I would like to stack the two &...
2
votes
1
answer
57
views
Establishing AAs protonation - which software to choose?
I need to establish protonation of aminoacids composing the studied protein. To my knowledge there are H++ and Pdb2pqr servers that are able to do that. However, they yield quite much different ...
2
votes
1
answer
162
views
How to calculate ionic radius in Gaussian?
I am working with lithium sulfate and am wondering how to figure out the ionic radius of the Li atoms. I did a DFT/B3LYP/G-311G ++ (d,p) calculation.
6
votes
1
answer
483
views
How accurate is DFTB+ in reproducing DFT result?
I am trying to study the shifts in the Fermi level of a graphene - ssDNA hydrid system, depending on the identity of the ssDNA. The structures are obtained from a classical MD simulation, and cannot ...
3
votes
1
answer
95
views
How do you make sure the reaction is diffusion limited (Gaussian)?
I am investigating the mechanism of hydrogen transfer between a hydroperoxyl radical and an aromatic hydroxyl group donating it. The problem is that I cannot find transition state structure. So far, I'...
4
votes
0
answers
64
views
What kinds of calculations have to execute for obtaining the velocity of the chemical precipitation of some salt on a thread from a water-salt mix?
Some salt ($\ce{NaCl}$ or another one) of $m$ grams is dissolved in some water of $V$ liters. After mixing up, the substance looks transparent or to put it more exact as a chemical solution. Then a ...
2
votes
0
answers
71
views
Determination of fluorescence lifetime by using oscillator strength?
$$k_r=2.88\times10^{-9}n^2\langle\tilde{v}\rangle_\text{av}^2 \int\varepsilon_{\tilde{v}}d\tilde{v}$$
$$f=4.32\times10^{-9}\int\varepsilon_{\tilde{v}}d\tilde{v}$$
I plan to use the above equations to ...
2
votes
2
answers
118
views
What is the minimum number of atoms required to run a Monte Carlo simulation meaningfully?
What is the minimum number of atoms required to run a Monte Carlo simulation meaningfully?
5
votes
1
answer
266
views
Error in UCCSD(T) Calculation in PySCF for S atom?
I am trying to calculate the UCCSD energy and UCCSD(T) correction for the S atom using PySCF. Using sto-3g basis set, I found that the UCCSD calculation would be completed successfully, but the ...
5
votes
0
answers
69
views
Loading Gromacs Trajectories using ASE (Atomic-Silumation-Environment)
I am trying to study some trajectories generated using Gromacs. I have been given tpr and pdb files of the simulation. What might be the most effective way to load them in ASE? Or is there way to ...
3
votes
2
answers
57
views
How to remove f functions from a cc-pVQZ basis set in CFOUR format?
I need to perform a calculation on a BN system with using cc-pVQZ stripped of the f functions. However, I don't really know how to do this and I have not been able to find a revised one already. I ...
3
votes
0
answers
41
views
Quantum ESPRESSO vs AMS [closed]
Both Quantum ESPRESSO and AMS (Amsterdam Modeling Suite) can be used to perform simulations of the materials and their properties. What are the differences between AMS and Quantum ESPRESSO? Is it ...
3
votes
0
answers
24
views
What are the general principles of setting parellelization options in VASP?
I am doing the computation of a perovskite system of 20 atoms in a supercell.
I know that options like NPAR and KPAR, and the numbers of cores used for computations can greatly affect the efficiency, ...
8
votes
1
answer
281
views
What should be the criteria of convergence over ENCUT?
For example, should a computation be considered to converge when the differences in energy per atom are smaller than a certain value?
2
votes
0
answers
28
views
Is it a good idea to always read charge and orbitals if I am working on similar systems?
For example, if I am investigating of the elastic properties of a crystal by doing a lot of computations of crystals applied with small strain.
Also, is it a good idea to read charge and orbitals when ...
8
votes
1
answer
339
views
Does PySCF use chemist's or physicist's notation for 1- and 2-electron integrals?
I am trying to evaluate 1- and 2-electron integrals using PySCF trough this code:
Hij = mol.intor('int1e_kin') + mol.intor('int1e_nuc')
Vijkl = mol.intor('int2e')
...
1
vote
1
answer
108
views
Why are my HPC calculations not running in SLURM?
Question
I submitted 2 jobs for the first time in a HPC to test it. It's a very light calculation that takes almost 5 minutes or less in my desktop. But it has been almost 1 hour since I submitted the ...
4
votes
1
answer
174
views
Smooth overlap of atomic positions (SOAP) output for Gaussian Process regression in sklearn
I want to use the smooth overlap of atomic positions (SOAP) as a descriptor to represent the atomic environment of a specific atom to predict chemical shifts. I have generated averaged SOAPs for ...
1
vote
1
answer
211
views
MPI error in Quantum ESPRESSO during optimization of particularly large molecules
I have been recently trying to optimize a particularly large system Mo2C in QE using the pw.x executable and encountering the following MPI errors:
...
4
votes
0
answers
41
views
How to calculate charge transfer (CT) in NBO Analysis using Gaussian 09?
I have a coronene structure, functionalized with hydroxyl groups.
I want to calculate the charge transfer (CT) between the coronene and the hydroxyl group, using NBO analysis as implemented in ...
4
votes
1
answer
128
views
Question on CPU timings and/or ratio resp. DLPNO-CCSD(T) and wB97X-2-D for protein of ca. 1M Da
Suppose one has a protein of about 1 million Da, and performs a DLPNO-CCSD(T)/def2-TZVPPD/AutoAux/AutoAux calculation, and then performs a wB97X-2-D/def2-TZVPPD/AutoAux/AutoAux calculation, both via ...
2
votes
0
answers
59
views
DFTB+ parameters for additional atoms
I wonder if anyone here knows about additional atomic parameters and interactions sets for atoms not currently present in the DFTB parameter sets. The parameter sets are located at DFTB parameters. ...