Messages

1

General Questions and Answers / Re: Why am I getting a small bandgap in density of states of graphene?

« on: March 24, 2017, 14:03 »

The k-mesh size is too coarse. You need to increase it for the geometry optimization and DOS calculations.

"bands_above_fermi_level=None," in the DOS calculations should not be "None".

2

General Questions and Answers / Re: GGA+U bandstructure calculation of monolayer mos2

« on: March 24, 2017, 08:21 »

The difference may be not related to the ATK version.

Please don't misuse DFT+U. Even some U value (in most case they may be very large) would give the band gap close to the experimental data, the physics meaning of this U may lose. To get a better prediction on the band gap, the DFT+U is not the best solution.

My suggest is to first try the use of linear response method in Quantum espresso or ABINT codes to estimate the U value for MoS2. The obtained U value would be more physics meaningful even such U values couldn't give the band gap to match the experimental results.

3

General Questions and Answers / Re: GGA+U bandstructure calculation of monolayer mos2

« on: March 23, 2017, 02:03 »

Please refer to the details of LDA+U (or GGA+U) methodology. The implementation of counting (calculating) occupation of local orbitals would also the affect the calculation results even if the same U is applied on the same orbital in the two different codes.

4

General Questions and Answers / Re: What's wrong with bader analysis in my calculation?

« on: March 22, 2017, 01:45 »

The ATK-DFT or VASP can generate the spin-up and spin-down charge densities. You need to store these charge densities separately into two different files.  And then use the same reference (ELF or all-electron charge density) to get the local volumes (basin) of every atom by using critic2 . In this way, you may be able to the local charge of spin-up and spin-down components for every atom. 

5

General Questions and Answers / Re: What's wrong with bader analysis in my calculation?

« on: March 21, 2017, 10:04 »

It cannot be done.

I recommend the use of critic2 to perform the Bader charge analysis even if the all-electron density is provided by the ATK-DFT code.
The critic2 (http://schooner.chem.dal.ca/wiki/Critic2) provides the superposition all-electron charge of atoms.  One can also use the ELF(electron localization function) as reference to determine the basin, local maximum and so on.

6

General Questions and Answers / Re: Optimization

« on: March 21, 2017, 09:54 »

Yes, you are right.

7

General Questions and Answers / Re: Optimization

« on: March 19, 2017, 06:39 »

1. Look at the manual :
http://docs.quantumwise.com/manuals/Types/OptimizeGeometry/OptimizeGeometry.html#optimizegeometry-f

With the default values for
"target_stress" , "constraints", and "optimizer_method" in OptimizeGeometry().
2. No.
3.  No.

8

General Questions and Answers / Re: charged molecule between electrodes

« on: March 18, 2017, 06:36 »

The center region or scatter region contains the conductor (i.e., the molecule in your case) and the buffer layers. The latter ones have the same components with those of electrodes. The spin setup for buffer layers should be same as those in the electrodes. The spin setup of molecule could be different to the spin setup of electrodes.

9

General Questions and Answers / Re: how do we calculate spin transport in finite bias?

« on: March 17, 2017, 01:51 »

"what's more, should we set initial-spin again when calculating transport under finite bias?"
Sure, you should set the initial spin again in the calculation of finite bias.

The code one is preferable.

10

General Questions and Answers / Re: how do we calculate spin transport in finite bias?

« on: March 16, 2017, 09:53 »

These three ways mentioned in your post mean just the different starting points (in particular, the charge density) for the anti-parallel spin configurations at finite bias. The different starting points may affect the convergence speed (e.g., the computation time).  In principles, all of them should give a nearly same result for a given finite bias, the little difference may appear because of the numerical noise. The most safe way for the convergence problem is the third one. 

11

General Questions and Answers / Re: Bloch state velocity in silicene = 0 !!!

« on: March 14, 2017, 10:38 »

Please show the band structure of your calculated system. From the band dispersion in such plot, one may judge whether your result is reasonable or not.

12

General Questions and Answers / Re: electron density

« on: March 13, 2017, 02:09 »

Hi,
how to calculate the electron density of any 2d material (example: MoS2) at the edge of conduction band?
Thanks in advance....

Please look at:
http://docs.quantumwise.com/manuals/Types/BlochState/BlochState.html#blochstate-c

The wavefunction at a given k-point and eigenvalue  can be accessed.

13

General Questions and Answers / Re: How to calculate the enthalpy of reaction in ATK?

« on: March 9, 2017, 03:41 »

I just check the atk2014 version. This version doesn't support the constraints on atoms for dynamic matrix calculations. Please use the latest version of ATK.
http://www.quantumwise.com/documents/manuals/ATK-2014/ReferenceManual/index.html/ref.dynamicalmatrixparameters.html

14

General Questions and Answers / Re: Calculating Quantum Capacitance

« on: March 8, 2017, 00:43 »

Yes, you are right.   
For example, the following papers estimated the quantum  capacitance from density of states:
J. Phys. Chem. C 2015, 119, 6464−6470
J. Phys. Chem. C 2014, 118, 4 − 15

15

General Questions and Answers / Re: How to calculate the enthalpy of reaction in ATK?

« on: March 6, 2017, 03:42 »

When you do the calculations of  force constants and dynamic matrix, you can fix the atoms of substrate. In such case, only the atoms of adsorbate molecule will be made by small displacement to calculate the vibration frequencies.
see:
http://docs.quantumwise.com/manuals/Types/DynamicalMatrix/DynamicalMatrix.html#dynamicalmatrix-c

constraints(),
The list of constrained atoms.