Questions tagged [mathematical-modeling]
Questions related to mathematical modeling of materials systems.
108
questions
5
votes
1
answer
46
views
Correct expressions for one-center Gaussian radial two-electron integrals
I am working on a new atomic solver library, for which I need expressions for the one-center two-electron integrals. Pitzer has furnished the necessary expressions for Slater-type and Gaussian-type ...
1
vote
0
answers
20
views
Log_linearization [migrated]
Can any one please suggest me how to log-linearize the utility maximization of the Ricardian households in a closed economy DSGE model. The equation is bellow :
5
votes
1
answer
134
views
Evaluating analytic gradients for overlap integrals
I have been working on a toy HF code for some time. From this reference I am trying to work out an analytic expression for the first derivatives of overlap integrals.
The original expression for ...
3
votes
0
answers
26
views
Convex hull calculations with WIEN2k
How does one do the calculation of a convex hull with the WIEN2k code for the formation energy of alloys?
4
votes
0
answers
39
views
Gauge constraint in the definition of the Z2 invariant
Cross-posted at Physics.SE.
In Fu and Kane's paper from 2006, the authors define the $\mathbb{Z}_2$ invariant for time-reversal invariant topological insulators as an obstruction to Stoke's theorem,
$$...
5
votes
2
answers
151
views
How is the equation for the position of a virtual site derived?
Cross posted on Math SE
I am trying to understand virtual sites in MD simulations, and I came across this configuration:
Here, coordinate $\mathbf{s}$ represents the virtual site, which is formed by ...
10
votes
1
answer
356
views
Inconsistencies in a famous point group table
The original resource
"Tables for Group Theory" By Peter Atkins, Mark Child and Courtenay Phillips (1970). I have been using this 2008 version, in which Pg. 37 has the following table:
Some ...
3
votes
1
answer
50
views
Is the following description talking about the simulation box or the boundary condition?
Is the following description talking about the simulation box or the boundary condition?
To prevent the argon atoms in the gas phase from running off to
infinity, we enclose the system in a "...
5
votes
0
answers
101
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-...
3
votes
1
answer
77
views
Cubic spline interpolation through a series of images
Cubic spline interpolation is used in a lot of methods like Nudged Elastic Band (NEB), or String or Growing String Method (GSM) either to reparameterize the path, or to estimate the position of the ...
2
votes
1
answer
55
views
What is the meaning of a "free energy per unit surface area" vs "pressure" plot?
I don't understand this graph. From the graph, does it mean that at high pressures the oxide layer is stable and this stability reduces at low pressures and temperatures?
4
votes
0
answers
89
views
Why are some point groups for single atoms listed as Oh in the CCCBDB?
Does anyone know why the CCCBDB database lists the tin cation and antimony atom as having the Oh point group (please see the screenshot below from CCCBDB when I look at "All molecules sorted by ...
4
votes
0
answers
46
views
Why, if the potential is different from the Coulomb one, but has spherical symmetry, the eigenvalues of the system are non-degenerate? [closed]
Question duplicated from here https://mathematica.stackexchange.com/questions/287433/why-if-the-potential-is-different-from-the-coulomb-one-but-has-spherical-symme and here https://physics....
4
votes
0
answers
144
views
Help with Mathematica for optimizing a Gaussian basis set for a hydrogen atom in a magnetic field
This is a follow-up to my previous question: Optimal Gaussian basis set for hydrogen atom in magnetic field
Brief description of the problem
I would like to find the ground and excited states of the ...
2
votes
1
answer
41
views
Why is strain rate used instead of absolute strain for modelling stress relaxation of a viscoelastic material using the Maxwell model?
I don't quite fully understand the derivation of the Maxwell model for prediction of stress relaxation for viscoelstic materials. usually the governing equation can be described in terms of strain ...
4
votes
0
answers
74
views
Predict if a group will be acceptor/donor of H bonds [closed]
Context:
I am a student of cheminformatics, and I am trying to solve a problem in predicting whether atom can be an acceptor/donor of hydrogen bonds, using Python programming language and RDKit ...
3
votes
1
answer
76
views
Molecular integrals in spherical harmonics?
I have been trying to write a basic HF code that can compute the MOs and total energy of a molecule. I followed the following reference for overlap and kinetic integrals. During a conversation with a ...
5
votes
0
answers
85
views
How can one determine the time-step for imaginary time evolution?
The real-time evolution of a quantum system described by the state $\psi(t)$ is formally given by,
$$ \psi(t+\Delta t) = \exp^{- i \hat{H} \Delta t / \hbar} \psi(t)\tag{1}$$
For dynamical simulations ...
4
votes
1
answer
108
views
Inverse Matrix Calculation
I am working on the charge transport calculation, based on the non-equilibrium green function method, for the Hamiltonian with the non-orthogonal basis. Suppose $H$ and $S$ are the Kohn-Sham ...
3
votes
1
answer
71
views
Estimate the time between different residues?
Imagine now, I have 5 structures, and I want to estimate the time when the 5 structures are closer to each other at the same time. I thought I could extract the distance between the different ...
7
votes
0
answers
50
views
Database (or table) of functions obeying symmetry properties of space groups
I'm looking for a table or database containing mathematical functions that belong to crystallographic space groups by having the same symmetry elements.
Imagine a 2D square planar lattice. It has a 4-...
6
votes
1
answer
76
views
What is the X in Almlöf and Taylor's "Unified treatment of energy derivatives?"
I have been studying the possible methods for basis set optimizations. One notable paper is "Energy-optimized GTO basis sets for LCAO Calculations. A Gradient Approach" by Knut Faegri Jr. ...
5
votes
1
answer
198
views
Constrained optimisation on a hypersphere
I am currently trying to implement a GS2 (Gonzalez-Schlegel second order) IRC algorithm in a python code. I am following the original paper ref(1).
The main problem is in the constrained optimisation ...
2
votes
1
answer
96
views
Understanding the complexity of geminal-based wavefunctions
Cross-posted on Reddit.
I have been reading through [1] to get a better understanding of geminal-based methods. Some short passages are included below:
The occupation of each orbital in the expansion ...
3
votes
0
answers
177
views
Understanding derivation of geminal based methods [closed]
Cross-posted on Reddit.
I have been reading through [1] to get a better understanding of geminal based methods. This is one of two questions that I'm asking about that paper on this site, the other ...
5
votes
1
answer
624
views
Numerical Evaluation of Hessian?
I am trying to understand the construction of Hessian matrix using numerical differentiation (NUMFREQ option in ORCA for example). I understand that Hessian is a square matrix containing second ...
5
votes
1
answer
105
views
How to Derive the Kong-Chakrabarty Mixing Rules
Background
In the world of atomistic modeling with classical force fields, one is often given a force field defined by like interactions (e.g. argon-argon interactions). If one is working with a ...
3
votes
1
answer
82
views
Decay rate of DSD-PBEPBE-D3BJ
I gave up on the wB97X-2(TQZ) functional, which would have an exact decay rate of -1/r, re: my previous question and decided to use DSD-PBEPBE-D3BJ, which is available on vanilla Psi4. I then looked ...
2
votes
0
answers
57
views
how to incorporate size polydispersity in colloidal suspensions [closed]
Basically, I'm working on spherical colloids and investigating the phase behavior of these monodisperse colloids using Monte Carlo simulations (Fortran 77 coding). Now I want to introduce size ...
1
vote
0
answers
55
views
Question on possible cancellations in the computation of Dyson orbitals from FCC/FCI wavefunctions [closed]
Dyson orbitals are mathematically meaningful, in the sense that they, unlike the DFT orbitals for instance, are actual observables of the wavefunction, instead of being "the other way around"...
1
vote
0
answers
50
views
Question on known expression (see post) of exact XC functional in analytic (but not closed) form [closed]
This paper analytically solves the "inverse DFT problem" of mapping the ground-state density to the exact XC functional. And the exact density, i.e. the FCC/FCI density, can be expressed as ...
2
votes
1
answer
91
views
Is there an analogue of the K-S theorem for double hybrids
In my previous question, I was basically asking whether the results of the double hybrid using the exact XC functional are the same as those of just the exact XC functional.
Even that sentence is hard ...
3
votes
0
answers
32
views
Is this method of obtaining only the exact "first excited state" density correct, under select conditions that makes it mathematically non-ambiguous? [closed]
(For context, see my other question here; this has been disproven (albeit contestedly), and is not even well-defined in the first place, so I'm going to ask a slightly different (and well-defined) ...
3
votes
1
answer
159
views
Is this method of obtaining only the exact "other" densities, and no other properties w.l.o.g, correct?
In Kohn-Sham DFT, the exact density of the ground-state wavefunction is given by summing the squares of the filled K-S orbitals.
Suppose now that one only needs the exact density (and nothing else w.l....
7
votes
1
answer
127
views
Is an "exact" double hybrid density the same as the "exact" DFT density?
Double hybrid approximate functionals have "unoccupied" Kohn-Sham orbitals in their formulations due to their MP2 component; however, the "exact functional" depends only on the &...
3
votes
1
answer
2k
views
Why is the Alphafold PAE (predicted aligned error) not symmetric?
Cross-posted on AI Stack Exchange.
We are running alphafold2 multimer on Google Colab to predict the association between two proteins.
It generally works fine, but we get an asymmetric PAE plot for a ...
5
votes
2
answers
299
views
How to apply FIRE to many atoms where P=F·v seems to be a vector rather than a scalar?
Introduction/Preamble
@SusiLehtola's answer to Basics of numerical energy minimization techniques used in molecular dynamics? mentions conjugate gradients, BFGS for energy minimization, Metropolis ...
5
votes
2
answers
274
views
Do the cc/pc/def2 basis sets mathematically converge to the CBS limit, assuming exact CI/DFT?
The cc-, pc- and def2- basis sets are often described as "systematic", in the sense that the results of these basis sets at different cardinal numbers (i.e. the number of ζ's) can be ...
12
votes
3
answers
3k
views
Artificial intelligence is a hot topic, but should I pursue it if I'm interested in Matter Modeling?
Since my previous question seems to have a large number of views, I would like to ask my central question.
I am a student now and should decide my major to research soon. AI/ML is a very hot research ...
1
vote
0
answers
104
views
What kind of chemistry background to get to pursue computational chemistry from an IT background?
What kind of chemistry background to get to pursue computational chemistry from an IT background?
I have no background in chemistry, and there might be others like me. But I know about computer models....
5
votes
1
answer
100
views
Cross-post: Matrix elements <n,k|x|n',k'> for Bloch states
Cross posted at Physics.SE
I believe this is just elementary QM, but I'm getting awfully confused. The question is drawn from this paper on Wannier-Stark localization (but is self-contained).
Let:
\...
4
votes
1
answer
182
views
What is the curve fitting model used for the "Qubit™ Protein Assay Kit"?
This is the example of Qubit™ Protein Assay curve.
It relates the concentration (x-axis) with fluorescence (y-axis).
What is the curve-fitting method it use for this? Clearly this is not simple linear ...
8
votes
2
answers
1k
views
What does it mean when the first order correction energy is 0?
Suppose I have the following Hamiltonian to start
$$ H_0 = \begin{pmatrix} 0 & 0 & 0 & 0\\ 0 & 0 & 2 & 0\\ 0 & 2 & 0 & 0\\ 0 & 0 & 0 & 0 \end{pmatrix} $$...
4
votes
1
answer
165
views
What are the precise forms for the SU(2) and rotation matrices in VASP?
Some time ago, I conducted a discussion with Dr. Gui-Bin Liu on topological materials here. One relevant thing he said was:
"The format of trace.txt file is only a necessary condition to be used ...
6
votes
1
answer
165
views
Understanding derivation of discretized FIRE algorithm
Introduction/Preamble
@SusiLehtola's answer to Basics of numerical energy minimization techniques used in molecular dynamics? mentions conjugate gradients, BFGS for energy minimization, Metropolis ...
6
votes
0
answers
57
views
Basics of configuration initialization strategies used in molecular dynamics? [duplicate]
This answer to Basics of numerical energy minimization techniques used in molecular dynamics? mentions
conjugate gradients, BFGS for energy minimization
and
Metropolis Monte Carlo (and the) the ...
10
votes
1
answer
303
views
Basics of numerical energy minimization techniques used in molecular dynamics?
The question below describes my plan to make a basic molecular dynamics calculation using a Python script rather than a canned, self-contained program.
There seems to be three parts:
a model of the ...
6
votes
0
answers
471
views
How to get a Wannier function for a tight-binding model numerically? [closed]
I have a question about construction of a Wannier function for a tight-binding model. Let's say we consider the tight-binding model of a 1D chain with two atoms (site A and B in a unit cell). In k-...
17
votes
1
answer
1k
views
Do we know for sure that all atomic and molecular wavefunctions decay exponentially as r goes to infinity?
Slater type orbitals (STO) are considered to be more accurate than gaussian type orbitals (GTO) for atomic and molecular QM calculations because - among other reasons - they decay with $e^{-\alpha r}$ ...
4
votes
0
answers
119
views
Can superconductivity be faked by using precise mathematics in a simulated circuit & specific materials in an actual circuit? [closed]
A type of electrical amplifier is proposed in which only magnetism and it’s analogous current exists for the most part. There is no voltage and no wattage (to speak of) in which electricity has to ...