I did a docking job and got a score of -5.69 kcal/mol and a value for MM-GBSA of -60 kcal/mol. I should analyze these results for my thesis but do not know how to comment on them. Does anyone have any insight about what this number might mean? Also on interactions, I wonder what can be said? A...
A few years ago there was a significant difference between CPUs and GPUs for performing calculations. It was quite clear when to go for a CPU or a GPU. Today there are GPUs with a very large number of cores and threads. I was away from the subject for a while. Is there still such a marked diff...
Does anyone know how to export a calculated Raman spectrum (my calculation was performed in Q-Chem) as a data-file, which can be used to plot the profile of a Gaussian-shape spectrum (in origin, excel, etc.), rathe than centroid data points?
I am looking for real data that concludes that implicit solvent models are definitively better than explicit solvent models for solvation free energy calculations. The motivation is a lovely reviewer who has said this. Evidence that they are not statistically different is also valid.
I am a theoretical physicist, so I am not very familiar with chemical names. The material I am interested in is $\ce{\alpha-NaMnO2}$. Would that be pronounced alpha sodium manganese dioxide? Is it important to say the $\alpha$? If it's relevant, I'm talking about Phys. Rev. Lett. 2020, 124 (19),...
The binding interaction energy during the adsorption of system B on system A can give an idea if the interaction is favorable (negative value) or unfavorable (positive value). Also, the absolute value represents the strength of the interaction. To determine the binding energy, the procedure is t...
I just saw I am a beginner in Density Functional Theory. What are some resources that could help me to learn the basics? and it reminded me that few years ago I'd asked How can I calculate the charge distribution of a water molecule? in Chemistry SE. There is some excellent advice there, but at ...
I am a DFT user and at some point in the future, I would like to write my own DFT code in Python to help gain a deeper understanding of DFT. As mentioned in a previous answer people have written their own DFT codes to understand more deeply how the theory and algorithms work. Where do I start? ...
In general, the CIF file with the crystal structure uses fractional coordinates for setting the atoms positions. In this type of coordinates, the atom positions are relative to the cell parameters. Another way to setting up the atom positions is to use Cartesian coordinates. In this type of coo...
I’m doing CCSD(T) calculations (for both transition and ground state geometries) on a series of small molecules in Gaussian09. Several similar calculations terminated without any problems but one TS calculation and one ground state calculation gave me this error after about 20 minutes. HESSIAN D...
The "curse of dimensionality" is an ubiquitous issue arise in both electronic structure and quantum molecular dynamics in computational chemistry, which refers to the exponential scaling of the computational cost for the systems of interest as the the number of their degrees of freedom. This...
So I am aware that many quantum chemical codes are able to add external fields to molecular calculations (for example to calculate polarizabilities). I also know that in some (eg. Orca) you have some control over the applied field, for example the strength. I am also aware that for example in Mol...
A lot of codes can do Quantum Mechanics/ Molecular Mechanics (QM/MM), and I am not naively asking if there is a program that can do QM/MM. What I am interested in is a program that is purely QM/MM meaning that it is designed and optimized specifically for this purpose. QM/MM is not an after-thou...
It is possible to use different types of thermostats in CP2K by specifying them in the corresponding section of the input as mentioned here. But there is no mention of specifying the barostat type in the manual. And I couldn't find any mention of the type of barostat being used in CP2K for NPT si...
A Brillouin zone is defined as a Wigner~Secitz primitive cell in the reciprocal lattice [1,2]. The construction of the first Brillouin zone for two different 2D lattices are shown below: The construction procedure for 3D lattices is basically the same as for 2D lattices. In the image below, ...
I've been doing some calculations for a water molecule, with a big uncontracted basis set: For oxygen I am using aug-cc-pV9Z supplemented with the "tight" functions from cc-pCV7Z. The s-type exponents are listed below, with all numbers up to 0.04456 constituting the s-type exponents of aug-cc-pV9...
I am not able to understand the math of this. After taking the Slater determinant as a wavefunction in Hartree-Fock, does the procedure to find wavefunctions remain the same as in the Hartree method i.e. the self-consistent method, or does it change?
In Kohn-Sham DFT calculations, Van der Waals corrections are often implemented in the structure optimization calculations because the typical functionals such as LDA and GGA are found to not treat long-range dispersion forces faithfully. But even in the context of Van der Waals corrections, there...
Does anyone have any knowledge or direction for creating a homemade Beowolf cluster for materials modeling?I would like to be able to run my own "hobby" simulations at home, so a low energy, low capital investment is what I seek. Raspberry Pi and ODroid use small, cheap, low energy single board c...
This might be a stupid question but to what extent will the initial configuration of a bulk phase geometry optimization calculation affect the final geometry? Most places say to start with experimental data but if there isn't any then is there a next best way of starting? Would populating your i...
Introduction It is know that the entropy of a system can be calculated using the formula (1): \begin{equation} \label{eq:entropy_Boltz} S = {k_B}\ln \Omega\tag1 \end{equation} where $k_B$ is the Boltzmann constant, $\Omega$ is the number of microstates of the system and $S$ is its entropy. Thi...
I'm new to DFT. Does the DFT take into account the spin-orbit interaction? On the one hand, this is a relativistic effect, perhaps DFT does not take it into account. However, on the other hand, I read that DFT is an extremely accurate method and the error of this method lies only in the inaccura...
As per @user1271772 's suggestion I am asking this question here again. I am trying to understand SCF cycle by trying to code up solved example from Quantum Chemistry by Levine (page 443, 5th edition). The problem is stated as: Do an SCF calculation for the helium-atom ground state using b...
This questions is not directly related to Matter Modeling, but it is related with the way we share/publish our research. Let suppose that someone deposited a manuscript in a preprint server. Then, he/she send the same manuscript to a Journal and it is accepted. In many journals you have to sign ...
In 2013 it was stated that v8.8 of DYNAMO was the version mostly used for the last 10-12 years. However I am struggling to find the code repository. It was developed by the LAMMPS developers and is used to mix eam potentials. DYNAMO was available on sourceforge.net, but now it seems to have be...
For systems with periodic boundary condition , plane wave basis is usually adopted. While for molecular systems, gaussian basis set is normally adopted. For periodic systems, the gaussian basis is transformed from real lattice space to momentum space (or k space) through Bloch function. My questi...
In other works and from the paper: Scaling Up Electronic Structure Calculations on Quantum Computers: The Frozen Natural Orbital Based Method of Increments. I have seen that the corrected CCSD correlation energy ($E_{c}^{CCSD}(corrected)$) is calculated from a CCSD computation of the truncated...
It is very well known that Kohn-Sham DFT underestimates bandgap. To get an accurate estimate of the bandgap, people often turn to Hybrid functionals (if they don't want to perform the actual quasiparticle calculation). I understand how the Hybrid functional, say for example, HSE works. It is an a...
How can an MM-GBSA(molecular mechanics energies combined with generalized Born and surface area continuum solvation) calculation be done from an induced fit docking (IFD) output? I performed an IFD and now I would like to run an MM-GBSA calculation, but when I try to start it, I get a warning mes...
In some science fields there are efforts to provide customized distributions, so you can get a productive environment with minimum set-up overhead. For example, in high energy particle physics, we have Scientific Linux, in mathematics we have MathLibre, and in cybersecurity research we have Kali ...
Kohn-Sham DFT appears to be so popular even though it is strictly a ground-state method - all calculations are done at 0 K. How then, is it so popular when describing materials that have real-life applications (possibly at room temperatures, or much higher temperatures)? Or am I mistakenly deemin...
Hartree-Fock density is free of self-interaction but lacks electron correlation effects, while the density from KS-DFT (using an xc functional or potential, both which are explicitly density-dependent and gradient-dependent) suffers from self-interaction errors but does incorporate electron corre...
I am studying a material which was first discovered in the early 1960s. It was reported then that it crystallized in the hexagonal $P6_3/mmc$ space group. Recently the material was synthesized again and found that it posses $P6_3mc$ symmetry. My own DFT calculations and another published work fin...
I am trying to input various charge states of various transition metal ion species into DFT calculations. The program requires an input of multiplicity = 2*(total spin) + 1. I have an undergraduate-level understanding of chemistry, up to the point where I understand orbital filling (1s,2s,2p,.....
When we calculate the band structure of certain material, we only have to calculate the value along the high symmetry point which enclose the Irreducible Brillouin Zone. Why the information lie in the IBZ is enough for us to know the information of entire BZ ? How to find the IBZ given a reciproc...
I find the experiments to be relatively easy to explain to my young cousins (who are planning to take scientific reserach as a career option). I find the theoretical works, especially the ones which require highly involved computational tools (read DFT) to be a bit difficult to explain to these 1...
I am using VASP package to model a material that has fractional site occupancies. An appropriate supercell (cif) was generated using supercell package and converted to POSCAR using VESTA. Then the structural relaxation was performed. Now I have the relaxed structure (CONTCAR) and I want to analyz...
I have a series of molecules calculated (DFT PBE0-D3/def2-TZVP performed in Q-Chem) and I'm ready to plot PES (Potential Energy Surface). I had to increase the grid to get rid of small imaginary frequencies in frequency calculations. All calculations were performed with SMD(water) (Geometry Optim...
I have a 9 x 9 x 9 periodic crystal structure in .xyz format and I would like to change it to .cif file format. Does anyone know how to convert the file? (e.g. 9 x 9 x 9 to unit cell then cif)
I was wondering if DFT (or specifically, LDA ?) can be exact in the limit of the homogeneous electron gas? In that case, shouldn't the self-interaction error perfectly cancel out? I realise that such a material might not exist, but in theory, is this true? I am aware that nothing is perfectly ex...
(This question is originally posted on physics stackexchange, but someone suggested me to post on this site, so there you go) I'm compiling the mathematical expression of SCAN (Strongly Constrained and Appropriately Normed) functionals' constraints, but apparently they are not very obvious from t...
Which quantum chemistry programs can compute natural bond orbitals (NBOs)? Which visualisation tools can I use subsequently to plot pictures of the NBOs?
Regarding the calculation of binding energy of molecular absorption on graphene, for which this has been discussed before, the process of calculating $E_{AB}$ is fairly straight forward, i.e., it is the single point energy of the relaxed complex AB. However, the calculation of $E_A$ and $E_B$ can...
Is it useful to energy minimize a docked pose before scoring it? I am speaking in the context of protein-ligand systems. Does the energy obtained after minimization correlate better with experimental affinity? Does someone have some published literature on that?
I am working on modelling of a material that has partial/mixed site occupancies. For example consider the tetragonal (P4mm) compound $\ce{(Ba_{0.67}Sr_{0.33})TiO3}$. Following is the atomic arrangement (visualised in VESTA) when looked along the Z axis. One of the ways to model such structures i...
I am looking for a free software able to draw 3d micelle configuration like the figure attached. Any suggestion?
I've found a page dedicated to quantum chemistry packages benchmarks, on GitHub, qmspeedtest. But most results there are several years old, and so probably outdated. Is there some place where we can find comparisons like these, but updated often, or at least more recently? I specified quantum che...
Some time ago I've read this news article, Mathematicians Discover the Perfect Way to Multiply, reporting a discovery published in 2019, where Harvey and Hoeven[1] found a algorithm able to execute multiplication in $N \log N$ steps. Compare with the $N^2$ we are used to while doing multiplicati...
In Kohn's Nobel Lecture, he claimed that: "In general the many-electron wave function $\Psi(\mathbf{r}_1,\ldots,\mathbf{r}_N)$ for a system of $N$ electrons is not a legitimate scientific concept, when $N\geq N_0$ where $N_0\simeq10^3$." He explains this in two ways: the first is that $\Psi$ cann...
I used VASP to do the band structure calculation on Graphene. Originally I would expect that the Dirac cone can be observed in the K point. If we watched qualitatively, we can see this behavior. But when I zoom in the region where the conduction band touches with the valences band, a bandgap can...
next day → last day (1407 days later) »