## Abstract

We propose a quantum critical detector (QCD) to amplify weak input signals. Our detector exploits a first-order discontinuous quantum-phase-transition and exhibits giant sensitivity (*χ* ∝ *N*^{2}) when biased at the critical point. We propose a model consisting of spins with long-range interactions coupled to a bosonic mode to describe the time-dynamics in the QCD. We numerically demonstrate dynamical features of the first order (discontinuous) quantum phase transition such as time-dependent quantum gain in a system with 80 interacting spins. We also show the linear scaling with the spin number *N* in both the quantum gain and the corresponding signal-to-quantum noise ratio during the time evolution of the device. Our work shows that engineering first order discontinuous quantum phase transitions can lead to a device application for metrology, weak signal amplification, and single photon detection.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

## 1. Introduction

Weak signal detection and amplification is a key process in metrology [1, 2], imaging [3], LIDAR [4], quantum communications and quantum computing [5]. There exist two widely used quantum amplification schemes for detecting weak signals. In the first kind, the weak input signal is directly amplified to generate a large output signal, as in quantum linear amplifiers [6–9]. It is widely utilized in circuit QED for parametric amplification in Josephson junction amplifiers [10]. Another scheme of amplification exists where the weak input perturbation functions as a control signal of an optimally biased device near the critical point. The practical realizations include single-electron transistors [11] and single-photon detectors [12]. This latter class can be classified as critically biased amplifiers. An important distinction of critically biased amplifiers from quantum linear amplifiers is that the input and output information carriers can be fundamentally different (eg: input photons and output electrons).

The amplification mechanism of a large class of weak signal detectors is based on thermodynamic (classical) phase transitions, which exhibits ultra-high sensitivity at the critical point, such as the superconducting single-photon detector [13] and the bubble chamber [14]. During the amplification, the incident weak signal triggers a phase transition of the critically biased detector to generate a large output signal. The goal of this paper is to propose a class of biased detectors with an amplification scheme that exploits quantum criticality in a first-order quantum phase transition (QPT). Our proposed device can be considered as an engineered quantum analog of widely utilized detectors which exploit naturally occurring thermodynamic phase transitions. We emphasize that the recent interest in engineering quantum phase transitions using qubit systems in ion traps [15], cold atoms [16] and circuit QED [17] forms an excellent platform to realize our proposed device. As QPTs happen at zero temperature, our proposed QCD may have higher signal-to-noise ratio and lower dark counting rate than detectors utilizing thermodynamic phase transitions.

A QPT at zero temperature describes an abrupt change in the ground state of a many-body system [18,19]. Most of the QPTs discovered in physical systems are of second-order [20–25] and they have also been proposed as a resource for metrology. However, no significant change in the number of excitations or energy transfer between subsystems occurs during a second-order QPT. To obtain an observable change in the output, a large parameter variation is required which can not be induced by a weak input signal (eg: single photon). Thus, high quantum gain can not be obtained using traditional second-order QPTs limiting their applicability for weak signal amplification. On the contrary, in a first-order (discontinuous) QPT, even a very small parameter variation at the critical point can lead to a significantly large change in the values of physical observables. Therefore, a universal model exhibiting first-order QPT with a well-understood microscopic mechanism can provide a natural platform for quantum amplification [6], quantum metrology [1,2], and lead to new types of single-photon detectors.

We emphasize that time dynamics near the critical point fundamentally determines whether QPTs can be a practical resource for amplification and detection. Specifically, critical sensitivity in QPT results only from the transition between the ground states of two different phases but can not be realized by a unitary adiabatic evolution. Thus, a dynamical detection event which necessarily involves excited states of the detector may not connect these two ground states leading to completely sub-optimal critical behavior. Even though few first-order QPT models have already been found [26–29], the dynamics of these QPTs around the critical point have not been revealed. The application of these first-order QPTs is therefore an open problem, as practical detection events and amplifications are fundamentally dynamical processes. In this paper, we overcome this important challenge related to time dynamics near the quantum critical point and also present a fundamental advancement in the design of weak signal detectors.

We introduce a first-order QPT model, composed of a bosonic mode and a spin ensemble with long-range interaction. We numerically show that there exists a critical point in this system where the sensitivity *χ* diverges with *N*^{2}-scaling (N the spin number). This scaling is much faster than previous first-order phase transitions [30,31] and provides extraordinary high sensitivity for weak signal detection. This first-order QPT and the giant sensitivity in our model fundamentally originate from the competition between two phases with long-range spin order. Exploiting this unique criticality, we propose a quantum critical detector (QCD) utilizing a weak input signal triggered first-order QPT [see the schematic in Fig. 1(a)]. We emphasize that the time-dynamics near a phase transition is a formidable challenge and we overcome this using the unique features of our proposed model. Via direct numerical evaluations, we demonstrate the time-dynamical features of a QCD consisting of 80 spins interacting with a bosonic mode. We show the existence of a first-order QPT in our QCD which sheds light on the microscopic origin of time-dependent quantum gain *g*(*t*). Finally, we show the linear scaling in both the maximum quantum gain and the corresponding signal-to-quantum noise ratio (SQNR) during the time-evolution of our detector revealing high figures of merit. Our work motivates quantum devices with engineered discontinuous phase transitions for single photon detection.

## 2. Detection using discontinuous quantum phase transition

The key element of a QCD is the first-order QPT based quantum amplification. Here, we introduce a specific first-order QPT model composed of a bosonic mode and an interacting spin ensemble. The input to the device can be a weak signal (eg: single photon) which either perturbs the coupling between the bosonic mode and the spins or the couping between the individual spins. The output of the detector i.e. the amplified signal is the macroscopic population of the bosonic mode (see Fig. 1). The Hamiltonian of the system is given by

*d̂*(

*d̂*

^{†}) denotes the output bosonic mode. Its frequency has been taken as the unit of energy

*ω*

_{0}= 1 and all the other parameters in the Hamiltonian have been rescaled by

*ω*

_{0}. The operators ${\widehat{\sigma}}_{j}^{\alpha}$ (

*α*=

*x*,

*y*,

*z*) are the Pauli matrices of the

*j*th spin. A magnetic field is applied along the

*z*-direction inducing a energy splitting

*∊*between spin states |↑〉

*and |↓〉*

_{j}*. The spin-boson coupling is along*

_{j}*x*-direction with homogeneous coupling strength

*λ*similar to the Dicke model [32]. The last term describes the all-to-all homogeneous coupling between the spins along the

*y*-direction akin to the Lipkin-Meshkov-Glick (LMG) model [21–23]. Here, we call this model as Dicke-LMGy model and emphasize the critical distinction from the Dicke-Ising model which only has nearest neighbour spin interactions [26]. The amplification process in the QCD is triggered by the weak input signal induced variation in the spin-boson coupling

*λ*or equivalently the spin-spin coupling

*J*.

To characterize the quantum phases and the corresponding QPTs in our Dicke-LMGy model,
we introduce two new magnetic order parameters (OPs): ${\zeta}_{M,x}={\u3008{\widehat{S}}_{x}^{2}\u3009}_{0}/{N}^{2}$ and ${\zeta}_{M,y}={\u3008{\widehat{S}}_{y}^{2}\u3009}_{0}/{N}^{2}$ (${\widehat{S}}_{\alpha}={\sum}_{j}{\widehat{\sigma}}_{j}^{\alpha}/2$ and 〈· · ·
〉_{0} means averaging on the ground state) characterizing
the magnetic fluctuations in the spins along *x* and
*y* axes, respectively. Such OPs can be probed
experimentally through spin noise spectroscopy [33]. Note, we do not choose the traditional
magnetic OP *M _{z}* =
〈

*Ŝ*〉

_{z}_{0}/

*N*[34], due to its incapability of characterizing the first-order QPT in our model [35]. The superradiant OP

*ζ*= 〈

_{S}*d̂*

^{†}

*d̂*〉

_{0}/

*N*is utilized to characterize the macroscopic excitation in the bosonic mode and functions as the output observable of our QCD. We emphasize that the fundamental competition between two ferromagnetic phases gives rise to a first order phase transition.

There exist three quantum phases in our model: paramagnetic-normal (PN) phase, ferromagnetic-normal (FN) phase, and ferromagnetic-superradiant (FS) phase as shown in Fig. 1(b) (more information about the quantum phases and ground-state wave functions can be obtained in Appendix A and B). In the *λJ*-plane, the the blue line determined by the critical spin-spin coupling strength *J*_{c,II} ≡ *∊*/2, the green line determined by the critical spin-boson coupling ${\lambda}_{c,\text{II}}\equiv \sqrt{\u220a}/2$, and the red line *J* = 2*λ*^{2} give the phase boundaries. We emphasize that the FN-FS boundary displays a first order QPT which is necessary for our QCD. There exists a unique triple-point (*λ*_{c,II}, *J*_{c,II}) determined by the energy splitting *∊*. The numerical demonstration of the phase diagram is given in Fig. 2 with the superradiant OP *ζ _{S}* and the magnetic OP

*ζ*

_{M,y}in panels (a) and (b), respectively. The OP

*ζ*

_{M,x}not shown here behaves similarly to

*ζ*. We exploit mean field theory [35,36] to verify this phase diagram obtained by our direct numerical calculation. We emphasize that the numerical approach we have introduced holds a fundamental advantage for dynamical amplification and noise calculations.

_{S}The most striking property of our tractable model is that the QPT between the FN phase
and the ferromagnetic-superradiant (FS) phase is of first order making QCD
a highly sensitive device. This first order transition occurs in the
strong coupling regime for large spin-spin coupling (*J*
> *J*_{c,II}) and
spin-boson coupling (*λ* >
*λ*_{c,II}). As shown by
the black, blue, and gray lines in Fig.
3(a) and 3(b), discontinuous changes are observed in both the
superradiant OP *ζ _{S}* and the magnetic OP

*ζ*

_{M,y}. Note that the overlapping red, pink and green curves correspond to low spin-spin coupling giving rise only to a second order transition between the PN and FS phase, not suitable for QCD. Once the spin-spin coupling increases, there is no paramagnetic phase and there exists only a ferromagnetic phase for all spin-boson coupling strengths. This ferromagnetism is evident by studying the magnetic order parameter in Fig. 3(b) (black, blue and gray curves). We emphasize that the ferromagnetic behavior shown in Fig. 3(b) is only along the y direction and the exact opposite trend occurs for the x direction (not shown). During the FN to FS transition which is first order, the energy prestored in the spins transfers to the bosonic mode coupled with a change in the spin fluctuations from the

*y*-direction to the

*x*-direction. The first-order QPT results from the competition between the FS phase – that arises from strong spin-boson coupling along the

*x*-axis and the FN phase – caused by large spin-spin coupling along the

*y*-axis. We also predict that this type of first-order QPT should also exist in the Ising XY-model [27]. Another important characteristic is that the first-order phase transition point is sensitive to the spin-spin coupling (bias) while the second order ones in Fig. 2(a) is not.

We now explicitly show that the FN-FS phase transition is of first order in Fig. 3(c) as required for QCD. Increasing the spin number, the phase transition shows critical scaling behavior. Here, the sensitivity of the system at the critical point is determined by first-order derivative of the superradiant OP *ζ _{S}*,

*N*is for consistency with the magnetic susceptibility. In Fig. 3(d), we plot the maximum of the sensitivity function (

*χ*

_{max}) at the first-order critical point ${\lambda}_{c,\text{I}}\equiv \sqrt{J/2}$ vs spin number (

*N*). In the thermodynamic limit

*N*→ ∞,

*χ*

_{max}diverges with speed ∝

*N*

^{2}, which is much faster than the $\sqrt{N}$-scaling obtained in the previous first-order dissipative transition [31] or the linear

*N*scaling in the first-order thermodynamic phase transition predicted by Imry [37]. There are two main reasons resulting in the sensitivity differences between our device and the one using a thermodynamic phase transition of the Dicke-Ising model in Ref. [30]. First, we use the first-order QPT at zero temperature instead of a first-order thermodynamic phase transition at a finite temperature as in [30]. ]. Second, the 1/

*N*

^{2}-scaling in sensitivity is obtained for small parameter variation at the critical point in the QPT. This is different from the 1/

*N*scaling in sensing the temperature changing in [30]. We predict that higher sensitivity on the estimation of system parameters, especially the coupling strength (i.e., the spin-boson coupling

*λ*or the spin-spin coupling

*J*), can also be obtained via utilizing the first-order QPT in the Dicke-Ising model. This

*N*

^{2}scaling, arising from competing phases, may be used to enhance the sensitivity in quantum metrology [38].

## 3. Dynamical quantum amplification

To utilize quantum criticality as a resource of quantum amplification, one has to study the dynamical behavior around the critical point. The challenge in exploiting QPTs for quantum amplification arises fundamentally from the fact that ground states of two different phases can not be connected via an adiabatic operation [39], as the gap of the system vanishes at the critical point in the thermodynamic limit. For example, starting from the ground state of the FN phase, the system can alternatively evolve to an arbitrary excited state instead of going to the ground state of the FS phase thereby completely negating critical amplification. To demonstrate the dynamical nonlinear amplification, here we show the dynamics of the first-order QPT in our model with 80 spins via direct numerical time-evolution. A QCD is demonstrated with the large quantum gain around the critical point after a QPT triggered by a time-varying parameter. We also show that a linear scaling with the spin number *N* in the quantum gain and SQNR of the QCD is obtained, instead of the *N*^{2} in the sensitivity of the first-order QPT.

The full measurement in our QCD is split into two main processes: transduction (absorption) and amplification. After the transduction, the excitations or energy in the input signal will be transferred into the detector. We model this process by considering a temporal variation in the spin-boson coupling. A physical mechanism for this process can be the radiation pressure of the photon falling on the walls of a cavity. An equivalent model can be a temporal variation in the spin-spin coupling. We note that time-dependent theory of phase transitions as well as single photon interaction with a system biased near a phase transition is an open challenge. To overcome this, we use the following approach to simulate the time dynamics in our system. We fix the trajectory of a parameter in the Hamiltonian and directly evaluate the time-dependent dynamics of the spins as well as the bosonic mode. This is numerically tractable due to the reduced Hilbert space in a system with all-to-all spin coupling.

To depict the weak signal interaction with our detector, we introduce a time-dependent variation of the spin-boson coupling strength *λ*(*t*) = *λ*_{0} +Δ*λ*× *P _{e}*(

*t*) and solve the entire time dynamics near the first order phase transition. The key amplification process in our QCD occurs from first-order QPT process triggered by this time-dependent parameter

*λ*(

*t*). The spin-boson coupling

*λ*

_{0}is optimally tuned to be very close to the critical point

*λ*

_{c,I}and the detector is prepared on the ground state of the system, which functions as the “quantum” bias. Thus, even a very small parameter variation Δ

*λ*(amplitude) can trigger a QPT and leads to efficient amplification. As mentioned before, amplification is defined as the enhanced number of bosons in the cavity mode. The time-dependent envelope

*P*(

_{e}*t*) of the coupling

*λ*(

*t*) is determined by the transduction process. One example of such a weak input signal is a single-photon. The absorption process is probabilistic [40] and results in a wave form of the quantum excitation in a detector (more information about

*P*(

_{e}*t*) can be found in Appendix C). The full dynamics of the whole system is governed by the time-dependent Hamiltonian

*H*(

*t*) in Eq. (1) by replacing

*λ*with

*λ*(

*t*). To characterize the detection sensitivity, we define the quantum gain of the amplification in our QCD as

*d̂*

^{†}(

*t*)

*d̂*(

*t*)〉 is the mean value of the time-dependent operator in the Heisenberg picture on the ground state of initial Hamiltonian

*H*(0) with

*λ*=

*λ*

_{0}.

The dynamical amplification in the QCD is demonstrated by the time-dependent gain as a function of the bias spin-boson coupling *λ*_{0} and time in Fig. 4. We see that the efficient amplification can only be obtained if the system is optimally biased close to the critical point. In the right subplot computed at fixed time *t* = 40, we explicitly show the high gain of our QCD around the critical point. Here, we show the highly sensitive nature of a time-dependent first-order QPT system to the initial bias. Similar to the enhanced decay of the Loschmidt echo by the criticality in a second-order QPT [19], the enhanced quantum gain in our QCD is a universal characteristic of criticality in a first-order QPT, which can be tested in experiments.

To show the high figures of merit of our QCD, we present the scaling of the quantum amplification with the spin number *N* in Fig. 5. The maximum gain is shown to be linearly proportional to *N*. This result should be contrasted with the amplification resulting from the second order phase transition with *J* ≤ *J*_{c,II} (see pink triangle line and the green circle line). The latter is much smaller than the one from the first-order QPT with *J* > *J*_{c,II} (see the blue diamond line). In the right subplot, we show the slope of *g*_{max} with the spin number for different spin-spin coupling. There exists a “phase transition” phenomenon in the slope at the same critical point *J*_{c,II} of the transition from second-order to first-order QPTs. To characterize the quantum noise in our QCD, we define the SQNR as [41]

*d̂*

^{†}(

*t*)

*d̂*(

*t*)]

^{2}〉 = 〈[

*d̂*

^{†}(

*t*)

*d̂*(

*t*)]

^{2}〉 − 〈

*d̂*

^{†}(

*t*)

*d̂*(

*t*)〉

^{2}is the variance of the bosonic excitation number operator. The corresponding SQNR for the three lines are displayed in the left subplot. It shows amplification based on first-order QPT has much higher SQNR than that of second-order QPTs. Similar to the quantum gain (the re-scaled excitation number in the output bosonic mode), the SQNR also increases linearly with the spin number, which is consistent with the SQNR of a coherent state as the final output state. Our proposed detector exploits quantum criticality of interacting spins for amplification. Hence, it is necessary for the amplification time and phase transition to complete before the decoherence of spins occurs. This is possible in principle as evidenced by recent demonstrations of quantum phase transitions in engineered ion trap [15], Rydberg atom [16], and circuit QED [17] systems.

## 4. Conclusion

We have proposed to engineer first-order (discontinuous) QPTs for weak-signal detection. Our proposed device is a quantum analog of widely used detectors which exploit thermodynamic phase transitions. One of the key differences between the thermodynamic and quantum phase transition based detectors is that the thermodynamic process requires the weak signal to change the temperature of the macroscopic system. In contrast, the quantum phase transition can work at constant low temperature and the weak signal can couple to other parameters/dynamical degrees of freedom in the system. As an example, we explicitly demonstrated the dynamical amplification by utilizing the high-sensitivity at the first-order critical point in the Dicke-LMG model we introduced. In recent experiments, second-order QPTs have been demonstrated with Bose-Einstein condensates [42], trapped ions [15], cold atoms [16], and superconducting qubits [17]. However, the quantum critical detector requires engineered first order phase transitions which is still an open challenge. We believe these platforms provide a starting point for physical realization of our proposed device where single photon perturbations can trigger a phase transition in systems with long-range interactions. Our work also motivates large scale macroscopic modeling of devices with a single quantum of energy propagating through it.

## A. Defining quantum phases in the detector

By defining the collective angular momentum operators of the *N* spins

*N*/2 [43]. As the total angular momentum of the spin ensemble is conserved, we can perform the numerical simulation within a subspace spanned by the

*N*+ 1 Dicke states [32]. This makes our model become tractable in numerical simulation.

There exist three quantum phases in our model: paramagnetic-normal (PN) phase, ferromagnetic-normal (FN) phase, and ferromagnetic-superradiant (FS) phase. The phase diagram of our model can be obtained via direct numerical evaluation. By exploiting a mean-field theory [36], we verify our numerical simulation and also obtain the exact analytical wave function of the ground states [35].

## B. Ground-State Wave Function of the detector

To reveal the underlying microscopic mechanism of our quantum critical detector (QCD), we need to understand the fundamental changes in the ground state of the system during the quantum phase transitions (QPTs). As we known, the Husimi *Q*-function is a quantum analog of the classical distribution function in the phase space, which provide a powerful tool to vividly show the ground-state wave function of the system in each phase.

The bosonic Husimi *Q* is defined as [44],

*ρ*is the density matrix of the ground state, |

_{g}*α*〉 is an arbitrary bosonic coherent state, and Tr

_{spin}[· · · ] means tracing off the degrees of freedom of the spins. For the spin degree of freedom, the Bloch sphere is usually utilized to characterize an arbitrary state of spin-1/2 particle. But for a spin-

*n*particle, a [(2

*n*+ 1)

^{2}− 1]-dimensional-sphere is required. It is impossible to show such a high-dimensional sphere. To overcome this issue, we introduce the spin Husimi

*Q*-function [45],

*N*is the spin number, Tr

_{boson}[· · · ] means tracing off the degrees of freedom of the bosonic mode, and |

*θ*,

*ϕ*〉 is a coherent spin state [46,47]. To let the angle

*θ*be the exact same polar angle of a spherical coordinate, we redefine the coherent spin state as,

*N*/2,

*N*/2〉 is the Dicke state with all the spins on the up state.

Our numerical approach allows us to directly calculate the ground states and the corresponding Husimi *Q*-functions of the system. In the paramagnetic-normal (PN) phase, the bosonic mode is on the vacuum state. As shown in Fig. 6(a), the bosonic *Q*-function has only one peak located at the origin. In this phase, the spins are on the Dicke state |*N*/2, −*N*/2〉, i.e., all the spins are polarized along negative *z*-axis. As shown in Fig. 7(a), the corresponding spin *Q*-function is a cigar-like structure lying along negative *z*-axis. From this quasi-probability distribution function, we can see that the mean magnetization *M _{z}* = 〈

*Ŝ*〉

_{z}_{0}/

*N*(〈· · · 〉

_{0}means averaging on the ground state) is a finite negative value, but the two magnetic order parameters ${\zeta}_{M,x}={\u3008{\widehat{S}}_{x}^{2}\u3009}_{0}/{N}^{2}$ and ${\zeta}_{M,y}={\u3008{\widehat{S}}_{y}^{2}\u3009}_{0}/{N}^{2}$ are very small and will go to zero in the limit

*N*→ ∞.

The ferromagnetic-superradiant (FS) phase has two degenerate ground states |*α*_{0}〉 ⊗ |*θ*_{0}, *π*〉 and |− *α*_{0}〉 ⊗ |*θ*_{0}, 0〉. Here, |*α*_{0}〉 is a bosonic coherent state and |*θ*_{0}, *ϕ*_{0}〉 (*ϕ*_{0} = 0, *π*) is a coherent spin state. The value of *α*_{0} and *θ*_{0} can be determined by the mean-field theory in the thermodynamic limit *N* → ∞ [35]. When the system adiabatically goes to the FS phase, the system can be on an arbitrary superposition of these two degenerate states. Thus, the ensemble mean value of the displacement of the bosonic mode 〈*d̂*^{†} + *d̂*〉_{0}, the magnetization along *x*-direction 〈*Ŝ _{x}*〉

_{0}, and the magnetization along

*y*-direction 〈

*Ŝ*〉

_{y}_{0}are all zero. But the excitation number reflected in the superradiant order parameter

*ζ*= 〈

_{S}*d̂*

^{†}

*d̂*〉

_{0}/

*N*and the magnetic noise characterized by the magnetic order parameter

*ζ*

_{M,x}are finite. In the numerical simulation, we choose the ground state as a symmetric superposition of these two degenerate states. In Fig. 6(b), we show the

*Q*-function of the bosonic mode. The two separated peaks on the real axis indicate the macroscopic excitation in the bosonic mode. In Fig. 7(b), we also see the rotation of the spins in the

*xz*-plane as the strong spin-boson coupling

*λ*is along

*x*-axis. These two branches correspond to the two degenerate states.

The ferromagnetic-normal (FN) phase also has two degenerate states |0〉 ⊗ |*θ*_{0}, *π*/2〉 and |0〉 ⊗ |*θ*_{0}, 3*π*/2〉. When the system adiabatically goes to the FN phase, the bosonic mode is always on the vacuum state as shown by the corresponding *Q*-function in Fig. 6(c). But the spins can be on an arbitrary superposition of these two coherent spin states |*θ*_{0}, *ϕ*_{0}〉 (*ϕ*_{0} = *π*/2, 3*π*/3). Thus, the ensemble mean of the magnetization in *xy*-plane is still zero. But the magnetic noise characterized by the magnetic order parameter *ζ*_{M,x} is finite. The *Q*-function of the spin in the FN phase is displayed in Fig. 7(c). The strong spin-spin coupling *J* along *y*-axis leads to the rotation of the spins in *yz*-plane. Here, we can also see that, when QPT from the FN phase to the FS phase occurs, the energy prestored in the spins transfers to the bosonic mode contributing to the macroscopic excitation in the bosonic mode. Also, a fundamental change in the spin noise from the *y*-direction to the *x*-direction can be observed.

The QPTs of PN↔FN and PN↔FS are of second order. Only the transition FN↔FS is a first order one. We emphasize that second order phase transitions which are continuous phase transitions do not possess the necessary features for quantum critical detection. In stark contrast, first order discontinuous phase transitions exhibit a giant response when a weak perturbation is applied making them an ideal resource for quantum critical detection.

## C. Analogy with practical single-photon detectors

To clearly explain the motivation of our work, we show the analogy of our proposed detector with the practical single-photon detectors. As we mentioned in the main text, there are two main amplification schemes: quantum linear amplifiers and critically biased amplifiers. We are focusing on the second amplification scheme, where the weak input signal does not get amplified directly. Instead, it functions as a control of an optimally biased critical system, which is significantly different from the first one, such as quantum linear amplifiers. In these critically biased amplifiers, the input and output information carriers can be fundamentally different (eg: input photons and output electrons) and the corresponding gain is defined as the ratio of the outputs with and without the input control signal. In this section, we show the analogy of our proposed QCD with the practical critical detectors to explain the motivation of our work more clearly. Finally, we show that first-order QPT-based devices can pave a way for new types of weak signal detector.

We first take the superconducting nanowire single-photon detector (SSNPD) as an example to explain the critical amplification scheme explicitly. The SSNPD is the best available and widely used near-infrared single-photon detector with ≥ 95% quantum efficiency [48], < 3 picoseconds timing jitter [49] and < 1 dark count per hour [50]. The input is a single-photon pulse—an extremely weak quantum signal. The current in the superconducting nanowire is biased very close to the critical current, thus even a single-photon pulse can break the superconductivity. The output signal is the voltage difference between the two ends of the superconducting nanowire. In the transduction (absorption) process, the incident single-photon pulse will general one resonantly excited electron. As the center frequency of the pulse is much larger than the energy gap of the superconductor, this highly excited electron will break hundreds of Cooper pairs, reduce the local density of the superconducting electrons, and finally triggers a phase transition from a superconductor to a normal metal. Before the absorption of the photon, the voltage difference is extremely small. After the absorption of the photon, the superconducting nanowire becomes a normal metal. An observable output voltage pulse will be generated to realize the amplification. The amplification process in the SNSPD can be modeled as a time-varying local temperature induced thermodynamic phase transition. Similarly, we describe the amplification process of our QCD as the time-dependent system parameter variation induced first-order QPT. During the amplification process, the exact dynamical change in the superconducting electron density and the effective time-dependent local temperature can be obtained by numerical simulation.

This type of critically biased amplifiers has been extensively used in practical measurements even outside the context of SNSPDs, like the photomultiplier tube, the single-photon avalanche diode, single-electron transistor, etc. However, our proposed QCD is the first quantum analog of the classical critical detectors. Similarly, we also need to bias the detector very close to the critical point. After the absorption of the input weak signal, a time-dependent variation in the system parameter (such as the spin-boson coupling or spin-spin interaction strength) is induced to trigger a first-order dynamical QPT. For a specific measurement process, we also need to calculate the exact form of time-dependent change in system parameters. Without loss of generality, we assume the time-varying spin-boson coupling change is proportional to the signal absorption probability *P _{e}*(

*t*) as shown in the subplot of Fig. 4 in the main text. The absorption probability

*P*(

_{e}*t*) can be obtained by calculating the master equation of a quantized time-dependent pulse interacting with a quantum system as shown in [40,51]. We show that, if the QCD is biased close to the critical point, a large amplification factor (the quantum gain) can be obtained.

During the dynamical amplification, the system parameter (coupling strength) is varied across the phase boundary time-dependently by the incident pulse. It is widely debated and is an open question whether QPTs retain their criticality during a dynamical process. The discontinuous change of the observable and the *N*^{2} sensitivity found in the first-order QPT only occurs if the system goes from the ground state in one phase to the ground state of another phase. However, in a dynamical process, the system will not go to the ground state of the other phase but evolves to some complicated excited state. Therefore, we numerically studied the dynamics of the system around the critical point. We showed for the first time that the high sensitivity also exists in a dynamical process and thus explicitly demonstrated the dynamical amplification.

## Funding

DARPA DETECT ARO award (W911NF-18-1-0074).

## References

**1. **V. Giovannetti, S. Lloyd, and L. Maccone, “Quantum-enhanced measurements: beating the standard quantum limit,” Science **306**, 1330–1336 (2004). [CrossRef] [PubMed]

**2. **V. Giovannetti, S. Lloyd, and L. Maccone, “Quantum metrology,” Phys. Rev. Lett. **96**, 010401 (2006). [CrossRef] [PubMed]

**3. **L. Lugiato, A. Gatti, and E. Brambilla, “Quantum imaging,” J. Opt. B: Quantum semiclassical optics **4**, S176 (2002). [CrossRef]

**4. **A. Swatantran, H. Tang, T. Barrett, P. DeCola, and R. Dubayah, “Rapid, high-resolution forest structure and terrain mapping over large areas using single photon lidar,” Sci. Reports **6**, 28277 (2016). [CrossRef]

**5. **M. A. Nielsen and I. Chuang, “Quantum computation and quantum information,” (2002).

**6. **C. M. Caves, “Quantum limits on noise in linear amplifiers,” Phys. Rev. D **26**, 1817 (1982). [CrossRef]

**7. **W. H. Louisell, A. Yariv, and A. E. Siegman, “Quantum fluctuations and noise in parametric processes. I.,” Phys. Rev. **124**, 1646–1654 (1961). [CrossRef]

**8. **B. R. Mollow and R. J. Glauber, “Quantum theory of parametric amplification. i,” Phys. Rev. **160**, 1076–1096 (1967). [CrossRef]

**9. **U. Gavish, B. Yurke, and Y. Imry, “Generalized constraints on quantum amplification,” Phys. Rev. Lett. **93**, 250601 (2004). [CrossRef]

**10. **A. Roy and M. Devoret, “Introduction to parametric amplification of quantum signals with josephson circuits,” Comptes Rendus Physique **17**, 740–755 (2016). [CrossRef]

**11. **M. H. Devoret and R. J. Schoelkopf, “Amplifying quantum signals with the single-electron transistor,” Nature **406**, 1039 (2000). [CrossRef] [PubMed]

**12. **M. Eisaman, J. Fan, A. Migdall, and S. V. Polyakov, “Invited review article: Single-photon sources and detectors,” Rev. scientific instruments **82**, 071101 (2011). [CrossRef]

**13. **G. Gol’Tsman, O. Okunev, G. Chulkova, A. Lipatov, A. Semenov, K. Smirnov, B. Voronov, A. Dzardanov, C. Williams, and R. Sobolewski, “Picosecond superconducting single-photon optical detector,” Appl. Phys. Lett. **79**, 705–707 (2001). [CrossRef]

**14. **D. A. Glaser, “Some effects of ionizing radiation on the formation of bubbles in liquids,” Phys. Rev. **87**, 665 (1952). [CrossRef]

**15. **J. Zhang, G. Pagano, P. W. Hess, A. Kyprianidis, P. Becker, H. Kaplan, A. V. Gorshkov, Z.-X. Gong, and C. Monroe, “Observation of a many-body dynamical phase transition with a 53-qubit quantum simulator,” Nature **551**, 601 (2017). [CrossRef] [PubMed]

**16. **H. Bernien, S. Schwartz, A. Keesling, H. Levine, A. Omran, H. Pichler, S. Choi, A. S. Zibrov, M. Endres, and M. Greiner *et al.*, “Probing many-body dynamics on a 51-atom quantum simulator,” Nature **551**, 579 (2017). [CrossRef] [PubMed]

**17. **R. Harris, Y. Sato, A. J. Berkley, M. Reis, F. Altomare, M. H. Amin, K. Boothby, P. Bunyk, C. Deng, C. Enderud, S. Huang, E. Hoskinson, M. W. Johnson, E. Ladizinsky, N. Ladizinsky, T. Lanting, R. Li, T. Medina, R. Molavi, R. Neufeld, T. Oh, I. Pavlov, I. Perminov, G. Poulin-Lamarre, C. Rich, A. Smirnov, L. Swenson, N. Tsai, M. Volkmann, J. Whittaker, and J. Yao, “Phase transitions in a programmable quantum spin glass simulator,” Science **361**, 162–165 (2018). [CrossRef] [PubMed]

**18. **S. Sachdev, *Quantum phase transitions* (Wiley Online Library, 2007).

**19. **H. T. Quan, Z. Song, X. F. Liu, P. Zanardi, and C. P. Sun, “Decay of loschmidt echo enhanced by quantum criticality,” Phys. Rev. Lett. **96**, 140604 (2006). [CrossRef] [PubMed]

**20. **E. Lieb, T. Schultz, and D. Mattis, “Two soluble models of an antiferromagnetic chain,” Annals Phys. **16**, 407–466 (1961). [CrossRef]

**21. **H. J. Lipkin, N. Meshkov, and A. Glick, “Validity of many-body approximation methods for a solvable model:(i). exact solutions and perturbation theory,” Nucl. Phys. **62**, 188–198 (1965). [CrossRef]

**22. **N. Meshkov, A. Glick, and H. Lipkin, “Validity of many-body approximation methods for a solvable model:(ii). linearization procedures,” Nucl. Phys. **62**, 199–210 (1965). [CrossRef]

**23. **A. Glick, H. Lipkin, and N. Meshkov, “Validity of many-body approximation methods for a solvable model:(iii). diagram summations,” Nucl. Phys. **62**, 211–224 (1965). [CrossRef]

**24. **K. Hepp and E. H. Lieb, “On the superradiant phase transition for molecules in a quantized radiation field: the dicke maser model,” Annals Phys. **76**, 360–404 (1973). [CrossRef]

**25. **Y. K. Wang and F. Hioe, “Phase transition in the dicke model of superradiance,” Phys. Rev. A **7**, 831 (1973). [CrossRef]

**26. **C. F. Lee and N. F. Johnson, “First-order superradiant phase transitions in a multiqubit cavity system,” Phys. review letters **93**, 083001 (2004). [CrossRef]

**27. **A. A. Ovchinnikov, D. V. Dmitriev, V. Y. Krivnov, and V. O. Cheranovskii, “Antiferromagnetic ising chain in a mixed transverse and longitudinal magnetic field,” Phys. Rev. B **68**, 214406 (2003). [CrossRef]

**28. **J. Vidal, R. Mosseri, and J. Dukelsky, “Entanglement in a first-order quantum phase transition,” Phys. Rev. A **69**, 054101 (2004). [CrossRef]

**29. **L. Del Re, M. Fabrizio, and E. Tosatti, “Nonequilibrium and nonhomogeneous phenomena around a first-order quantum phase transition,” Phys. Rev. B **93**, 125131 (2016). [CrossRef]

**30. **S. Gammelmark and K. Mølmer, “Phase transitions and heisenberg limited metrology in an ising chain interacting with a single-mode cavity field,” New J. Phys. **13**, 053035 (2011). [CrossRef]

**31. **M. Raghunandan, J. Wrachtrup, and H. Weimer, “High-density quantum sensing with dissipative first order transitions,” Phys. Rev. Lett. **120**, 150501 (2018). [CrossRef] [PubMed]

**32. **R. H. Dicke, “Coherence in spontaneous radiation processes,” Phys. Rev. **93**, 99 (1954). [CrossRef]

**33. **V. S. Zapasskii, “Spin-noise spectroscopy: from proof of principle to applications,” Adv. Opt. Photonics **5**, 131–168 (2013). [CrossRef]

**34. **P. Pfeuty, “The one-dimensional ising model with a transverse field,” ANNALS Phys. **57**, 79–90 (1970). [CrossRef]

**35. **L.-P. Yang and Z. Jacob, Engineering Quantum Phase Transitions for Weak Signal Detection (2018). Under preparation.

**36. **Y. Zhang, L. Yu, J.-Q. Liang, G. Chen, S. Jia, and F. Nori, “Quantum phases in circuit qed with a superconducting qubit array,” Sci. reports **4**, 4083 (2014). [CrossRef]

**37. **Y. Imry, “Finite-size rounding of a first-order phase transition,” Phys. Rev. B **21**, 2042 (1980). [CrossRef]

**38. **M. Skotiniotis, P. Sekatski, and W. Dür, “Quantum metrology for the ising hamiltonian with transverse magnetic field,” New J. Phys. **17**, 073032 (2015). [CrossRef]

**39. **J. Dziarmaga, “Dynamics of a quantum phase transition and relaxation to a steady state,” Adv. Phys. **59**, 1063–1189 (2010). [CrossRef]

**40. **L.-P. Yang, H. X. Tang, and Z. Jacob, “Concept of quantum timing jitter and non-markovian limits in single-photon detection,” Phys. Rev. A **97**, 013833 (2018). [CrossRef]

**41. **H. Yuen, “States that give the maximum signal-to-quantum noise ratio for a fixed energy,” Phys. Lett. A **56**, 105–106 (1976). [CrossRef]

**42. **K. Baumann, C. Guerlin, F. Brennecke, and T. Esslinger, “Dicke quantum phase transition with a superfluid gas in an optical cavity,” Nature **464**, 1301 (2010). [CrossRef] [PubMed]

**43. **L.-P. Yang, Y. Li, and C. Sun, “Franck-condon effect in central spin system,” The Eur. Phys. J. D **66**, 300 (2012). [CrossRef]

**44. **K. Husimi, “Some formal properties of the density matrix,” Proc. Physico-Mathematical Soc. Jpn. 3rd Ser. **22**, 264–314 (1940).

**45. **C. T. Lee, “q representation of the atomic coherent states and the origin of fluctuations in superfluorescence,” Phys. Rev. A **30**, 3308–3310 (1984). [CrossRef]

**46. **J. Radcliffe, “Some properties of coherent spin states,” J. Phys. A: Gen. Phys. **4**, 313 (1971). [CrossRef]

**47. **F. Arecchi, E. Courtens, R. Gilmore, and H. Thomas, “Atomic coherent states in quantum optics,” Phys. Rev. A **6**, 2211 (1972). [CrossRef]

**48. **A. E. Lita, A. J. Miller, and S. W. Nam, “Counting near-infrared single-photons with 95% efficiency,” Opt. Express **16**, 3032–3040 (2008). [CrossRef] [PubMed]

**49. **B. Korzh, Q. Zhao, S. Frasca, J. Allmaras, T. Autry, E. Bersin, M. Colangelo, G. Crouch, A. Dane, and T. Gerrits *et al.*, “Demonstrating sub-3 ps temporal resolution in a superconducting nanowire single-photon detector,” arXiv preprint arXiv:1804.06839 (2018).

**50. **C. Schuck, W. H. Pernice, and H. X. Tang, “Waveguide integrated low noise nbtin nanowire single-photon detectors with milli-hz dark count rate,” Sci. Reports **3**, 1893 (2013). [CrossRef]

**51. **B. Q. Baragiola, R. L. Cook, A. M. Brańczyk, and J. Combes, “N-photon wave packets interacting with an arbitrary quantum system,” Phys. Rev. A **86**, 013811 (2012). [CrossRef]