def one_and_two_body_interaction(p, q, a, b) -> cirq.OP_TREE: yield CRxxyy(self.hamiltonian.one_body[p, q].real * time).on( control_qubit, a, b) yield CRyxxy(self.hamiltonian.one_body[p, q].imag * time).on( control_qubit, a, b) yield rot111(-2 * self.hamiltonian.two_body[p, q] * time).on( control_qubit, a, b)
def trotter_step( self, qubits: Sequence[cirq.Qid], time: float, control_qubit: Optional[cirq.Qid]=None ) -> cirq.OP_TREE: if not isinstance(control_qubit, cirq.Qid): raise NotImplementedError('Control qudit must be specified.') n_qubits = len(qubits) # Change to the basis in which the one-body term is diagonal yield bogoliubov_transform( qubits, self.one_body_basis_change_matrix.T.conj()) # Simulate the one-body terms. for p in range(n_qubits): yield rot11(rads= -self.one_body_energies[p] * time ).on(control_qubit, qubits[p]) # Simulate each singular vector of the two-body terms. prior_basis_matrix = self.one_body_basis_change_matrix for j in range(len(self.eigenvalues)): # Get the two-body coefficients and basis change matrix. two_body_coefficients = self.scaled_density_density_matrices[j] basis_change_matrix = self.basis_change_matrices[j] # Merge previous basis change matrix with the inverse of the # current one merged_basis_change_matrix = numpy.dot(prior_basis_matrix, basis_change_matrix.T.conj()) yield bogoliubov_transform(qubits, merged_basis_change_matrix) # Simulate the off-diagonal two-body terms. yield swap_network( qubits, lambda p, q, a, b: rot111( -2 * two_body_coefficients[p, q] * time).on( control_qubit, a, b)) qubits = qubits[::-1] # Simulate the diagonal two-body terms. yield (rot11(rads= -two_body_coefficients[k, k] * time).on( control_qubit, qubits[k]) for k in range(n_qubits)) # Update prior basis change matrix. prior_basis_matrix = basis_change_matrix # Undo final basis transformation. yield bogoliubov_transform(qubits, prior_basis_matrix) # Apply phase from constant term yield cirq.rz(rads= -self.hamiltonian.constant * time).on(control_qubit)
def one_and_two_body_interaction_reverse_order(p, q, a, b) -> cirq.OP_TREE: yield rot111(-self.hamiltonian.two_body[p, q] * time).on( cast(cirq.Qid, control_qubit), a, b) yield CRyxxy(0.5 * self.hamiltonian.one_body[p, q].imag * time).on( cast(cirq.Qid, control_qubit), a, b) yield CRxxyy(0.5 * self.hamiltonian.one_body[p, q].real * time).on( cast(cirq.Qid, control_qubit), a, b)
def two_body_interaction(p, q, a, b) -> cirq.OP_TREE: yield rot111(-2 * self.hamiltonian.two_body[p, q] * time).on( control_qubit, a, b)