| //
|
| // Copyright (c) 2026 BEL ESPRIT D ACCORD TRUST HOLDINGS INC
|
| // All rights reserved.
|
|
|
| // BlackHoleEngine/Quantinuum.qs
|
| // Hardware-specific: Quantinuum H1/H2 (all-to-all, native ZZ)
|
| namespace BlackHoleEngine.Quantinuum {
|
| open Microsoft.Quantum.Intrinsic;
|
| open Microsoft.Quantum.Canon;
|
| open BlackHoleEngine.Primitives;
|
|
|
| operation QuantinuumXX(q1 : Qubit, q2 : Qubit, theta : Double) : Unit is Adj + Ctl {
|
| within {
|
| H(q1); H(q2);
|
| } apply {
|
| Exp(, theta, );
|
| }
|
| }
|
|
|
| operation QuantinuumHaydenPreskill(
|
| qubits : Qubit,
|
| scrambleTime : Double,
|
| couplings : Double,
|
| fields : Double
|
| ) : Unit is Adj + Ctl {
|
| let n = Length(qubits);
|
| let dt = scrambleTime / 5.0;
|
| let steps = 5;
|
|
|
| for _ in 1..steps {
|
| for i in 0..n-1 {
|
| Rz(2.0 * fields * dt, qubits);
|
| }
|
|
|
| for layer in GetXXLayers(n, couplings) {
|
| for (i, j) in layer {
|
| let J = couplings;
|
| if (AbsD(J) > 1e-10) {
|
| QuantinuumXX(qubits, qubits, 2.0 * J * dt);
|
| }
|
| }
|
| }
|
| }
|
| }
|
|
|
| operation QuantinuumLindbladStep(
|
| ancilla : Qubit,
|
| qubits : Qubit,
|
| viol : ViolatingAssignment,
|
| gamma_dt : Double
|
| ) : Unit {
|
| let sqrt_gamma = Sqrt(gamma_dt);
|
| let theta = 2.0 * Arcsin(sqrt_gamma);
|
|
|
| Ry(theta, ancilla);
|
|
|
| let (va, vb, vc) = viol;
|
| if (va) { X(qubits); }
|
| if (vb) { X(qubits); }
|
| if (vc) { X(qubits); }
|
|
|
| // Quantinuum has high-fidelity CCZ
|
| CCZ(ancilla, qubits, qubits);
|
| CCZ(ancilla, qubits, qubits);
|
| CCZ(ancilla, qubits, qubits);
|
|
|
| if (va) { X(qubits); }
|
| if (vb) { X(qubits); }
|
| if (vc) { X(qubits); }
|
|
|
| let outcome = M(ancilla);
|
| Reset(ancilla);
|
| }
|
| }
|
|
|
