Nikita Astrakhantsev
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
Nikita Astrakhantsev has a total of 56 co-authors. Here are the top co-authors by number of shared publications:
- Victor Valerjevich Braguta
Joint Institute for Nuclear Research
16 shared publications - Andrey Yuryevich Kotov
Jülich Supercomputing Centre, Forschungszentrum Jülich, Jülich, Germany
12 shared publications - Alexander Aleksandrovich Nikolaev
Department of Physics, College of Science, Swansea University, Swansea SA2 8PP, United Kingdom
6 shared publications - Artem A. Roenko
Joint Institute for Nuclear Research
3 shared publications - Ernst-Michael Ilgenfritz
Bogoliubov Laboratory of Theoretical Physics Joint Institute for Nuclear Research 141980 Dubna Russia
3 shared publications - A. Yu. Kotov
Bogoliubov Laboratory of Theoretical Physics
2 shared publications
Latest Publications
- Open Access
- 11 February 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (3)
Electrical conductivity of the quark-gluon plasma in the presence of strong magnetic fields
We compute the electrical conductivity of the strongly interacting medium in the presence of strong magnetic background fields, e B = 4 , 9 GeV 2 , and for different values of the temperature, both in the confined and in the deconfined quark-gluon plasma (QGP) phase. The conductivity is obtained from the Euclidean lattice time correlator of the electrical current, computed on gauge configurations sampled from Monte-Carlo simulations of an improved staggered discretization of N f = 2 + 1 QCD. We perform the inverse Laplace transform of the correlator adopting a recently proposed version of the standard Backus-Gilbert procedure for the inversion. The results obtained in the QGP phase show a sizable enhancement of the conductivity in the direction parallel to the magnetic field, as well as a suppression in the direction orthogonal to it. The enhancement can interpreted through various different mechanisms, typically involving the chiral anomaly, among which the chiral magnetic effect (CME): for this reason we try to interpret our observations based on a CME-inspired , we extract the relaxation time of this process, extrapolate it to the continuum limit and compare it to previous results, finding it lower than expected in the explored range of temperatures.
- 11 February 2025
- Published by American Physical Society (APS) in Physical Review D
- Vol. 111 (3)
Open AccessElectrical conductivity of the quark-gluon plasma in the presence of strong magnetic fields
We compute the electrical conductivity of the strongly interacting medium in the presence of strong magnetic background fields, e B = 4 , 9 GeV 2 , and for different values of the temperature, both in the confined and in the deconfined quark-gluon plasma (QGP) phase. The conductivity is obtained from the Euclidean lattice time correlator of the electrical current, computed on gauge configurations sampled from Monte-Carlo simulations of an improved staggered discretization of N f = 2 + 1 QCD. We perform the inverse Laplace transform of the correlator adopting a recently proposed version of the standard Backus-Gilbert procedure for the inversion. The results obtained in the QGP phase show a sizable enhancement of the conductivity in the direction parallel to the magnetic field, as well as a suppression in the direction orthogonal to it. The enhancement can interpreted through various different mechanisms, typically involving the chiral anomaly, among which the chiral magnetic effect (CME): for this reason we try to interpret our observations based on a CME-inspired , we extract the relaxation time of this process, extrapolate it to the continuum limit and compare it to previous results, finding it lower than expected in the explored range of temperatures.
- Open Access
- 22 May 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 109 (9) , 094511
QCD equation of state at nonzero baryon density in an external magnetic field
This paper is devoted to the study of quantum chromodynamics (QCD) equation of state in external magnetic field and nonzero baryon density. Our study is carried out by means of lattice simulation with dynamical staggered quarks at the physical masses. The simulation is conducted at imaginary baryon chemical potential what allowed us to overcome the sign problem. We expand the pressure in the baryon imaginary chemical potential and study three leading nonzero coefficients in this expansion. These coefficients were calculated for the following values of magnetic field: , 0.6, with the lattice sizes , , . Using these data we take continuum limit for the coefficients. Our results indicate considerable enhancement of the expansion coefficients by the magnetic field.- 22 May 2024
- Published by American Physical Society (APS) in Physical Review D
- Vol. 109 (9) , 094511
Open AccessQCD equation of state at nonzero baryon density in an external magnetic field
This paper is devoted to the study of quantum chromodynamics (QCD) equation of state in external magnetic field and nonzero baryon density. Our study is carried out by means of lattice simulation with dynamical staggered quarks at the physical masses. The simulation is conducted at imaginary baryon chemical potential what allowed us to overcome the sign problem. We expand the pressure in the baryon imaginary chemical potential and study three leading nonzero coefficients in this expansion. These coefficients were calculated for the following values of magnetic field: , 0.6, with the lattice sizes , , . Using these data we take continuum limit for the coefficients. Our results indicate considerable enhancement of the expansion coefficients by the magnetic field. - 11 March 2024
- Published by American Physical Society (APS) in Physical Review A
- Vol. 109 (3) , 032408
Challenges of variational quantum optimization with measurement shot noise
Quantum enhanced optimization of classical cost functions is a central theme of quantum computing due to its high potential value in science and technology. The variational quantum eigensolver (VQE) and the quantum approximate optimization algorithm (QAOA) are popular variational approaches that are considered the most viable solutions in the noisy-intermediate scale quantum (NISQ) era. Here, we study the scaling of the quantum resources, defined as the required number of circuit repetitions, to reach a fixed success probability as the problem size increases, focusing on the role played by measurement shot noise, which is unavoidable in realistic implementations. Simple and reproducible problem instances are addressed, namely, the ferromagnetic and disordered Ising chains. Our results show that: (1) VQE with the standard heuristic Ansatz scales comparably to direct brute-force search when energy-based optimizers are employed. The performance improves at most quadratically using a gradient-based optimizer. (2) When the parameters are optimized from random guesses, also the scaling of QAOA implies problematically long absolute runtimes for large problem sizes. (3) QAOA becomes practical when supplemented with a physically inspired initialization of the parameters. Our results suggest that hybrid quantum-classical algorithms should possibly avoid a brute force classical outer loop, but focus on smart parameters initialization.
- 11 March 2024
- Published by American Physical Society (APS) in Physical Review A
- Vol. 109 (3) , 032408
Challenges of variational quantum optimization with measurement shot noise
Quantum enhanced optimization of classical cost functions is a central theme of quantum computing due to its high potential value in science and technology. The variational quantum eigensolver (VQE) and the quantum approximate optimization algorithm (QAOA) are popular variational approaches that are considered the most viable solutions in the noisy-intermediate scale quantum (NISQ) era. Here, we study the scaling of the quantum resources, defined as the required number of circuit repetitions, to reach a fixed success probability as the problem size increases, focusing on the role played by measurement shot noise, which is unavoidable in realistic implementations. Simple and reproducible problem instances are addressed, namely, the ferromagnetic and disordered Ising chains. Our results show that: (1) VQE with the standard heuristic Ansatz scales comparably to direct brute-force search when energy-based optimizers are employed. The performance improves at most quadratically using a gradient-based optimizer. (2) When the parameters are optimized from random guesses, also the scaling of QAOA implies problematically long absolute runtimes for large problem sizes. (3) QAOA becomes practical when supplemented with a physically inspired initialization of the parameters. Our results suggest that hybrid quantum-classical algorithms should possibly avoid a brute force classical outer loop, but focus on smart parameters initialization.
- Open Access
- 7 December 2023
- Published by American Physical Society (APS) in Physical Review Research
- Vol. 5 (4) , 043214
Understanding symmetry breaking in twisted bilayer graphene from cluster constraints
Twisted bilayer graphene is an exciting platform for exploring correlated quantum phases, extremely tunable with respect to both the single-particle bands and the interaction profile of electrons. Here, we investigate the phase diagram of twisted bilayer graphene as described by an extended Hubbard model on the honeycomb lattice with two fermionic orbitals (valleys) per site. Besides the special extended cluster interaction , we incorporate the effect of gating through an on-site Hubbard-interaction . Within quantum Monte Carlo, we find valence-bond solid, Nel-valley antiferromagnetic or charge-density wave phases. Further, we elucidate the competition of these phases by noticing that the cluster interaction induces an exotic constraint on the Hilbert space, which we dub the cluster rule, in analogy to the famous pyrochlore spin-ice rule. Formulating the perturbative Hamiltonian by projecting into the cluster-rule manifold, we perform exact diagonalization and construct the fixed-point states of the observed phases. Finally, we compute the local electron density patterns as signatures distinguishing these phases, which could be observed with scanning tunneling microscopy. Our paper capitalizes on the notion of cluster constraints in the extended Hubbard model of twisted bilayer graphene and suggests a scheme towards realization of several symmetry-breaking insulating phases in a twisted-bilayer graphene sheet.- 7 December 2023
- Published by American Physical Society (APS) in Physical Review Research
- Vol. 5 (4) , 043214
Open AccessUnderstanding symmetry breaking in twisted bilayer graphene from cluster constraints
Twisted bilayer graphene is an exciting platform for exploring correlated quantum phases, extremely tunable with respect to both the single-particle bands and the interaction profile of electrons. Here, we investigate the phase diagram of twisted bilayer graphene as described by an extended Hubbard model on the honeycomb lattice with two fermionic orbitals (valleys) per site. Besides the special extended cluster interaction , we incorporate the effect of gating through an on-site Hubbard-interaction . Within quantum Monte Carlo, we find valence-bond solid, Nel-valley antiferromagnetic or charge-density wave phases. Further, we elucidate the competition of these phases by noticing that the cluster interaction induces an exotic constraint on the Hilbert space, which we dub the cluster rule, in analogy to the famous pyrochlore spin-ice rule. Formulating the perturbative Hamiltonian by projecting into the cluster-rule manifold, we perform exact diagonalization and construct the fixed-point states of the observed phases. Finally, we compute the local electron density patterns as signatures distinguishing these phases, which could be observed with scanning tunneling microscopy. Our paper capitalizes on the notion of cluster constraints in the extended Hubbard model of twisted bilayer graphene and suggests a scheme towards realization of several symmetry-breaking insulating phases in a twisted-bilayer graphene sheet. - Open Access
- 28 September 2023
- Published by American Physical Society (APS) in Physical Review Research
- Vol. 5 (3) , 033225
Phenomenological theory of variational quantum ground-state preparation
The variational approach is a cornerstone of computational physics, considering both conventional and quantum computing computational platforms. The variational quantum eigensolver algorithm aims to prepare the ground state of a Hamiltonian exploiting parametrized quantum circuits that may offer an advantage compared to classical trial states used, for instance, in quantum Monte Carlo or tensor network calculations. While, traditionally, the main focus has been on developing better trial circuits, we show that the algorithm's success, if optimized within stochastic gradient descent (SGD) or quantum natural gradient descent (QNGD), crucially depends on other parameters such as the learning rate, the number of measurements to estimate the gradient components, and the Hamiltonian gap . Within the standard SGD or QNGD, we first observe the existence of a finite value below which the optimization is impossible, and the energy variance resembles the behavior of the specific heat in second-order phase transitions. Second, when is above such threshold level, and learning is possible, we develop a phenomenological model that relates the fidelity of the state preparation with the optimization hyperparameters and . More specifically, we observe that the computational resources scale as , and we propose a symmetry enhancement of the variational ansatz as a way to increase the closing gap. We test our understanding on several instances of two-dimensional frustrated quantum magnets, which are believed to be the most promising candidates for near-term quantum advantage through variational quantum simulations.- 28 September 2023
- Published by American Physical Society (APS) in Physical Review Research
- Vol. 5 (3) , 033225
Open AccessPhenomenological theory of variational quantum ground-state preparation
The variational approach is a cornerstone of computational physics, considering both conventional and quantum computing computational platforms. The variational quantum eigensolver algorithm aims to prepare the ground state of a Hamiltonian exploiting parametrized quantum circuits that may offer an advantage compared to classical trial states used, for instance, in quantum Monte Carlo or tensor network calculations. While, traditionally, the main focus has been on developing better trial circuits, we show that the algorithm's success, if optimized within stochastic gradient descent (SGD) or quantum natural gradient descent (QNGD), crucially depends on other parameters such as the learning rate, the number of measurements to estimate the gradient components, and the Hamiltonian gap . Within the standard SGD or QNGD, we first observe the existence of a finite value below which the optimization is impossible, and the energy variance resembles the behavior of the specific heat in second-order phase transitions. Second, when is above such threshold level, and learning is possible, we develop a phenomenological model that relates the fidelity of the state preparation with the optimization hyperparameters and . More specifically, we observe that the computational resources scale as , and we propose a symmetry enhancement of the variational ansatz as a way to increase the closing gap. We test our understanding on several instances of two-dimensional frustrated quantum magnets, which are believed to be the most promising candidates for near-term quantum advantage through variational quantum simulations. - Open Access
- 14 September 2023
- Published by American Physical Society (APS) in Physical Review Research
- Vol. 5 (3) , 033187
Time evolution of uniform sequential circuits
Simulating time evolution of generic quantum many-body systems using classical numerical approaches has an exponentially growing cost either with evolution time or with the system size. In this work we present a polynomially scaling hybrid quantum-classical algorithm for time evolving a one-dimensional uniform system in the thermodynamic limit. This algorithm uses a layered uniform sequential quantum circuit as a variational Ansatz to represent infinite translation-invariant quantum states. We show numerically that this Ansatz requires a number of parameters polynomial in the simulation time for a given accuracy. Furthermore, this favorable scaling of the Ansatz is maintained during our variational evolution algorithm. All steps of the hybrid optimization are designed with near-term digital quantum computers in mind. After benchmarking the evolution algorithm on a classical computer, we demonstrate the measurement of observables of this uniform state using a finite number of qubits on a cloud-based quantum processing unit. With more efficient tensor contraction schemes, this algorithm may also offer improvements as a classical numerical algorithm.
- 14 September 2023
- Published by American Physical Society (APS) in Physical Review Research
- Vol. 5 (3) , 033187
Open AccessTime evolution of uniform sequential circuits
Simulating time evolution of generic quantum many-body systems using classical numerical approaches has an exponentially growing cost either with evolution time or with the system size. In this work we present a polynomially scaling hybrid quantum-classical algorithm for time evolving a one-dimensional uniform system in the thermodynamic limit. This algorithm uses a layered uniform sequential quantum circuit as a variational Ansatz to represent infinite translation-invariant quantum states. We show numerically that this Ansatz requires a number of parameters polynomial in the simulation time for a given accuracy. Furthermore, this favorable scaling of the Ansatz is maintained during our variational evolution algorithm. All steps of the hybrid optimization are designed with near-term digital quantum computers in mind. After benchmarking the evolution algorithm on a classical computer, we demonstrate the measurement of observables of this uniform state using a finite number of qubits on a cloud-based quantum processing unit. With more efficient tensor contraction schemes, this algorithm may also offer improvements as a classical numerical algorithm.
- Open Access
Understanding Symmetry Breaking in Twisted Bilayer Graphene from Cluster Constraints
Twisted bilayer graphene is an exciting platform for exploring correlated quantum phases, extremely tunable with respect to both the single-particle bands and the interaction profile of electrons. Here, we investigate the phase diagram of twisted bilayer graphene as described by an extended Hubbard model on the honeycomb lattice with two fermionic orbitals (valleys) per site. Besides the special extended {\it cluster interaction} $Q$, we incorporate the effect of gating through an onsite Hubbard-interaction $U$. Within Quantum Monte Carlo (QMC), we find valence-bond-solid, Néel-valley antiferromagnetic or charge-density wave phases. Further, we elucidate the competition of these phases by noticing that the cluster interaction induces an exotic constraint on the Hilbert space, which we dub {\it the cluster rule}, in analogy to the famous pyrochlore spin-ice rule. Formulating the perturbative Hamiltonian by projecting into the cluster-rule manifold, we perform exact diagonalization and construct the fixed-point states of the observed phases. Finally, we compute the local electron density patterns as signatures distinguishing these phases, which could be observed with scanning tunneling microscopy. Our work capitalizes on the notion of cluster constraints in the extended Hubbard model of twisted bilayer graphene, and suggests a scheme towards realization of several symmetry-breaking insulating phases in a twisted-bilayer graphene sheet.Open AccessUnderstanding Symmetry Breaking in Twisted Bilayer Graphene from Cluster Constraints
Twisted bilayer graphene is an exciting platform for exploring correlated quantum phases, extremely tunable with respect to both the single-particle bands and the interaction profile of electrons. Here, we investigate the phase diagram of twisted bilayer graphene as described by an extended Hubbard model on the honeycomb lattice with two fermionic orbitals (valleys) per site. Besides the special extended {\it cluster interaction} $Q$, we incorporate the effect of gating through an onsite Hubbard-interaction $U$. Within Quantum Monte Carlo (QMC), we find valence-bond-solid, Néel-valley antiferromagnetic or charge-density wave phases. Further, we elucidate the competition of these phases by noticing that the cluster interaction induces an exotic constraint on the Hilbert space, which we dub {\it the cluster rule}, in analogy to the famous pyrochlore spin-ice rule. Formulating the perturbative Hamiltonian by projecting into the cluster-rule manifold, we perform exact diagonalization and construct the fixed-point states of the observed phases. Finally, we compute the local electron density patterns as signatures distinguishing these phases, which could be observed with scanning tunneling microscopy. Our work capitalizes on the notion of cluster constraints in the extended Hubbard model of twisted bilayer graphene, and suggests a scheme towards realization of several symmetry-breaking insulating phases in a twisted-bilayer graphene sheet. - Open Access
Time Evolution of Uniform Sequential Circuits
Simulating time evolution of generic quantum many-body systems using classical numerical approaches has an exponentially growing cost either with evolution time or with the system size. In this work, we present a polynomially scaling hybrid quantum-classical algorithm for time evolving a one-dimensional uniform system in the thermodynamic limit. This algorithm uses a layered uniform sequential quantum circuit as a variational ansatz to represent infinite translation-invariant quantum states. We show numerically that this ansatz requires a number of parameters polynomial in the simulation time for a given accuracy. Furthermore, this favourable scaling of the ansatz is maintained during our variational evolution algorithm. All steps of the hybrid optimization are designed with near-term digital quantum computers in mind. After benchmarking the evolution algorithm on a classical computer, we demonstrate the measurement of observables of this uniform state using a finite number of qubits on a cloud-based quantum processing unit. With more efficient tensor contraction schemes, this algorithm may also offer improvements as a classical numerical algorithm.
Open AccessTime Evolution of Uniform Sequential Circuits
Simulating time evolution of generic quantum many-body systems using classical numerical approaches has an exponentially growing cost either with evolution time or with the system size. In this work, we present a polynomially scaling hybrid quantum-classical algorithm for time evolving a one-dimensional uniform system in the thermodynamic limit. This algorithm uses a layered uniform sequential quantum circuit as a variational ansatz to represent infinite translation-invariant quantum states. We show numerically that this ansatz requires a number of parameters polynomial in the simulation time for a given accuracy. Furthermore, this favourable scaling of the ansatz is maintained during our variational evolution algorithm. All steps of the hybrid optimization are designed with near-term digital quantum computers in mind. After benchmarking the evolution algorithm on a classical computer, we demonstrate the measurement of observables of this uniform state using a finite number of qubits on a cloud-based quantum processing unit. With more efficient tensor contraction schemes, this algorithm may also offer improvements as a classical numerical algorithm.
- 29 June 2023
- Published by American Physical Society (APS) in Physical Review B
- Vol. 107 (24) , 245145
Interacting topological quantum chemistry of Mott atomic limits
Topological quantum chemistry (TQC) is a successful framework for identifying (noninteracting) topological materials. Based on the symmetry eigenvalues of Bloch eigenstates at maximal momenta, which are attainable from first principles calculations, a band structure can either be classified as an atomic limit, in other words adiabatically connected to independent electronic orbitals on the respective crystal lattice, or it is topological. For interacting systems, there is no single-particle band structure and hence, the TQC machinery grinds to a halt. We develop a framework analogous to TQC, but employing -particle Green's function to classify interacting systems. Fundamentally, we define a class of interacting reference states that generalize the notion of atomic limits, which we call Mott atomic limits, and are symmetry protected topological states. Our formalism allows to fully classify these reference states (with ), which can themselves represent symmetry protected topological states. We present a comprehensive classification of such states in one dimension and provide numerical results on model systems. With this, we establish Mott atomic limit states as a generalization of the atomic limits to interacting systems.- 29 June 2023
- Published by American Physical Society (APS) in Physical Review B
- Vol. 107 (24) , 245145
Interacting topological quantum chemistry of Mott atomic limits
Topological quantum chemistry (TQC) is a successful framework for identifying (noninteracting) topological materials. Based on the symmetry eigenvalues of Bloch eigenstates at maximal momenta, which are attainable from first principles calculations, a band structure can either be classified as an atomic limit, in other words adiabatically connected to independent electronic orbitals on the respective crystal lattice, or it is topological. For interacting systems, there is no single-particle band structure and hence, the TQC machinery grinds to a halt. We develop a framework analogous to TQC, but employing -particle Green's function to classify interacting systems. Fundamentally, we define a class of interacting reference states that generalize the notion of atomic limits, which we call Mott atomic limits, and are symmetry protected topological states. Our formalism allows to fully classify these reference states (with ), which can themselves represent symmetry protected topological states. We present a comprehensive classification of such states in one dimension and provide numerical results on model systems. With this, we establish Mott atomic limit states as a generalization of the atomic limits to interacting systems.- Open Access
Equation of State of dense QCD in external magnetic field
In this proceeding we present our first results of the study of the QCD Equation of State at non-zero baryon density and in external magnetic field. We focused on the first three non-vanishing expansion coefficients of pressure in chemical potential and their dependence on magnetic field. The study is carried out within lattice simulations with $N_f=2+1$ dynamical quarks with physical quark masses. To overcome the sign problem, the simulations are carried out at imaginary baryon chemical potential. Our results suggest that external magnetic field considerably enhances the expansion coefficients and modifies their dependence on temperature.Open AccessEquation of State of dense QCD in external magnetic field
In this proceeding we present our first results of the study of the QCD Equation of State at non-zero baryon density and in external magnetic field. We focused on the first three non-vanishing expansion coefficients of pressure in chemical potential and their dependence on magnetic field. The study is carried out within lattice simulations with $N_f=2+1$ dynamical quarks with physical quark masses. To overcome the sign problem, the simulations are carried out at imaginary baryon chemical potential. Our results suggest that external magnetic field considerably enhances the expansion coefficients and modifies their dependence on temperature.