Questions tagged [computational-chemistry]

The branch of chemistry dealing with the use of mathematical models to solve problems in chemistry.

Filter by
Sorted by
Tagged with
2 votes
2 answers
146 views

Starting configuration for Molecular Dynamics

I have to build a bilayer using a structure obtained in a previous MD simulation. In the structure, some chains of the molecules cross the simulation box as shown below. I added water layers to the ...
Ema 's user avatar
  • 21
0 votes
0 answers
17 views

pp.x command in quantum espresso

I am trying to generate the .pp and .cube file using pp.x -i (input file name) command. But, the .cube file is not getting generated and the output.log file also shows an error message. this is my pp....
Bilal Ahmad's user avatar
4 votes
2 answers
364 views

Choice between DAE or ODE formulation for chemical systems

Consider a simple ODE system describing the evolution of two chemical species undergoing the reaction $A = B$ : $$ \frac{dn_A}{dt} = - k * n_A $$ $$ \frac{dn_B}{dt} = k * n_A $$ We can discretize ...
Anon_Chem's user avatar
1 vote
0 answers
33 views

Algorithm for constructing a supercell using an arbitrary integer matrix [closed]

Software like Diamond and VESTA allow users to build a supercell by extending the structure along the basis vectors of the latttice. However, I'd like to do something more general and build a ...
brainandforce's user avatar
0 votes
0 answers
28 views

Calculating radial distribution function in a cylindrical pore

I have tried calculating radial distribution functions in a cylindrical pore geometry. The pore has z-axis periodicity and it does not have x and y periodicity. But I got some troubles on that. The ...
Yeongkyu Lee's user avatar
0 votes
0 answers
126 views

Optimising transition state geometries using gaussian 16

I'm attempting to carry out transition state optimisations in Gaussian and a lot of my calculations are coming up with errors, the main one being no analytic 2nd derivatives for this method. How can ...
Beth's user avatar
  • 11
-1 votes
1 answer
65 views

Does it make any sense to acquire some sort of knowledge about manufacturing or engineering for computational design optimization? [closed]

Does it make any sense to acquire some sort of knowledge about manufacturing or engineering for computational design optimization? I've viewed computational science to be an abstract science ...
mavavilj's user avatar
  • 427
1 vote
1 answer
261 views

How to perform FFT from plane-wave basis function coefficients to real space?

I have a 3D grid in real space of grid spacing $L$ and say 21 grid points in each direction, containing e.g. a charge distribution. This is stored as a numpy array of shape ...
Protocola's user avatar
1 vote
0 answers
41 views

Calculating the species mass consumption from implicit reaction-term in diffusion-reaction equation

The 1D diffusion equation with a chemical source term has the following form: $$\frac{\partial Y}{\partial t} = D \frac{\partial^2 Y}{\partial x^2} - k Y,$$ where $Y$ is the molar concentration of the ...
user306429's user avatar
0 votes
1 answer
764 views

Who can explain the minimum image convention in molecular dynamic simulations?

How to choose the cutoff radius so that the atoms do not interact with its periodic image? Especially when simulating macromolecules (proteins).
Muhriddin Mahkamov's user avatar
2 votes
1 answer
105 views

Quasi-Newton Method with a Transformed Hessian

I've recently came across an implementation of the BFGS algorithm but it has an additional step where the Hessian is transformed after the each update. This transformation is done so that certain ...
Unskilled's user avatar
  • 123
0 votes
0 answers
70 views

Examples of kinetic modeling with optimization techniques in Python

We are looking for studies or tools that implemented kinetic modeling with parameter estimation differential evolution or similar optimization techniques in Python. We are trying to understand what ...
Jean Carlo Machado's user avatar
6 votes
4 answers
994 views

How important is learning hardware/architecture for scientific computing?

Apologies if this is a bit of a soft, unclear, or opinion-based question. I'm a relatively new PhD student in a (computational) quantum chemistry group. My group develops and maintains a few software ...
tmph's user avatar
  • 337
0 votes
0 answers
43 views

Reference request: Textbook similar in structure to "Computational Materials Science: An Introduction" by June Lee for Computational Biology

In the Computational Materials Science: An Introduction by June Lee, he discusses molecular dynamics and density functional theory with examples from LAMMPS and QuantumEspresso, and explains LAMMPS ...
Josh's user avatar
  • 1
0 votes
0 answers
57 views

Moller Plesset Theory- MP2 correction

I have a molecule LiBH2 with C2V point group. I used 3-21G basis set and obtained RHF coefficients and symmetry of orbitals. I wish to find the set of orbitals (a,b,i,j)( where a and b are occupied ...
Tapasvi Bhatt's user avatar
2 votes
1 answer
239 views

Applications of Julia in Chemistry and Molecular Physics?

I was wondering if there are any Theoretical & Computational Chemistry (MM, QM) codes or publications out there that are based primarily on the Julia programming language?
epalos's user avatar
  • 121
1 vote
0 answers
32 views

Why is scoring a separate problem even after docking is done?

With regards to protein-ligand interactions, we talk about two concepts, docking and scoring. I understand that docking means to find the best orientation of a ligand at the active site of a protein, ...
Mahathi Vempati's user avatar
1 vote
0 answers
110 views

Tight binding model calculation with Extended Huckel Approximation

I've been reading Albright's Orbital Interactions in Chemistry. In the chapter on solids, he provided a general approach to find the band structure of a solid state system Now if we are to model a ...
Macrophage's user avatar
5 votes
0 answers
99 views

Evolutionary dynamics in vascularised tumors, PDE-ODE coupled system

I have to solve the following PDE-ODE system $$ \displaystyle{\partial_{t} n = \bigl[a(s) - b(s)(y - h(s))^{2} - d\int_{\mathbb{R}} n \, dy \; + \; \beta \, \partial^2_{yy} n \;} \\\\ \...
AleManara's user avatar
0 votes
0 answers
180 views

How to handle system of chemical reactions for a batch reactor SciPy solver

I have a system of chemical reactions where the rate equations represent a batch reactor model. The model is a system of ODEs which is solved with the SciPy ...
wigging's user avatar
  • 325
1 vote
1 answer
2k views

Getting started with Computational Chemistry

I´m now a chemistry grad student and I feel the need to get involved with computational chemistry and coding in the chemical field (in general). I have a very simple question: What is the best way to ...
TheChemist's user avatar
5 votes
3 answers
668 views

Maximize a function of an orthogonal matrix

I'm trying to write up a small code that, given a set of normal vibrational modes for a molecule, will convert them to localized vibrational modes. To do this I'm following the procedure from J. Chem. ...
Tyberius's user avatar
  • 1,023
0 votes
1 answer
130 views

Software to simulate molten salt flow and thermodynamic operations

I was curious if there was any software (preferably in C++, Java, and/or python) that could be used to simulate the following: Heat capacity of a fluid Heat transfer through a liquid and a solid ...
t1r3d's user avatar
  • 1
40 votes
23 answers
10k views

Good examples of "two is easy, three is hard" in computational sciences

I recently encountered a formulation of the meta-phenomenon: "two is easy, three is hard" (phrased this way by Federico Poloni), which can be described, as follows: When a certain problem is ...
Anton Menshov's user avatar
  • 8,592
2 votes
2 answers
656 views

What equation should I fit this set of data points to?

I have done an experiment Estimation of silver nitrate by potentiometric titration with standard KCl solution. A plot of $\dfrac{\Delta E}{\Delta V}$ versus Volume of KCl solution gives the ...
ntk47's user avatar
  • 91
1 vote
1 answer
65 views

What are the things I should keep in mind before doing an analysis of my gromacs simulation?

I did a liquid argon simulation at 100k. I forgot during analysis that I need to accommodate for the Periodic Boundary Conditions before doing any analysis which included distance. What are other such ...
fireball.1's user avatar
3 votes
1 answer
200 views

Non-linear flux interface condition - variational formulation

Context: I am working on implementing this paper and I am struggling to come up with a variational formulation for the Butler-Volmer interface conditions. To simplify my question I consider the ...
Bernd Schwarzenbacher's user avatar
2 votes
0 answers
92 views

Calculate activation energy of inter molecular displacement on a nano particle

I'm simulating a nano particle and wish to calculate the activation energy of the reaction where an atom is displaced from one point on the lattice to a neighboring void( might be at higher energy ...
fireball.1's user avatar
1 vote
1 answer
227 views

Oil/Water interface simulation using GROMACS

I would like to perform a decane/water(tip4p) interface simulation using GROMACS to analyse the fluctuation of the interface over a long period of time (50ps - 1ns). Though I have a rough idea on ...
Ramanathan Varadharajan's user avatar
1 vote
0 answers
65 views

Parallelizing molecular simulation with full configuration energy

First, let just me stress that I'm not a an expert in computation chemistry, so now the problem: We have GCMC molecular simulation, in the Grand Canonical ensemble, using the standard metropolis ...
dual0's user avatar
  • 11
0 votes
1 answer
60 views

Forcefield for studying a two phase system, consisting of noble gases and gaseous organic components

I am preparing to do a MD-simulation of a system which consists of two phases. The gaseous phase consists of $H_{2}$, $Ar$, $He$, $N_{2}$, $CH_{4}$ and $CO_{2}$. In this system, $p=0.355 ~bar$ and $T\...
Bdrs's user avatar
  • 151
0 votes
1 answer
130 views

List of Tersoff potentials?

I would like to know some GNU libraries or lists of downloadable Tersoff potentials. I couldn't find that many of them, and I also would like to know if there is a reason as to why this is the case.
CoffeeIsLife's user avatar
5 votes
1 answer
5k views

Solve rate equations with different reaction orders using SciPy ode

The image below illustrates the kinetic scheme I am trying to model. My first example focuses on the conversion of the wood-oil to non-volatiles and volatiles. The rate equations are as follows: $$ ...
wigging's user avatar
  • 325
2 votes
0 answers
397 views

What are alternatives to openbabel to convert CML file to image?

I recently started to perform natural resonance theory (NRT) analysis with NBO6. This program may automatically create a chemical markup language (cml) file with the treated structures. Unfortunately ...
Martin - マーチン's user avatar
0 votes
3 answers
173 views

Quantum Chemical Calculations is there a book for which method to use with what problem?

Does anyone know of a book that will outline which quantum chemical methods are appropriate for what problems? I am trying to make informed choices before I start using computational resources. It is ...
user2267258's user avatar
5 votes
1 answer
1k views

What determines the usual chemistry textbook plots of atom orbitals?

In elementary chemistry textbooks you often have pictures like the following one: Are there any conventions how to get them? I am not sure, but I guess that it are contour plots with only one iso-...
student's user avatar
  • 231
2 votes
1 answer
930 views

Is the Finite Volume Method conservative when the source term depends on the variable?

I trying to do a simulation where there are two quantities, $\Delta\theta(x,t) = \theta(x,t) - \theta_{o}$ and $\Delta\nu(x,t) = \nu(x,t) - \nu_{o}$. These quantities are chemical concentrations ...
user avatar
8 votes
0 answers
562 views

DIIS method to accelerate SCF convergence for stretched geometries

I am implementing from scratch an Hartree-Fock calculation in the STO-3G basis set to perform Born-Oppenheimer molecular dynamics. I have a Restricted Hartree-Fock procedure that can reproduce very ...
user avatar
-1 votes
1 answer
1k views

Gaussian and external basis set. How to write the input correctly [closed]

I have to use the LANL2DZpd basis set. It can be obtained from EMSL Basis Set Exchange. But trying to calculate simple Me2Se molecule : ...
XuMuK's user avatar
  • 101
6 votes
2 answers
2k views

Borrow computational power from machines around me

I'm a computational chemist working with little computational power and dealing with increasingly demanding chemical systems. My work machines are part of a big network (mixed with windows and linux) ...
HCSthe2nd's user avatar
  • 163
1 vote
0 answers
97 views

Time dependent self-consistent equations

I am facing the following problem. I need to solve numerically a set of coupled equations $$i\frac{d}{dt}f_{n}^{(i)}(t) = \left[U\cdot n(n-1) + \mu\cdot n\right]f_{n}^{(i)}(t) - \sqrt{n+1}\Phi_i^{*}\...
WoofDoggy's user avatar
  • 203
4 votes
1 answer
425 views

OpenFOAM precipitation/crystallisation reaction solver

This question is probably related to chemical engineers that are around this forum. I am looking into writing a solver for a precipitation reaction (struvite to be precise): A+B+C+D <=> E+F A, B, ...
user3800363's user avatar
12 votes
0 answers
137 views

Are there any standardized file formats for point group character tables?

Character tables are an important tool for symmetry analysis in many computational chemistry software packages. Are there any standardized file formats for point group character tables? This may seem ...
jvtrudel's user avatar
  • 221
1 vote
1 answer
86 views

Shifting points in fortran

Hello I am trying to shift points which have been previously generated in the square area. I am having a trouble with some additional conditions how they should have been shifted. ...
Beginner in fort's user avatar
1 vote
1 answer
32 views

Computational complexity of chemical dynamics for biological studies [closed]

How many CPU cycle is required to simulate a complete human body from it's very initial stem cell using classical algorithms and also is it possible to use similar algorithms for simulating stem cells ...
user16908's user avatar
3 votes
1 answer
2k views

Computer Build for Scientific Computing

I am currently a .NET software developer(SQL Server, ASP.NET, C#, MVC & Web Forms). In my spare time I'm researching different areas of environmental science. E.g.(Hydrology, Ecology, Atmospheric ...
IEnjoyEatingVegetables's user avatar
3 votes
1 answer
81 views

Is there any numerical application whose performance heavily depends on the division operation?

I am an undergraduate student majoring in computer science. Recently, I am interested in the division operation, which is not directly supported by some architectures. While some architectures ...
konjac's user avatar
  • 131
2 votes
1 answer
334 views

Global optimization methods in computational chemistry

I'm looking for a current and comprehensive overview (like a review article) of global optimization methods and their application in computational chemistry. Mostly I'm interested in geometry ...
tmartin's user avatar
  • 209
1 vote
0 answers
451 views

Center of Mass for molecule in Gaussian09

Gaussian09 prints coordinates where the local coordinate system as center at $C = (\Sigma n_i*X_i)/(\Sigma n_i)$ where $n_i$ is the Mass Numer of the atom. However, the Mass Number is not actually ...
drjrm3's user avatar
  • 2,119
2 votes
1 answer
864 views

Installing GPAW (Python DFT package) on a computing cluster

I'm in the middle of installing GPAW on a cluster but I'm running into problems building the "gpaw-python" interpreter, its job is to allow users to execute parallelised GPAW calculations. Following ...
Josh's user avatar
  • 253