Questions tagged [computational-chemistry]

For questions about computational chemistry.

Filter by
Sorted by
Tagged with
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 ...
leeman's user avatar
  • 141
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?
Rubi Agrawal's user avatar
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 ...
user avatar
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 ...
Don Aborah's user avatar
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 ...
sunfishstanford's user avatar
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 ...
diamond999's user avatar
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 ...
isolated matrix's user avatar
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 ...
Imane Yamari's user avatar
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 ...
agaitaarino's user avatar
  • 1,501
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 ...
asquith's user avatar
  • 41
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. ...
Poccia's user avatar
  • 53
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 ...
Kevin Freddo's user avatar
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-...
D_C's user avatar
  • 113
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," &...
hamza's user avatar
  • 21
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 ...
babak ab's user avatar
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 ...
anupama's user avatar
  • 121
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 ...
Cro's user avatar
  • 51
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 ...
dav267's user avatar
  • 133
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 ...
jessss's user avatar
  • 21
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 ...
mle's user avatar
  • 733
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 ...
mle's user avatar
  • 733
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 ...
Kevin Freddo's user avatar
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 ...
Kevin Freddo's user avatar
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 ...
Kevin Freddo's user avatar
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 ...
Duy-Minh Hoang's user avatar
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 ...
Solarflare0's user avatar
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 ...
kal92's user avatar
  • 11
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
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 &...
Salvios's user avatar
  • 63
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 ...
farmaceut's user avatar
  • 536
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.
Shelby Franke's user avatar
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 ...
Hemanth Haridas's user avatar
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'...
farmaceut's user avatar
  • 536
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 ...
SFriendly's user avatar
  • 877
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 ...
Yasir Ulutas's user avatar
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?
user366312's user avatar
  • 2,352
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 ...
Hemanth Haridas's user avatar
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 ...
Saha_1994's user avatar
  • 969
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 ...
farmaceut's user avatar
  • 536
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 ...
Zeina 's user avatar
  • 31
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, ...
蕭力諶's user avatar
  • 401
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?
蕭力諶's user avatar
  • 401
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 ...
蕭力諶's user avatar
  • 401
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') ...
Rodrigo Sánchez Martínez's user avatar
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 ...
Camilla's user avatar
  • 1,957
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 ...
C_Swann22's user avatar
  • 573
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: ...
Avrana's user avatar
  • 79
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 ...
Bruno's user avatar
  • 233
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 ...
Kanghun Kim's user avatar
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. ...
bhopshang's user avatar
  • 233

1
2 3 4 5 6