S. Li was supported by the U.S. Department of Energy under contract DE–AC02–76SF00515, and P. Lowdon by the Agence Na-tionale de la Recherche under project No. ANR–18–ERC1–0002. P.J. Silva acknowledges support from FCT under the contracts SFRH/BPD/40998/2007 and SFRH/BPD/109971/2015, and both O. Oliveira and P.J. Silva acknowledge support from FCT under the contracts UID/FIS/04564/2016 and CERN/FIS-COM/0029/2017. This work was granted access to the HPC resources of the PDC Cen-ter for High Performance Computing at the KTH Royal Institute of Technology, Sweden, made available within the Distributed Eu-ropean Computing Initiative by the PRACE-2IP, receiving fund-ing from the European Community’s Seventh Framework Pro-gramme (FP7/2007-2013) under grand agreement no. RI-283493. The use of Lindgren has been provided under DECI-9 project COIM-BRALATT. We acknowledge that the results of this research have been achieved using the PRACE-3IP project (FP7 RI312763) re-source Sisu based in Finland at CSC. The use of Sisu has been provided under DECI-12 project COIMBRALATT2. We also acknowl-edge the Laboratory for Advanced Computing at the University of Coimbra (http://www.uc .pt /lca) for providing access to the HPC re-source Navigator.