diff --git a/.github/workflows/Python.yml b/.github/workflows/Python.yml
index 4c199b81c..b3d714f92 100644
--- a/.github/workflows/Python.yml
+++ b/.github/workflows/Python.yml
@@ -9,7 +9,7 @@ jobs:
 
     strategy:
       matrix:
-        python-version: ['3.7', '3.8', '3.9'] #, '3.10']
+        python-version: ['3.8', '3.9'] #, '3.10']
         os: [ubuntu-latest, macos-latest]
 
     steps:
diff --git a/LICENSE b/LICENSE
index b4318ac92..8f9664407 100644
--- a/LICENSE
+++ b/LICENSE
@@ -1,4 +1,4 @@
-Copyright 2011-2021 Max-Planck-Institut für Eisenforschung GmbH
+Copyright 2011-2022 Max-Planck-Institut für Eisenforschung GmbH
 
 DAMASK is free software: you can redistribute it and/or modify
 it under the terms of the GNU Affero General Public License as published by
diff --git a/PRIVATE b/PRIVATE
index 2495e15d0..5598ec489 160000
--- a/PRIVATE
+++ b/PRIVATE
@@ -1 +1 @@
-Subproject commit 2495e15d064dea61e0faf078f817e32a637d6eda
+Subproject commit 5598ec4892b0921fccf63e8570f9fcd8e0365716
diff --git a/python/damask/VERSION b/python/damask/VERSION
index 663e61aee..750043c40 100644
--- a/python/damask/VERSION
+++ b/python/damask/VERSION
@@ -1 +1 @@
-v3.0.0-alpha5-495-g814d89532
+v3.0.0-alpha5-518-g4fa97b9a3
diff --git a/python/setup.py b/python/setup.py
index 1740b54a5..abb8053a9 100644
--- a/python/setup.py
+++ b/python/setup.py
@@ -16,7 +16,7 @@ setuptools.setup(
     url='https://damask.mpie.de',
     packages=setuptools.find_packages(),
     include_package_data=True,
-    python_requires = '>=3.7',
+    python_requires = '>=3.8',
     install_requires = [
         'pandas>=0.24',                                                                             # requires numpy
         'numpy>=1.17',                                                                              # needed for default_rng
@@ -30,7 +30,7 @@ setuptools.setup(
         'Intended Audience :: Science/Research',
         'Topic :: Scientific/Engineering',
         'Programming Language :: Python :: 3',
-        'License :: OSI Approved :: GNU General Public License v3 or later (GPLv3+)',
+        'License :: OSI Approved :: GNU Affero General Public License v3 or later (AGPLv3+)',
         'Operating System :: OS Independent',
     ],
 )
diff --git a/src/phase_mechanical_eigen_thermalexpansion.f90 b/src/phase_mechanical_eigen_thermalexpansion.f90
index 087cf7d3d..70bf84cd8 100644
--- a/src/phase_mechanical_eigen_thermalexpansion.f90
+++ b/src/phase_mechanical_eigen_thermalexpansion.f90
@@ -99,11 +99,11 @@ module subroutine thermalexpansion_LiAndItsTangent(Li, dLi_dTstar, ph,me)
   associate(prm => param(kinematics_thermal_expansion_instance(ph)))
     Li = dot_T * ( &
                   prm%A(1:3,1:3,1) &                                                                ! constant  coefficient
-                + prm%A(1:3,1:3,2)*(T - prm%T_ref)**1 &                                             ! linear    coefficient
+                + prm%A(1:3,1:3,2)*(T - prm%T_ref) &                                                ! linear    coefficient
                 + prm%A(1:3,1:3,3)*(T - prm%T_ref)**2 &                                             ! quadratic coefficient
                 ) / &
          (1.0_pReal &
-               + prm%A(1:3,1:3,1)*(T - prm%T_ref)**1 / 1.0_pReal &
+               + prm%A(1:3,1:3,1)*(T - prm%T_ref)    / 1.0_pReal &
                + prm%A(1:3,1:3,2)*(T - prm%T_ref)**2 / 2.0_pReal &
                + prm%A(1:3,1:3,3)*(T - prm%T_ref)**3 / 3.0_pReal &
          )
diff --git a/src/phase_mechanical_elastic.f90 b/src/phase_mechanical_elastic.f90
index bb530747b..8972367d5 100644
--- a/src/phase_mechanical_elastic.f90
+++ b/src/phase_mechanical_elastic.f90
@@ -101,32 +101,32 @@ pure module function elastic_C66(ph,en) result(C66)
     T = thermal_T(ph,en)
 
     C66(1,1) = prm%C_11(1) &
-             + prm%C_11(2)*(T - prm%T_ref)**1 &
+             + prm%C_11(2)*(T - prm%T_ref) &
              + prm%C_11(3)*(T - prm%T_ref)**2
 
     C66(1,2) = prm%C_12(1) &
-             + prm%C_12(2)*(T - prm%T_ref)**1 &
+             + prm%C_12(2)*(T - prm%T_ref) &
              + prm%C_12(3)*(T - prm%T_ref)**2
 
     C66(4,4) = prm%C_44(1) &
-             + prm%C_44(2)*(T - prm%T_ref)**1 &
+             + prm%C_44(2)*(T - prm%T_ref) &
              + prm%C_44(3)*(T - prm%T_ref)**2
 
 
     if (any(phase_lattice(ph) == ['hP','tI'])) then
       C66(1,3) = prm%C_13(1) &
-               + prm%C_13(2)*(T - prm%T_ref)**1 &
+               + prm%C_13(2)*(T - prm%T_ref) &
                + prm%C_13(3)*(T - prm%T_ref)**2
 
       C66(3,3) = prm%C_33(1) &
-               + prm%C_33(2)*(T - prm%T_ref)**1 &
+               + prm%C_33(2)*(T - prm%T_ref) &
                + prm%C_33(3)*(T - prm%T_ref)**2
 
     end if
 
     if (phase_lattice(ph) == 'tI') then
       C66(6,6) = prm%C_66(1) &
-               + prm%C_66(2)*(T - prm%T_ref)**1 &
+               + prm%C_66(2)*(T - prm%T_ref) &
                + prm%C_66(3)*(T - prm%T_ref)**2
     end if