46
General Questions and Answers / Questions about nlread
« on: January 1, 2024, 05:17 »
Hello, and happy new year everyone.
Recently, I had a Deviceconfiguration that successfully converged with PLDOS. It took more than 10 hours
I wanted to reuse it with transmission spectrum, but since it consumed a lot of time, I was planning to redo it with 'nlread function'.
Here is my code
# DeviceConfiguration
Device_configuration = nlread('Si80_111_5e19_4040.hdf5',DeviceConfiguration_0)[0]
#Choose the first configuration in the list
Device_configuration = Device_configuration_list[0]
# %% Set Calculator
Device_configuration.setCalculator(calculator)
Device_configuration.update()
nlsave('Si80_111_5e19_5050trans.hdf5', Device_configuration)
# %% TransmissionSpectrum
kpoints = MonkhorstPackGrid(
na=50,
nb=50,
nc=1,
symmetries=[
([[ 1., 0., 0.],
[ 0., 1., 0.],
[ 0., 0., 1.]], [ 0., 0., 0.]),
([[-1., 0., 0.],
[ 0.,-1., 0.],
[ 0., 0.,-1.]], [ 0., 0., 0.]),
],
force_timereversal=True,
k_point_shift=None,
shift_to_gamma=[True, True, True]
)
self_energy_calculator = RecursionSelfEnergy(
storage_strategy=NoStorage()
)
device_transmission_method = GreensFunction(
processes_per_contour_point=1,
)
transmission_spectrum = TransmissionSpectrum(
configuration=Device_configuration,
energies=numpy.linspace(0.0, 0.0, 1)*eV,
kpoints=kpoints,
kpoints_weights=None,
self_energy_calculator=self_energy_calculator,
energy_zero_parameter=AverageFermiLevel,
infinitesimal=1e-06*eV,
device_transmission_method=device_transmission_method,
enforce_zero_in_band_gap=True
)
nlsave('Si80_111_5e19_5050trans.hdf5', transmission_spectrum)
However when I ran the calculation nothing appeared on the log file, and the calculation did not run.
Could somebody please help me with how to use 'nlread' properly?
Thank you.
Recently, I had a Deviceconfiguration that successfully converged with PLDOS. It took more than 10 hours
I wanted to reuse it with transmission spectrum, but since it consumed a lot of time, I was planning to redo it with 'nlread function'.
Here is my code
# DeviceConfiguration
Device_configuration = nlread('Si80_111_5e19_4040.hdf5',DeviceConfiguration_0)[0]
#Choose the first configuration in the list
Device_configuration = Device_configuration_list[0]
# %% Set Calculator
Device_configuration.setCalculator(calculator)
Device_configuration.update()
nlsave('Si80_111_5e19_5050trans.hdf5', Device_configuration)
# %% TransmissionSpectrum
kpoints = MonkhorstPackGrid(
na=50,
nb=50,
nc=1,
symmetries=[
([[ 1., 0., 0.],
[ 0., 1., 0.],
[ 0., 0., 1.]], [ 0., 0., 0.]),
([[-1., 0., 0.],
[ 0.,-1., 0.],
[ 0., 0.,-1.]], [ 0., 0., 0.]),
],
force_timereversal=True,
k_point_shift=None,
shift_to_gamma=[True, True, True]
)
self_energy_calculator = RecursionSelfEnergy(
storage_strategy=NoStorage()
)
device_transmission_method = GreensFunction(
processes_per_contour_point=1,
)
transmission_spectrum = TransmissionSpectrum(
configuration=Device_configuration,
energies=numpy.linspace(0.0, 0.0, 1)*eV,
kpoints=kpoints,
kpoints_weights=None,
self_energy_calculator=self_energy_calculator,
energy_zero_parameter=AverageFermiLevel,
infinitesimal=1e-06*eV,
device_transmission_method=device_transmission_method,
enforce_zero_in_band_gap=True
)
nlsave('Si80_111_5e19_5050trans.hdf5', transmission_spectrum)
However when I ran the calculation nothing appeared on the log file, and the calculation did not run.
Could somebody please help me with how to use 'nlread' properly?
Thank you.
