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 - k.cerda

Pages: [1]
1
General Questions and Answers / Effective Potential unit
« on: January 10, 2022, 13:26 »
Hi

Effective potential unit is displayed as eV in ATK. Is this right or should it be volts as the classical potential unit?

Thanks

2
Hello

I want to share simulation results in a public data storage server so that the people reading my paper will be able to view simulation output files if they need to. However, how can I share the hdf5 output of QuantumATK? Is this a generic format so that the people who do not have QuantumATK will be able to view these hdf5's in any other software? Or is there any way to convert these hdf5 outputs to any other generic file format, if so how? Thanks.

3
General Questions and Answers / Re: freeEnergy()
« on: November 2, 2021, 23:42 »
Dear shuresko

Thanks a lot for the helpful answer. Best regards.

4
General Questions and Answers / freeEnergy()
« on: November 1, 2021, 23:48 »
Hi

The potential energy calculated by TotalEnergy in ATK is always negative. But freeEnergy() calculated from phonons is positive. Why is this difference? Thanks.

5
Dear shuresko

Thanks for your detailed explanations. I'll try optimization again...

6
Dear shuresko

Thanks again for your reply. I attached the script I used for the dynamic matrix calculation and the IdealGasThermoChemistry object. Looking forward for your valuable ideas. Thanks again.

7
Dear Shuresko

Thanks a lot for your reply. However the IdealGasThermo object gives the following error when applied:
Quote
+------------------------------------------------------------------------------+
| Ideal Gas Thermo Chemistry Report (at standard pressure 1 bar)               |
|   Configuration:      C398 H579                                 |
|   Total energy:      -6224.1 eV                                              |
|   Zero point energy:  345.31 eV                                              |
| ---------------------------------------------------------------------------- |
zipdir\NL\Analysis\IdealGasThermoChemistry.py:473: RuntimeWarning: invalid value encountered in log
| temperature (K)  entropy (meV/K)  free energy (eV)                           |
|     0.00          0.000000        -5878.83940665                             |
|   100.00               nan               nan                                 |
|   200.00               nan               nan                                 |
|   250.00               nan               nan                                 |
|   298.15               nan               nan                                 |
|   300.00               nan               nan                                 |
|   350.00               nan               nan                                 |
|   400.00               nan               nan                                 |
|   450.00               nan               nan                                 |
|   500.00               nan               nan                                 |
|   600.00               nan               nan                                 |
|   700.00               nan               nan                                 |
|   800.00               nan               nan                                 |
|   900.00               nan               nan                                 |
|  1000.00               nan               nan                                 |
|  1100.00               nan               nan                                 |
|  1200.00               nan               nan                                 |
|  1300.00               nan               nan                                 |
|  1400.00               nan               nan                                 |
|  1500.00               nan               nan                                 |
|  1600.00               nan               nan                                 |
|  1700.00               nan               nan                                 |
|  1800.00               nan               nan                                 |
|  1900.00               nan               nan                                 |
|  2000.00               nan               nan                                 |
|  2100.00               nan               nan                                 |
|  2200.00               nan               nan                                 |
|  2300.00               nan               nan                                 |
|  2400.00               nan               nan                                 |
|  2500.00               nan               nan                                 |
|  2600.00               nan               nan                                 |
|  2700.00               nan               nan                                 |
|  2800.00               nan               nan                                 |
|  2900.00               nan               nan                                 |
|  3000.00               nan               nan                                 |
|  3100.00               nan               nan                                 |
|  3200.00               nan               nan                                 |
|  3300.00               nan               nan                                 |
|  3400.00               nan               nan                                 |
|  3500.00               nan               nan                                 |
|  3600.00               nan               nan                                 |
|  3700.00               nan               nan                                 |
|  3800.00               nan               nan                                 |
|  3900.00               nan               nan                                 |
|  4000.00               nan               nan                                 |
|  4100.00               nan               nan                                 |
|  4200.00               nan               nan                                 |
|  4300.00               nan               nan                                 |
|  4400.00               nan               nan                                 |
|  4500.00               nan               nan                                 |
|  4600.00               nan               nan                                 |
|  4700.00               nan               nan                                 |
|  4800.00               nan               nan                                 |
|  4900.00               nan               nan                                 |
|  5000.00               nan               nan                                 |
|  5100.00               nan               nan                                 |
|  5200.00               nan               nan                                 |
|  5300.00               nan               nan                                 |
|  5400.00               nan               nan                                 |
|  5500.00               nan               nan                                 |
|  5600.00               nan               nan                                 |
|  5700.00               nan               nan                                 |
|  5800.00               nan               nan                                 |
|  5900.00               nan               nan                         

8
Anyone here?

9
General Questions and Answers / Gibbs free energy calculation
« on: October 7, 2021, 13:27 »
Hello

The total energy of a molecule is negative when calculated using TotalEnergy() function: https://docs.quantumatk.com/manual/Types/TotalEnergy/TotalEnergy.html

When I calculate the free (vibrational) energy of the same molecule it is a positive number using the freeEnergy() function: https://docs.quantumatk.com/manual/Types/PhononDensityOfStates/PhononDensityOfStates.html

I want to calculate the Gibbs free energy therefore should I add the total energy obtained from TotalEnergy() and the vibrational free energy obtained from phonon_density_of_states.freeEnergy()? E.g. if the TotalEnergy() gives -20eV and freeEnergy() gives +1eV, is it true that Gibbs free energy of the molecule is -20+1=-19eV?

Thanks a lot.

10
Thanks a lot for your reply.

11
Hello

The total energy gives the total energy  ('Total free energy' at T<> 0 K) and the ground state energy (i.e., at T=0 K). The difference is σS. Is it possible to calculate the total free energy at a broadening other than T=0, e.g. at σ1 using a correction term such as σ1S? Thanks.

12
General Questions and Answers / Re: possible bug?
« on: September 14, 2021, 13:11 »
Thanks for your reply.

13
General Questions and Answers / possible bug?
« on: August 29, 2021, 07:03 »
I've used counterpoise LCAO calculator to simulate a big structure and took over 36 hrs. I now have the hdf5 file of the converged calculation.

When I use the existing converged hdf5 file to calculate analyses such as electron density, potential and many more, it works as expected, uses the converged data to extract these quantities very fastly.

But when I try to calculate the total energy of the system, ATK doesn't use the converged data to calculate the total energy, instead it re-starts the counterpoise LCAO SCF loops again which takes lots of hours to complete.

Is this a bug or is there a setting to calculate the total energy from converged data of counterpoise corrected LCAO calculation?

Edit: The following total energy calculation triggers a restart of the counterpoise LCAO calculation:

Code
# -------------------------------------------------------------
# Analysis from File
# -------------------------------------------------------------
path = '/home/kcerda/20082021/converged.hdf5'
configuration = nlread(path, object_id='MoleculeConfiguration_0')[0]

# -------------------------------------------------------------
# Total Energy
# -------------------------------------------------------------
total_energy = TotalEnergy(configuration)
nlsave('/home/kcerda/20082021/convergedTE.hdf5', total_energy)
nlprint(total_energy)

whereas the following calculation uses the converged data and calculates the potential instantly without a restart of SCF loops:

Code
# -------------------------------------------------------------
# Analysis from File
# -------------------------------------------------------------
path = '/home/kcerda/20082021/converged.hdf5'
configuration = nlread(path, object_id='MoleculeConfiguration_0')[0]

# -------------------------------------------------------------
# Electrostatic Potential
# -------------------------------------------------------------
electrostatic_potential = ElectrostaticPotential(configuration)
nlsave('/home/kcerda/20082021/convergedPOT.hdf5', electrostatic_potential)

14
Hi,

I'm just learning to use QuantumATK for my thesis work. I need to simulate organic molecules which contain 5500-15.000 atoms. Therefore DFT does not seem to be feasible due to the large number of atoms. I turned to Semi-Empirical method and it also does not seem to be capable of converging with such number of atoms.  Today I was reading the QuantumATK manual pages and came across the following paragraph:

Quote
"Note

The Hückel parameters have been fitted for non-self-consistent calculations. To use the parameters in self-consistent calculations, the self-consistent onsite shifts must be compensated by a reverse shift of the vacuum_levels."

on page: https://docs.quantumatk.com/manual/SE.html.

Combining this expression with the fact that the no SCF checkbox is checked by default in Semi-Empirical Calculator, can we say that the results of  non-self-consistent calculations with Semi Empirical Calculator in QuantumATK are valid and accurate? In other words, are the results of no SCF Semi-Empirical calculations valid? I understand that self-consistent field is no doubt more accurate but given the very large number of atoms, can the results of non-self-consistent Semi Empirical package be considered valid for a thesis or publication work?

Thanks,
K. Cerda

Pages: [1]