Topics

1

General Questions and Answers / Mobility workflow completes Electron–Phonon Coupling but no mobility output

« on: July 6, 2026, 06:56 »

I am trying to calculate the carrier mobility of K₃BiI₉ using the Mobility (Full) workflow. Although most of the workflow appears to execute successfully, I am unable to obtain the final mobility results and would appreciate your guidance. The following intermediate calculations appear to be completed successfully:

✔ Self-consistent field (SCF)
✔ Band structure
✔ Dynamical Matrix
✔ Hamiltonian Derivatives
✔ Phonon Band Structure
✔ Electron–Phonon Coupling (at least up to lifetime calculation)

The only missing output is the final mobility result.
Is the Mobility (Full) workflow fully supported for Hybrid (HSE) + Noncollinear SOC calculations in QuantumATK X-2025.06?
Could the HDF5 file become corrupted even when the Job Manager reports that the task has finished successfully?
Is there any known issue in QuantumATK X-2025.06 that could prevent the mobility object from being written after the electron–phonon coupling calculation?
Is there any additional debug option or verbose logging that I can enable to determine where the workflow stops after the lifetime calculation?
Would you recommend performing the mobility calculation using PBE/PBE+SOC instead of HSE+SOC, or should the current setup work correctly?

2

General Questions and Answers / Absorption data in Optical Spectrum Analyzer

« on: February 26, 2026, 12:22 »

The absorption data in the optical spectrum analyzer displays a graph of absorption (cm-1) vs wavelength (nm) upto 1722nm but the exported data of the same figure provides wavelength values from 380 to 779nm only.  Also the units in the data and the visualizer are different. I have manually converted the absoption values but the data for wavelength after 779nm is not given. Why is it so and how to fix it?

3

General Questions and Answers / Transport properties for Optimized Structures

« on: August 14, 2025, 14:35 »

Hello everyone,

I’m facing an issue when trying to convert a bulk structure into a device for transport property calculations. I can successfully perform the bulk-to-device conversion for the original (unoptimized) structures. I want to apply pressure on the material and calculate the transport properties which is an available feature in the "Optimize Geometry" block. However, after geometry optimization, I’m unable to convert the optimized structure into a device.

Could you please explain why this happens and how I can convert the optimized structure into a device?

Note:     Ori.png: Original structure>Surface Cleave (Thickness=5)>Device from Bulk
Optgeom1.png: Optimized Structure>Surface Cleave
Optgeom2.png: Surface Cleave (of Optimized structure)>Device from Bulk

Additionally, I’d like to know how to calculate the transport properties of amorphous materials, and what would be an maximum number of atoms to use for such calculations.

Thanks in advance for your help!

4

Links to Resources and Publications / Paper from 2024 ("Device from bulk" inactive in GUI)

« on: August 11, 2025, 14:12 »

Hi all,

I’m trying to calculate the thermal properties of Bi₂Te₃ (space group R3m) following the methodology from a paper (https://doi.org/10.1016/j.scriptamat.2024.115991). I was able to successfully compute the density of states (DOS) for the bulk structure.

However, I’m having trouble creating a device configuration from the bulk structure. The “Device from Bulk” option in the GUI is grayed out (inactive), so I can’t define electrodes.

Could anyone explain why this option might be unavailable, and how I should proceed to set up a device model from the bulk structure?

Thanks in advance!

5

General Questions and Answers / CALCULATED CARRIER CONCENTRATION

« on: August 4, 2025, 14:06 »

Greetings,
I am currently facing an issue while attempting to simulate carrier concentration—specifically, the function "carrier_conc" is not recognized. My objective is to compute carrier concentration as a function of chemical potential.

Could you please guide me on how to correctly perform these simulations?

Here is the code I am using:


# %% BulkConfiguration_0 (1)

# Set up lattice
vector_a = [3.950731, 0.0, 0.0]*Angstrom
vector_b = [0.0, 3.950731, 0.0]*Angstrom
vector_c = [0.0, 0.0, 3.950731]*Angstrom
lattice = UnitCell(vector_a, vector_b, vector_c)

# Define elements
elements = [Titanium, Tin, Oxygen, Oxygen, Oxygen]

# Define coordinates
fractional_coordinates = [[ 0.5,  0.5,  0.5],
                          [ 0. ,  0. ,  0. ],
                          [ 0. ,  0.5,  0.5],
                          [ 0.5,  0.5,  0. ],
                          [ 0.5,  0. ,  0.5]]

# Set up configuration
bulkconfiguration_0_1 = BulkConfiguration(
    bravais_lattice=lattice,
    elements=elements,
    fractional_coordinates=fractional_coordinates
    )

bulkconfiguration_0_1_name = "bulkconfiguration_0_1"


# %% Set LCAOCalculator

# %% LCAOCalculator

k_point_sampling = KpointDensity(
    density_a=4.0 * Angstrom, density_b=4.0 * Angstrom, density_c=4.0 * Angstrom
)

numerical_accuracy_parameters = NumericalAccuracyParameters(
    k_point_sampling=k_point_sampling
)

calculator = LCAOCalculator(
    numerical_accuracy_parameters=numerical_accuracy_parameters,
    checkpoint_handler=NoCheckpointHandler,
)


# %% Set Calculator

bulkconfiguration_0_1.setCalculator(calculator)

bulkconfiguration_0_1.update()

nlsave('BulkConfiguration_0_1_results.hdf5', bulkconfiguration_0_1)


# %% DensityOfStates

density_of_states = DensityOfStates(configuration=bulkconfiguration_0_1)
carrier_conc=calculateCarrierConcentration(contribution=Electron,fermi_shift=None,temperature=None)
nlsave('BulkConfiguration_0_1_results.hdf5', density_of_states)
nlsave('BulkConfiguration_0_1_results.hdf5', carrier_conc)


Thank you in advance for your support.

6

General Questions and Answers / WYCKOFF POSITIONS FROM CRYSTAL STRUCTURE PREDICTION SCRIPTOR

« on: July 31, 2025, 09:59 »

Dear Experts,
I need your assistance regarding the output of a crystal structure prediction script. I want to create a structure for a material and have been following the crystal structure prediction tutorial. However, I noticed that the output does not include structural details such as the Wyckoff positions or lattice parameters (a, b, c, and angles α, β, γ).

How can I extract these details from the predicted structures?

7

General Questions and Answers / Boltzman Transport Calculations

« on: July 22, 2025, 12:52 »

I am currently working on the following research paper: https://www.mdpi.com/2073-4352/12/9/1317, and I am seeking guidance on using Boltzmann transport theory to calculate the Seebeck coefficient and thermal conductivity.

Specifically, I am having difficulty understanding how to compute the Seebeck and Hall coefficients as functions of chemical potential. I would greatly appreciate your insights or any relevant references that could help clarify this process.

8

General Questions and Answers / Optimized geometry configuration is not visible

« on: April 28, 2025, 14:51 »

Hi everyone, I'm facing an error after running a Molecular Dynamics simulation that involves the following steps- Heating, quenching and geometry optimization. The simulation finishes successfully. However, when I try to open the optimized geometry in the Builder, the screen appears blank and no structure is visible. But i get to see the configuration in the data preview but if we try to open it in the viewer/Builder the blank screen without configuration is viewing. What could be the possible reasons for this behavior? How can I resolve or debug this issue?

9

General Questions and Answers / Documentation example is not working

« on: April 24, 2025, 10:02 »

https://docs.quantumatk.com/_downloads/230cae9ab632421ed113e69e59fb14aa/TiO2_Matsui.py
shows error
Traceback (most recent call last):
  File "/home/qtk/quantumatk/W-2024.09/bin/../atkpython/bin/atkpython", line 8, in <module>
    sys.exit(__run_atkpython())
             ^^^^^^^^^^^^^^^^^
  File "zipdir/ATKExecutables/atkwrappers/__init__.py", line 803, in __run_atkpython
  File "Tio2.py", line 32, in <module>
    potentialSet.addParticleType(ParticleType(symbol='O', mass=15.9994*atomic_mass_unit, charge=-1.098))
                                 ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "build/atkpython/lib/python3.11/site-packages/tremolox/potentials/particletype.py", line 90, in __init__
  File "build/atkpython/lib/python3.11/site-packages/tremolox/potentials/potentialbase.py", line 335, in setParameter
  File "build/atkpython/lib/python3.11/site-packages/tremolox/potentials/potentialbase.py", line 155, in __setitem__
  File "build/atkpython/lib/python3.11/site-packages/tremolox/potentials/checks.py", line 139, in apply
TypeError: The parameter charge = -1.098 must have a unit that can be converted to e.