Show Posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.


Messages - Jess Wellendorff

Pages: 1 [2] 3 4 ... 63
16
The brillouin zone points in the repeated cell are not high-symmetry points in the same way as for the minimal cell, so the band paths are different. No surprise there. I suggest you look are the DOS instead; if you normalize properly (e.g. on a per-atom basis), and use well-converged k-point grids for DOS, the DOS should be similar no matter the repetition. 

17
In addition to using Packmol to construct the box with randomly positioned H2O molecules, a simple Force-Field based geometry optimization (minimization of both forces on atoms and stress on the supercell) might give a good starting point for DFT.

Flow:
1. use Packmol to construct the box with randomly positioned H2O molecules (perhaps aim for experimental density of molecules).
2. change cell type from UnitCell to SimpleCubic.
3. use ForceField calculator and OptimizeGeometry to relax forces and stress (take care to select a force field that actually works for water).

Example script attached. Runs in 2 mins on a laptop. Resulting water structure looks pretty OK to me, especially with the expected directional non-covalent bonding between H and O on different molecules.

If you believe finite-temperate molecular dynamics is needed to get better or more realistic water structure, that is also an option, but less simple than flow given above.

18
Dear Jie Jiang,
we are looking into this. Get back to you soon.
Regards,
Jess

19
General Questions and Answers / Re: SOGGA.PBE
« on: January 2, 2020, 09:57 »
1) spin-orbit coupling breaks time-reversal symmetry, so k-point reduction from time-reversal must be switched off.
2) looks like with LCAO uou can also use OMX for SOGGA. Not sure if OMX is available for PlaneWave also.

We have no plane-wave tutorial specifically for Q-2019.12 release. But we do have this one: https://docs.quantumatk.com/tutorials/pw_intro/pw_intro.html

20
General Questions and Answers / Re: Measuring lattice constant
« on: January 2, 2020, 09:52 »
It will be printed in the log file. But most easy to select the optimizwed BulkConfiguration on the Lab Floor, and use Text Representation widget to see info on lattice constants, atomic positions, etc.

21
Dear Bertha,
the fuzz factor simply determines how close two atoms must be in order for the GUI to draw a bond between them, so it is safe to tune it a bit according to chemical intuition. As for the k-point sampling along direction A, the equilibrium lattice constant should converge to a constant value as k_A increases, so you should probably check that.

22
The DFT-LCAO engine does not support partil charge analysis. This is also clear if you search for PertialCharges in the online documentation:
https://docs.quantumatk.com/manual/Types/PartialCharges/PartialCharges.html
"... an attached classical calculator"

So you need to define a ForceField calculator and attach it to the bulk configuration.

24
General Questions and Answers / Re: Band structure Energy
« on: May 20, 2019, 14:16 »
It is very simple: The band structure plot is Energy vs. k-point position in the BZ, so the x-axis is a route in the BZ. The y-axis is Energy referenced to the Fermi level, so -2 eV means the energy is 2 eV below the Fermi level. Just a standard band structure plot.

26
General Questions and Answers / Re: Bulk MoS2 band gap
« on: January 21, 2019, 10:27 »
0.89 eV is what you would expect for bulk MoS2 if using a LDA or GGA density functional, there are many examples in literature, see e.g. https://file.scirp.org/pdf/Graphene_2014102014384642.pdf

27
Yes, non-STD means in this case non-interacting.
No, improved screening towards the electrodes are accomplished by (1) increasing the central region size and/or (2) doping the device.
The STD method relies in the dynamical matrix, so the simulation engine must provide reliable forces. Tight-binding models are in general not suited for this. On the other hand, classical force fields, as implement in the ATK-ForceField calculator, are idel for the purpose (i.e. for computing the dynamical matrix), provided that a good force field exists for the system of interest. I believe this is the case for a graphene-based device.

28
The problem has nothing to do with the QuantumATK software. It is most likely because you do not set the QE relaxation parameter "cell_dofree", which by default is 'all' in quantum espresso. See http://www.quantum-espresso.org/Doc/INPUT_PW.html. If you have questions regarding such QE settings, I suggest you post them on the QE users mailing list: users@lists.quantum-espresso.org

Your VNL-ATK 2015 did what you asked it to do (produced a simple QE input file), and it is now your own responsibility to adapt it correctly to your specific needs.

Kind regards,
Jess Wellendorff

29
The band structure of a device electrode does not change under finite-bias NEGF calculations, only the Fermi level is shifted because of the applied bias. In fact, the electrode electronic structure is calculated before the NEGF calculations starts, and is fixed for the entire NEGF calculations, also under finite bias (except for a possible shift of the Fermi level).

I suggest you use the ProjectedLocalDensityOfStates analysis to investigate the effect of biasing your device. See e.g. this tutorial:
https://docs.quantumwise.com/tutorials/inas_p-i-n_junction/inas_p-i-n_junction.html#band-gap-and-band-edges. Perhaps consider also consider computing the Voltage Drop: https://docs.quantumwise.com/tutorials/atk_transport_calculations/atk_transport_calculations.html#voltage-drop

30
Are you calculating the effective mass along the correct/expected direction? (in the M-K direction in reciprocal space)

Pages: 1 [2] 3 4 ... 63