Wolfram Function Repository
Instantuse addon functions for the Wolfram Language
Function Repository Resource:
A quantum cellular automaton model that evolves the tensor product of a collection of initial qubits using arbitrary compositions of unitary operators for a finite number of steps
ResourceFunction["QuantumTensorAutomaton"][matrices,initial,t] evolves initial for t steps using a composition of unitary matrices given by matrices. 
Apply the operator CNOT 10 times to an initial state:
In[1]:= 

Out[2]= 

Visualize the state vectors for the preceding:
In[3]:= 

Out[3]= 

Compute the probabilities (norms squared) for the preceding:
In[4]:= 

Out[4]= 

Plot the probabilities (norms squared) against the steps (time) for the preceding:
In[5]:= 

Out[5]= 

Compose multiple operators:
In[6]:= 

Out[6]= 

Random unitary operators will also work:
In[7]:= 

Out[7]= 

This work is licensed under a Creative Commons Attribution 4.0 International License