Francisco José Peña
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Research Interests
Francisco José Peña has a total of 59 co-authors. Here are the top co-authors by number of shared publications:
- Patricio Vargas
Universidad Tecnica Federico Santa Maria
38 shared publications - Oscar Andrés Negrete
Departamento de Física, Universidad Técnica Federico Santa María, Valparaíso 2340000, Chile
15 shared publications - Bastian Castorene
Pontificia Universidad Católica de Valparaíso
9 shared publications - Eugenio E. Vogel
Universidad de La Frontera
7 shared publications - Nathan M. Myers
Pacific Northwest National Laboratory
6 shared publications - Michel Aguilera
Instituto de Física, Pontificia Universidad Católica de Valparaíso, Casilla 4950, Valparaíso 2373223, Chile
5 shared publications
Latest Publications
- 7 October 2026
- Published by American Physical Society (APS) in Physical Review E
- Vol. 114 (4) , L042101
Thermodynamic signatures of spectral compression in weakly non-Hermitian Dirac fermions
Non-Hermitian band descriptions capture how loss, gain, and environmental coupling reshape quantum matter, yet most experimental probes remain wave based or dynamical. Here we develop an equilibrium charge-response route to spectral compression in Dirac materials. In the real-spectrum regime, invariance of the particle-number operator under the similarity transformation makes the quantum capacitance exactly that of the Hermitian partner, so the whole non-Hermitian content is carried by the Petermann factor. In a minimal nonreciprocal graphene model, hopping imbalance suppresses the Dirac velocity, enhancing the low-energy density of states and the capacitance as an exceptional point is approached. At charge neutrality the capacitance stays linear in temperature, with its slope enhanced relative to the Hermitian value. A perpendicular magnetic field recasts the same compression as a collapse of the Landau-level ladder, drawing more levels into the thermal window. In two dimensions, the capacitance enhancement relative to the Hermitian partner coincides exactly with the Petermann factor, even though the two have entirely different origins. This relation provides a controlled starting point for extending equilibrium thermodynamic probes to non-Hermitian electronic systems with interactions, disorder, and reservoir coupling.
- 7 October 2026
- Published by American Physical Society (APS) in Physical Review E
- Vol. 114 (4) , L042101
Thermodynamic signatures of spectral compression in weakly non-Hermitian Dirac fermions
Non-Hermitian band descriptions capture how loss, gain, and environmental coupling reshape quantum matter, yet most experimental probes remain wave based or dynamical. Here we develop an equilibrium charge-response route to spectral compression in Dirac materials. In the real-spectrum regime, invariance of the particle-number operator under the similarity transformation makes the quantum capacitance exactly that of the Hermitian partner, so the whole non-Hermitian content is carried by the Petermann factor. In a minimal nonreciprocal graphene model, hopping imbalance suppresses the Dirac velocity, enhancing the low-energy density of states and the capacitance as an exceptional point is approached. At charge neutrality the capacitance stays linear in temperature, with its slope enhanced relative to the Hermitian value. A perpendicular magnetic field recasts the same compression as a collapse of the Landau-level ladder, drawing more levels into the thermal window. In two dimensions, the capacitance enhancement relative to the Hermitian partner coincides exactly with the Petermann factor, even though the two have entirely different origins. This relation provides a controlled starting point for extending equilibrium thermodynamic probes to non-Hermitian electronic systems with interactions, disorder, and reservoir coupling.
- Open Access
- 17 September 2026
- Published by Wiley in Advanced Quantum Technologies
- Vol. 9 (9) , e70470
Fundamental Work Scaling and Non‐Extensivity in Critical Engines
We present a general analytical framework for critical two‐isothermal engines that emerge operationally from quasi‐static quantum Stirling cycles across ground‐state level crossings (GLCs) in the low‐temperature regime, where reversible heat exchange is governed by the structural entropy change while the internal energy remains constant. Within this ideal‐reservoir, low‐temperature equilibrium model, the Primarch Formula gives an exact quasistatic ensemble‐average expression linking extracted work and efficiency directly to macroscopic ground‐state degeneracies. In this reversible setting, the engines reach Carnot efficiency without a classical regenerator. For the perturbative population pattern analyzed here, thermal excitations reduce the ideal work and efficiency. Validated against exact numerical simulations of generalized N ‐th spin‐s Heisenberg models with nontrivial interactions, the framework is applied to the one‐dimensional antiferromagnetic Ising model, revealing a profound connection to number theory. Governed by Fibonacci‐Lucas and parity‐dependent critical degeneracies, the engine exhibits three distinct work‐scaling regimes: persistent logarithmic non‐extensivity, recovery of extensivity in the thermodynamic limit, and asymptotic extensivity with a logarithmic finite‐size correction. In all three regimes, Carnot efficiency is attained within the reversible equilibrium model.
- 17 September 2026
- Published by Wiley in Advanced Quantum Technologies
- Vol. 9 (9) , e70470
Open AccessFundamental Work Scaling and Non‐Extensivity in Critical Engines
We present a general analytical framework for critical two‐isothermal engines that emerge operationally from quasi‐static quantum Stirling cycles across ground‐state level crossings (GLCs) in the low‐temperature regime, where reversible heat exchange is governed by the structural entropy change while the internal energy remains constant. Within this ideal‐reservoir, low‐temperature equilibrium model, the Primarch Formula gives an exact quasistatic ensemble‐average expression linking extracted work and efficiency directly to macroscopic ground‐state degeneracies. In this reversible setting, the engines reach Carnot efficiency without a classical regenerator. For the perturbative population pattern analyzed here, thermal excitations reduce the ideal work and efficiency. Validated against exact numerical simulations of generalized N ‐th spin‐s Heisenberg models with nontrivial interactions, the framework is applied to the one‐dimensional antiferromagnetic Ising model, revealing a profound connection to number theory. Governed by Fibonacci‐Lucas and parity‐dependent critical degeneracies, the engine exhibits three distinct work‐scaling regimes: persistent logarithmic non‐extensivity, recovery of extensivity in the thermodynamic limit, and asymptotic extensivity with a logarithmic finite‐size correction. In all three regimes, Carnot efficiency is attained within the reversible equilibrium model.
- Open Access
Exact Combinatorial Density of States for the Critical 1D Ising Model
This work presents an exact microcanonical combinatorial analysis of the one-dimensional antiferromagnetic Ising model. At the primary ground-state level crossing , degeneracies follow the Fibonacci and Lucas sequences for open chains and periodic rings, respectively. We extend this framework to the complete excitation spectrum, demonstrating that the density of states is constructed from topological defects governed by linear Diophantine equations and p-fold Fibonacci convolutions. Open boundaries act as fractional defects, densifying the chain spectrum into energy steps of , whereas the closed ring remains quantized in units of . Notably, this exact topological counting exposes non-trivial spectral gaps near the fully polarized limit, strictly forbidding the penultimate macroscopic energy levels in both topologies. Using the transfer-matrix formalism, we derive exact closed-form expressions for the critical degeneracies at all energy levels. These results provide a rigorous analytical foundation for extracting exact residual entropies and exposing the intrinsic number-theoretic architecture of quantum critical manifolds.Open AccessExact Combinatorial Density of States for the Critical 1D Ising Model
This work presents an exact microcanonical combinatorial analysis of the one-dimensional antiferromagnetic Ising model. At the primary ground-state level crossing , degeneracies follow the Fibonacci and Lucas sequences for open chains and periodic rings, respectively. We extend this framework to the complete excitation spectrum, demonstrating that the density of states is constructed from topological defects governed by linear Diophantine equations and p-fold Fibonacci convolutions. Open boundaries act as fractional defects, densifying the chain spectrum into energy steps of , whereas the closed ring remains quantized in units of . Notably, this exact topological counting exposes non-trivial spectral gaps near the fully polarized limit, strictly forbidding the penultimate macroscopic energy levels in both topologies. Using the transfer-matrix formalism, we derive exact closed-form expressions for the critical degeneracies at all energy levels. These results provide a rigorous analytical foundation for extracting exact residual entropies and exposing the intrinsic number-theoretic architecture of quantum critical manifolds. - 17 April 2026
- Published by American Physical Society (APS) in Physical Review E
- Vol. 113 (4) , 044126
Reaching maximum efficiency in quantum Stirling engines using multilayer graphene
In this work, quantum Stirling engines based on monolayer, AB-stacked bilayer, and ABC-stacked trilayer graphene under perpendicular magnetic fields are analyzed. Performance maps of the useful work (ηW) reveal a robust optimum at low magnetic fields and moderately low temperatures, with all stackings capable of reaching Carnot efficiency under suitable configurations. The AB bilayer achieves this across the broadest parameter window while sustaining finite work, the monolayer exhibits highly constrained regimes, and the trilayer shows smoother trends with sizable ηW. These results identify multilayer graphene, particularly the AB bilayer, as a promising platform for efficient Stirling engines, while also highlighting the versatility of the monolayer in realizing all four operational regimes of the Stirling cycle.
- 17 April 2026
- Published by American Physical Society (APS) in Physical Review E
- Vol. 113 (4) , 044126
Reaching maximum efficiency in quantum Stirling engines using multilayer graphene
In this work, quantum Stirling engines based on monolayer, AB-stacked bilayer, and ABC-stacked trilayer graphene under perpendicular magnetic fields are analyzed. Performance maps of the useful work (ηW) reveal a robust optimum at low magnetic fields and moderately low temperatures, with all stackings capable of reaching Carnot efficiency under suitable configurations. The AB bilayer achieves this across the broadest parameter window while sustaining finite work, the monolayer exhibits highly constrained regimes, and the trilayer shows smoother trends with sizable ηW. These results identify multilayer graphene, particularly the AB bilayer, as a promising platform for efficient Stirling engines, while also highlighting the versatility of the monolayer in realizing all four operational regimes of the Stirling cycle.
- Open Access
Finite-Size Thermodynamics of the Two-Dimensional Dipolar Q-Clock Model
We present a fully controlled thermodynamic study of the two-dimensional dipolar Q-state clock model on small square lattices with free boundaries, combining exhaustive state enumeration with noise-free evaluation of canonical observables. We resolve the complete energy spectra and degeneracies for the Ising case () on lattices of size , and for clock symmetries on a lattice, tracking how the competition between exchange and long-range dipolar interactions reorganizes the low-energy manifold as the ratio is varied. Beyond a finite-size characterization, we identify several qualitatively new thermodynamic signatures induced solely by dipolar anisotropy. First, we demonstrate that ground-state level crossings generated by long-range interactions appear as exact zeros of the specific heat in the limit , establishing an unambiguous correspondence between microscopic spectral rearrangements and macroscopic caloric response. Second, we show that the shape of the associated Schottky-like anomalies encodes detailed information about the degeneracy structure of the competing low-energy states: odd lattices () display strongly asymmetric peaks due to unbalanced multiplicities, whereas the even lattice () exhibits three critical values of accompanied by nearly symmetric anomalies, reflecting paired degeneracies and revealing lattice parity as a key organizing principle. Third, we uncover a symmetry-driven crossover with increasing Q: while the and models retain sharp dipolar-induced critical points and pronounced low-temperature structure, for , the energy landscape becomes sufficiently smooth to suppress ground-state crossings altogether, yielding purely thermal specific-heat maxima. Altogether, our results provide a unified, size- and symmetry-resolved picture of how long-range anisotropy, lattice parity, and discrete rotational symmetry shape the thermodynamics of mesoscopic magnetic systems. We show that dipolar interactions alone are sufficient to generate nontrivial critical-like caloric behavior in clusters as small as , establishing exact finite-size benchmarks directly relevant for van der Waals nanomagnets, artificial spin-ice arrays, and dipolar-coupled nanomagnetic structures.Open AccessFinite-Size Thermodynamics of the Two-Dimensional Dipolar Q-Clock Model
We present a fully controlled thermodynamic study of the two-dimensional dipolar Q-state clock model on small square lattices with free boundaries, combining exhaustive state enumeration with noise-free evaluation of canonical observables. We resolve the complete energy spectra and degeneracies for the Ising case () on lattices of size , and for clock symmetries on a lattice, tracking how the competition between exchange and long-range dipolar interactions reorganizes the low-energy manifold as the ratio is varied. Beyond a finite-size characterization, we identify several qualitatively new thermodynamic signatures induced solely by dipolar anisotropy. First, we demonstrate that ground-state level crossings generated by long-range interactions appear as exact zeros of the specific heat in the limit , establishing an unambiguous correspondence between microscopic spectral rearrangements and macroscopic caloric response. Second, we show that the shape of the associated Schottky-like anomalies encodes detailed information about the degeneracy structure of the competing low-energy states: odd lattices () display strongly asymmetric peaks due to unbalanced multiplicities, whereas the even lattice () exhibits three critical values of accompanied by nearly symmetric anomalies, reflecting paired degeneracies and revealing lattice parity as a key organizing principle. Third, we uncover a symmetry-driven crossover with increasing Q: while the and models retain sharp dipolar-induced critical points and pronounced low-temperature structure, for , the energy landscape becomes sufficiently smooth to suppress ground-state crossings altogether, yielding purely thermal specific-heat maxima. Altogether, our results provide a unified, size- and symmetry-resolved picture of how long-range anisotropy, lattice parity, and discrete rotational symmetry shape the thermodynamics of mesoscopic magnetic systems. We show that dipolar interactions alone are sufficient to generate nontrivial critical-like caloric behavior in clusters as small as , establishing exact finite-size benchmarks directly relevant for van der Waals nanomagnets, artificial spin-ice arrays, and dipolar-coupled nanomagnetic structures. - Open Access
- 12 January 2026
- Published by MDPI AG in Preprints.org
Finite-Size Thermodynamics of the Two-Dimensional Dipolar Q-Clock Model
We present a fully controlled thermodynamic study of the two-dimensional dipolar $Q$-state clock model on small square lattices with free boundaries, combining exhaustive state enumeration with noise-free evaluation of canonical observables. We resolve the complete energy spectra and degeneracies $\{E_n,c_n\}$ for the Ising case ($Q=2$) on lattices of size $L=3,4,5$, and for clock symmetries $Q=4,6,8$ on a $3\times3$ lattice, tracking how the competition between exchange and long-range dipolar interactions reorganizes the low-energy manifold as the ratio $\alpha = D/J$ is varied. Beyond a finite-size characterization, we identify several qualitatively new thermodynamic signatures induced solely by dipolar anisotropy. First, we demonstrate that ground-state level crossings generated by long-range interactions appear as exact zeros of the specific heat in the limit $C(T \rightarrow 0,\alpha)$, establishing an unambiguous correspondence between microscopic spectral rearrangements and macroscopic caloric response. Second, we show that the shape of the associated Schottky-like anomalies encodes detailed information about the degeneracy structure of the competing low-energy states: odd lattices ($L = 3,5$) display strongly asymmetric peaks due to unbalanced multiplicities, whereas the even lattice ($L = 4$) exhibits three critical values of $\alpha$ accompanied by nearly symmetric anomalies, reflecting paired degeneracies and revealing lattice parity as a key organizing principle. Third, we uncover a symmetry-driven crossover with increasing $Q$: while the $Q=2$ and $Q=4$ models retain sharp dipolar-induced critical points and pronounced low-temperature structure, for $Q \ge 6$ the energy landscape becomes sufficiently smooth to suppress ground-state crossings altogether, yielding purely thermal specific-heat maxima. Altogether, our results provide a unified, size- and symmetry-resolved picture of how long-range anisotropy, lattice parity, and discrete rotational symmetry shape the thermodynamics of mesoscopic magnetic systems. We show that dipolar interactions alone are sufficient to generate nontrivial critical-like caloric behavior in clusters as small as $3\times3$, establishing exact finite-size benchmarks directly relevant for van der Waals nanomagnets, artificial spin-ice arrays, and dipolar-coupled nanomagnetic structures.- 12 January 2026
- Published by MDPI AG in Preprints.org
Open AccessFinite-Size Thermodynamics of the Two-Dimensional Dipolar Q-Clock Model
We present a fully controlled thermodynamic study of the two-dimensional dipolar $Q$-state clock model on small square lattices with free boundaries, combining exhaustive state enumeration with noise-free evaluation of canonical observables. We resolve the complete energy spectra and degeneracies $\{E_n,c_n\}$ for the Ising case ($Q=2$) on lattices of size $L=3,4,5$, and for clock symmetries $Q=4,6,8$ on a $3\times3$ lattice, tracking how the competition between exchange and long-range dipolar interactions reorganizes the low-energy manifold as the ratio $\alpha = D/J$ is varied. Beyond a finite-size characterization, we identify several qualitatively new thermodynamic signatures induced solely by dipolar anisotropy. First, we demonstrate that ground-state level crossings generated by long-range interactions appear as exact zeros of the specific heat in the limit $C(T \rightarrow 0,\alpha)$, establishing an unambiguous correspondence between microscopic spectral rearrangements and macroscopic caloric response. Second, we show that the shape of the associated Schottky-like anomalies encodes detailed information about the degeneracy structure of the competing low-energy states: odd lattices ($L = 3,5$) display strongly asymmetric peaks due to unbalanced multiplicities, whereas the even lattice ($L = 4$) exhibits three critical values of $\alpha$ accompanied by nearly symmetric anomalies, reflecting paired degeneracies and revealing lattice parity as a key organizing principle. Third, we uncover a symmetry-driven crossover with increasing $Q$: while the $Q=2$ and $Q=4$ models retain sharp dipolar-induced critical points and pronounced low-temperature structure, for $Q \ge 6$ the energy landscape becomes sufficiently smooth to suppress ground-state crossings altogether, yielding purely thermal specific-heat maxima. Altogether, our results provide a unified, size- and symmetry-resolved picture of how long-range anisotropy, lattice parity, and discrete rotational symmetry shape the thermodynamics of mesoscopic magnetic systems. We show that dipolar interactions alone are sufficient to generate nontrivial critical-like caloric behavior in clusters as small as $3\times3$, establishing exact finite-size benchmarks directly relevant for van der Waals nanomagnets, artificial spin-ice arrays, and dipolar-coupled nanomagnetic structures. - 30 December 2025
- Published by American Physical Society (APS) in Physical Review E
- Vol. 112 (6) , 065508
Ratio between Seebeck coefficient and entropy per particle as a tool for elementary charge determination
In this work, we investigate the relationship between the Seebeck coefficient (S), and the differential entropy per particle (DEP, s), as a tool for characterizing charge carriers in two-dimensional systems. Using armchair silicene nanoribbons as a model platform, we analyze how both quantities and their ratio depend on chemical potential at room temperature. While the Seebeck coefficient captures transport properties through the energy dependence of the electronic transmission, the DEP is directly connected to the system's electronic entropy, offering a direct thermodynamic alternative for estimating S. We evaluate these transport-thermodynamic properties considering diverse ribbon widths, defining metallic and semiconducting regimes. We find both quantities S and s, are highly interconnected within the ribbon's band gap energy region, and their ratio s/S converges to the elementary charge e across that energy window, fulfilling the Kelvin formula S=s/e. On the contrary, s/S is undefined for gapless ribbons in the energy window of the first transmission channel. These results establish the ratio between the DEP and the Seebeck coefficient as a reliable and complementary probe for the determination of the elementary charge, and to identify the cleanness of electronic band gaps as s/S matches with e.
- 30 December 2025
- Published by American Physical Society (APS) in Physical Review E
- Vol. 112 (6) , 065508
Ratio between Seebeck coefficient and entropy per particle as a tool for elementary charge determination
In this work, we investigate the relationship between the Seebeck coefficient (S), and the differential entropy per particle (DEP, s), as a tool for characterizing charge carriers in two-dimensional systems. Using armchair silicene nanoribbons as a model platform, we analyze how both quantities and their ratio depend on chemical potential at room temperature. While the Seebeck coefficient captures transport properties through the energy dependence of the electronic transmission, the DEP is directly connected to the system's electronic entropy, offering a direct thermodynamic alternative for estimating S. We evaluate these transport-thermodynamic properties considering diverse ribbon widths, defining metallic and semiconducting regimes. We find both quantities S and s, are highly interconnected within the ribbon's band gap energy region, and their ratio s/S converges to the elementary charge e across that energy window, fulfilling the Kelvin formula S=s/e. On the contrary, s/S is undefined for gapless ribbons in the energy window of the first transmission channel. These results establish the ratio between the DEP and the Seebeck coefficient as a reliable and complementary probe for the determination of the elementary charge, and to identify the cleanness of electronic band gaps as s/S matches with e.
- 1 July 2025
- Published by Wiley in Advanced Quantum Technologies
- Vol. 8 (10)
Quantum Level‐Crossing Induced by Anisotropy in Spin‐1 Heisenberg Dimers: Applications to Quantum Stirling Engines
This work explores the thermodynamic performance of a quantum Stirling heat engine implemented with an anisotropic spin‐1 Heisenberg dimer as the working medium. Using the Hamiltonian of the system, the interplay of anisotropy, magnetic field, and exchange interactions and their influence is analyzed on the energy spectrum and the quantum level crossing. These results reveal that double‐degenerate point (DDP) and a triple‐degenerate point (TDP) play pivotal roles in shaping the operational regimes and efficiency of the quantum Stirling engine. At those points, the Carnot efficiency reaches higher work output and enhanced stability, making it a robust candidate for optimal thermodynamic performance. These findings highlight the potential of anisotropic spin systems as viable platforms for quantum heat engines and contribute to advancing the field of quantum thermodynamics.
- 1 July 2025
- Published by Wiley in Advanced Quantum Technologies
- Vol. 8 (10)
Quantum Level‐Crossing Induced by Anisotropy in Spin‐1 Heisenberg Dimers: Applications to Quantum Stirling Engines
This work explores the thermodynamic performance of a quantum Stirling heat engine implemented with an anisotropic spin‐1 Heisenberg dimer as the working medium. Using the Hamiltonian of the system, the interplay of anisotropy, magnetic field, and exchange interactions and their influence is analyzed on the energy spectrum and the quantum level crossing. These results reveal that double‐degenerate point (DDP) and a triple‐degenerate point (TDP) play pivotal roles in shaping the operational regimes and efficiency of the quantum Stirling engine. At those points, the Carnot efficiency reaches higher work output and enhanced stability, making it a robust candidate for optimal thermodynamic performance. These findings highlight the potential of anisotropic spin systems as viable platforms for quantum heat engines and contribute to advancing the field of quantum thermodynamics.
- Open Access
Entropy Alternatives for Equilibrium and Out-of-Equilibrium Systems
We introduce a novel entropy-related function, non-repeatability, designed to capture dynamical behaviors in complex systems. Its normalized form, mutability, has been previously applied in statistical physics as a dynamical entropy measure associated with any observable stored in a sequential file. We now extend this concept by calculating the sorted mutability for the same data file previously ordered by increasing or decreasing value. To present the scope and advantages of these quantities, we analyze two distinct systems: (a) Monte Carlo simulations of magnetic moments on a square lattice, and (b) seismic time series from the United States Geological Survey catalog. Both systems are well established in the literature, serving as robust benchmarks. Shannon entropy is employed as a reference point to assess the similarities and differences with the proposed measures. A key distinction lies in the sensitivity of non-repeatability and mutability to the temporal ordering of data, which contrasts with traditional entropy definitions. Moreover, sorted mutability reveals additional insights into the critical behavior of the systems under study.
Open AccessEntropy Alternatives for Equilibrium and Out-of-Equilibrium Systems
We introduce a novel entropy-related function, non-repeatability, designed to capture dynamical behaviors in complex systems. Its normalized form, mutability, has been previously applied in statistical physics as a dynamical entropy measure associated with any observable stored in a sequential file. We now extend this concept by calculating the sorted mutability for the same data file previously ordered by increasing or decreasing value. To present the scope and advantages of these quantities, we analyze two distinct systems: (a) Monte Carlo simulations of magnetic moments on a square lattice, and (b) seismic time series from the United States Geological Survey catalog. Both systems are well established in the literature, serving as robust benchmarks. Shannon entropy is employed as a reference point to assess the similarities and differences with the proposed measures. A key distinction lies in the sensitivity of non-repeatability and mutability to the temporal ordering of data, which contrasts with traditional entropy definitions. Moreover, sorted mutability reveals additional insights into the critical behavior of the systems under study.
- 14 March 2025
- Published by American Physical Society (APS) in Physical Review E
- Vol. 111 (3) , 034118
Entropy, entanglement, and susceptibility of three qubits near quantum criticality
In this work, we investigate a system of three entangled qubits within the XXX model, subjected to an external magnetic field in the z direction and incorporating an anisotropy term along the y axis. We explore the system's thermodynamics by calculating its magnetic susceptibility and analyzing how this quantity encodes information about entanglement. By deriving rigorous bounds for the susceptibility, we demonstrate that violations of these bounds serve as an entanglement witness. Our results show that anisotropy enhances entanglement, extending the temperature range over which it persists. Furthermore, when the system is traced over the degrees of freedom of two qubits, the reduced entropy of the remaining single qubit corresponds to that of a system in a thermal bath with temperature T>0K.
- 14 March 2025
- Published by American Physical Society (APS) in Physical Review E
- Vol. 111 (3) , 034118
Entropy, entanglement, and susceptibility of three qubits near quantum criticality
In this work, we investigate a system of three entangled qubits within the XXX model, subjected to an external magnetic field in the z direction and incorporating an anisotropy term along the y axis. We explore the system's thermodynamics by calculating its magnetic susceptibility and analyzing how this quantity encodes information about entanglement. By deriving rigorous bounds for the susceptibility, we demonstrate that violations of these bounds serve as an entanglement witness. Our results show that anisotropy enhances entanglement, extending the temperature range over which it persists. Furthermore, when the system is traced over the degrees of freedom of two qubits, the reduced entropy of the remaining single qubit corresponds to that of a system in a thermal bath with temperature T>0K.