diff --git a/README.md b/README.md
index ebc0c51..a3ba034 100644
--- a/README.md
+++ b/README.md
@@ -252,3 +252,7 @@ with the following keys.
 | `tec_i`        | Amperes         | TEC output current feedback derived from `i_tec`     |
 | `tec_u_meas`   | Volts           | Measurement of the voltage across the TEC            |
 | `pid_output`   | Amperes         | PID control output                                   |
+
+## PID Tuning
+
+The thermostat implements a PID control loop for each of the TEC channels, more details on setting up the PID control loop can be found [here](./doc/PID%20tuning.md).
diff --git a/doc/PID tuning.md b/doc/PID tuning.md
new file mode 100644
index 0000000..fbc15bd
--- /dev/null
+++ b/doc/PID tuning.md	
@@ -0,0 +1,81 @@
+# PID Tuning
+
+## Note on hardware setup
+
+When using a TEC module with the Thermostat, the Thermostat expects the thermal load (where the thermistor is connected) to heat up when a positive current flow from the TEC + terminal, through the TEC, to the TEC - terminal, and cool down when the current flows in the reverse direction.
+
+The heat sinking side of the TEC module should be thermally bonded to a large heat-sinking thermal mass to ensure maximum temperature stability, a large optical table had provided good results in tests.
+
+The thermal load under control should be well insulated from the surrounding for maximum stability, closed cell foam had been tested showing good results.
+
+## Real time plot
+
+When tuning Thermostat PID parameters, it is helpful to view the temperature, PID output and other data in the form of a real time graph.
+
+To use the Python real-time plotting utility, run
+
+```shell
+python pytec/plot.py
+```
+
+![default view](./assets/default%20view.png)
+
+## Temperature Setpoints and Thermal Load
+
+A PID controller with the same set of PID parameters may not work identically across all temperatures, especially when comparing the performance of a TEC module cooling a load versus heating a load. This is due to self ohmic heating of the TEC module aiding efficiency when  heating, but harming efficiency when cooling.
+
+When a PID loop is expected to operate the TEC in both heating and cooling modes, it is important to verify the loop performance in both modes.
+
+For systems expected to operate at a narrow range of temperatures, it is a good idea to tune the PID loop at the temperature region of interest.
+
+The same is also true for controlling loads that are expected to produce heat, e.g. laser cooling blocks. Testing the loop performance across varying amount of thermal load is needed to ensure stability in operation.
+
+## Manual Tuning
+
+Below are some general guidelines for manually tuning PID loops. Note that every system is different, and some of the values mentioned below may not apply to all systems.
+
+1. To start the manual tuning process, set the kp, ki and kd parameters to 0.
+
+2. Begin by increasing kp until the temperature begins to oscillate. Offset between the target temperature and the actual temperature can be ignored for now.
+
+3. Reduce kp by 30%, increase ki until the offset between target and actual temperature is eliminated.
+
+4. Increase kd until the maximum allowable amount of overshoot is observed.
+
+5. Some tweaking will be needed to obtain the desired result, especially when trying to balance between minimizing overshoot and maximizing response speed.
+
+## Auto Tuning
+
+A PID auto tuning utility is provided in the Pytec library. The auto tuning utility drives the the load to a controlled oscillation, observes the ultimate gain and oscillation period and calculates a set of PID parameters.
+
+To run the auto tuning utility, run
+
+```shell
+python pytec/autotune.py
+```
+
+After some time, the auto tuning utility will output the auto tuning results, below is a sample output
+
+```shell
+Ku: 0.7553203471147422
+Pu: 75.93899999999977
+rule: ziegler-nichols
+kp: 0.45319220826884526
+Ki: 0.011935690706194357
+Kd: 4.301870387965967
+rule: tyreus-luyben
+kp: 0.3432930977636503
+Ki: 0.0020549280832497956
+Kd: 4.137825730504864
+        .
+        .   
+        .
+```
+
+At the end of the test, the ultimate gain `Ku`, oscillation period `Pu` and a few sets of recommended PID parameters are calculated and displayed.  
+
+Multiple suggested sets of PID parameters based on different calculation rules are displayed. While all sets are expected to work, the different sets trade off response time with overshoot differently, and testing is needed to see which set works best for the system on hand.
+
+With a well designed and constructed setup, the PID parameters calculated by the auto tune utility can provide mK level control stability with no manual tweaking.
+
+![tec3pessen](./assets/tec3pessen.png)
diff --git a/doc/assets/default view.png b/doc/assets/default view.png
new file mode 100644
index 0000000..6a77c6e
Binary files /dev/null and b/doc/assets/default view.png differ
diff --git a/doc/assets/tec3pessen.png b/doc/assets/tec3pessen.png
new file mode 100644
index 0000000..bddb90c
Binary files /dev/null and b/doc/assets/tec3pessen.png differ