Papers

The variational quantum eigensolver (VQE) is a hybrid quantum-classical algorithm designed for noisy intermediate-scale quantum (NISQ) hardware to estimate eigenvalues of many-body Hamiltonians. Unlike fully quantum approaches such as quantum phase estimation (QPE), VQE trades deep coherent circuits for repeated state preparation, measurement, and classical optimization, making it more compatible with limited qubit counts and finite coherence times. Recent developments have focused on reducing quantum resource requirements while retaining chemically meaningful wavefunction structure. In this paper, we examine recent progress in VQE methods for molecular electronic structure with an emphasis on three themes: (i) strategies for circuit and ansatz complexity reduction, including adaptive and selectively screened approaches, (ii) chemically motivated workflows that combine VQE with orbital optimization, fragmentation, and localized active-space ideas to better address strong correlation, and (iii) extensions of VQE to excited-state calculations. Throughout, we emphasize the tradeoffs among parameter count, gate depth, symmetry preservation, measurement overhead, and classical preprocessing, and discuss where these approaches may become most useful for chemically challenging active spaces. We also highlight benchmarking considerations for assessing both accuracy and resource requirements, and conclude with a perspective on regimes in which VQE may offer the greatest long-term value, particularly multireference active spaces and low-lying excited-state manifolds.

The training of Quantum Neural Networks (QNNs) is hindered by the high computational cost of gradient estimation and the barren plateau problem, where optimization landscapes become intractably flat. To address these challenges, we introduce Weighted Stochastic Block Descent (WSBD), a novel optimizer with a dynamic, parameter-wise freezing strategy. WSBD intelligently focuses computational resources by identifying and temporarily freezing less influential parameters based on a gradient-derived importance score. This approach significantly reduces the number of forward passes required per training step and helps navigate the optimization landscape more effectively. Unlike pruning or layer-wise freezing, WSBD maintains full expressive capacity while adapting throughout training. Our extensive evaluation shows that WSBD converges on average 63.9% faster than Adam for the popular ground-state-energy problem, an advantage that grows with QNN size. We provide a formal convergence proof for WSBD and show that parameter-wise freezing outperforms traditional layer-wise approaches in QNNs. Project page: https://github.com/Damrl-lab/WSBD-Stochastic-Freezing-Optimizer.

Entangling interactions between distant qubits can be mediated via an additional degree of freedom. In conventional trapped-ion schemes, realizing a well-defined, coherent gate typically requires spectrally addressing a specific bus mode. As the ion number increases, the coupling to each individual motional mode becomes weaker, so gates on large ion strings mediated by a single mode are necessarily slow. Moreover, addressing a large number of modes demands complex driving schemes, and the fundamentally perturbative character of these approaches imposes constraints on achievable gate speed and fidelity. Here, we introduce a scheme for entangling trapped-ion qubits using a time-dependent magnetic-field gradient, in which all axial motional modes participate in mediating the interaction and the gate construction is nonperturbative. The framework can be used to implement both multi-qubit gates and two-qubit gates between arbitrary pairs in a linear ion string. Through several explicit examples, we highlight the advantages over existing magnetic-gradient schemes and show how gates on multiple ion pairs can be carried out simultaneously.

We construct a new nonlinear deformed Schrödinger structure using a nonlinear derivative operator which depends on a parameter $q$. This operator recovers Newton derivative when $q \rightarrow 1$. Using this operator we propose a deformed Lagrangian which gives us a deformed nonlinear Schrödinger equation with a nonlinear kinetic energy term and a standard potential $V(\vec{x})$. We analytically solve the nonlinear deformed Schrödinger equation for $V(\vec{x}) = 0$ and $q \simeq1$. This model has a continuity equation, the energy is conserved, as well as the momentum and also interacts with electromagnetic field. Planck relation remains valid and in all steps we easily recover the undeformed quantities when the deformation parameter goes to 1. Finally, we numerically solve the equation of motion for the free particle in any spatial dimension, which shows a solitonic pattern when the space is equal to one for particular values of $q$.

Besides well-known conditions of locality or factorisability, deriving the Bell inequalities requires assuming that the distribution of hidden variables and Alice's and Bob's measurement settings be independent of each other. We show that (analogously to violations of locality due to action at a distance) certain violations of this Measurement Independence assumption can be associated with a notion of signalling in principle, thus making them also testable in principle, and spell out the appropriate conditions. Accordingly, we show that by imposing no-signalling one can prove a version of Bell's theorem that does not require the assumption of Measurement Independence. We discuss the "Schulman model" as an example, as well as lessons for "experimental metaphysics".

It is known that two local unitaries (LU) equivalent states possess the same amount of entanglement and can be used to perform the same tasks in quantum information theory (QIT). For a protocol for a task in QIT, we call a protocol LU invariant if two LU-equivalent states are either both suitable for the protocol or neither is. So far, no one has discussed whether a protocol for a task in QIT is LU invariant. In [Phys. Rev. A, 74, 062320 (2006)], Agrawal and Pati proposed the perfect teleportation protocol (PTP) and the protocol for superdense coding to transmit 2-bit classical information by sending one qubit (PSDC-2) and 3-bit classical information by sending two qubits (PSDC-3). In this paper, we show that PTP and PSDC-2 are LU invariant. That is, two LU equivalent states are suitable for PTP and PSDC-2 or neither of them is. We show that PSDC-3 is not LU invariant. We also indicate that the teleportation proposed in <cite>Nielsen</cite> is not LU invariant. We give a necessary and sufficient condition for a state to be suitable for PTP, PSDC-2, and PSDC-3, respectively. Via the LU invariance of PTP and PSDC-2, we prove that a state is suitable for PTP and PSDC-2 if and only if it has 1 ebit of shared entanglement, respectively and find all genuine entangled states and separable states which are suitable for PTP and PSDC-2, respectively. So far, no one has indicated that PTP and PSDC-2 do not require genuine entanglement. Agrawal and Pati suggested to study if there are subclasses of W SLOCC class which are suitable for PSDC-3. So far, it still remains an unsolved question. We show that any state of the SLOCC class W is not suitable for PSDC-3.

Integer linear programming (ILP) remains computationally challenging due to its NP-complete nature despite its central role in scheduling, logistics, and design optimization. We introduce a fully quantum Metropolis-Hastings algorithm for ILP that implements a coherent random walk over the discrete feasible region using only reversible quantum circuits, without quantum-RAM assumptions or classical pre/post-processing. Each walk step is a unitary update that prepares coherent candidate moves, evaluates the objective and constraints reversibly – including a constraint-satisfaction counter to enforce feasibility – and encodes Metropolis acceptance amplitudes via a low-overhead linearized rule. At the logical level, the construction uses $\mathcal{O}(n\log_2 N)$ qubits to represent $n$ integer variables over the interval $[-N,\,N-1]$, and the Toffoli-equivalent cost per Metropolis step grows linearly with the total logical qubit count. Using explicit ripple-carry adder constructions, we support linear objectives and mixed equality/inequality constraints. Numerical circuit-level simulations on a broad ensemble of randomly generated instances validate the predicted linear resource scaling and exhibit progressive thermalization toward low-cost feasible solutions under the annealing schedule. Overall, the method provides a coherent, resource-characterized baseline for fully quantum constraint programming and a foundation for incorporating additional quantum speedups in combinatorial optimization.

We present a theoretical study of quantum nonreciprocity induced via a Dzyaloshinskii-Moriya interaction (DMI) in an otherwise achiral, waveguide quantum electrodynamics. Using the full quantum master equation and input-output formalism for two-level systems coupled to a one-dimensional waveguide and driven by a coherent field, we show that an engineered DMI enables strong nonreciprocity in an otherwise reciprocal system, with tunable behavior governed by driving strength, detunings, and phase of the DMI. Using it not only demonstrates nonreciprocal transmission but also demonstrates nonreciprocal quantum entanglement and photon bunching. The system can end up in a pure state as certain decohering channels do not participate. The pure state leads to power-independent perfect transparency. Conditions are derived and depend on the propagation phase, the relative detuning of the two qubits, and the exchange interaction. At these pure-state points, the steady-state entanglement is reciprocal and admits a closed-form expression; away from them, phase control generates strong entanglement nonreciprocity. The DMI also reshapes photon statistics, redistributing two-photon correlations and shifting superbunching from transmission (no DMI) to reflection at finite DMI. These results establish DMI as a versatile resource for engineering nonreciprocity, transparency, entanglement, and photon correlations in waveguide QED, enabling isolators, routers, and superbunching light sources without requiring chiral waveguides.

Krylov subspace methods quantify operator growth in quantum many-body systems through Lanczos coefficients that encode how operators spread under time evolution. Although these diagnostics were originally motivated by questions of chaos and integrability, quadratic fermionic Hamiltonians are often expected to exhibit trivial Lanczos structure. Here we show that, in the long-range Kitaev chain, Lanczos coefficients generated from local boundary operators sharply diagnose whether the lowest excitation gap is controlled by boundary-localized or bulk-extended modes. We introduce the $Krylov$ $staggering$ $parameter$ for the Lanczos coefficients. In the short-range Kitaev chain with balanced hopping and pairing, we derive analytically for arbitrary system size (valid in the thermodynamic limit) and show that this quantity is exactly constant and its sign cleanly distinguishes the topological phase with Majorana edge modes from the trivial phase. Away from that limit, long-range couplings and pairing-hopping imbalance deform the simple flat structure and analytical control is lost, nevertheless, we show that the sign pattern of the diagnostic still tracks whether the lowest excitation gap is controlled by boundary modes or by bulk excitations. These results are enabled by an exact single-particle operator Lanczos algorithm, as derived in this work, which reduces the recursion from exponentially large operator space to a finite-dimensional linear problem and achieves machine precision for chains of hundreds of sites. Krylov diagnostics thus emerge as practical probes of boundary-versus-bulk low-energy physics in topological superconductors with broken $U(1)$ symmetry and algebraically decaying couplings.

Maxwell's demon elucidates the value of information in thermodynamics, using measurement and feedback: he evolves an equilibrated gas into a nonequilibrium state, from which one might extract work. The demon can evolve the system farther from equilibrium, on average, if the particles obey Bose-Einstein statistics than if they are distinguishable. We experimentally support this decade-and-a-half-old prediction by comparing indistinguishable with distinguishable photons. We use a fully programmable linear-optics platform, whose local photon statistics were shown recently to behave thermally. Our demon nondestructively measures the number of photons in a subset of the modes. Guided by the outcome, he conditionally interchanges the measured and unmeasured modes. This interchange creates a positive temperature difference between a mode in a particular subset and a mode in the other. The temperature difference is greater, on average, if the photons are indistinguishable. This result bolsters a long-standing prediction about the interplay among thermodynamics, information, and quantum particle statistics. Additionally, this work suggests a thermodynamic means of weakly validating boson-sampling platforms.

Resolving signals with closely spaced frequencies is central to applications in communications, spectroscopy and sensing. Recent results have shown that quantum sensing protocols can exhibit superresolution, the ability to discriminate between spectral lines with arbitrarily small frequency separation. Here, we revisit this problem from the perspective of quantum control theory, utilizing the filter function formalism to derive general, analytic conditions on quantum control protocols for achieving superresolution. Building on these conditions, we develop an optimal control framework, the utility of which is demonstrated through numerical identification of superresolution control protocols in the presence of realistic, experimentally-relevant constraints. We further extend our results to entangled initial states and assess their potential advantage. Our approach is broadly applicable to a wide variety of quantum sensing platforms, and it provides a systematic path to discover novel protocols that surpass conventional resolution limits in these systems.

In this work, we present a quantum algorithm for direct first-principles simulation of electron-nuclear dynamics on a first-quantized real-space grid. Our algorithm achieves best-in-class efficiency for block-encoding the pre-Born-Oppenheimer molecular Hamiltonian by harnessing the linear scaling of swap networks for implementing the quadratic number of particle interactions, while using a novel alternating sign implementation of the Coulomb interaction that exploits highly optimized arithmetic routines. We benchmark our approach for a series of scientifically and industrially relevant chemical reactions. We demonstrate over an order-of-magnitude reduction in costs compared to previous state-of-the-art for the $\rm NH_3+BF_3$ reaction, achieving a Toffoli cost of $8.7\times10^{9}$ per femtosecond using $1362$ logical qubits (system + ancillas). Our results significantly lower the resources required for fault-tolerant simulations of photochemical reactions, while providing a suite of algorithmic primitives that are expected to serve as foundational building blocks for a broader class of quantum algorithms.

We identify fundamental limitations in machine learning by demonstrating that non-trivial mixed-state phases of matter are computationally hard to learn. Focusing on unsupervised learning of distributions, we show that autoregressive neural networks fail to learn global properties of distributions characterized by locally indistinguishable (LI) states. We demonstrate that conditional mutual information (CMI) is a useful diagnostic for LI: we show that for classical distributions, long-range CMI of a state implies a spatially LI partner. By introducing a restricted statistical query model, we prove that nontrivial phases with long-range CMI, such as strong-to-weak spontaneous symmetry breaking phases, are hard to learn. We validate our claims by using recurrent, convolutional, and Transformer neural networks to learn the syndrome and physical distributions of toric/surface code under bit flip noise. Our findings suggest hardness of learning as a diagnostic tool for detecting mixed-state phases and transitions and error-correction thresholds, and they suggest CMI and more generally “non-local Gibbsness” as metrics for how hard a distribution is to learn.

In seminal work (arxiv:quant-ph/9707021) Alexei Kitaev proposed topological quantum computing (arXiv:cond-mat/0010440, arxiv:quant-ph/9707021, arXiv:quant-ph/0001108, arXiv:0707.1889), whereby logic gates of a quantum computer are conducted by creating, braiding and fusing anyonic particles on a two-dimensional plane. Furthermore, he showed the proposal is inherently robust to local perturbations (arXiv:cond-mat/0010440, arxiv:quant-ph/9707021, arXiv:1001.0344, arXiv:1001.4363) when anyons are created as quasiparticle excitations of a topologically ordered lattice model prepared at zero temperature. Over the decades following this proposal there have been considerable technological developments towards the construction of a fault-tolerant quantum computer. Rather than maintaining some target ground state at zero temperature, a modern approach is to actively correct the errors a target state experiences, where we use noisy quantum circuit elements to identify and subsequently correct for deviations from the ideal state. We present an error-correction scheme that enables us to carry out robust universal quantum computation by braiding anyons. We show that our scheme can be carried out on a suitably large device with an arbitrarily small failure rate assuming circuit elements are below some threshold level of local noise. The error-corrected scheme we have developed therefore enables us to carry out fault-tolerant topological quantum computation using modern quantum hardware that is now under development.

Quasi-probability decompositions (QPDs) have proven essential in many quantum algorithms and protocols – one replaces a “difficult” quantum circuit with an ensemble of “easier” circuit variants whose weighted outcomes reproduce any target observable. This, however, inevitably yields an increased configuration variance beyond Born-rule shot noise. We develop a broad framework for accounting for and reducing this variance and prove that stratified sampling – under ideal proportional allocation – results in an unbiased estimator with a variance that is never worse than naïve sampling (with equality only in degenerate cases). Furthermore, we provide a classical dynamic programme to enable stratification on arbitrary product-form QPDs. Numerical simulations of typical QPDs, such as Probabilistic Error Cancellation (PEC) and Probabilistic Angle Interpolation (PAI), demonstrate constant-factor reductions in overall variance (up to $\sim 60$–$80\%$ in an oracle model) and robust $\sim 10\%$ savings in the pessimistic single-shot regime. Our results can be applied immediately to reduce the net sampling cost of practically relevant QPDs that are commonly used in near term and early fault-tolerant algorithms without requiring additional quantum resources.

Bell states and quantum teleportation play crucial roles in quantum information and computation. However, a comprehensive theoretical study of both topics remains to be carried out. This work aims to investigate key algebraic properties of generalized Bell states and explore the topological features of quantum teleportation. First, the basis theorem and the basis group are introduced to show that the extension of a generalized Bell basis by a unitary matrix still forms an orthonormal basis. Then, a twist operator is defined to establish a connection between a generalized multi-qubit Bell state and a tensor product of two-qubit Bell states. In addition, the Temperley–Lieb algebra, the braid group relations, and the Yang–Baxter equation are employed to provide a topological description of generalized Bell states and quantum teleportation. The results demonstrate that our approach not only offers a clear illustration of relevant quantum information protocols but also reveals the topological nature of quantum entanglement and teleportation.

Transmon qubits are a cornerstone of modern superconducting quantum computing platforms. Temporal fluctuations of energy relaxation in these qubits are widely attributed to microscopic two-level systems (TLSs) in device dielectrics and interfaces, yet isolating individual defects typically relies on tuning the qubit or the TLS into resonance. We demonstrate a novel spectroscopy method for fixed-frequency transmons based on multilevel relaxation: repeated preparation of the second excited state and simultaneous $T_1$ extraction of the first and second excited states reveals characteristic correlations in the decay rates of adjacent transitions. From these correlations we identify one or more dominant TLSs and reconstruct their frequency drift over time. Remarkably, we find that TLSs detuned by $\gtrsim 100\,\mathrm{MHz}$ from the qubit transition can still significantly influence relaxation. The proposed method provides a powerful tool for TLS spectroscopy without the need to tune the transmon frequency, either via a flux-tunable inductor or AC-Stark shifts.

This work investigates the Floquet dynamics of electrons and excitons (particle-hole pairs) in a Dirac material referred to as Kekulé-distorted graphene. Specifically, we examine the role played by a high frequency driving electromagnetic field on the tunneling and blocking by a potential barrier on both the charged single particles as well as the neutral composite particles. We demonstrate that the small effective masses of the electron and hole for the energy spectrum of this Kekulé distorted graphene leads to practically almost perfect transmission across a symmetric potential barrier for any angle of incidence of impinging excitons. However, this unexpected Klein paradox for excitons does not hold for the single-particle electrons. The reduced total transmission of electron due to Kekulé distortion is more suppressed due to irradiation. Additionally, we calculate and investigate the exciton binding energy since the quantum tunneling of a bound electron-hole pair across a potential barrier is governed by its mass measured in the center of mass and binding energy of the composite pair. Thus, irradiation with circularly polarized light fundamentally modifies exciton formation, coherence and transport properties, thereby producing unusual topological behaviors. These behaviors are unlike conventional Dirac materials. Possible technical applications of the results arising from our investigation include valleytronics due to the folding of the valleys, thereby making intervalley coupling feasible. Other practical applications include optoelectronics due to Floquet tuning of energy spectrum and transport properties.

Reciprocity is a fundamental symmetry present in many natural phenomena and engineered systems. Distinct situations where this symmetry is broken are typically grouped under the umbrella term "nonreciprocity", colloquially defined by: the action of A on B $\neq$ the action of B on A. In this review, we elucidate what nonreciprocity is by providing an introduction to its most salient classes: nonvariational dynamics, violations of Newton's third law, broken detailed balance, nonreciprocal responses and nonreciprocity of arbitrary linear operators. Next, we point out where to find these manifestations of non-reciprocity, from ensembles of particles with field mediated interactions to synthetic neural networks and open quantum systems. Given this breadth of contexts and the lack of an all-encompassing definition, it makes it all the more intriguing that some general conclusions can be gathered, when distinct definitions of nonreciprocity overlap. We explore what these universal consequences are with a special emphasis on collective phenomena that arise in nonreciprocal many-body systems. The topics covered include nonreciprocal phase transitions and non-normal amplification of noise and perturbations. We conclude with some open questions.

Identifying where quantum models may offer practical benefits in near term quantum machine learning (QML) requires moving beyond isolated algorithmic proposals toward systematic and empirical exploration across models, datasets, and hardware constraints. We introduce MerLin, an open-source framework designed as a discovery engine for photonic and hybrid quantum machine learning. MerLin integrates optimized strong simulation of linear optical circuits into standard PyTorch and scikit learn workflows, enabling end-to-end differentiable training of quantum layers. MerLin is designed around systematic benchmarking and reproducibility. As an initial contribution, we reproduce eighteen state-of-the-art photonic and hybrid QML works spanning kernel methods, reservoir computing, convolutional and recurrent architectures, generative models, and modern training paradigms. These reproductions are released as reusable, modular experiments that can be directly extended and adapted, establishing a shared experimental baseline consistent with empirical benchmarking methodologies widely adopted in modern artificial intelligence. By embedding photonic quantum models within established machine learning ecosystems, MerLin allows practitioners to leverage existing tooling for ablation studies, cross-modality comparisons, and hybrid classical-quantum workflows. The framework already implements hardware-aware features, allowing tests on available quantum hardware while enabling exploration beyond its current capabilities, positioning MerLin as a forward-looking co-design tool linking algorithms, benchmarks, and hardware.

The quantum Mpemba effect is a counterintuitive phenomenon in which a state initially farther from equilibrium relaxes more rapidly than one that starts nearer to equilibrium. In the context of finite-dimensional quantum batteries interacting with an environment, we introduce the notion of an ergotropic Mpemba crossing (EMC), defined by the intersection of ergotropy trajectories during the dynamics. For qubit batteries subjected to amplitude damping noise, we derive a condition for the occurrence of EMC in terms of the relative coherence of the initial states and fully characterize the region of state space that exhibits EMC with respect to a fixed reference state. Interestingly, our analysis reveals that under anisotropic Pauli noise, the emergence of EMC is jointly governed by the coherence and the energy of the initial states. To elucidate the physical origin of EMC, we decompose ergotropy into coherent and incoherent contributions and show that, in qubit systems, the coherent component plays a crucial role for EMC, an observation that strikingly does not extend to three-level batteries, providing the dimensional role in observing EMC. Further, by extending our analysis to non-Markovian environments, we demonstrate that, unlike the Markovian case, non-Markovian dynamics can give rise to multiple Mpemba crossings, with the total number of crossings always being odd. Also, we propose a set-up showing that EMC can be utilized to suppress the ergotropic loss, occurred due to dissipation. Moreover, analyzing the connection between the EMC and the conventional state Mpemba effect reveals that, for qubits, an EMC necessarily entails a state Mpemba crossing, while this correspondence breaks down for qutrits, where EMCs may arise without any state Mpemba crossing.

Instantaneous quantum polynomial quantum circuit Born machines (IQP-QCBMs) have been proposed as quantum generative models with a classically tractable training objective based on the maximum mean discrepancy (MMD) and a potential quantum advantage motivated by sampling-complexity arguments, making them an exciting model worth deeper investigation. While recent works have further proven the universality of a (slightly generalized) model, the next immediate question pertains to its trainability, i.e., whether it suffers from the exponentially vanishing loss gradients, known as the barren plateau issue, preventing effective use, and how regimes of trainability overlap with regimes of possible quantum advantage. Here, we provide significant strides in these directions. To study the trainability at initialization, we analytically derive closed-form expressions for the variances of the partial derivatives of the MMD loss function and provide general upper and lower bounds. With uniform initialization, we show that barren plateaus depend on the generator set and the spectrum of the chosen kernel. We identify regimes in which low-weight-biased kernels avoid exponential gradient suppression in structured topologies. Also, we prove that a small-variance Gaussian initialization ensures polynomial scaling for the gradient under mild conditions. As for the potential quantum advantage, we further argue, based on previous complexity-theoretic arguments, that sparse IQP families can output a probability distribution family that is classically intractable, and that this distribution remains trainable at initialization at least at lower-weight frequencies.

A fundamental challenge in photonics-based deterministic quantum information processing is to realize key transformations on time scales shorter than those of detrimental decoherence and loss mechanisms. This challenge has been addressed through device-focused approaches that aim to increase nonlinear interactions relative to decoherence rates. In this work, we adopt a complementary architecture-focused approach by proposing a recirculating quantum photonic network (RQPN) that minimizes the duration of quantum information processing tasks, thereby reducing the requirements on nonlinear interaction rates. The RQPN consists of a network of all-to-all connected nonlinear cavities with dynamically controlled waveguide couplings, and it processes information by capturing a photonic input state, recirculating photons between the cavities, and releasing a photonic output state. We demonstrate the RQPN's architectural advantage through two examples: first, we show that processing all qubits simultaneously yields faster operations than single- and two-qubit decompositions of the three-qubit Toffoli gate. Second, we demonstrate implementations of a measurement-free correction for single-photon loss, achieving up to seven-fold speedups and significantly improved hardware efficiency relative to state-of-the-art architecture proposals. Our work shows that a single hardware-efficient recirculating architecture substantially reduces the temporal overhead of multi-qubit gates and quantum error correction, thereby lowering the barrier to experimental realizations of deterministic photonic quantum information processing.

We generalize Frenkel's integral formula for traces of operators to operators. The resulting formula holds for bounded self-adjoint positive operators and $p$-Schatten class of compact positive operators.

Finding optimal solutions of conflicting objectives is a daily matter in many industrial applications, with multi-objective optimization trying to find the best solutions to them. The advent of quantum computing has led to researchers wondering if the promised exponential advantage can be obtained for these problems by variational quantum multi-objective optimization (QMOO) algorithm. Here, we improve it by introducing a Pareto Archive and dominated solutions substitution, clearly improving in hyper-volume convergence at additional quantum and classical cost. We propose the use of RMNK-landscapes as a unifying testbed for benchmarking QMOO, as it is common in classical multi-objective field. By devising a generic classical-to-quantum mapping of these landscapes, we perform a numerical hyperparameter tuning of QMOO, significantly enhancing its performance. Finally, we compare QMOO against well-known classical solvers for multi-objective tasks, NSGA-II/III, showing comparable results in small instances. Our results demonstrate that QMOO, when carefully tuned for the task at hand, might be advantageous on harder problems than its classical counterparts.

Recent advances in quantum computing and the increasing availability of quantum hardware have substantially enhanced the practical relevance of quantum approaches to discrete optimization. Among these, the Quadratic Unconstrained Binary Optimization (QUBO) formulation provides a unifying modeling framework for a broad class of $\mathbf{NP}$-hard problems and is naturally suited to quantum computing and quantum-inspired algorithms. Location science, network design, and logistics represent core application domains of discrete optimization, combining high practical impact with substantial computational challenges. In this work, we develop QUBO formulations for several fundamental problems in these domains, including a nonlinear integer formulation of the Discrete Ordered Median Problem (DOMP). Beyond their modeling relevance, these QUBO formulations serve as representative benchmark problems for assessing quantum algorithms and quantum hardware. We further derive a tight bound for the penalty parameter ensuring equivalence between the QUBO formulation and its underlying integer program. Finally, we conduct a comprehensive computational study using QAOA, WS-QAOA, and classical heuristics for QUBO instances of the $p$-Median Problem and the Fixed-Charge Facility Location Problem (FCFLP), and introduce two effective warm-start strategies for WS-QAOA based on its linear programming relaxation.

Open quantum systems are often described by a Lindblad master equation, which relies on a set of approximations, most importantly the rotating-wave approximation which is only valid for weak damping. In the Lindblad setting, dissipative processes are described through jump operators, distinguishing between absorption and emission of photons. This enables the simple identification of emitted photons which provides a straightforward way to obtain the radiation statistics. Outside the rotating-wave limit, the Lindblad approach does not work. Open quantum systems can then be described by, e.g., the quantum Langevin equation. However, in this framework the number of emitted photons is not easily accessible. In this work, we point out how to obtain the photon counting statistics from a quantum Langevin equation and provide an expression for the photon current operator, for arbitrary systems coupled to linear environments. As an example, we employ the method to study the radiation statistics of a damped harmonic oscillator at finite temperature beyond the rotating-wave approximation. We show that even outside the rotating-wave limit, the most important contribution to the radiation statistics can be captured by an effective Lindblad equation, thus extending the range of possible applications of the Lindblad framework.

Quantum Electrodynamics predicts that the vacuum must behave as a nonlinear optical medium: the vacuum optical index should increase when it is stressed by intense electromagnetic fields. The DeLLight (Deflection of Light by Light) project aims to measure it by using intense and ultra-short laser pulses. The experiment uses a Sagnac interferometer to amplify the tiny deflection signal of a low-intensity probe pulse crossing the vacuum refractive-index gradient produced by an external high-intensity pump pulse. The measurement of the amplified signal by a CCD camera requires a high spatial resolution, which is limited by the ultimate quantum noise of the CCD. However, interferometric phase noise induced by the mechanical vibrations of the interferometer is also amplified and degrades spatial resolution. To overcome this, we propose a new method named High-Frequency Phase Noise Suppression (HFPNS), based on the addition of a delayed replica (5 ns) of the probe pulse. The delayed pulse, which is not affected by the pump but is subject to the same vibration noise, enables offline subtraction of correlated phase noise. In this work, we present an experimental proof-of-concept on a prototype interferometer operating with a limited amplification factor ($\mathcal{A}\simeq25$), about 10 times smaller than the required value of the final experiment. We have succeeded in reducing phase noise by a factor of 40, resulting in a residual noise level 2.3 times higher than the expected quantum noise. The residual noise is linked to delay-line instabilities and incident beam pointing fluctuations present during these tests. This result validates HFPNS as a robust method for future quantum-noise-limited interferometric measurements of vacuum optical nonlinearity, though additional stabilization and higher interferometric amplification are still needed.

Quantum technologies are powered by platforms to generate complex non-classical states of matter or light to realize applications. We investigate the non-classical properties of high-harmonic generation in semiconductors, an emerging photonic platform. Measuring the click statistics of three double-digit orders, we evaluate witness operators to certify the non-classicality of the generated states. We show that higher-order harmonics driven by a coherent laser are squeezed and entangled. The properties of the emission are well retrieved with an entangled Gaussian state model, obtained by numerical state optimization to multiple observables. Additionally, we perform inter-order heralded measurements to engineer the quantum state of the emission. The heralded states have distinct properties, showing sub-Poissonian photon statistics. Further, we witness the generation of a quantum non-Gaussian state, a resource highly relevant for quantum information. With this, we establish high-harmonic generation as a platform for generating quantum optical resources.

Controlling ensembles of weakly coupled spins typically requires computationally expensive multispin optimisations. We present a compact framework that enables control of weakly coupled spin systems (of any spin), but using RF pulses optimised for a single spin-1/2. We do this by explicitly creating a GRAPE pulse with fixed 'active' evolution times using single spin-1/2 methods, and pulsing on one spin at a time. By enforcing this form uniformly across offsets ('band-schematic' pulses),chemical shift and scalar coupling evolution of the entire system can be precisely controlled. We demonstrate the approach by constructing band-schematic pulses and a continuously irradiated joint INEPT (JINEPT) that achieves band-selective transfer $I_z \rightarrow 2I_zS_z$. The framework is implemented in the software Seedless, which both rapidly generates such pulses and analyses the schematic form of arbitrary pulses, enabling robust multi-spin control, without multi-spin optimisation.