I Introduction and motivation
Since, Feynman 1 proposed the idea of quantum computing and stated that quantum computers can simulate quantum mechanical systems exponentially faster, outstanding advancement has been made in simulation of quantum systems. In the last decade, the quantum simulation of closed and open systems has got overwhelming response among research communities. It promises powerful applications in the field of high energy physics, quantum chemistry and condensed matter, which are intractable on classical computers. The field of quantum computing is concern with the behavior and nature of energy at the quantum level to improve the efficiency of computations. The main aim of running quantum algorithms on quantum computers to solve the various computational tasks more efficiently and in less time as compared to existing classical ones. The quantum principle of superposition and entanglement is the backbone of quantum algorithms, which allows us to perform operations in large Hilbert spaces exponentially.
Tensor network states are a new language for quantum manybody systems based on pure quantum mechanical phenomena ’entanglement’ 2
. Tensor network states are classified on the basis of dimensions along which the tensors are crossed. Thus, it manages the exponentially growing Hilbert space by restricting the entanglement between two parts of the quantum system
4. Although, the dimension of Hilbert space increases exponentially with the increase in size of the system, there exists some quantum manybody systems whose implementations are realizable by classical computers. The density matrix renormalization group (DMRG) has been applied in lattice field theory to study particle physics with Monte Carlo simulations 33; 34. A matrix product state (MPS) (onedimensional system) is an example of tractable quantum systems numerically. In MPS, low entangled states are efficiently represented, which is not with large dimensions tensor network states. It has some basic properties such as dense nature, finite correlation, translational invariance and onedimensional area law 3.In the last decade, the simulations of time dependent quantum states have been widely used in different physical systems with MPS 43; 44; 45; 46. There are several softwares consisting high performance libraries and parameters build upon tensor network theory such as C++ library (ZKCM_QC) 37 for multiprecision, QCMaquis 38 for optimization purpose, EvoMPS 36 for simulating timedependent (real or imaginary) onedimensional manyparticles, DMRG++ 39 an open source implementation of DMRG algorithm, iTensor 40 a C++ library for executing tensor network algorithms, Uni10 41 an opensource and free C++ library for the construction of tensor network algorithms, QuTiP 20 open source python framework, OSMPSopensource MPS 42, Symbolic C++ quantum simulation 47 and many more.
In the past few years, the research has been increased in simulating experiments on IBM (International Business Machines Corporation) Q Experience platform. IBM has given access to real quantum computers and simulators, which allows researchers to develop, test and implement their experiments 5. Through IBM, we can examined the simulation results on classical computer and analysis on available quantum hardware to get a feel of the quantum system. It is becoming far more widespread and efficiently used to simulate several computational problems quickly 6. Currently, the platform has been extensively used to implement several experiments such as quantum tunneling simulation 7, Ising model simulation hard problems 9, quantum algorithms 10, quantum error correction 11
, quantum machine learning
14, benchmarking the quantum gates 28 to name a few.Recently, matrix product state (MPS) provides a stepping stone to various advancements in field of condensed matter and quantum computational theory: unsupervised learning using MPS
31, quantum dynamics 17; 19, open source MPS simulation 20, simulating quantum computation 22, simulation of open quantum systems 18, quantum finite state machines of MPS 1623and neural network representation
24 and MPS based efficient, productive applications for high performance computers 29 and can be employed in various emerging technologies such as optical computing, quantum cryptography, image recognition and dynamic quantum clustering 30. MPS provides an efficient approximation of realistic local Hamiltonians and can be generated by tensors sequentially. We focused on the maximum entangled MPS (GHZ and W states). Zizzi 32 has shown that it is possible to preserve the probabilities by executing a reversible quantum measurement on closed systems. In previous paper, we have efficiently simulated MPS with a broader quantum computational theory and investigated their relationship with quantum finitestate machine (QFSM) using unitary criteria 16. In this paper, we have implemented simulation on customized real IBM quantum computer and simulator. Further, the probability distributions among qubits is investigated. The paper is organized as follows: Section 2 is devoted to the family of matrix product state. In Section 3, simulation results are presented. Finally, Section 4 is the conclusion.Ii Matrix product state
Matrix product state is complete. Basically, it concede the extent of entanglement in bond dimensions. In fact, any pure quantum state can be described by substituting the coefficients e.g. rankN tensor by Nrank 3 tensors and 2rank by 2 tensors. In MPS, a pure quantum state is represented as:
(1) 
where are complex square matrices, d is dimension, represents the indices i.e. {0, 1} for qubits and Tr() denotes trace of matrices 17. Figure 1 shows the MPS as onedimensional array of tensors and an instance of finite system of 5 sites 3. The GHZ state and W state can be represented using MPS:

GHZ state: A GreenbergerHorneZeilinger (GHZ) 26 state of Nspins 1/2 is represented as a
(2) It is defined as maximum entangled state, consisting some nontrivial entanglement features 27. Pan et al. 48 presented schemes for preparing GHZ state remotely. The 4qubit GHZ state is represented as
. Its matrix product state and unitary matrix (
U) are represented as:(3) 
W state: The nqubit W state is represented as :
(4) The W state refers to the superposition of pure entangled states with same coefficients 27. It is different from above GHZ state. It represents a multipartite entanglement, where one of the qubits is in up state , while others are in down state . The 4qubit W state is represented as 17. The matrix representation is given as:
(5) (6)
Iii Simulation results
Before, we proceed to the simulation of matrix product state, it is useful to define the notion of some operators used over single qubit in quantum computation. The controlled NOT operator (CNOT) has two inputs as well as outputs 17. The flip operation is performed over the target qubit when the first qubit is 1. The quantum circuit for CNOT operator is shown in Fig 3 and its matrix representation is given in Eq. (10). The other highly used quantum operators are flip operator (X): flips the qubit, identity operator (I): it does nothing i.e. outputs the qubit as it is, Hadamard operator (H): it generates superposition of states with equal probability in computational basis. Their matrices are given in Eq. (79) respectively. In quantum circuits, such operators can be represented using quantum gates. The universal set of quantum gates are shown in Fig 2.

Flip operator (X):
(7) 
Identity operator (I):
(8) 
Hadamard operator (H):
(9) 
CNOT operator:
(10)
iii.1 GHZ state
Experimentally, we utilized the IBM Quantum Experience (IBMQE), a universal customized (2qubit, 3qubit and 4qubit) quantum computer to perform our simulation. The following are the steps to test experiment on IBMQE:

Construct the quantum circuit for the computational problem through Python programming language and specify it using graphical interface.

Execute the quantum circuit on the IBM simulator and test whether it produces works correctly or not.

The implementation of the quantum circuit is performed on the hardware processor for given number of shots. During each shot, the specified quantum processor is assigned to implement the quantum circuit.

Finally, the measurement is performed and the probability of each process is calculated.
Here, we started with the maximally entangled 2qubit GHZ state, which is represented as . Initially, we have taken two quantum as well as classical registers. Both quantum registers are initialized to zero. The Hadamard gate and CNOT gate is used to create an entanglement between both qubits, followed by measurement of quantum state q[0] and q[1] to classical registers c[0] and c[1] respectively. Correspondingly, the quantum circuit for 2qubit GHZ state is shown in Fig. 4. Further, the probability distribution is determined as 53.3906 % for state and 49.6094 % for respectively i.e. on summation it satisfies unitary property. The histogram of probability distribution for is shown in Fig. 5.
Next, we have considered 3qubit GHZ state interpreted as . To simulate the on quantum computer, we have taken three quantum registers to carry out computation and three classical registers for measurement purpose at the end. The Hadamard and CNOT operation between first and second quantum qubit is similar to . Then, the second CNOT operation is performed between second and third quantum qubit i.e. flips the q[3] whenever the q[2] is in excited state. The final results can be restirved from classical registers. It creates the 3qubit entangled GHZ state and its equivalent quantum circuit is given in Fig 6. Further, the probability distribution is calculated, which comes out to be same i.e. for both and respectively.
Further, the complex 4qubit GHZ quantum state is taken. A 4qubit quantum circuit for is designed to check whether it is capable of showing quantum behavior on implementing repeatedly. Following the above procedure, the Hadamard gate and CNOT gate produces an entanglement between q[3] and q[4] as shown in Fig 8. Finally, the measurement is readout after one shot and yields the occurrence of state with 50.2930 % and with 49.7070 % respectively, shown in Fig 9. It can be easily checked that after every shot on quantum simulator, the summation of probabilities of occurrence of both quantum states is comes out to be 1. Thus, it confirms the probability distribution of four spins1/2 GHZ state into single state i.e. highly entangled quantum state.
iii.2 W state
The GHZ state and W state represent different types of entanglement, where W state is said to be less entangled as compared to GHZ. On measurement, GHZ state is collapse into mixture of other states. But, in case of W state, it leaves the bipartite entanglements on measuring one of its subsystems.
In this subsection, we have considered 3 and 4qubit W state and their equivalent quantum circuits are designed and tested on quantum simulator. The 3qubit W state is considered. Here, we have used flip gate (X) to flip the qubit and Pauli matrix () to rotate the spin in ydirection with given theta (). After performing the flip operation over the third qubit, barrier is assigned between second and third to avoid the optimization of consecutive gates in a quantum circuit, as shown in Fig 10. Then, CNOT gate is applied successively between , and , to make one of qubits in excited state and other in ground state same time.
Finally, the equivalent quantum circuit for is tested on IBMQE and results are plotted in Fig. 11. The probability distribution of states is examined as approx 33 % for each and , and 34 % for . It agrees the superimposed 3qubit W state into singlet state. Furthermore, we constructed the quantum circuit for complicated 4qubit W state . Following the above computational procedure of , we have used combination of Hadamard gates before and after the CNOT operation. Thus, it execute the CNOT operation in reverse direction using the Kronecker product () i.e. equal to () matrix. The constructed quantum circuit for is shown in Fig 12. Finally, after simulation, the measurement is retrieved from classical registers and the histogram of probability distribution among the states is plotted in Fig 13. It shows the occurrence of state and with approx 25 % each, 24 % of and state with 26 % approximately. Thus, it can be easily checked that nqubit GHZ and W state can be efficiently simulated on a quantum computer after constructing complicated quantum circuits.
Iv Conclusion
In conclusion, we have demonstrated the two highly entangled matrix product states (GHZ and W). The main focus of this paper is towards the simulation on a quantum computer. We have designed quantum circuits for (2, 3 and 4qubit) GHZ state and (3 and 4qubit) W state. Further, the constructed quantum circuits have been simulated on a quantum computer and their probability distribution among the quantum states is investigated. To conclude, we have noticed that any matrix product state can be efficiently simulated on a quantum computer. This allows us to employ these circuits for efficient and productive quantum computational theoretical and experimental applications such as quantum machine learning, MPS based quantum algorithms and can be used in condensed matter physics for quantum state tomography in future.
Acknowledgments
Amandeep Singh Bhatia was supported by Maulana Azad National Fellowship (MANF), funded by Ministry of Minority Affairs, Government of India.
Comments
There are no comments yet.