5. Infrared and Raman spectra using DFT+U+V#
Using different functionals might be needed in different real case scenario. Here, we demonstrate how the workflows in the package can be employed to predict the infrared and Raman spectra using both DFT and DFT+U+V.
We will make use of LiCoO2 for the demonstration.
Let’s get started!
Show cell content
from local_module import load_temp_profile
# If you download this file, you can run it with your own profile.
# Put these lines instead:
# from aiida import load_profile
# load_profile()
data = load_temp_profile(
name="functionals-tutorial",
add_computer=True,
add_pw_code=True,
add_sssp=True,
)
5.1. Creating the structure#
Let’s define the LiCoO2 structure.
from aiida.orm import StructureData
a, b, c, d = 1.4060463552647, 0.81178124180108, 4.6012019181836, 1.6235624832021
cell = [[a, -b, c], [0.0, d, c], [-a, -b, c]]
sites = [
['Co', 'Co', (0, 0, 0)],
['O', 'O', (0, 0, 3.6020728736387)],
['O', 'O', (0, 0, 10.201532881212)],
['Li', 'Li', (0, 0, 6.9018028772754)],
]
structure = StructureData(cell=cell)
for site in sites:
structure.append_atom(position=site[2], symbols=site[0], name=site[1])
structure.store()
<StructureData: uuid: a489219a-22da-4135-879d-9483bd1380e3 (pk: 87)>
5.2. Spectra using PBEsol#
Let’s first import the WorkChain and run it using ‘simple’ PBEsol!
As usual, we use the get_builder_from_protocol to get a prefilled builder with all inputs. These inputs should be considered not as converged parameters, but as a good starting point. You may also need to tweak some parameters, e.g. add magnetization etc., depending on your case.
from aiida.plugins import WorkflowFactory
from aiida.engine import run_get_node
IRamanSpectraWorkChain = WorkflowFactory("vibroscopy.spectra.iraman")
# To make the calculations run faster
scf_overrides = {'scf':{'pw':{'parameters':{'SYSTEM':{'ecutwfc':30, 'ecutrho':30*8}}}}}
builder = IRamanSpectraWorkChain.get_builder_from_protocol(
code=data.pw_code,
structure=structure,
protocol="fast",
overrides={'dielectric':scf_overrides, 'phonon':scf_overrides}
)
results, calc = run_get_node(builder)
Show cell output
08/25/2023 06:55:26 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [106|IRamanSpectraWorkChain|run_spectra]: submitting `HarmonicWorkChain` <PK=108>
08/25/2023 06:55:26 PM <799519> aiida.engine.processes.functions: [WARNING] function `generate_preprocess_data` has invalid type hints: unsupported operand type(s) for |: 'AbstractNodeMeta' and 'NoneType'
08/25/2023 06:55:26 PM <799519> aiida.engine.processes.functions: [WARNING] function `generate_phonopy_data` has invalid type hints: unsupported operand type(s) for |: 'AbstractNodeMeta' and 'NoneType'
08/25/2023 06:55:27 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [108|HarmonicWorkChain|run_phonon]: submitting `PhononWorkChain` <PK=114>
08/25/2023 06:55:27 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [108|HarmonicWorkChain|run_dielectric]: submitting `DielectricWorkChain` <PK=118>
08/25/2023 06:55:28 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_base_scf]: launching base scf PwBaseWorkChain<143>
08/25/2023 06:55:29 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [143|PwBaseWorkChain|run_process]: launching PwCalculation<147> iteration #1
08/25/2023 06:55:29 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [114|PhononWorkChain|run_base_supercell]: launching base supercell scf PwBaseWorkChain<152>
08/25/2023 06:55:29 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [152|PwBaseWorkChain|run_process]: launching PwCalculation<155> iteration #1
08/25/2023 06:55:47 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [152|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 06:55:47 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [152|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 06:55:48 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [143|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 06:55:48 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [143|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 06:55:49 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [114|PhononWorkChain|run_forces]: submitting `PwBaseWorkChain` <PK=172> with supercell n.o 1
08/25/2023 06:55:50 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [114|PhononWorkChain|run_forces]: submitting `PwBaseWorkChain` <PK=174> with supercell n.o 2
08/25/2023 06:55:51 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [114|PhononWorkChain|run_forces]: submitting `PwBaseWorkChain` <PK=176> with supercell n.o 3
08/25/2023 06:55:52 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [114|PhononWorkChain|run_forces]: submitting `PwBaseWorkChain` <PK=178> with supercell n.o 4
08/25/2023 06:55:54 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [172|PwBaseWorkChain|run_process]: launching PwCalculation<181> iteration #1
08/25/2023 06:55:54 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [174|PwBaseWorkChain|run_process]: launching PwCalculation<184> iteration #1
08/25/2023 06:55:54 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [176|PwBaseWorkChain|run_process]: launching PwCalculation<187> iteration #1
08/25/2023 06:55:54 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [178|PwBaseWorkChain|run_process]: launching PwCalculation<190> iteration #1
08/25/2023 06:55:55 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_nscf]: launching base scf PwBaseWorkChain<193>
08/25/2023 06:55:55 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [193|PwBaseWorkChain|run_process]: launching PwCalculation<196> iteration #1
08/25/2023 06:56:52 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [193|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 06:56:52 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [193|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 06:56:55 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_null_field_scfs]: launching PwBaseWorkChain<216> with null electric field 0
08/25/2023 06:56:55 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_null_field_scfs]: launching PwBaseWorkChain<217> with null electric field 1
08/25/2023 06:56:55 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_null_field_scfs]: launching PwBaseWorkChain<218> with null electric field 2
08/25/2023 06:56:55 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_null_field_scfs]: launching PwBaseWorkChain<219> with null electric field 3
08/25/2023 06:56:55 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_null_field_scfs]: launching PwBaseWorkChain<220> with null electric field 4
08/25/2023 06:56:57 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [216|PwBaseWorkChain|run_process]: launching PwCalculation<229> iteration #1
08/25/2023 06:56:57 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [217|PwBaseWorkChain|run_process]: launching PwCalculation<232> iteration #1
08/25/2023 06:56:57 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [218|PwBaseWorkChain|run_process]: launching PwCalculation<235> iteration #1
08/25/2023 06:56:57 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [219|PwBaseWorkChain|run_process]: launching PwCalculation<238> iteration #1
08/25/2023 06:56:58 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [220|PwBaseWorkChain|run_process]: launching PwCalculation<241> iteration #1
08/25/2023 06:56:58 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [178|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 06:56:58 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [178|PwBaseWorkChain|on_terminated]: cleaned remote folders of calculations: 190
08/25/2023 06:56:58 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [174|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 06:56:59 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [174|PwBaseWorkChain|on_terminated]: cleaned remote folders of calculations: 184
08/25/2023 06:57:13 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [176|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 06:57:14 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [176|PwBaseWorkChain|on_terminated]: cleaned remote folders of calculations: 187
08/25/2023 06:57:15 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [172|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 06:57:15 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [172|PwBaseWorkChain|on_terminated]: cleaned remote folders of calculations: 181
08/25/2023 06:57:25 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [114|PhononWorkChain|on_terminated]: cleaned remote folders of calculations: 155 181 184 187 190
08/25/2023 06:59:52 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [219|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 06:59:52 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [219|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 06:59:52 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [217|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 06:59:52 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [217|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 06:59:59 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [220|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 06:59:59 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [220|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:00:02 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [216|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:00:02 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [216|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:00:15 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [218|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:00:16 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [218|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:00:17 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<280> with electric field index 0 and sign 1.0 iteration #0
08/25/2023 07:00:18 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<283> with electric field index 1 and sign 1.0 iteration #0
08/25/2023 07:00:19 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<286> with electric field index 2 and sign 1.0 iteration #0
08/25/2023 07:00:20 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<289> with electric field index 3 and sign 1.0 iteration #0
08/25/2023 07:00:21 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<292> with electric field index 4 and sign 1.0 iteration #0
08/25/2023 07:00:22 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<295> with electric field index 5 and sign 1.0 iteration #0
08/25/2023 07:00:24 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [280|PwBaseWorkChain|run_process]: launching PwCalculation<298> iteration #1
08/25/2023 07:00:24 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [283|PwBaseWorkChain|run_process]: launching PwCalculation<301> iteration #1
08/25/2023 07:00:24 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [286|PwBaseWorkChain|run_process]: launching PwCalculation<304> iteration #1
08/25/2023 07:00:24 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [289|PwBaseWorkChain|run_process]: launching PwCalculation<307> iteration #1
08/25/2023 07:00:25 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [292|PwBaseWorkChain|run_process]: launching PwCalculation<310> iteration #1
08/25/2023 07:00:25 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [295|PwBaseWorkChain|run_process]: launching PwCalculation<313> iteration #1
08/25/2023 07:06:49 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [283|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:06:50 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [283|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:07:47 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [280|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:07:48 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [280|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:07:56 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [295|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:07:56 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [295|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:08:12 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [292|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:08:12 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [292|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:08:16 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [289|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:08:16 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [289|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:08:59 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [286|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:08:59 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [286|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:09:01 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<346> with electric field index 0 and sign 1.0 iteration #1
08/25/2023 07:09:02 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<349> with electric field index 1 and sign 1.0 iteration #1
08/25/2023 07:09:03 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<352> with electric field index 2 and sign 1.0 iteration #1
08/25/2023 07:09:04 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<355> with electric field index 3 and sign 1.0 iteration #1
08/25/2023 07:09:05 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<358> with electric field index 4 and sign 1.0 iteration #1
08/25/2023 07:09:07 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<361> with electric field index 5 and sign 1.0 iteration #1
08/25/2023 07:09:08 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [346|PwBaseWorkChain|run_process]: launching PwCalculation<364> iteration #1
08/25/2023 07:09:08 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [349|PwBaseWorkChain|run_process]: launching PwCalculation<367> iteration #1
08/25/2023 07:09:08 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [352|PwBaseWorkChain|run_process]: launching PwCalculation<370> iteration #1
08/25/2023 07:09:09 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [355|PwBaseWorkChain|run_process]: launching PwCalculation<373> iteration #1
08/25/2023 07:09:09 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [358|PwBaseWorkChain|run_process]: launching PwCalculation<376> iteration #1
08/25/2023 07:09:09 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [361|PwBaseWorkChain|run_process]: launching PwCalculation<379> iteration #1
08/25/2023 07:16:20 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [349|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:16:20 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [349|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:16:33 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [361|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:16:33 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [361|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:17:42 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [346|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:17:42 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [346|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:17:42 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [358|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:17:42 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [358|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:17:45 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [355|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:17:45 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [355|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:18:43 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [352|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:18:43 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [352|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:18:46 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|run_numerical_derivatives]: launching NumericalDerivativesWorkChain<425> for computing numerical derivatives.
/opt/conda/envs/aiida-vibroscopy-demo/lib/python3.9/site-packages/phonopy/structure/symmetry.py:618: UserWarning: Symmetry of Born effective charge is largely broken. The difference is:
[[[ 1.74716261e+00 1.36724849e-02 6.43411053e-02]
[ 7.23837435e-03 -1.74736367e+00 1.28601784e+00]
[ 5.14728843e-02 -3.13662889e-02 -8.73967565e-13]]
[[-4.89696131e-01 -7.23837435e-03 -1.64874082e-02]
[-2.41279145e-03 4.89495065e-01 -5.84699794e-01]
[-1.84980678e-02 1.04554296e-02 3.41060513e-13]]
[[-4.89696131e-01 -7.23837435e-03 -1.64874082e-02]
[-2.41279145e-03 4.89495065e-01 -5.84699794e-01]
[-1.84980678e-02 1.04554296e-02 3.41060513e-13]]
[[-7.68574610e-01 8.04263817e-04 -3.13662889e-02]
[-3.21705527e-03 7.68373544e-01 -1.17422517e-01]
[-1.52810125e-02 1.04554296e-02 3.41060513e-13]]]
warnings.warn("\n".join(lines))
/opt/conda/envs/aiida-vibroscopy-demo/lib/python3.9/site-packages/phonopy/structure/symmetry.py:618: UserWarning: Symmetry of Born effective charge is largely broken. The difference is:
[[[ 3.81040090e+00 3.05620250e-02 1.25465155e-01]
[-3.21705527e-03 -3.81080303e+00 7.07752157e-01]
[ 1.41550432e-01 -1.09379879e-01 4.02131905e-04]]
[[-1.14949406e+00 -1.60852763e-02 -3.13662889e-02]
[ 4.02131908e-03 1.14909193e+00 -3.47441968e-01]
[-4.82558290e-02 4.26259823e-02 4.02131910e-04]]
[[-1.14949406e+00 -1.60852763e-02 -3.13662889e-02]
[ 4.02131908e-03 1.14909193e+00 -3.47441968e-01]
[-4.82558290e-02 4.26259823e-02 4.02131910e-04]]
[[-1.51141278e+00 3.21705527e-03 -6.27325777e-02]
[-6.43411053e-03 1.51101065e+00 -1.28682206e-02]
[-4.34302461e-02 2.57364421e-02 4.02131909e-04]]]
warnings.warn("\n".join(lines))
/opt/conda/envs/aiida-vibroscopy-demo/lib/python3.9/site-packages/phonopy/structure/symmetry.py:618: UserWarning: Symmetry of Born effective charge is largely broken. The difference is:
[[[ 4.49814699e+00 3.61918718e-02 1.45839839e-01]
[-6.70219847e-03 -4.49861615e+00 5.14996929e-01]
[ 1.71576281e-01 -1.35384409e-01 5.36175874e-04]]
[[-1.36942670e+00 -1.90342437e-02 -3.63259157e-02]
[ 6.16602260e-03 1.36895755e+00 -2.68356026e-01]
[-5.81750827e-02 5.33494998e-02 5.36175879e-04]]
[[-1.36942670e+00 -1.90342437e-02 -3.63259157e-02]
[ 6.16602260e-03 1.36895755e+00 -2.68356026e-01]
[-5.81750827e-02 5.33494998e-02 5.36175879e-04]]
[[-1.75902550e+00 4.02131908e-03 -7.31880073e-02]
[-7.50646229e-03 1.75855635e+00 2.19832115e-02]
[-5.28133240e-02 3.08301130e-02 5.36175879e-04]]]
warnings.warn("\n".join(lines))
/opt/conda/envs/aiida-vibroscopy-demo/lib/python3.9/site-packages/aiida_vibroscopy/calculations/symmetry.py:163: UserWarning: Symmetry of dChi/dr tensors is largely broken.
The max difference is:
733.3581532308866
warnings.warn('\n'.join(lines))
/opt/conda/envs/aiida-vibroscopy-demo/lib/python3.9/site-packages/aiida_vibroscopy/calculations/symmetry.py:163: UserWarning: Symmetry of dChi/dr tensors is largely broken.
The max difference is:
1997.5518823010948
warnings.warn('\n'.join(lines))
/opt/conda/envs/aiida-vibroscopy-demo/lib/python3.9/site-packages/aiida_vibroscopy/calculations/symmetry.py:163: UserWarning: Symmetry of dChi/dr tensors is largely broken.
The max difference is:
2747.682091653869
warnings.warn('\n'.join(lines))
08/25/2023 07:18:56 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [118|DielectricWorkChain|on_terminated]: cleaned remote folders of calculations: 147 196 229 232 235 238 241 298 301 304 307 310 313 364 367 370 373 376 379
08/25/2023 07:18:59 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [108|HarmonicWorkChain|on_terminated]: cleaned remote folders of calculations: 155 181 184 187 190 147 196 229 232 235 238 241 298 301 304 307 310 313 364 367 370 373 376 379
5.3. Powder Raman spectra using PBEsol#
Let’s plot the powder Raman spectra computed using PBEsol.
from aiida_vibroscopy.utils.plotting import get_spectra_plot
vibro = calc.outputs.vibrational_data.numerical_accuracy_4
polarized_intensities, unpolarized_intensities, frequencies_pbesol, labels = vibro.run_powder_raman_intensities(frequency_laser=532, temperature=300)
total_intensities_pbesol = polarized_intensities + unpolarized_intensities
plt = get_spectra_plot(frequencies_pbesol, total_intensities_pbesol)
plt.show()
5.4. Spectra using PBEsol+U+V#
5.4.1. Defining the HubbardStructureData#
The Hubbard input should be specified alongside with the structural data, in a unique data called HubbardStructureData. Refer to the aiida-quantumespresso for more insights. What you need to know is that you can run a full Hubbard corrected calculations by just specifying the HubbardStructureData as an input for the structure; this will work in any workchain of the package!
Attention
If you want to compute self-consistently the Hubbard parameters in automated fashion, please have a look at our aiida-quantumespresso-hp plugin, fully automated and parallelisable! The SelfConsistentHubbardWorkChain will deliver for you a fully relaxed structure with fully self-consistent ab-initio Hubbard parameters!
Note
For this demonstration we will not compute the Hubbard parameters ab-initio. These ones that we will use were computed beforehand in a separate job run in order to demonstrate quickly the properties we are interested in.
from aiida_quantumespresso.common.hubbard import Hubbard
from aiida_quantumespresso.data.hubbard_structure import HubbardStructureData
hubbard_list = [
(0, '3d', 0, '3d', 7.2362, (0, 0, 0), 'V'),
(0, '3d', 2, '2p', 0.2999, (-1, 0, -1), 'V'),
(0, '3d', 1, '2p', 0.2999, (0, 0, -1), 'V'),
(0, '3d', 1, '2p', 0.2999, (-1, 0, 0), 'V'),
(0, '3d', 2, '2p', 0.2999, (0, -1, -1), 'V'),
(0, '3d', 2, '2p', 0.2999, (-1, -1, 0), 'V'),
(0, '3d', 1, '2p', 0.2999, (0, -1, 0), 'V')
]
a, b, c, d = 1.4023265572925, 0.80963361535606, 4.6735772081287, 1.6192672305121
cell = [[a, -b, c], [0.0, d, c], [-a, -b, c]]
sites = [
['Co', 'Co', (0, 0, 0)],
['O', 'O', (0, 0, 3.6492980821051)],
['O', 'O', (0, 0, 10.371433542381)],
['Li', 'Li', (0, 0, 7.010365812143)],
]
hubbard_structure = HubbardStructureData(cell=cell, sites=sites)
hubbard_structure.hubbard = Hubbard.from_list(hubbard_list)
hubbard_structure.store()
<HubbardStructureData: uuid: 87553170-f713-4065-a8b0-effccfdc70fb (pk: 455)>
5.4.2. Running the IRamanSpectraWorkChain with the HubbardStructureData#
Let’s run the workchain for IR and Raman using the just defined structure with Hubbard parameters.
builder = IRamanSpectraWorkChain.get_builder_from_protocol(
code=data.pw_code,
structure=hubbard_structure,
protocol="fast",
overrides={'dielectric':scf_overrides, 'phonon':scf_overrides}
)
results_hubabrd, calc_hubbard = run_get_node(builder)
Show cell output
08/25/2023 07:19:02 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [474|IRamanSpectraWorkChain|run_spectra]: submitting `HarmonicWorkChain` <PK=476>
08/25/2023 07:19:02 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [476|HarmonicWorkChain|run_phonon]: submitting `PhononWorkChain` <PK=482>
08/25/2023 07:19:02 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [476|HarmonicWorkChain|run_dielectric]: submitting `DielectricWorkChain` <PK=486>
08/25/2023 07:19:04 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_base_scf]: launching base scf PwBaseWorkChain<511>
08/25/2023 07:19:04 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [511|PwBaseWorkChain|run_process]: launching PwCalculation<515> iteration #1
08/25/2023 07:19:04 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [482|PhononWorkChain|run_base_supercell]: launching base supercell scf PwBaseWorkChain<522>
08/25/2023 07:19:05 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [522|PwBaseWorkChain|run_process]: launching PwCalculation<525> iteration #1
08/25/2023 07:19:21 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [522|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:19:21 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [522|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:19:23 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [482|PhononWorkChain|run_forces]: submitting `PwBaseWorkChain` <PK=538> with supercell n.o 1
08/25/2023 07:19:24 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [482|PhononWorkChain|run_forces]: submitting `PwBaseWorkChain` <PK=540> with supercell n.o 2
08/25/2023 07:19:25 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [482|PhononWorkChain|run_forces]: submitting `PwBaseWorkChain` <PK=542> with supercell n.o 3
08/25/2023 07:19:27 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [482|PhononWorkChain|run_forces]: submitting `PwBaseWorkChain` <PK=544> with supercell n.o 4
08/25/2023 07:19:28 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [538|PwBaseWorkChain|run_process]: launching PwCalculation<547> iteration #1
08/25/2023 07:19:28 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [540|PwBaseWorkChain|run_process]: launching PwCalculation<550> iteration #1
08/25/2023 07:19:28 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [542|PwBaseWorkChain|run_process]: launching PwCalculation<553> iteration #1
08/25/2023 07:19:29 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [544|PwBaseWorkChain|run_process]: launching PwCalculation<556> iteration #1
08/25/2023 07:19:34 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [511|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:19:34 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [511|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:19:38 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_nscf]: launching base scf PwBaseWorkChain<567>
08/25/2023 07:19:39 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [567|PwBaseWorkChain|run_process]: launching PwCalculation<570> iteration #1
08/25/2023 07:20:41 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [544|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:20:42 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [544|PwBaseWorkChain|on_terminated]: cleaned remote folders of calculations: 556
08/25/2023 07:20:42 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [567|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:20:42 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [567|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:20:44 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_null_field_scfs]: launching PwBaseWorkChain<590> with null electric field 0
08/25/2023 07:20:44 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_null_field_scfs]: launching PwBaseWorkChain<591> with null electric field 1
08/25/2023 07:20:44 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_null_field_scfs]: launching PwBaseWorkChain<592> with null electric field 2
08/25/2023 07:20:44 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_null_field_scfs]: launching PwBaseWorkChain<593> with null electric field 3
08/25/2023 07:20:45 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_null_field_scfs]: launching PwBaseWorkChain<594> with null electric field 4
08/25/2023 07:20:46 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [590|PwBaseWorkChain|run_process]: launching PwCalculation<603> iteration #1
08/25/2023 07:20:46 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [591|PwBaseWorkChain|run_process]: launching PwCalculation<606> iteration #1
08/25/2023 07:20:47 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [592|PwBaseWorkChain|run_process]: launching PwCalculation<609> iteration #1
08/25/2023 07:20:47 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [593|PwBaseWorkChain|run_process]: launching PwCalculation<612> iteration #1
08/25/2023 07:20:47 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [594|PwBaseWorkChain|run_process]: launching PwCalculation<615> iteration #1
08/25/2023 07:20:47 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [538|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:20:48 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [538|PwBaseWorkChain|on_terminated]: cleaned remote folders of calculations: 547
08/25/2023 07:20:48 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [540|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:20:48 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [540|PwBaseWorkChain|on_terminated]: cleaned remote folders of calculations: 550
08/25/2023 07:20:55 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [542|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:20:55 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [542|PwBaseWorkChain|on_terminated]: cleaned remote folders of calculations: 553
08/25/2023 07:21:03 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [482|PhononWorkChain|on_terminated]: cleaned remote folders of calculations: 525 547 550 553 556
08/25/2023 07:24:40 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [591|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:24:40 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [591|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:24:45 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [593|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:24:45 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [593|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:24:46 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [594|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:24:46 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [594|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:24:54 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [592|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:24:55 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [592|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:25:02 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [590|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:25:03 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [590|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:25:03 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<650> with electric field index 0 and sign 1.0 iteration #0
08/25/2023 07:25:04 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<653> with electric field index 1 and sign 1.0 iteration #0
08/25/2023 07:25:05 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<656> with electric field index 2 and sign 1.0 iteration #0
08/25/2023 07:25:06 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<659> with electric field index 3 and sign 1.0 iteration #0
08/25/2023 07:25:07 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<662> with electric field index 4 and sign 1.0 iteration #0
08/25/2023 07:25:09 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<665> with electric field index 5 and sign 1.0 iteration #0
08/25/2023 07:25:10 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [650|PwBaseWorkChain|run_process]: launching PwCalculation<668> iteration #1
08/25/2023 07:25:10 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [653|PwBaseWorkChain|run_process]: launching PwCalculation<671> iteration #1
08/25/2023 07:25:10 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [656|PwBaseWorkChain|run_process]: launching PwCalculation<674> iteration #1
08/25/2023 07:25:11 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [659|PwBaseWorkChain|run_process]: launching PwCalculation<677> iteration #1
08/25/2023 07:25:11 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [662|PwBaseWorkChain|run_process]: launching PwCalculation<680> iteration #1
08/25/2023 07:25:11 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [665|PwBaseWorkChain|run_process]: launching PwCalculation<683> iteration #1
08/25/2023 07:32:47 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [653|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:32:47 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [653|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:33:32 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [665|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:33:32 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [665|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:33:36 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [662|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:33:37 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [662|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:33:46 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [650|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:33:46 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [650|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:33:52 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [659|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:33:52 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [659|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:34:57 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [656|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:34:58 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [656|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:34:59 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<716> with electric field index 0 and sign 1.0 iteration #1
08/25/2023 07:35:00 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<719> with electric field index 1 and sign 1.0 iteration #1
08/25/2023 07:35:01 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<722> with electric field index 2 and sign 1.0 iteration #1
08/25/2023 07:35:02 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<725> with electric field index 3 and sign 1.0 iteration #1
08/25/2023 07:35:03 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<728> with electric field index 4 and sign 1.0 iteration #1
08/25/2023 07:35:05 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_electric_field_scfs]: launching PwBaseWorkChain<731> with electric field index 5 and sign 1.0 iteration #1
08/25/2023 07:35:06 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [716|PwBaseWorkChain|run_process]: launching PwCalculation<734> iteration #1
08/25/2023 07:35:06 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [719|PwBaseWorkChain|run_process]: launching PwCalculation<737> iteration #1
08/25/2023 07:35:06 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [722|PwBaseWorkChain|run_process]: launching PwCalculation<740> iteration #1
08/25/2023 07:35:07 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [725|PwBaseWorkChain|run_process]: launching PwCalculation<743> iteration #1
08/25/2023 07:35:07 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [728|PwBaseWorkChain|run_process]: launching PwCalculation<746> iteration #1
08/25/2023 07:35:07 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [731|PwBaseWorkChain|run_process]: launching PwCalculation<749> iteration #1
08/25/2023 07:44:23 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [719|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:44:24 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [719|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:44:36 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [731|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:44:37 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [731|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:44:48 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [716|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:44:49 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [716|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:45:11 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [728|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:45:11 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [728|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:45:16 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [725|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:45:16 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [725|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:46:12 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [722|PwBaseWorkChain|results]: work chain completed after 1 iterations
08/25/2023 07:46:12 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [722|PwBaseWorkChain|on_terminated]: remote folders will not be cleaned
08/25/2023 07:46:16 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|run_numerical_derivatives]: launching NumericalDerivativesWorkChain<795> for computing numerical derivatives.
/opt/conda/envs/aiida-vibroscopy-demo/lib/python3.9/site-packages/phonopy/structure/symmetry.py:618: UserWarning: Symmetry of Born effective charge is largely broken. The difference is:
[[[-1.77461208e-01 5.17945883e-04 -6.21535060e-03]
[-2.33075647e-03 1.77331722e-01 -5.40994476e-01]
[-2.04588624e-02 5.17945883e-03 6.47432355e-05]]
[[ 8.78565704e-02 -1.29486471e-04 2.20127000e-03]
[ 3.88459412e-04 -8.79860569e-02 1.55254279e-01]
[ 7.12175589e-03 -1.55383765e-03 6.47432353e-05]]
[[ 8.78565704e-02 -1.29486471e-04 2.20127000e-03]
[ 3.88459412e-04 -8.79860569e-02 1.55254279e-01]
[ 7.12175589e-03 -1.55383765e-03 6.47432353e-05]]
[[ 1.48909441e-03 -5.17945883e-04 1.55383765e-03]
[ 1.29486471e-03 -1.61858088e-03 2.30485918e-01]
[ 6.73329648e-03 -1.81281059e-03 6.47432354e-05]]]
warnings.warn("\n".join(lines))
/opt/conda/envs/aiida-vibroscopy-demo/lib/python3.9/site-packages/phonopy/structure/symmetry.py:618: UserWarning: Symmetry of Born effective charge is largely broken. The difference is:
[[[ 2.35082688e-01 1.55383765e-03 -1.24307012e-02]
[-5.17945883e-03 -2.35212174e-01 -1.33630039e-01]
[-3.62562118e-02 1.50204306e-02 1.29486471e-04]]
[[-4.33132245e-02 -2.58972941e-04 4.14356706e-03]
[ 1.29486471e-03 4.31837380e-02 -6.99226908e-03]
[ 1.13948094e-02 -5.43843177e-03 1.29486471e-04]]
[[-4.33132245e-02 -2.58972941e-04 4.14356706e-03]
[ 1.29486471e-03 4.31837380e-02 -6.99226908e-03]
[ 1.13948094e-02 -5.43843177e-03 1.29486471e-04]]
[[-1.48974185e-01 -1.03589177e-03 3.62562118e-03]
[ 2.58972941e-03 1.48844698e-01 1.48132523e-01]
[ 1.45024847e-02 -4.14356706e-03 1.29486471e-04]]]
warnings.warn("\n".join(lines))
/opt/conda/envs/aiida-vibroscopy-demo/lib/python3.9/site-packages/phonopy/structure/symmetry.py:618: UserWarning: Symmetry of Born effective charge is largely broken. The difference is:
[[[ 3.72597320e-01 1.89913490e-03 -1.45024847e-02]
[-6.12902628e-03 -3.72726806e-01 2.15810692e-03]
[-4.15219950e-02 1.83007545e-02 1.51067549e-04]]
[[-8.70364894e-02 -3.02135098e-04 4.79099942e-03]
[ 1.59699981e-03 8.69070029e-02 -6.10744517e-02]
[ 1.28191606e-02 -6.73329648e-03 1.51067549e-04]]
[[-8.70364894e-02 -3.02135098e-04 4.79099942e-03]
[ 1.59699981e-03 8.69070029e-02 -6.10744517e-02]
[ 1.28191606e-02 -6.73329648e-03 1.51067549e-04]]
[[-1.99128611e-01 -1.20854039e-03 4.31621569e-03]
[ 3.02135098e-03 1.98999124e-01 1.20681391e-01]
[ 1.70922141e-02 -4.92048589e-03 1.51067549e-04]]]
warnings.warn("\n".join(lines))
/opt/conda/envs/aiida-vibroscopy-demo/lib/python3.9/site-packages/aiida_vibroscopy/calculations/symmetry.py:163: UserWarning: Symmetry of dChi/dr tensors is largely broken.
The max difference is:
64.69056804204402
warnings.warn('\n'.join(lines))
/opt/conda/envs/aiida-vibroscopy-demo/lib/python3.9/site-packages/aiida_vibroscopy/calculations/symmetry.py:163: UserWarning: Symmetry of dChi/dr tensors is largely broken.
The max difference is:
320.8443495419144
warnings.warn('\n'.join(lines))
/opt/conda/envs/aiida-vibroscopy-demo/lib/python3.9/site-packages/aiida_vibroscopy/calculations/symmetry.py:163: UserWarning: Symmetry of dChi/dr tensors is largely broken.
The max difference is:
421.79293757428945
warnings.warn('\n'.join(lines))
08/25/2023 07:46:28 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [486|DielectricWorkChain|on_terminated]: cleaned remote folders of calculations: 515 570 603 606 609 612 615 668 671 674 677 680 683 734 737 740 743 746 749
08/25/2023 07:46:30 PM <799519> aiida.orm.nodes.process.workflow.workchain.WorkChainNode: [REPORT] [476|HarmonicWorkChain|on_terminated]: cleaned remote folders of calculations: 525 547 550 553 556 515 570 603 606 609 612 615 668 671 674 677 680 683 734 737 740 743 746 749
5.5. Powder Raman spectra using PBEsol+U+V#
Let’s plot the powder Raman spectra computed using PBEsol+U+V.
vibro = calc_hubbard.outputs.vibrational_data.numerical_accuracy_4
polarized_intensities, unpolarized_intensities, frequencies_hubbard, labels = vibro.run_powder_raman_intensities(frequency_laser=532, temperature=300)
total_intensities_hubbard = polarized_intensities + unpolarized_intensities
plt = get_spectra_plot(frequencies_hubbard, total_intensities_hubbard)
plt.show()
5.6. Comparison with experiments#
We compare the above results with the experimental results from M. Inaba et al., Chem. Lett., 24, 889 (1995).
Show cell source
import numpy as np
import matplotlib.pyplot as plt
from matplotlib.ticker import AutoMinorLocator
from aiida_vibroscopy.utils.broadenings import multilorentz
broadening = 10.0
# ----- Canvas
_, ax = plt.subplots()
# ----- Experiments
converter = lambda s: float(s.decode("UTF-8").replace(",", "."))
x_range, y_range = np.loadtxt("./_static/lico2_raman_exp.txt", unpack=True, converters={0:converter, 1:converter})
y_range = (y_range-y_range.min())/(y_range.max()-y_range.min())
ax.plot(x_range, y_range+2, linewidth=2.0, color='black', linestyle='--', label="Exp. M. Inaba (1995)")
# ----- PBEsol+U+V
y_range = multilorentz(x_range, frequencies_hubbard, total_intensities_hubbard, broadening)
y_range /= y_range.max()
ax.plot(x_range, y_range+1, linewidth=2.0, linestyle='-', label="PBEsol+U+V")
# ----- PBEsol
y_range = multilorentz(x_range, frequencies_pbesol, total_intensities_pbesol, broadening)
y_range /= y_range.max()
ax.plot(x_range, y_range+0, linewidth=2.0, linestyle='-', label="PBEsol")
# ----- Ticks and Labels
ax.set_yticklabels('')
ax.set_yticks([])
ax.tick_params(axis='both', which='both', direction='in')
ax.set_ylabel('Intensity (arb. units)') # Add a y-label to the axes.
ax.set_xlabel('Wavenumber (cm$^{-1}$)')
ax.xaxis.set_minor_locator(AutoMinorLocator())
ax.yaxis.set_minor_locator(AutoMinorLocator())
# ----- Legend
ax.legend(bbox_to_anchor=(0,1,1,0.1), loc='lower left', frameon=False)
plt.show()
Warning
We are here quite far from concluding something, i.e. which functional is the best. We used unconverged paramaters to fit the time of a tutorial. We should make sure to converge cutoff, k points and electronic charge convergence.
Moreover, temperature can affect comparison with experiments. In general, though, we can say that thermal expansion would lead to renormalization of the frequencies, usually reducing them. Thus, it is better to overestimate in general the phonon frequencies - although note that this is not a necessary condition.