Papers

The ability to switch and program the spectral response of quantum modes via deterministically located plasmonic nanoantennae presents opportunities for wide spectrum of applications from biosensors to quantum computing. Due to its topology, toroidal nanoantenna (TNA) focuses immense amount of three-dimensional (3D) local electric field by toroidal moment while allowing pre and post positioning around quantum emitters (QEs). Here, within local-response finite difference time domain (FDTD) simulations, we demonstrate high-contrast spectral switching of the radiative decay channel of a dipolar QE coupled to a TNA by introducing effective Lorentzian quantum objects (QOs). At optimized TNA geometries, Fano interference between the broadband plasmonic continuum and narrow quantum transitions of QOs suppresses both radiative and non-radiative decay channels near 850 nm, yielding an observable full switching that traps energy within the hybrid mode instead of re-emitting it. To show the promises of the concept, we further demonstrate systems with multiple QOs where spectral degeneracy enhances the transparency bandwidth, while detuning generates distinct minima, enabling individually addressable spectral responses. These results establish plasmonic TNAs as promising architectures for spectral detection and individual mode switching of single- or multi-QO configurations and empowers the user for the implementation of photonic processing of continuous photon sources.

We study the van der Waals-London force, which is typically associated with fluctuating electric dipoles in atoms, in a mesoscopic circuit consisting of two inductively coupled superconducting loops. We investigate the inductive van der Waals-London interaction using both semiclassical and quantum electrodynamic (QED) approaches. The semiclassical model predicts a repulsive interaction due to anticorrelated current fluctuations. In contrast, the QED framework, which incorporates virtual photon exchange, reveals a predominantly attractive force. A key contribution comes from a state-independent two-photon exchange, which is absent in the semiclassical description and undetectable by spectroscopy. Our study introduces a theoretical framework for exploring the van der Waals force between individual artificial atoms via controlled mesoscopic circuits.

We propose an experimentally accessible scheme for realizing tunable nonclassical light in cavity-coupled reconfigurable atomic arrays. By coherently controlling the collective interference phase, the system switches from single-photon blockade to high-purity multiphoton bundle emission, unveiling a hierarchical structure of photon correlations dictated by atom-number parity and cavity detuning. The scaling of photon population identifies the transition between superradiant and subradiant regimes, while parity- and phase-dependent spin correlations elucidate the microscopic interference processes enabling coherent multiphoton generation. This work establishes a unified framework connecting cooperative atomic interactions to controllable nonclassical photon statistics and introduces a distinct interference-enabled mechanism that provides a practical route toward high-fidelity multiphoton sources in scalable cavity QEDs.

We consider a problem of geometric phase generation in a system of two interacting bosons confined in a narrow ring potential with a localized defect. Geometric phase emerges from variation of parameters of the defect. Particle interaction is taken into account within a framework of the Lieb-Liniger model. The energy spectrum is evaluated and its dependence on the parameters of the problem is described. It is shown that the interaction leads to increase of the geometric phase for a given contour of variations. The work is motivated by earlier proposed ideas of quantum gyroscope and quantum accelerometer based on atomic Bose-Einstein condensates.

A laser beam is often modelled by a pure coherent state. In fact its state is mixed, even if it has coherent-state photon-number statistics (Poissonian), because the phase must vary. We consider such an ideal laser beam, with phase diffusion rate $\ell$, equal to its (Lorentzian) spectral width. We show that the beam entropy is extensive, with an entropy flow of $\dot{S} = \kB\sqrt{\dot{N}\ell}$, where $\dot{N}$ is the number flow. We give an intuitive explanation for this remarkably simple result, and compare it to the entropy flow of a unidirectional thermal beam.

Maintaining stable entanglement distribution through perturbed fiber links is essential for practical quantum-optics experiments, yet it remains challenging because of polarization fluctuations and phase or temporal-delay variations. We demonstrate stable entangled-photon transmission using a cross-aligned polarization-maintaining fiber (CAPMF) structure composed of two polarization-maintaining fiber sections with mutually orthogonal principal axes. The CAPMF configuration passively compensates polarization fluctuations without real-time active polarization control. We theoretically analyze the CAPMF structure and experimentally verify its stabilization performance under external mechanical perturbations. In the experiment, the single-mode fiber configuration yields an average visibility of $0.7655$ and a CHSH value of $S=1.7714$, whereas the CAPMF configuration maintains an average visibility of $0.9843$ and a CHSH value of $S=2.6838$. These results show that CAPMF offers a simple and robust architecture for stabilizing fiber-interface sections in practical entanglement-distribution systems.

The exceptional point has presented considerably interesting and counterintuitive phenomena associated with nonreciprocity, precision measurement, and topological dynamics. The Liouvillian exceptional point (LEP), involving the interplay of energy loss and decoherence inherently relevant to quantum jumps, has recently drawn much attention due to capability to fully capture quantum system dynamics and naturally facilitate non-Hermitian quantum investigations. It was also predicted that quantum jumps could give rise to third-order LEPs in two-level quantum systems for its high dimensional Liouvillian superoperator, which, however, has never been experimentally confirmed until now. Here we report the first observation of the third-order LEPs emerging from quantum jumps in an ultracold two-level trapped-ion system. Moreover, by combining decay with dephasing processes, we present the first experimental exploration of LEPs involving combinatorial effect of decay and dephasing. In particular, due to non-commutativity between the Lindblad superoperators governing LEPs for decay and dephasing, we witness the movement of LEPs driven by the competition between decay and dephasing occurring in an open quantum system. This unique feature of non-Hermitian quantum systems paves a new avenue for modifying nonreciprocity, enhancing precision measurement, and manipulating topological dynamics by tuning the LEPs.

The determination of the quantum chromatic number of graphs has attracted considerable attention recently. However, there are few families of graphs whose quantum chromatic numbers are determined. A notable exception is the family of orthogonality graphs, whose quantum chromatic numbers are fully determined. In this paper, we extend these results by determining the exact quantum chromatic number of several subgraphs of the orthogonality graphs. Using the technique of combinatorial designs, we also determine the quantum chromatic number of the distance-2 Hamming graph, whose edges consist of binary vectors of Hamming distance 2, for infinitely many length.

Decision diagram (DD)-based quantum circuit simulators represent quantum states and gates using DDs, enabling memory-efficient and fast simulations for some quantum circuits like Shor. Although it is known that DD size and processing time vary depending on the variable order in classical circuits, there has not been much research on the variable order under quantum circuit simulation. One existing study pointed out that dynamic reordering worsens the simulation time and numerical accuracy, and there is no comprehensive research on static orders in the context of quantum circuit simulation. Therefore, this paper proposes a scoring-based heuristic method for determining a static variable order that enables efficient DD-based quantum circuit simulation. When applied to benchmark circuits, the default original variable orders resulted in slow simulations, whereas the proposed method achieved speedups of up to 150x. Furthermore, the proposed order completed the simulation of Shor's 1011 factorization in 5 hours on a single-core laptop, although it was not completed within two days previously.

A Mn dopant in a III-V semiconductor produces a highly-entangled, coherent triplet ground state not fully captured by single-determinant theories of electron structure. We directly construct an analytic form for its ground-state wavefunction, finding surprising spin-charge correlations not revealed by semiclassical calculations. Spin-correlated circulating currents associated with the dopant yield remarkably large magnetic fringe fields of $\sim$1$\,\mu$T at distances of $\sim 10$~nm from Mn in GaAs, potentially detectable by NV-diamond magnetometry while the dopant spin coherently precesses.

In this article, a new formula for computing Cayley's first hyperdeterminant in terms of the Levi-Civita symbol is given. It is then shown that this formula can be used to compute the hyperdeterminant of symmetric tensors in polynomial time with respect to their order (assuming fixed side length). Applications to quantifying the entanglement of states of bosonic quantum systems are then discussed. Additionally, in order to obtain the fast calculation of the hyperdeterminant on symmetric tensors, generalized elimination and duplication matrices are defined and their explicit formulas are derived.

Barren plateaus, which means the training gradients become extremely small, pose a major challenge in optimizing parameterized quantum circuits, often making the learning process impractically slow or stall. This work shows why using neural networks to generate quantum circuit parameters helps overcome this difficulty. We introduce a geometric viewpoint that describes how the parameters produced by neural networks evolve during training. Our analysis shows that these parameters follow smooth and efficient paths that avoid the flat regions in the training that cause barren plateaus. This provides a computational explanation for the improved trainability observed in recent neural network-assisted quantum learning methods. Overall, our findings bridge ideas from quantum machine learning and computational optimization, offering new insight into the structure of quantum models and guiding future approaches for designing more trainable quantum circuits or parameter initialization.

The projected rapid growth of battery cell production over the next decade demands advanced diagnostic tools for quality control, ageing prediction, and recycling. Most existing techniques lack the spatial and temporal resolution required to capture internal electrochemical processes non-invasively. Here, we present magnetic imaging of current densities in battery cells, a sensitive quantum-magnetometry method that uses optically pumped magnetometers (OPMs) to perform real-time imaging of internal dynamics in open-circuit configuration. We demonstrate this approach for monitoring relaxation processes in 6000 mA h lithium-ion cells following pulsed discharges across a range of pulse durations and currents as well as states of charge. The measurement results are benchmarked against superconducting-quantum-interference-device (SQUID) magnetometry and validated with three-dimensional finite element simulations. Equivalent circuit models are employed to interpret the relaxation profiles, revealing spatially resolved features and transient magnetic-field signatures that are inaccessible with complementary non-invasive techniques such as electrochemical impedance spectroscopy (EIS). This work establishes OPM-based magnetic imaging of battery current density as a powerful diagnostic tool with potential impact on cell development, manufacturing quality assurance, and second-life assessment.

We present a variational quantum algorithm that solves the one-dimensional diffusion problem with a space-dependent diffusion constant $D(x)$. This problem is relevant for the exchange of hydroxide ions across a two-layer membrane in an alkaline electrolyzer, where the concentration of OH$^-$ ion determines the chemical stability for longer time periods. We use $16$ to $64$ grid points across the membrane, resulting from $n=4$ to $6$ data qubits for the ideal statevector and shot-based quantum simulations implemented using Qiskit. For these qubit numbers, the depth of the parametric quantum circuit has been chosen to ensure sufficient expressibility. The state preparation requires particular attention since the diffusivity $D$ is piecewise constant in the different layers with discontinuities at the interface. Furthermore, we compare different classical optimization schemes with respect to their convergence in the VQA method. We demonstrate the applicability of the quantum algorithm to a problem with non-trivial boundary conditions and jump conditions of the diffusion constant and outline possible extensions of the proof-of-concept application case of quantum computing. Our simulations show that pronounced hydroxide ion concentration gradients, and thus chemical instabilities, can occur only when the ratio of diffusivity in both layers of the membrane exceeds approximately 50.

The phenomenon of Entangled Two-Photon Absorption (ETPA) presents a persistent controversy in the literature, evidenced by a wide disparity in the reported values for the $\sigma_E$ cross-sections. Much of this discrepancy is attributed to the difficulty in discriminating ETPA from various background processes that can mimic its signal, such as linear absorption or scattering. Given this need to certify the presence of ETPA unequivocally, this work introduces a key strategy to isolate the ETPA contribution through its spectral signature. This involves modeling the molecule as a two-photon notch filter and inducing a controlled asymmetric overlap with the joint spectral intensity (JSI) of the incident photons. This asymmetry is used to generate a measurable distortion in the transmitted JSI, and, complementarily, as a reduction in the visibility of the Hong-Ou-Mandel (HOM) dip. To ensure the experimental feasibility of this technique, a comprehensive analysis of the experimental conditions required for its detection is presented, establishing the absorption efficiency limits that must be overcome given the constraints of entangled-photon-pair flux and detector noise. A preliminary experiment addressing the ETPA detection limits using the standard RhB dye is presented.

We investigate quantum features and non-classical nature of two-spin-$1/2$ NMR systems at thermal equilibrium under external magnetic fields. More specifically, using suitable quantifiers, we analyze quantum coherence, mixedness, and entanglement in NMR systems and examine their features within the system. We derive closed-form analytical expressions for the quantum elements and show how they depend on temperature and magnetic field strength. We demonstrate that at zero temperature, the system exhibits a quantum critical point, characterized by non-analytic behavior in the measures of coherence, and a sharp peak in mixedness. Moreover, we analyze the entanglement of the system using a suitable entanglement witness. This provides an experimentally friendly setting for testing entanglement in NMR systems. In other words, the witness links the entanglement in the system to quantum observables, making it directly provable in NMR experiments. We establish a connection between quantum information quantifiers and experimentally accessible NMR spectra of the system, enabling the quantification of entanglement, coherence, and mixedness through NMR signal processing.

Quantum many-body (QMB) systems are generally computationally hard: the computing resources necessary to simulate them exactly can often exceed the existing computation resources by orders of magnitude. For this reason, Richard Feynman proposed the concept of a quantum simulator: quantum systems engineered to obey a prescribed evolution equation and repeating the experiment multiple times. Experimentally, however, as we explain below, the vast majority of observables describing the system are inaccessible. Thus, while Feynman's idea addresses the problem of simulating quantum dynamics, it leaves unsolved the equally fundamental problem of inferring the underlying physics from the limited observables accessible in experiments. Indeed, many complex phenomena associated with QMB systems remain elusive. Perhaps, the most important example is identifying phase transitions in QMB systems when no simple order-parameter exists, which poses major challenges to this day. Complicating the problem further is the fact that, in most cases, it is impossible to learn from numerical simulations, as the underlying systems are often too large to be computable, and small QMB can show strong finite size effects, masking the presence of the transition. Here, we present an unsupervised machine learning approach to study QMB experiments, specifically aimed at detecting phase transitions and crossovers directly from raw experimental measurements. We demonstrate our methodology on systems undergoing Many-Body Localization cross-over and Mott-to-Superfluid phase-transition, showing that it reveals collective phenomena from the very partial experimental data and without any model-specific prior knowledge of the system. This approach offers a general and scalable route for data-driven discovery of emergent phenomena in complex quantum many-body systems.

We investigate quantum entanglement in two-spin-1/2 NMR systems at thermal equilibrium under external magnetic fields. We derive closed-form analytical expressions for the entanglement of the system and show how the entanglement depends on temperature and magnetic field strength, resulting in a threshold temperature beyond which entanglement vanishes. We demonstrate that at zero temperature, the system exhibits a quantum critical point, characterized by non-analytic behavior in the measure of entanglement. We further develop analytical criterion for level crossing, which serves as a condition for identifying quantum critical points in both homonuclear and heteronuclear systems, and apply it to multiple settings to analyze their quantum critical points. We establish a direct link between the quantum entanglement quantifier and experimentally accessible NMR observables, enabling entanglement to be quantified through NMR signal processing. This provides a practical framework for characterizing quantum correlations using standard NMR experiments. These findings provide insights into the thermal control of quantum features, with implications for quantum-enhanced NMR, low-temperature spectroscopy, and emerging quantum technologies.

The iconic problem of photon modes in a spherical cavity has been discussed in the literature; however, conflicting results have been reported \cite{Heitler,Davydov_QuantumMechanics}. For this reason, the solution of this problem is worked out in detail here, starting with the Maxwell equations and applying boundary conditions at the surface of the bounding perfect conductor. Contrary to the treatments in the literature \cite{Heitler,Davydov_QuantumMechanics}, the allowed frequencies for photons in the sphere are given by two different conditions, one for electric and one for magnetic multipole photons. After establishing the modes and their allowed frequencies, we write down the second-quantized vector potential in the spherical geometry. Based on these spherical modes, bipartite photon entanglement is investigated showing that there are in-principle 40 different types of entangled photon states. Finally, we include some appendices about photon plane-wave and spherical-wave helicity states, helicity spherical harmonic vectors, and rotation of the helicity states and 3-d vectors using the Wigner $D$-matrix.

Emission and absorption spectra of molecular films are impacted by low-frequency molecular vibrations. These vibrations define the linewidths of the absorption and emission spectral peaks, as well as the Stokes shift. In cavities that use a molecular film as an active medium, low-frequency molecular vibrations facilitate the thermalization of light, enabling the formation of Bose-Einstein condensation. In this work, I employ perturbation theory for Lindblad superoperators and derive the contribution of the low-frequency molecular vibrations to the thermalization rate of light in a weak coupling regime between light and matter. The derived thermalization rate applies for any cavity design but depends on the local microscopic properties of low-frequency molecular vibrations. I provide an estimation for the thermalization rate, which requires only knowledge of the macroscopic parameters of the system: light-matter interaction strength, resonant frequencies of the cavity and excitons, number of molecules in the illuminated area, and the linewidth temperature dependence of the 0-0 peak in the emission spectra of standalone molecular film.

We provide a systematic evaluation of the sample-based quantum diagonalization (SQD) method for electronic structure based on the W4-11 thermochemistry dataset, comprising 124 total atomization, 83 bond dissociation, 20 isomerization, 505 heavy-atom transfer, and 13 nucleophilic substitution processes, covering diverse bonding situations and reaction mechanisms. This is the largest study assessing the accuracy and precision of a quantum-hybrid algorithm on a digital quantum device across a variety of molecular systems and chemical reactions, using 16.85 hours on the superconducting quantum processor ibm_rensselaer and 724.22 node hours on the supercomputer AiMOS. To ensure a fair comparison, our study employs commensurate resource allocation for both classical and quantum simulations. Although SQD exhibits large statistical deviations from ground-state reference energies, energy extrapolations yield CCSD-level accuracy. While bond-breaking reactions show a systematic improvement as computational resources increase, nucleophilic substitution or heavy atom transfer reactions do not. The limitations quantified in this manuscript indicate opportunities for improvement in SQD-based algorithms. This work provides a benchmark and community resource for exploring new quantum algorithms and devices, supported by an online benchmark challenge and an open-source Python library for direct comparison.

We use the Gaussian Phase-Space Representation to solve the real-time dynamic of interacting fermions in 1D, 2D, and 3D systems. The method is exact up to a spiking point, which represents a limit on the practical simulation time. The spiking can be delayed, and the practical simulation time extended, by adjusting the gauges of the representation, resulting in different equivalent stochastic differential equations. Here, we work on the so-called diffusion gauge and propose an algorithm to find efficiently new implementations of the noise terms. Compared with our initial results [F. Rousse et al. 2024, J. Phys. A: Math. Theor. 57, 015303], the new method achieves a significantly longer practical simulation time and can be applied to significantly larger systems.

Quantum information scrambling has emerged as a powerful tool for studying the dynamics of chaotic quantum many-body systems, assessing benchmarking protocols, and even investigating exotic black hole models. During quantum information scrambling, localized quantum information disperses across the entire system, hiding from the observers who can only access part of it. On the other side, we have a fundamental cryptographic primitive called secret-sharing schemes, where a dealer shares a quantum secret with a group of parties, such that any subset of parties above a specific threshold size can reconstruct it. In this paper, we demonstrate that two protocols, Haar scrambling, which serves as a baseline for other scrambling techniques, and quantum secret sharing, are indeed equivalent. The scheme one gets out of this is of a ramp secret-sharing nature rather than a threshold scheme. We expect this because of the inherent randomness present in Haar unitaries. Moreover, by varying the purity of the initial states, we demonstrate that it is possible to obtain all possible types of ramp secret-sharing schemes. Furthermore, utilizing complexity theoretic arguments, we argue that the protocol can be implemented efficiently, i.e., in polynomial time. Finally, we explore potential applications, ranging from security implications for distributed quantum networks to cryptography in the Noisy Intermediate-Scale Quantum (NISQ) era, and draw insights for fundamental physics, such as many-body physics and quantum black holes.

Quantum thermodynamics has emerged as a central field for understanding how energy conversion processes occur in microscopic systems. In these systems, effects such as coherence, entanglement, and non-Markovianity play key roles. In this thesis, we explore different ways to describe quantum thermodynamics, using two main approaches: one based on entropy and the other on ergotropy. First, we introduce a generalized approach to quantify non-Markovianity through the breakdown of monotonicity in thermodynamic functions. In this context, the entropy-based heat flow serves as a practical tool to witness and measure quantum memory in unital maps that do not reverse the sign of the internal energy. Next, we analyze the dynamics of ergotropy in open qubits under both Markovian and non-Markovian evolutions. We identify phenomena such as freezing and sudden death of ergotropy, and we establish an analytical relation between the change in ergotropy and the environment-induced work . This provides a clear physical interpretation for the additional term in the first law in the entropy-based formulation. Finally, we propose an ergotropy-based thermodynamic formulation, in which heat is reinterpreted in terms of the change of the passive state associated with the density operator governing the quantum dynamics. This approach allows one to measure non-Markovianity of unital maps more generally and accurately than the entropy-based heat flow. This advantage comes from the direct link between heat and von Neumann entropy, a property ensured by invariance under passive transformations. Moreover, the out-of-equilibrium temperature naturally remains non-negative, similarly to equilibrium thermodynamics.

Characterizing non-Markovian quantum dynamics is currently hindered by the self-inconsistency and high computational complexity of existing quantum comb tomography (QCT) methods. In this work, we propose a self-consistent framework that unifies the quantum comb, instrument set, and initial states into a single geometric entity, termed as the Comb-Instrument-State (CIS) set. We demonstrate that the CIS set naturally resides on a product Stiefel manifold, allowing the tomography problem to be solved via efficient unconstrained Riemannian optimization while automatically preserving physical constraints. Numerical simulations confirm that our approach is computationally scalable and robust against gate definition errors, significantly outperforming conventional isometry-based QCT methods. Our work indicates the potential to efficiently learn quantum comb with fewer computational resources.

In an era where data underpins decision-making across science, politics, and economics, ensuring high data quality is of paramount importance. Conventional computing algorithms for enhancing data quality, including anomaly detection, demand substantial computational resources, lengthy processing times, and extensive training datasets. This work aims to explore the potential advantages of quantum computing for enhancing data quality, with a particular focus on detection. We begin by examining quantum techniques that could replace key subroutines in conventional anomaly detection frameworks to mitigate their computational intensity. We then provide practical demonstrations of quantum-based anomaly detection methods, highlighting their capabilities. We present a technical implementation for detecting volatility regime changes in stock market data using quantum reservoir computing, which is a special type of quantum machine learning model. The experimental results indicate that quantum-based embeddings are a competitive alternative to classical ones in this particular example. Finally, we identify unresolved challenges and limitations in applying quantum computing to data quality tasks. Our findings open up new avenues for innovative research and commercial applications that aim to advance data quality through quantum technologies.

Multiqubit gates that involve three or more qubits are usually thought to be of little significance for fault-tolerant quantum error correction because single gate faults can lead to errors of high Pauli weight. However, recent works have shown that multiqubit gates can be beneficial for measurement-free fault-tolerant quantum error correction and for fault-tolerant stabilizer readout in unrotated surface codes. In this work, we investigate multiqubit Rydberg gates that are useful for fault-tolerant quantum error correction in single-species neutral-atom platforms and can be implemented with global laser pulses that do not individually address atomic sites. We develop an open-source Python package to generate analytical, few-parameter pulses that implement the desired gates while minimizing gate errors due to Rydberg-state decay. The tool also allows us to identify parameter-optimal pulses, characterized by a minimal parameter count for the pulse ansatz. Measurement-free quantum error correction protocols require CCZ gates, which we analyze for atoms arranged in symmetric and asymmetric configurations. We investigate the performance of these schemes for various single-, two-, and three-qubit gate error rates, showing that break-even performance of measurement-free QEC is within reach of current hardware. Moreover, we study Floquet quantum error correction protocols that comprise two-body stabilizer measurements. Those can be realized using global three-qubit gates, and we show that this can lead to a significant reduction in shuttling operations. Simulations with realistic circuit-level noise indicate that applying three-qubit gates for stabilizer measurements in Floquet codes can yield competitive logical qubit performance in experimentally relevant error regimes.

The emergence of quantum chaos in a system of trapped interacting bosons with externally impressed rotation is studied through spectral form factor (SFF) and power spectrum using exact diagonalization. Two distinct interaction regimes are considered: the moderate, when the interaction energy is small compared to the trap energy and the strong, when the interaction energy is comparable to the trap energy. In the moderate interaction regime, the SFF for the non-rotating case exhibits a dip-plateau structure with absence of linear ramp, indicating integrable behavior, while for the single-vortex state the SFF exhibits a discernible linear ramp consistent with pseudo-integrable behavior. In the strong interaction regime, the non-rotating case exhibits emergence of a linear ramp with small time span in SFF, indicating that the system has moved further towards chaotic regime but continues to be pseudo-integrable. For the single-vortex and the multi-vortex states in strong interaction regime, the span of the linear ramp in SFF increases progressively with rotation, indicating the system has moved into strong chaotic regime consistent with Gaussian orthogonal ensemble. The power spectrum results with exponent lying in the interval $1 \lesssim \alpha \lesssim 2$ are consistent with the findings of SFF. An understanding of the observed crossover from integrable to quantum chaos is presented in terms of the macroscopic occupation of a single-particle quantum state—the Bose-Einstein condensation—and its depletion driven by interaction and rotation.

This work studies the relationship between parametric amplification (or particle creation), adiabaticity and irreversibility in the non-quasi-static regime of a time-dependent quantum harmonic oscillator (TDHO) that evolves unitarily. We provide analytical results for the evolution of the TDHO valid for any functional value of the frequency, which enables us to monitor the behavior of the thermodynamical magnitudes in the non-quasi-static regime. In the latter, the largest modes of the energy eigenstates commonly undergo a process of spontaneous thermalization, where the concept of temperature naturally arises from the unitary evolution of the oscillator, i.e. without relation to any external source of temperature or thermal bath. As the evolution is unitary, this thermalization process can be reversible. We extend the standard definitions of quantum heat and work to account for the change in the populations of the energy levels and relate it to the diagonal entropy to obtain an effective temperature from the associated second principle of thermodynamics. The qualitative behavior of this temperature is similar to the temperature of the largest modes and is only significant in the non quasi-static region, returning to the original value $T\rightarrow 0$ afterwards. This effect is briefly analyzed in the case of a periodically driven quantum system, where the energy eigenstate mixing and the temperature grow and return to the original value periodically.

We introduce color Heisenberg-Lie (super)algebras graded by the abelian groups $Z_3^2$, $Z_2^p\times Z_3^2$ for $p=1,2,3$, and investigate the properties of their associated multi-particle quantum paraoscillators. In the Rittenberg-Wyler's color Lie (super)algebras framework the above abelian groups are the simplest ones which induce mixed brackets interpolating commutators and anticommutators. These mixed brackets allow to accommodate two types of parastatistics: one based on the permutation group (beyond bosons and fermions in any space dimension) and an anyonic parastatistics based on the braid group. In both such cases the two broad classes of paraparticles are given by parabosons and parafermions. Mixed-bracket parafermions are created by nilpotent operators; they satisfy a generalized Pauli exclusion principle leading to roots-of-unity truncations in their multi-particle energy spectrum (braided Majorana qubits and their Gentile-type parastatistics are recovered in this color Lie superalgebra setting). Mixed-bracket parabosons do not admit truncations of the spectrum; the minimal detectable signature of their parastatistics is encoded in the measurable probability density of two indistinguishable parabosonic oscillators in a given energy eigenstate.