Loading...
Technical Insight

Magazine Feature
This article was originally featured in the edition:
Issue 2 2024

Quantum PICs: Empowering designers with accurate simulations

News

Integrated photonics is an exciting potential route to quantum computation and communication. The Ansys Lumerical qINTERCONNECT solver can precisely simulate quantum PICs, enabling engineers to design these circuits without needing extensive quantum mechanics knowledge.

By Sebastian Gitt, Ruoshi Xu, and Ahsan Alam, R&D engineers at Ansys

Of the many promises of the quantum revolution, two of the most exciting are powerful quantum computing and secure quantum communication channels impervious to tampering. Researchers are making rapid progress in exploring integrated photonics as a potential route to both of these. However, crafting photonic circuits for these applications often demands a profound grasp of quantum mechanics and the nuanced behaviour of qubits – the elemental units of quantum information.

To expedite the advancement of this technology, there is therefore a pressing industry demand for a commercially accessible tool that can empower engineers who have limited quantum mechanics expertise to design and refine these circuits.Furthermore, as these quantum PICs grow larger and more sophisticated, it is increasingly critical to suppress errors due to circuit imperfections, which impede their reliability and scalability.

These flaws include manufacturing variances in circuit components, photons that are not perfectly indistinguishable, and losses in the circuit. Consequently, there is a growing need for simulation tools that can also model realistic quantum PICs.

For purely classical photonic chips, the simulation can be performed by commercially available tools like Ansys Lumerical INTERCONNECT, which can accurately model the evolution of a single photon in any PIC by calculating the S-parameter. However, for scenarios with more than one photon, non-classical interference between photons must be accounted for.


Figure 1. A circuit for a user-defined unitary matrix of any size can be constructed using Ansys Lumerical qINTERCONNECT. Using the Clements decomposition method, the circuit consists of a square mesh of beamsplitters, which in turn can be customised by the user to contain custom components or to be a different geometry altogether.

Enter the Lumerical qINTERCONNECT solver by Ansys, a cutting-edge solution poised to address this need. By facilitating precise simulation of quantum PICs, compatible with Process Design Kits (PDKs) from photonic foundries, qINTERCONNECT streamlines the design and simulation process for engineers, supporting them to navigate this intricate domain with confidence and ease.

The qINTERCONNECT solver uses the classical scattering matrix from INTERCONNECT to generate a quantum scattering matrix, which can then be used to calculate the evolution of a quantum state in a PIC [1]. By using the quantum state in the Fock (number) basis as input to the solver, qINTERCONNECT returns the output of the PIC as a density matrix. Additionally, the solver can report the fidelity of the final state with respect to an expected final state and the probability of success for measuring a given output.

Since qINTERCONNECT uses INTERCONNECT as a starting point to simulate the PIC, it benefits from the advanced features of INTERCONNECT. These include the ability to design with PDKs from different photonic foundries, to account for manufacturing variability in photonic components by using statistically enabled compact models, and advanced nonlinear model support for elements like nonlinear optical fibres.

While qINTERCONNECT works in the Fock basis, it supports more exotic inputs in terms of Fock states. Users can encode logical qubit states in terms of physical Fock states, as well as coherent states and squeezed states. Additionally, the output from non-classical photon sources, such as physics-based simulations of optical nonlinearities, can be used as inputs.

Two areas that have recently attracted increased interest are programmable unitary matrices, which have uses across many different areas within quantum photonics, and continuous-variable quantum key distribution (CV-QKD) for secure data communication. qINTERCONNECT, as we will demonstrate, can be used in both of these contexts.



Figure 2. (a) The probabilities of measuring up to three photons in the first channel of the output of the unitary matrix circuit for an input state of indistinguishable photons |1,1,1,0 >⟩ (red) and a state of distinguishable photons (blue) at six different interaction strengths, represented different times, as in [7]. (b) The probabilities of measuring any given output state across all four channels at the output of the unitary matrix circuit for the same indistinguishable state as in (a). (c) Same as (b) but coupling gap values are taken from a normal distribution with a mean value equal to the ideal coupling gap value of 200 nm and with a standard deviation of 10 nm.

Generating and simulating programmable unitary matrix circuits

Unitary operations are fundamental building blocks in both classical and quantum linear systems, since any linear operator can be decomposed into a product of unitary and diagonal operators [2]. In classical systems, they are commonly used in applications related to matrix multiplication and neural networks, whereas in quantum systems they can be used for boson sampling, graph similarity algorithms, and the simulation of quantum systems [3-5].

The main principle behind constructing a programmable unitary matrix lies in the fact that an N-dimensional unitary can be factorised into a product of 2-dimensional unitary matrices [6]. Hence, a collection of programmable 2-dimensional unitary matrices can be used to construct a programmable unitary matrix of any size.

Since Mach-Zehnder interferometers (MZIs) can be used for these 2-dimensional unitary matrices, integrated photonics is a particularly convenient platform for implementing programmable unitary matrices. A common implementation is the so-called Clements decomposition, which improves upon the design of the earlier proposed Reck decomposition by using a symmetrical design that reduces the optical depth of the circuit and is more robust to losses [7].

With Ansys Lumerical qINTERCONNECT, users can construct a PIC in INTERCONNECT representing their own user-defined unitary matrix of arbitrary size. Furthermore, custom circuit elements can be used in the construction of each MZI. This enables workflows where the designer starts with component design with Multiphysics tools like Ansys Lumerical FDTD, CHARGE, or MODE, then generates an accurate, statistically enabled compact model of the optimised element using CML Compiler, and finally creates the PIC in INTERCONNECT to simulate how quantum states will evolve through it.


Figure 3. An overview of the pilot-assisted CV-QKD circuit, from transmitter to receiver, being simulated in INTER-CONNECT.

Furthermore, the user has the option of using more complex MZIs as the fundamental building block in the unitary circuit, which can provide a greater robustness against manufacturing imperfections and losses, at the cost of increased circuit complexity [8]. Figure 1 depicts an overview of how the user constructs the unitary matrix circuit, as well as each hierarchical level of the circuit.

Once the INTERCONNECT project file containing the programmable unitary matrix circuit has been generated, it is ready for simulation in qINTERCONNECT, a quantum circuit solver that can calculate the frequency domain response of a PIC in the Fock, or photon number, basis. The solver takes the classical S-parameters given by INTERCONNECT and uses them to calculate a quantum S-matrix (QS-matrix), which relates input and output states defined in a Fock basis prescribed by the number of channels and a maximum total number of photons.

The QS-matrix specifies the quantum mechanical state of light exiting the circuit in terms of the non-classical input state. Consequently, the qINTERCONNECT solver enables users to specify an input state in the Fock basis, either as a state
in bra-ket notation or a density matrix, and the solver will yield the output state in terms of probability amplitudes for measuring discrete photon detection patterns in the output channels. Furthermore, if a target output state is supplied, the solver can also calculate the fidelity.

As an example, the qINTERCONNECT solver can be used to construct a photonic circuit to model the unitary evolution of a non-interacting Bose-Hubbard model, or “hopping” Hamiltonian with four sites. It can then simulate the probabilities of success for various measurement outcomes for three-photon input states of both distinguishable and indistinguishable photons.

A simulation was performed using an input state of |1,1,1,0>, corresponding to one photon in each of the first three channels and no photon in the final channel. Figure 2a and 2b show the measurement results for the output of the first channel using both distinguishable and indistinguishable photons, as well as the results for measuring all four channels using an input state consisting of indistinguishable photons, which agree with established results [5].

Using this workflow, it is possible to characterise the impact of manufacturing uncertainty on these measurement outcomes. For example, users can vary the coupler gap spacing in the directional coupler elements that make up each MZI in the unitary circuit, and then repeat the simulations, as illustrated in Figure 2c. Here, the coupling gap values are drawn from a distribution with a standard deviation of 10 nm, modelling the effects of fabrication variance. The impact this has on the measured probabilities of detecting various output states is clearly demonstrated.


Figure 4. (a) Calibration of the shot noise variance is accomplished by measuring the shot noise from the detector when the signal is disconnected and the LO is connected, and the electronic noise from when both the signal and LO are disconnected. (b) The power spectra for the shot noise and electronic noise.

Continuous-variable quantum key distribution

Ansys Lumerical tools can also be used to model circuits for continuous-variable quantum key distribution (CV-QKD). In these applications, keys are encoded in the quadrature components of strongly attenuated coherent states using quadrature phase shift keying (QPSK). Since the information in the quadrature components can be read out using homodyne detection, CV-QKD can be fully implemented using well-established telecom technology. This stands in contrast with discrete-variable QKD, which requires single-photon sources and detectors.

An eavesdropper intercepting the quantum signal will necessarily introduce additional noise into a CV-QKD system, which can be detected by the receiving party. If the noise level inherent in the system is minimised, the additional noise introduced due to the eavesdropper can be detected, and the protocol can be aborted until the channel is determined to be secure. Since the amount of noise tolerated in the system before the protocol is aborted is so low, it is important to be able to simulate the noise in the system accurately.

Given that the coherent states are weak, they cannot be used as an accurate phase reference to recover the encoded keys. However, a stronger pilot reference can be transmitted alongside the quantum signal for phase correction at the receiver side. To avoid crosstalk between the pilot and the signal, the pilot can be first shifted to a different frequency using suppressed carrier single-sideband modulation (SC-SSB), and then rotated to an orthogonal polarisation [9]. In addition, a pulse carver can be used to reduce the information that can be obtained from monitoring the transition between symbols. Figure 3 depicts an overview of the full CV-QKD circuit.

The most important figure of merit for a QKD system is the secure key rate, which, in general, is improved by minimising noise and losses in the system. Consequently, the goal of the simulation is to determine the transmittivity and total excess loss, the latter of which is attributable to the eavesdropper and must therefore be kept as low as possible.

Since the transmittivity and excess noise must be given in shot noise units, the first step is to perform a shot noise unit calibration, which consists of two simulations. First, the transmission line is disconnected from the receiver and a simulation is performed with only the local oscillator (LO) connected to the homodyne detectors to measure shot noise plus noise intrinsic to the detector. Next, a second simulation, with the LO disconnected, measures detector noise. The variances of both measured voltage values are subtracted to obtain the shot noise variance. These two steps are then repeated eight times to obtain an averaged value for the shot noise calibration. Figure 4 shows the power spectra for the shot noise and electrical background noise.


Figure 5. (a) Accumulated phase drift over time between the pilot and the LO. (b) Digital signal processing (DSP) is performed using the measured quadrature values of the pilot tone to correct to the phase of the quantum signal. (c) Comparison of simulated parameters in INTERCONNECT against experimentally determined values [9].

After the shot noise calibration is complete, a simulation is performed with the transmitter connected to the receiver and the measured quadrature values are normalised by the previously measured shot noise variance, which is necessary to relate the measured results to the language of Gaussian information [10]. As before, the results are averaged over eight simulations. Then, the pilot tone is used to correct the phase drift between the LO and the signal.

Next, a fraction of the key information is disclosed between the sender and the receiver, which allows the receiver to calculate the transmittivity and the total noise in the link. If the receiver is inaccessible to an eavesdropper, the receiver noise can, under a relaxed security assumption, be treated as trusted noise, and is no longer attributed to the eavesdropper. In this case, the total noise minus detector noise will be another important figure of merit (ξtot - ξdet). Comparing the simulated parameters to those obtained experimentally and reported in the literature, we can see consistent results in the transmittivity (T), signal-to-noise ratio (SNR), and total noise (ξtot), as shown in Figure 5.

Since the transmittivity is affected not only by the transmission channel losses but also by detector losses and any coupling losses between the transmission channel and the detector, its value will be strongly dependent on the details of any given experimental implementation. Then, assuming a reconciliation efficiency value of β = 0.97, a frame-error rate of FER = 0.05, and a fraction disclosed for parameter estimation of ν = 0.25, and using the parameters in Figure 5, the secret key rate is found to be 38.0 Mbit/s over 13 km [9].

A versatile tool for designing and refining quantum PICs
Thanks to Ansys Lumerical qINTERCONNECT, you no longer need to be an expert in quantum mechanics to design a functional and effective integrated quantum photonic circuit. By leveraging foundry PDKs and advanced statistical and nonlinear compact models for different building blocks in their PICs, designers can account for issues originating from fabrication variations, propagation losses, and indistinguishable photons.

Designers can investigate the evolution of quantum states propagating through their quantum PICs using various types of sources ranging from physical Fock states to logical qubit states, such as coherent states or squeezed states to even non-classical photon sources. This makes qINTERCONNECT a versatile tool for the design and evaluation of integrated quantum photonic circuits.







Lightwave Logic receives ECOC Innovation Award for Hybrid PIC/Optical Integration Platform
Coherent wins ECOC award for datacentre innovation
HyperLight announces $37 million funding round
Jabil expands silicon photonics capabilities
Ephos raises $8.5 million for glass-based photonic chips
Designing for manufacture: PAM-4 transmitters using segmented-electrode Mach-Zehnder modulators
OpenLight and Epiphany partner on PIC ecosystem
NewPhotonics and SoftBank team up on advanced photonics
POET and Mitsubishi collaborate on 3.2T optical engines
Integrated photonic platforms: The case for SiC
Integrating high-speed germanium modulators with silicon photonics and fast electronics
Lightium Secures $7 Million Seed Funding
Revolutionising optoelectronics with high-precision bonding
Fraunhofer IMS invites participation in PIC engineering runs
Advances in active alignment engines for efficient photonics device test and assembly
Aeva announces participation at IAA Transportation 2024
Sumitomo Electric announces participation in ECOC 2024
Quside receives NIST certification for quantum entropy source
DustPhotonics launches industry-first merchant 1.6T silicon photonics engine
Arelion and Ciena announce live 1.6T wave data transmission
DGIST leads joint original semiconductor research with the EU
POET Technologies reorganises engineering team
A silicon chip for 6G communications
South Dakota Mines wins $5 million from NSF for Quantum Materials Institute
HieFo indium phosphide fab resumes production
Coherent launches new lasers for silicon photonics transceivers
AlixLabs wins funding from PhotonHub Europe
Sandia National Labs and Arizona State University join forces
Perovskite waveguides for nonlinear photonics
A graphene-based infrared emitter
Atom interferometry performed with silicon photonics
A step towards combining the conventional and quantum internet

×
Search the news archive

To close this popup you can press escape or click the close icon.
Logo
x
Logo
×
Register - Step 1

You may choose to subscribe to the PIC Magazine, the PIC Newsletter, or both. You may also request additional information if required, before submitting your application.


Please subscribe me to:

 

You chose the industry type of "Other"

Please enter the industry that you work in:
Please enter the industry that you work in: