Concurrence of dynamical phase transitions at finite temperature in the fully connected transverse-field Ising model

We construct the finite-temperature dynamical phase diagram of the fully connected transverse-field Ising model from the vantage point of two disparate concepts of dynamical criticality. An analytical derivation of the classical dynamics and exact diagonalization simulations are used to study the dynamics after a quantum quench in the system prepared in a thermal equilibrium state. The different dynamical phases characterized by the type of non-analyticities that emerge in an appropriately defined Loschmidt-echo return rate directly correspond to the dynamical phases determined by the spontaneous breaking of $\mathbb{Z}_2$ symmetry in the long-time steady state. The dynamical phase diagram is qualitatively different depending on whether the initial thermal state is ferromagnetic or paramagnetic. Whereas the former leads to a dynamical phase diagram that can be directly related to its equilibrium counterpart, the latter gives rise to a divergent dynamical critical temperature at vanishing final transverse-field strength.Comment: journal article, 15 pages, 12 figures. Final versio

Fast logarithmic Fourier-Laplace transform of nonintegrable functions

We present an efficient and very flexible numerical fast Fourier-Laplace transform, that extends the logarithmic Fourier transform (LFT) introduced by Haines and Jones [Geophys. J. Int. 92(1):171 (1988)] for functions varying over many scales to nonintegrable functions. In particular, these include cases of the asymptotic form $f(\nu\to0)\sim\nu^a$ and $f(|\nu|\to\infty)\sim\nu^b$ with arbitrary real $a>b$. Furthermore, we prove that the numerical transform converges exponentially fast in the number of data points, provided that the function is analytic in a cone $|\Im{\nu}|<\theta|\Re{\nu}|$ with a finite opening angle $\theta$ around the real axis and satisfies $|f(\nu)f(1/\nu)|<\nu^c$ as $\nu\to 0$ with a positive constant $c$, which is the case for the class of functions with power-law tails. Based on these properties we derive ideal transformation parameters and discuss how the logarithmic Fourier transform can be applied to convolutions. The ability of the logarithmic Fourier transform to perform these operations on multiscale (non-integrable) functions with power-law tails with exponentially small errors makes it the method of choice for many physical applications, which we demonstrate on typical examples. These include benchmarks against known analytical results inaccessible to other numerical methods, as well as physical models near criticality.Comment: 14 pages, 8 figure

Implicit Peer Triplets in Gradient-Based Solution Algorithms for ODE Constrained Optimal Control

It is common practice to apply gradient-based optimization algorithms to numerically solve large-scale ODE constrained optimal control problems. Gradients of the objective function are most efficiently computed by approximate adjoint variables. High accuracy with moderate computing time can be achieved by such time integration methods that satisfy a sufficiently large number of adjoint order conditions and supply gradients with higher orders of consistency. In this paper, we upgrade our former implicit two-step Peer triplets constructed in [Algorithms, 15:310, 2022] to meet those new requirements. Since Peer methods use several stages of the same high stage order, a decisive advantage is their lack of order reduction as for semi-discretized PDE problems with boundary control. Additional order conditions for the control and certain positivity requirements now intensify the demands on the Peer triplet. We discuss the construction of 4-stage methods with order pairs (4,3) and (3,3) in detail and provide three Peer triplets of practical interest. We prove convergence for s-stage methods, for instance, order s for the state variables even if the adjoint method and the control satisfy the conditions for order s-1, only. Numerical tests show the expected order of convergence for the new Peer triplets.Comment: 47 pages, 5 figure

Dynamical Quantum Phase Transitions: A Geometric Picture

The Loschmidt echo (LE) is a purely quantum-mechanical quantity whose determination for large quantum many-body systems requires an exceptionally precise knowledge of all eigenstates and eigenenergies. One might therefore be tempted to dismiss the applicability of any approximations to the underlying time evolution as hopeless. However, using the fully connected transverse-field Ising model (FC-TFIM) as an example, we show that this indeed is not the case, and that a simple semiclassical approximation to systems well described by mean-field theory (MFT) is in fact in good quantitative agreement with the exact quantum-mechanical calculation. Beyond the potential to capture the entire dynamical phase diagram of these models, the method presented here also allows for an intuitive geometric interpretation of the fidelity return rate at any temperature, thereby connecting the order parameter dynamics and the Loschmidt echo in a common framework. Videos of the post-quench dynamics provided in the supplemental material visualize this new point of view.Comment: Accepted version. 7 pages with 4 Figures in main file. 3 pages including 2 Figures of supplemental material. 3 videos linked in the references of the main fil