diff --git a/src/DAMASK_abaqus_exp.f b/src/DAMASK_abaqus_exp.f
index 1fab2472d..dc755f2e8 100644
--- a/src/DAMASK_abaqus_exp.f
+++ b/src/DAMASK_abaqus_exp.f
@@ -37,6 +37,7 @@ subroutine DAMASK_interface_init
    dateAndTime                                                                                      ! type default integer
  call date_and_time(values = dateAndTime)
  write(6,'(/,a)') ' <<<+-  DAMASK_abaqus_exp  -+>>>'
+ write(6,'(/,a)')   ' Roters et. al, Computational Materials Science, 2018'
  write(6,'(/,a)')              ' Version: '//DAMASKVERSION
  write(6,'(a,2(i2.2,a),i4.4)') ' Date:    ',dateAndTime(3),'/',&
                                             dateAndTime(2),'/',&
diff --git a/src/DAMASK_abaqus_std.f b/src/DAMASK_abaqus_std.f
index d15682c58..cf60781ce 100644
--- a/src/DAMASK_abaqus_std.f
+++ b/src/DAMASK_abaqus_std.f
@@ -37,6 +37,7 @@ subroutine DAMASK_interface_init
    dateAndTime                                                                                      ! type default integer
  call date_and_time(values = dateAndTime)
  write(6,'(/,a)') ' <<<+-  DAMASK_abaqus_std  -+>>>'
+ write(6,'(/,a)')   ' Roters et. al, Computational Materials Science, 2018'
  write(6,'(/,a)')              ' Version: '//DAMASKVERSION
  write(6,'(a,2(i2.2,a),i4.4)') ' Date:    ',dateAndTime(3),'/',&
                                             dateAndTime(2),'/',&
diff --git a/src/DAMASK_marc.f90 b/src/DAMASK_marc.f90
index 27640e5f2..fe636dc52 100644
--- a/src/DAMASK_marc.f90
+++ b/src/DAMASK_marc.f90
@@ -54,6 +54,7 @@ subroutine DAMASK_interface_init
 
  call date_and_time(values = dateAndTime)
  write(6,'(/,a)') ' <<<+-  DAMASK_Marc  -+>>>'
+ write(6,'(/,a)')   ' Roters et. al, Computational Materials Science, 2018'
  write(6,'(/,a)')              ' Version: '//DAMASKVERSION
  write(6,'(a,2(i2.2,a),i4.4)') ' Date:    ',dateAndTime(3),'/',&
                                             dateAndTime(2),'/',&
diff --git a/src/DAMASK_spectral.f90 b/src/DAMASK_spectral.f90
index 068aebbc6..c68511cb7 100644
--- a/src/DAMASK_spectral.f90
+++ b/src/DAMASK_spectral.f90
@@ -161,6 +161,7 @@ program DAMASK_spectral
 ! init DAMASK (all modules)
  call CPFEM_initAll(el = 1_pInt, ip = 1_pInt)
  write(6,'(/,a)')   ' <<<+-  DAMASK_spectral init  -+>>>'
+ write(6,'(/,a)')   ' Roters et. al, Computational Materials Science, 2018'
  write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"
 
diff --git a/src/homogenization_RGC.f90 b/src/homogenization_RGC.f90
index e1f4f4cee..acf8f5a0f 100644
--- a/src/homogenization_RGC.f90
+++ b/src/homogenization_RGC.f90
@@ -116,6 +116,10 @@ subroutine homogenization_RGC_init(fileUnit)
    line = ''
  
  write(6,'(/,a)')   ' <<<+-  homogenization_'//HOMOGENIZATION_RGC_label//' init  -+>>>'
+ write(6,'(/,a)')   ' Tjahjanto et. al, International Journal of Material Forming, 2(1):939–942, 2009 '
+ write(6,'(/,a)')   ' doi: 10.1007/s12289-009-0619-1 '
+ write(6,'(/,a)')   ' Tjahjanto et. al, Modelling and Simulation in Materials Science and Engineering, 18:015006, 2010 '
+ write(6,'(/,a)')   ' doi: 10.1088/0965-0393/18/1/015006 '
  write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"
 
diff --git a/src/spectral_damage.f90 b/src/spectral_damage.f90
index 1ac3c4c73..03bb5766b 100644
--- a/src/spectral_damage.f90
+++ b/src/spectral_damage.f90
@@ -103,6 +103,8 @@ subroutine spectral_damage_init()
    SNESVISetVariableBounds
 
  write(6,'(/,a)') ' <<<+-  spectral_damage init  -+>>>'
+ write(6,'(/,a)') ' Shanthraj et. al, Handbook of Mechanics of Materials, volume in press, '
+ write(6,'(/,a)') ' chapter Spectral Solvers for Crystal Plasticity and Multi-Physics Simulations. Springer, 2018 '
  write(6,'(a15,a)')   ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"
  
diff --git a/src/spectral_interface.f90 b/src/spectral_interface.f90
index b45e8316c..ba644ea9c 100644
--- a/src/spectral_interface.f90
+++ b/src/spectral_interface.f90
@@ -115,6 +115,7 @@ subroutine DAMASK_interface_init()
 
  call date_and_time(values = dateAndTime)
  write(6,'(/,a)') ' <<<+-  DAMASK_spectral  -+>>>'
+ write(6,'(/,a)') ' Roters et. al, Computational Materials Science, 2018'
  write(6,'(/,a)')              ' Version: '//DAMASKVERSION
  write(6,'(a,2(i2.2,a),i4.4)') ' Date:    ',dateAndTime(3),'/',&
                                             dateAndTime(2),'/',&
diff --git a/src/spectral_mech_AL.f90 b/src/spectral_mech_AL.f90
index 27a1def0c..86e38c81c 100644
--- a/src/spectral_mech_AL.f90
+++ b/src/spectral_mech_AL.f90
@@ -147,7 +147,9 @@ subroutine AL_init
    SNESSetFromOptions
    
  write(6,'(/,a)') ' <<<+-  DAMASK_spectral_solverAL init  -+>>>'
- write(6,'(a15,a)')   ' Current time: ',IO_timeStamp()
+ write(6,'(/,a)') ' Shanthraj et. al, International Journal of Plasticity, 66:31–45, 2015 '
+ write(6,'(/,a)') ' doi: 10.1016/j.ijplas.2014.02.006 '
+ write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"
 
 !--------------------------------------------------------------------------------------------------
diff --git a/src/spectral_mech_Basic.f90 b/src/spectral_mech_Basic.f90
index 171eeacad..080874ec1 100644
--- a/src/spectral_mech_Basic.f90
+++ b/src/spectral_mech_Basic.f90
@@ -135,7 +135,9 @@ subroutine basicPETSc_init
    SNESSetFromOptions
    
  write(6,'(/,a)') ' <<<+-  DAMASK_spectral_solverBasicPETSc init  -+>>>'
- write(6,'(a15,a)')   ' Current time: ',IO_timeStamp()
+ write(6,'(/,a)') ' Shanthraj et. al, International Journal of Plasticity, 66:31–45, 2015 '
+ write(6,'(/,a)') ' doi: 10.1016/j.ijplas.2014.02.006 '
+ write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"
 
 !--------------------------------------------------------------------------------------------------
diff --git a/src/spectral_mech_Polarisation.f90 b/src/spectral_mech_Polarisation.f90
index acd713c70..cafddd611 100644
--- a/src/spectral_mech_Polarisation.f90
+++ b/src/spectral_mech_Polarisation.f90
@@ -145,7 +145,9 @@ subroutine Polarisation_init
    SNESSetFromOptions
    
  write(6,'(/,a)') ' <<<+-  DAMASK_spectral_solverPolarisation init  -+>>>'
- write(6,'(a15,a)')   ' Current time: ',IO_timeStamp()
+ write(6,'(/,a)') ' Shanthraj et. al, International Journal of Plasticity, 66:31–45, 2015 '
+ write(6,'(/,a)') ' doi: 10.1016/j.ijplas.2014.02.006 '
+ write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"
 
 !--------------------------------------------------------------------------------------------------
diff --git a/src/spectral_thermal.f90 b/src/spectral_thermal.f90
index f89184543..b021f0c48 100644
--- a/src/spectral_thermal.f90
+++ b/src/spectral_thermal.f90
@@ -103,7 +103,9 @@ subroutine spectral_thermal_init
    
  mainProcess: if (worldrank == 0_pInt) then
    write(6,'(/,a)') ' <<<+-  spectral_thermal init  -+>>>'
-   write(6,'(a15,a)')   ' Current time: ',IO_timeStamp()
+   write(6,'(/,a)') ' Shanthraj et. al, Handbook of Mechanics of Materials, volume in press, '
+   write(6,'(/,a)') ' chapter Spectral Solvers for Crystal Plasticity and Multi-Physics Simulations. Springer, 2018 '
+   write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"
  endif mainProcess
  
diff --git a/src/spectral_utilities.f90 b/src/spectral_utilities.f90
index e3383f3d1..c5497c613 100644
--- a/src/spectral_utilities.f90
+++ b/src/spectral_utilities.f90
@@ -215,7 +215,9 @@ subroutine utilities_init()
    vecSize = 3_C_INTPTR_T, &
    tensorSize = 9_C_INTPTR_T
 
- write(6,'(/,a)')   ' <<<+-  spectral_utilities init  -+>>>'
+ write(6,'(/,a)') ' <<<+-  spectral_utilities init  -+>>>'
+ write(6,'(/,a)') ' Eisenlohr et. al, International Journal of Plasticity, 46:37–53, 2013 '
+ write(6,'(/,a)') ' doi: 10.1016/j.ijplas.2012.09.012 '
  write(6,'(a15,a)') ' Current time: ',IO_timeStamp()
 #include "compilation_info.f90"