Newest Questions
4,805 questions
1
vote
0
answers
13
views
Why does the magnetic structure change after non-collinear AFM band calculation in SIESTA?
I'm performing a non-collinear antiferromagnetic (AFM) band structure calculation using SIESTA on a Mn-based compound. However, I noticed that after the calculation, the magnetic structure changed: ...
- 61
2
votes
1
answer
40
views
Getting the order of execution of functions in LAMMPS
I have written an MD code myself to model a system of interacting charged particles. I am using Nose-Hoover thermostat. To validate my results, I have also modeled the same system in LAMMPS taking ...
- 183
4
votes
1
answer
60
views
Why double hybrid functionals produce different results in different packages?
I'm calculating the single-point energy of the Nitrogen dimer using a double-hybrid functional called "mPW2PLYP". This is a widely known double hybrid functional available for use in ...
- 181
3
votes
0
answers
13
views
The best choice of the Extension Shell in the Partial Hessian Vibrational Analysis (PHVA)
I am interested in vibrational analysis (Hessian) for the QM part of a QM/MM system, with the MM part being frozen. This is (usually) done using the Partial Hessian Vibrational Analysis (PHVA). If I ...
- 703
1
vote
0
answers
14
views
does anyone has the Potential file for the U - Zr Alloy?
Does anyone have the interatomic potential file for the U-ZR alloys? If yes, please share with me
- 49
3
votes
1
answer
134
views
Different results from two different integral routines
I've been working on a Hartree-Fock code called Planck, where I've implemented overlap and kinetic integrals using three different methods: THO, OS, and Hermite, following Chapter 9 of Molecular ...
- 4,144
3
votes
0
answers
47
views
Is there a typo in this open shell CCSD formula from Scuseria?
In Scuseria (1991), the author provides an implementation of open-shell restricted Hartree-Fock coupled-cluster singles-doubles in the Appendix.
Equation A.5 is given as:
$$
\begin{equation}
v^{\...
- 31
5
votes
1
answer
107
views
What is the method to examine atomic bonds and hydrogen interactions after an optimization calculation?
I am currently working on a project to construct an energy profile for a reaction, which involves many optimization calculations. As a result, I need to verify numerous structures and ensure the ...
- 111
4
votes
0
answers
27
views
Calculating the IR spectrum of a system in an arbitrary vibrational state with propagation-based methods
In the context of nuclear quantum dynamics, there are two ways of computing an IR spectrum for the vibrational ground state $|\Psi_0 \rangle$:
energy differences and transition dipole moments, or
...
- 448
5
votes
1
answer
572
views
After optimizing all of my structures, I discovered that one of them needs tighter convergence criteria. Is it mandatory for all to be re-done?
I'm working on optimizing a transition state for proton transfer and C-N bond cleavage.
I successfully obtained a transition state with a negative frequency corresponding to the expected process. ...
- 111
1
vote
1
answer
36
views
How to Calculate Quantum capacitance from DOS(VASP-DFT) of Graphene 4x4 monolayer?
My query might be very basic but I am not able to figure out where I am going wrong in calculating Quantum capacitance for the Graphene 4x4 monolayer.I am using vasp to calculate Density of states.
My ...
- 111
1
vote
0
answers
10
views
How to consider LHE and RHE for a given system in tbtrans (SIESTA)?
I am working on Janus MXenes right now. I want to use TBTrans-3.2 to compute the systems' transport properties. Could you please clarify which part will be the scattering region (SR) and how to take ...
0
votes
0
answers
12
views
There are problems when calculating the projected transmission rate using tbtrans
I am calculating a four-electrode single-molecule device based on a graphene-naphthalene molecule-graphene structure, which has a total of 338 atoms. I have already calculated the transmission rate. ...
- 21
0
votes
0
answers
9
views
I encountered the problem of band crossings when using TBTrans to calculate the projected transmission
I'm calculating a four-electrode single-molecule device based on a graphene-naphthalene molecule-graphene structure, which consists of a total of 338 atoms. I've already calculated the transmission ...
- 21
2
votes
0
answers
12
views
Problem when performing the calculation of projected transmission using TBTrans
I am calculating a four-electrode single-molecule device based on graphene-naphthalene molecule-graphene, which has a total of 338 atoms. I have already calculated the transmission rate. When using ...
- 21
2
votes
0
answers
14
views
How may I make a heterostructure from two different crystal structures with different shapes?
I have two different crystal structure and I want to create a heterostructure but the problem is that the shape of their unit cell is different. I know how to combine them if their unit cell are same....
- 21
3
votes
0
answers
61
views
Using libmsym for symmetrizing Molecular Orbitals
I have written a Hartree-Fock code (mostly as a learning exercise, and to use it as a tool in my lessons) that utilizes the libmsym library (following the advice here) to identify the symmetry of ...
- 4,144
3
votes
0
answers
20
views
Symmetry error in slab relaxation calculation with quantum espresso pw.x
I´m doing ralaxation calculations in some slabs to compute adsortion energies of Li in several adsortion sites. In some outputs from QE pw.x I get the following error:
...
2
votes
0
answers
50
views
NEB convegence problem with LUA engine in siesta
I am having problems to converge NEB calculations with LUA engine within SIESTA. The NEB forces remain high after many days of calculations. Many of my NEB calculations so far never converge. I tried ...
- 123
3
votes
0
answers
20
views
How does CRYSTAL14 properties use fort.98 to generate 3D charge density from ECH3?
I’m using the properties program from CRYSTAL14, and I want to understand how the ECH3 function processes the fort.98 file to compute the 3D charge density grid (i.e., the values written in DENS_CUBE....
10
votes
1
answer
204
views
How to specify isotopes in PySCF?
I'm using Qiskit-Nature with PySCF driver to calculate PES of simple molecules. And my question is about PySCF.
If there is a way to specify isotopes of atoms in PySCF, how may I do it?
PS: I've ...
- 103
9
votes
1
answer
91
views
Recent Advances in Density Functional Theory (Mathematical Side)
I have been learning DFT with resources such as Engel and Dreizler's book, and it seems that higher level mathematics such as functional analysis appear frequently, especially in the more foundational ...
- 403
1
vote
0
answers
32
views
How to read LAMMPS trajectory file when the number of atoms is not same as that in topology file?
I'm running a LAMMPS simulation with 2197 atoms, when writing the trajectory file, I only write the information of atoms that enter a specific region in the simulation box. Because of this the number ...
- 199
7
votes
1
answer
63
views
How to find the symmetry unique atoms of a molecule with a specified point group by any software or tool?
Suppose, for cyclopropane, the highest point group is D3h. However, I want to define it to its subgroup C2v for my NwChem calculation. But for specifying a different point group NwChem requires the ...
- 71
2
votes
0
answers
34
views
Omega variant of B2PLYP functional in case of charge-separated (zwitterionic) ground state
I wonder if B2PLYP and wB2PLYP functionals are similar for the normal ground state but wB2PLYP might be better in the case of charge-separated (zwitterionic) ground state?
- 703
3
votes
1
answer
91
views
How do I input k-path in the bands.in file?
I use seek k-path website to find the suggested k-path and its coordinates for band structure calculation. Following is the input file (bands.in)
...
- 193
5
votes
1
answer
49
views
How do I find the updated cell parameters?
After performing relax-vc calculations for any system, I would like to find the cell parameters celldm(1), celldm(2), celldm(3), etc. using my output file. But I don't understand which of the values ...
- 193
5
votes
2
answers
220
views
How to Keep Spin Vectors Fixed in Non-Collinear Antiferromagnetic SIESTA Calculations？
When using SIESTA to perform structural relaxation on non-collinear antiferromagnetic systems, I would like to incorporate magnetic effects while ensuring that the spin directions of certain atoms ...
- 61
2
votes
0
answers
20
views
How can I access the original xlo/xhi box bounds from a LAMMPS dump file in MDAnalysis?
I'm working with a LAMMPS trajectory in MDAnalysis using the default LAMMPSDUMP format, like this:
...
- 199
4
votes
0
answers
39
views
Fourier transform expression of the coulomb interaction between the sheets(graphene)
Do you know any expression which can present the Fourier form of coulomb interaction between two graphene sheets separated by vacuum (as the figure above)?
Here S(w) is the transmission equation. ...
- 41