##################################
# README for the GW100.aiida file.
##################################

We provide an GW100.aiida zip file which contains the export of the AiiDA database for all the calculations. You can export all the nodes into your own AiiDA database (>v1.5) by typing:

```
verdi import GW100.aiida
```

You will find a group for each molecule, by typing "verdi group list -A":

GW100/cell_fcc-dnn=13A/<name_of_the_molecule>

In addition, a Dict node with the "type_dict" extra equal to "instructions", in which you can find the uuid 
of the YamboConvergence node, which collects all the YamboWorkflow (single Yambo calculation) nodes.

To access the instruction Dict, it is possible to use the QueryBuilder:

```
from aiida import orm

qb = QueryBuilder()
qb.append(Group, filters={
        'label': <GROUP_NAME>}, 
              tag ='g')
qb.append(
    orm.Dict, 
    filters={
        'extras.type_dict':'instructions',},
        with_group='g',
        tag = 'structure_imported',
)
```

To directly access the YamboConvergence:

```
from aiida_yambo.workflows.yamboconvergence import YamboConvergence
qb = QueryBuilder()
qb.append(Group, filters={
        'label':'GW100/cell_fcc-dnn=13A/{}'.format(mol)}, 
              tag ='g')
qb.append(
    orm.Dict, 
    filters={
        'extras.type_dict':'instructions',},
        with_group='g',
        tag = 'structure_imported',
)
```

## Access the convergence history for a system:

You can access all the convergence history by means of the pandas library:

```
import pandas as pd

YamboConvergence_node = load_node(<YamboConvergence node uuid>)
    
    df = pd.DataFrame(YamboConvergence_node.outputs.history.get_dict()) 

```
The final converged YamboWorkflow can be accessed by: 
    
```
    final_workflow_uuid = df[df['useful'] == True]['uuid'].values[0]
```

Calculations can be accessed by their corresponding YamboWorkflow caller, whose uuid is store in df['uuid'] column of the df DataFrame.
##########

# Evaluation of the extrapolation:

The extrapolation procedure (A/x + b)(C/y + d) --> homo or lumo = b*d can be reproduced for each molecule, once exported the GW100.aiida and installed the aiida-yambo plugin, by using the analyse_molecule.py script:


```
verdi run analyse_molecule.py --mol <name_of_the_molecule>
```