Quantum correlation of microwave two-mode squeezed state generated by nonlinearity of InP HEMT

This study significantly concentrates on cryogenic InP HEMT high-frequency circuit analysis using quantum theory to find how the transistor nonlinearity can affect the quantum correlation of the modes generated. Firstly, the total Hamiltonian of the circuit is derived, and the dynamic equation of the motion contributed is examined using the Heisenberg-Langevin equation. Using the nonlinear Hamiltonian, some components are attached to the intrinsic internal circuit of InP HEMT to address the circuit characteristics fully. The components attached are arisen due to the nonlinearity effects. As a result, the theoretical calculations show that the states generated in the circuit are mixed, and no pure state is produced. Accordingly, the modified circuit generates the two-mode squeezed thermal state, which means one can focus on calculating the Gaussian quantum discord to evaluate quantum correlation. It is also found that the nonlinearity factors (addressed as the nonlinear components in the circuit) can intensely influence the squeezed thermal state by which the quantum discord is changed. Finally, as the primary point, it is concluded that although it is possible to enhance the quantum correlation between modes by engineering the nonlinear components; however, attaining quantum discord greater than unity, entangled microwave photons, seems a challenging task since InP HEMT operates at 4.2 K.


Introduction
Quantum correlation as a fundamental issue is significantly represented in quantum applications such as quantum information [1][2][3][4][5][6][7] and quantum sensors [8][9][10][11][12].The entanglement has been synonymously applied in the same way with quantum correlation in most applications since the quantum system contains only pure states [5].In contrast, when the mixed states are generated by a quantum system, which covers most quantum systems, such as quantum radars [8,10], the term entanglement cannot be used over quantum correlation.It is because some separable mixed states can introduce "residual correlation" that cannot be fitted by any classical probability distributions [5][6][7].In other words, for a bipartite system with separable state ρAB expressed as ρAB = ∑Pi ρAi × ρBi, where ρAi and ρBi are density matrices of the subsystems; the states ρAi and ρBi may be physically non-distinguishable.Consequently, all information about the subsystems cannot be locally retrieved due to the nonorthogonality of the states [7].Some recently published studies have shown that classically correlated states might show the signature of quantumness [13][14].To evaluate the quantumness, the term "quantum discord" has been applied [1][2][3][4][5][6][7].The quantum discord can quantify the residual correlation or the signature of quantumness [5][6][7], which captures all quantum correlations in a bipartite state.However, several different quantifiers have been introduced to evaluate the nonclassicality of the correlation (quantum correlation); but the most popular one is quantum discord [15].Quantum discord is the difference between the quantum correlation within a quantum state and its classical correlations [5][6][7].Some studies show that quantum discord can be defied for both qubits and continuous variable systems [7,16].Therefore, studying the quantum discord for the continuous variables is valuable due to the critical applications of the continuous variables, such as in quantum computation and quantum communications [18].For some reason mentioned above, and since the circuit designed in this work interacts with the environment as any real quantum systems interact inevitably with the surrounding, we have to focus on the continuous variable and calculate the quantum discord for these variables.As the main point, the influence of the thermal noise generated by the circuit is studied on the quantum correlation.
Thus, concerning the mentioned points, the study significantly focuses on the "Gaussian quantum discord" and uses the close formula satisfied to all families of the Gaussian state.Accordingly, the Gaussian state includes the important class of the squeezed thermal state.This unique state is realized by applying twomode squeezing to a pair of single-mode thermal states [5,7].Recent studies have shown that the squeezed thermal state, generally the Gaussian state, can be decomposed into the EPR (Einstien-Podolsky-Rosen) state plus the local action of a "phase-sensitive Gaussian channel" [5].
As mentioned, this study significantly focuses on the two-mode squeezing thermally state, which is generated by the nonlinearity of the InP HEMT.It has to be initially shown that due to the nonlinearity of InP HEMT in the circuit, the states generated in the system are mixed.This means that the quantum discord as a quantifier can completely define the quantum correlation rather than the entanglement.Then, it will be verified that the introduced nonlinearity in the circuit generates the microwave two-mode squeezed thermal state.Finally, it introduces some critical factors (components) relating to nonlinearity to enhance the quantum correlation between modes.The mentioned components, as the circuit elements, are attached to the original internal circuit of InP HEMT to make the modified circuit; Then, the modified circuit is simulated using PySpice.The simulation in PySpice is performed to check some critical issues from a circuit analysis point of view.Consequently, it will show that the attached nonlinearity can add a challenging trade-off.
It is better to indicate that this study completes the latter similar work [17], in which, for simplicity, it was supposed that all states produced by the InP HEMT are pure states.Another difference between the present work with [17] is that this study concentrates on the microwave two-mode squeezed thermal state and the squeezing parameters that affect the quantum correlation.Additionally, in this study, the thermally excited photons in the Drain due to the related conductance are introduced at 450 K (Drain conductance noise temperature), making the circuit behave like an accurate model.
The present study initially defines the system and its crucial elements that can affect the primary goal.In the next step, the Hamiltonian of the system is derived, and using the nonlinear part of the Hamiltonian, the dynamics equation of the motion is introduced.In addition, it is shown that the circuit can generate a microwave two-mode squeezed thermal state by applying some nonlinear coefficients.In the other step, it is theoretically shown that the introduced circuit generates the mixed states.Finally, the theory related to quantum discord and its deriving is introduced and discusses by which parameters or factors one can manipulate the quantum discord.

System definition
The circuit studied in this work is schematically shown in Fig. 1.It contains InP HEMT as a nonlinear active element, an input and output matching network to match the input and output impedance, and a DC stabilization circuit.The generic equivalent circuits of the input and output matching network are shown in the inset figures.In the inset figures, it is shown in a usual way that any transmission line can be modeled with an equivalent lumped element circuit.As shown in the schematic, InP HEMT is biased via Vg and Vd to operate in the desired region, and the small signal RF input wave is applied to the circuit through an input capacitor (Cin).Lg1 and Ld1 are key elements in the stabilization circuit.As an essential element, InP HEMT's nonlinear equivalent circuit is attached as another inset figure.The attached inset tries to completely show all factors that could affect the operation of InP HEMT at cryogenic temperature; it is mainly for quantum applications in which the number of microwave photons is dramatically kept at a low level.For the reason that the thermally excited photons can affect the low photons quantum applications, the noise generated by resistances in the circuit is considered.This makes the circuit behave like a real one.
For instance, 4KTBnRs is the voltage-noise generated by Rs in the circuit, in which K, T, and Bn are, respectively, Boltzmann constant, operational temperature, and noise bandwidth.
One of the essential components in the inset circuit is ids, defined as the dependent current source controlled by the voltage.As can be seen from the related relationship, the amount of the current is controlled by gm, defined as the intrinsic transconductance of the circuit; that is additionally manipulated by gm2 and gm3 called the second and third nonlinearity factors (generally called nonlinearity factors).This paper focuses on the latter factors and their effect on the circuit, through which the circuit's resonance frequency is changed.In addition, the coupling signals between the resonators in the circuit can be strongly changed.
The term "coupling between the resonators" is used because, from a general view, the circuit illustrated in Fig. 1 can be supposed as the two separate oscillators oscillating in the gate and drain side of the transistor, which are coupled to each other through the InP HEMT nonlinear circuit.In other words, the nonlinearity in InP HEMT changes the features of the oscillators, including the resonance frequency and their impedance.
Fig. 1 the schematic of the circuit containing the input and output matching networks, stability network, and active elements (InP HEMT) and its internal circuit.

Hamiltonian of the circuit
To derive the classical Hamiltonian, one has to use Legendre transformation as H(φk,Qk) = ∑k (φk.Qk) -Lc, where Lc and Qk are, respectively, the system's Lagrangian and conjugate variables of the "coordinate variable" φk, calculated through Qk = ∂Lc/∂(∂φk/∂t) [17,20].The total Hamiltonian of the circuit introduced in Fig. 1 is given by:   ( ) ( ) where CA = Cin + C1 + Cgs + Cgd, CB = Cgd + C2, and Cc = Cgd.The total Hamiltonian of the system is written as Ht = HL + HN, where HL stands for linear Hamiltonian and HN contains the nonlinear terms.The nonlinearity is arisen due to the terms inside the curly bracket.This study looks to generate two-mode squeezed thermal states using the nonlinearity introduced by InP HEMT.For this reason, the focus is just laid on HN and its analysis, using which I. it will theoretically show that the generated states by the circuit presented in Fig. 1 are mixed, and II.there is the nonlinearity effect that can generate the two-mode squeezed thermal state.In contrast, HL is used to define the steady-state operational point, energy level, and also the linear part of Hamiltonian has some terms that can generate the coherent state.The definition of HL and its parameters have been introduced in Appendix A (Eq. A1 and A2).Consequently, the nonlinear Hamiltonian is given by: where gN2 = gm2 + 6gm3[∂φ1/∂t]DC, CM 2 = CB(CA+CN)-Cc To clarify the effects of the nonlinearity, some components created because of the InP HEMT nonlinearity are attached to the main circuit (inset figure in Fig. 1) and illustrated in Fig. 2.These elements are theoretically derived using nonlinear Hamiltonian expressed in Eq. 2. The red scribble-dotted line on the circuit contains a variable capacitor (CN) and controllable current depending on the nonlinearity factors of the InP HEMT.The circuit illustrated in Fig. 2 is an approximated circuit that may model the InP HEMT transistor behavior.
The present study uses this model to quantum mechanically analyze the InP HEMT nonlinearity effect on the quantum correlation that may occur between modes generated in the circuit at cryogenic temperature.
As clearly seen in Fig. 2, CN and ids_N directly affect the second oscillator resonance properties; nonetheless, the nonlinear effects are coupled to the first oscillator through Cgd.This means that the coupling properties of the oscillators are strongly affected by the nonlinear properties.For example, the cross-correlation between the coupled oscillator's modes is severely influenced by the nonlinearity factors defined in the circuit.In other words, the classicality correlated modes generated by the coupled oscillators can be affected in such a way as to show quantumness.This phenomenon strongly depends on nonlinearity and its impacts.
As an important quantifier, the quantum discord [5][6][7] is selected to evaluate the quantum correlation between the states generated in the circuit.The quantum discord rather than the entanglement is chosen to show the effects because the circuit discussed in this study shows mixed states.This point will be explored in the next section.In the complementary part of this study, the circuit shown in Fig. 2 is analyzed using PySpice to investigate the effect of the attached components on InP HEMT's main feature, such as DC characterization.Using PySpice, rather than any specific CAD simulator tools, gives some degree of freedom by which the designer can easily attach any elements or sub-elements to the circuit and make the desired pack.We also selected PySpice to work with because all theoretical simulations to calculate the mean photons number, quantum discord, the smaller Symplectic eigenvalue, and quantum mutual information were done in Python.
It should be noted that the general aim of the attached elements to the circuit is to study their effects on the performance of InP HEMT operating at 4.
where ak and ak + (k = 1, 2) are the annihilation and creation operators, respectively.For simplicity, a few constants are defined by which Eq. 3 is simplified as: i g a a a a i g a a a a i g a a a a a a where gN11, gN21, gN31, gN41, gN51, and gN61, are the related constants defined in Appendix A (Eq. A3).These constants are strongly dependent on the nonlinearity factor.The dynamics equation of motion of the circuit is calculated using the Heisenberg-Langevin equation; A simple way to achieve a stationary and robust calculation in continuous modes to calculate quantum correlation, i.e. entanglement [28][29] or quantum discord [9], is to select a constant point that the system is driven and works with.Since the interaction field is so strong, it is appropriate to focus on linearization and calculate the quantum fluctuation around the semi-classical constant point.For linearization, the oscillator modes are expressed as the summation of the stationary (constant) and fluctuating parts as a1 = A1+δa1 and a2 = A2+δa2, where the capital letter (A1 and A2) denotes the system's steady-state points, and δ indicates the fluctuation around the steady-state point.
The steady-state points related to the circuit are calculated and presented in Appendix A (Eq. A7).Thus, the linearized equations around the steady-state points are given by: where γa11, γa12, γa13, γa21, γa22, γa23, and γa24 are constant rates depending on the stationary points of the system.These rates can be complex numbers defined in Appendix A (Eq. A4).In Eq. 5, κ1, κ2, ω1, and ω2 are the first and second oscillators' decay rates and the contributed frequencies, respectively.Additionally, δain_1 and δain_2 are the input noises fluctuation.They obey the correlation function , in which kB and T are the Boltzmann's constant and operational temperature, respectively [9,17,20].N(ω) is the equilibrium mean of the thermal photon numbers at the frequency ω.Finally, capital δ expressed in δ(s-s') is Dirac's function.
In the following, Eq. 5 is transformed to the Fourier domain to simplify the algebra in the frequency domain to calculate the quantum discord, which is introduced as: where ∆1 and ∆2 are the oscillators detuning frequencies, which are calculated as ∆1 = ω -ω1 and ∆2 = ω -ω2, where ω is the RF incident frequency.Notably, Eq. 6 is used to calculate the < δa1 + δa1>, < δa2 + δa2>, and < δa1δa2> as the mean photon number of the first-, second-oscillator, and the phase-sensitive crosscorrelation [17][18][19][20].

Mixed states generation because of the nonlinearity of InP HEMT
This part shows that the state of the oscillators is dispersed, meaning that all states generated by the circuit become mixed states.One can utilize the first perturbation theory to calculate the oscillator's state changes.
The change of a typical state is calculated using the first perturbation theory by |j> (1) = ∑i≠j {(<i|HN|j>)/(Ei-Ej)}×|j> [20], where |j> is the pure state of the oscillators, and Ei and Ej are the contributed energies.|j> (1)   is the final state of the first oscillator that may differ from |j>; it depends on the system and the related Hamiltonian.The results of the calculation are presented in Eq. 7 as: ( ) In this equation, the data related to the constants used are given in Appendix A (A6).This equation shows that the final states related to the first and second oscillators |j1> (1) and |j2> (1) , are strongly influenced by the nonlinear (perturbation) Hamiltonian.For example, the state of the first oscillator is coupled to |j1±1> and | j±2> due to the nonlinearity effect.For instance, if one fixed the first oscillator state at |0> as a pure state, the final state of this oscillator will be found in the superposition state of |1> and |2>, meaning that HN causes the final states to be mixed.In other words, the state |0> is mixed with |1> and |2>.
In addition, it is necessary to calculate the energy of the oscillators using Ej = <ji|H0+HN|ji> [20], where ji = 1,2.The oscillators associated energies due to the total Hamiltonian are given by: From Eq. 8, it is clear that the first term is arisen due to the linear Hamiltonian, and the terms inside the curly bracket are generated because of the nonlinear Hamiltonian effect.It is called the perturbation effect on the energy levels of the oscillators.In this equation, A1 and A2 are the steady-state points (DC points) of LC1 and LC2, where the oscillators are designed to operate on.The DC points can be calculated using Heisenberg-Langevin equations in the steady-state [15,16]; these points are essentially affected by the circuit's DC bias and the linear Hamiltonian effect.The DC points related to the circuit are calculated and presented in Appendix A (Eq. A7).It is clear from the mentioned equation that the circuit's steady state points strongly depend on the circuit's generated noise (shown in Fig. 1), and the DC bias point.As an essential result of this research, it is found that the quantities such as g12 ' , g22 ' , and Cq1q2 handling the displacement in the circuit (linear Hamiltonian) strongly affect the steady-state points.This gives any engineer a degree of freedom to control and manipulate the steady-state point in which the circuit is established to be operated.

Two-mode squeezed thermal state
A squeezed-coherent state is created via the acting of the squeezed and displacement operators on the vacuum state defined as |α,ζ> = D(α)S(ζ)|0>, where |0> is the vacuum state [20,25].In this study, it can be easily shown that the coherent state can be generated by the linear part of the total Hamiltonian.In contrast, a squeezed state for each oscillator needs quadratic terms such as ai 2 and ai +2 in the exponent.Nonetheless, for two-mode squeezing, the Hamiltonian should have the terms like {a1a2 + a1 + a2 + } [20][21][22][23][24][25].Therefore, the two-mode squeezed state becomes analyzed by the evolution of exp[HNt/iħ], where HN is defined in Eq. 4.
Based on this definition, any quadratic terms like {a1a2 + a1 + a2 + } in the Hamiltonian may generate twomode squeezing.Consequently, using exp[HNt/iħ], the two-mode squeezing state of the nonlinear circuit studied is presented by: Eq. 9 shows that the coupled oscillator can generate two-mode squeezing through the transistor's nonlinearity.That means that the nonlinearity created by the transistor couples two oscillators so that the coupled two-mode become squeezed.The two-mode squeezing parameter is defined as ζ12 = -2gN11Im(A1) + 2gN41Re(A2) + 2gN61Re(A1), where Re{} and Im{} indicate the real and imaginary parts, respectively.It is apparent from the relationship that ζ12 can be a complex number, which means that the two-mode squeezing parameter contains a phase, which determines the angle of the quadrature.
It should be noted that "t" in the exponent (exp[Ht/jħ]) can be determined from t < min{1/κ1, 1/κ2}, where κ1 and κ2 are the first and second oscillators' decay rates.By selecting "t", the system is forced to generate squeezing before the resonator decaying [25].

Quantum discord
The calculation of the quantum discord initiates using two generalizations of the classical mutual information [1][2][3][4][5][6][7].The mutual information is used primarily to evaluate the total correlations between subsystems [15].In the first generalization, the quantum mutual information for two systems, A and B, is defined as I(ρAB) = S(ρA) + S(ρB) -S(ρAB), where S(ρA) = -Tr(ρAlog2ρA) is the von Neumann entropy of system A, and S(ρAB) is the conditional von Neumann entropy [5][6][7].In the following, for the system studied in this work, the first oscillator is shortly called A, and the second one is called B. The conditional entropy arises because the measurement process disturbs the state on which a physical system is set.In other words, the applied measurement on subsystem B may change the state of subsystem A.
The second generalization introduces the entropic quantity C(A|B), by which the classical correlation in the joint state ρAB is calculated.The classical correlation defines the maximum information about one subsystem depending on the measurement types applied to the other subsystem.The entropic quantity, by considering the generalization, is defined as C(ρAB) = S(ρA) -Smin(ρAB), where the only difference between mutual quantum information is Smin(ρAB).This term is the conditional minimized entropy of system A over all possible measurements on system B. That is generally described as positive operator valued measures (POVMs) [5][6][7].Thus, quantum discord is usually defined as D(ρAB) = I(ρAB) -C(ρAB), and by substituting from the above definitions, it becomes D(ρAB) = S(ρB) -S(ρAB) + Smin(A|B).
The output entropy associated with the heterodyne detection is equal to the average entropy of the output ensemble A [5].Since entropy is invariant under displacements, it may write S(A|hetB) = h(τ + η), where h(x) ≡ (x+0.5)log2(x+0.5)-(x-0.5)log2(x-0.5).In other words, the von Neumann entropy of an n-mode Gaussian state with CM expressed in Eq. 10 is calculated as S(VAB) = ∑i=1 N f(νi), where νi are the related Symplectic eigenvalues [6].However, there is a heterodyne detection for which it is optimal for minimization of the output entropy, so the Gaussian discord is optimal.The "Gaussian quantum discord" of a two-mode Gaussian state ρAB, assuming a two-mode squeezed thermal state in this study, can be defined as the quantum discord satisfying the conditional entropy restricted to the generalized Gaussian POVMs on B [6].As a result, in the system defined, the second oscillator entropy and the off-diagonal elements in the CM significantly affect the system's quantum discord.For this reason, in this article, we specially pay attention to the interferences between the oscillators and determine which quantities can affect the quantum discord.Perhaps the other critical factor that may be considered to enhance the quantum discord is h(τ + η), by which the classical correlation is decreased.As mentioned in [5][6][7], this critical factor strongly depends on the type of measurements.

Results and Discussions
In the following, the quantum correlation generated between the two-mode squeezing state by the circuit shown in Fig. 1 is discussed.It attempted to focus on the crucial parameters in the circuit (Fig. 2) that one can manipulate in such a way as to enhance the quantum correlation between modes.As shown in Fig. 2, CN and gN2 are two crucial factors, mainly the function of gm2 and gm3, that the study concentrates on them to manipulate the quantum correlation.Notably, it was theoretically shown the latter section that gN2 is the main factor by which the final states of the oscillators became mixed and also plays a significant role in generating the two-mode squeezing thermal state.
Using the information mentioned above, in the following, we attempt to show that the quantum correlation can be generated between the microwave photons created by the two coupled oscillators.As a quantifier, the quantum discord is calculated, a property almost all quantum states hold.However, other quantities, such as the quantum information and classical correlation generated between modes in the system, are analyzed.The author thinks that comparing the quantum discord, quantum information, and classical correlation gives a solid sense for anyone to find which parameter can strongly enhance the quantum correlation.In addition, it makes it possible to know what portion of the classical correlation may appear in the quantum correlation.This means that all classically correlated states don't retain signatures of quantumness.
Before discussing the simulation results, it is necessary to care about a more severe and crucial point, which is InP HEMT operating at 4.2 K; this makes an astonishing thermally exciting noise in the circuit.Also, the internal circuit effects cause the quantum discord amplitude to remain limited.Thus, it seems impossible to compare the results of this study with the quantum discord generated by a typical opto-mechanicalmicrowave system's quantum discord [9].It should be noted that the recent system works at 10 mK, whereas InP HEMT operates at 4.2 K; this means that the thermally excited photons dramatically affect the quantum correlation.
In this study, one of the aims is to design a circuit with an operating frequency (finc = f1+f2) of around 12 GHz, where finc, f1, and f2 are the RF incident frequency, the first and second oscillator's frequency, respectively.In the range of mentioned frequency and operating at 4.2 K, the mean thermally exciting photons were calculated as nth1 ≡ < δain_1 + δain_1 > or nth2 ≡ < δain_2 + δain_2 > ~ 12, which is so greater than the average number of photons generated by the first and second oscillators n1 ≡ < δa1 + δa1> or n2 ≡ < δa2 + δa2> ~ 0.2.With knowledge of the thermally excited photons as the inserting noise and also the oscillators' mean photon number, in the following, the Gaussian quantum discord of the two-modes Gaussian states generated by the coupled oscillators is calculated.All simulations in this work containing the theoretical simulation of quantum correlation and also the simulation of the nonlinearity effects on the InP HEMT DC characteristics in PySpice were performed using the data from Table .1.The simulation results for quantum discord, the smaller Symplectic eigenvalue, quantum mutual information, and classical correlation at 4.2 K are shown in Fig. 3.The result in Fig. 3a indicates that increasing nonlinearity factor, gN2, significantly enhances the quantum correlation.The graph shows that the quantum discord is strongly improved in the two different frequency ranges.It means that there is an avoided-level crossing in the graph (like a gap), where the quantum correlation between modes vanishes.
According to the demonstrated figure, in some regions, D(ρAB) fluctuates around zero.It has been shown that a state with zero quantum discord represents a classical probability, while a positive discord indicates quantumness even in a separable mixed state [6,7].Therefore, for 0<D(ρAB)<1, the states may be entangled or separable, but for D(ρAB) >1, the states are entangled.
Generally, using the results illustrated, this paper wants to answer a few important questions: 1.Is it possible to generate quantum correlation using InP HEMT nonlinearity operating at 4.2 K? 2. By what factor the nonlinearity generated by the InP HEMT can affect the quantum correlation?3. Is it possible to generate entangled microwave photons at 4.2 K? The Gaussian quantum A discord is calculated to answer the mentioned questions, and the measurement is performed on the second oscillator to minimize the classical correlation.Fig. 3a demonstrates that the amplitude of quantum A discord is less than unity 0< D(ρAB)<1, which means that the states generated by the oscillators may have either separable or entangled states.Consequently, if one applies a real InP HEMT operating at 4.2 K and considers all the noises and nonlinearity effects (Fig. 2), the generation of the entangled microwave photons seems impossible.This is mainly because more thermally excited photons are generated in microwave frequencies.Nonetheless, it is possible to create a partial quantum correlation between the two-mode Gaussian states (mixed states in this study).An inset figure is attached to Fig. 3a to show the critical points where the quantum discord is maximized.However, the interesting fact about this graph is the gap at which the quantum discord is minimized or completely disappeared.To discuss this odd behavior, the smaller Symplectic eigenvalue (đ-) of the partially transposed state is illustrated in Fig. 3b.For better illustration, the logarithmic scale of the quantity (log10(đ-)) is shown, at which two inset figures are attached.The minimum value for đis around 1.1 (log10(1.1)= 0.095).It has been shown that the Gaussian state is entangled if and only if đ-< 0.5 [7] (đ-< 1 [6]).As a result, the curve illustrated in Fig. 3b shows that the state generated by the circuit is not entangled because d-is always greater than 1.1 for each gN2.In addition, the study focuses on the smaller Symplectic eigenvalue because there is a pure consistency between the two figures illustrated in Fig. 3a and 3b.It can be found that the more quantum discord the circuit can generate, the less đthe circuit can produce.
However, we learned from the latter section that some other quantifiers could help to elucidate the quantum discord's behavior.In other words, the difference between quantum discord, quantum mutual information, and classical correlation was mathematically discussed in the latter section.Thus, to thoroughly analyze the quantum discord, the simulation results of the quantum information and classical correlation are illustrated in Fig. 3c and Fig. 3d.The comparison between the quantum mutual information and classical correlation reveals that at the gap, where the quantum discord is minimized, the quantum mutual information and classical correlation are maximized.This means that at the gap mentioned, the modes totally become separable, and there is no quantumness between them.To discuss more technically, one may consider the definition of the Gaussian quantum discord, one-way classical correlation, and quantum mutual information The same effect could be found in the quantum mutual information graph.x-axis and y-axis are annexed in the figures.The other graph in Fig. 4b shows the operating temperature effect on the quantum A discord.We suppose that if InP HEMT could operate at 1.2 K, what happens on the quantum A discord? Fig. 3b shows a significant enhancement in quantum discord amplitude.This directly contributes to the level of the thermally excited noise generations, which dramatically decreases as the temperature drop from 4.2 K to 1.2 K.The other interesting point is that by reducing the temperature to around 1.2 K, at lower gN2 the quantum discord amplitude is increased.In the following, we will study that increasing gN2 can dramatically change the InP HEMT DC characterization, by which, for instance, the power dissipation can be severely increased.In other words, manipulating gN2 to enhance the quantum discord adds a crucial trade-off to RF circuit engineering.So far, the study shows that by engineering some parameters related to the internal circuit of the InP HEMT specially focused on the nonlinearity factors, it was possible to enhance the quantum correlation between the two-mode Gaussian states.Initially, it was theoretically shown that the states generated by the InP HEMT nonlinear circuit became mixed.The results showed that the quantum discord is increased for gN2 between 0.8-1.2A/V 2 , meaning that the quantum correlation may be created between the states.Finally, arising a trade-off in the design was discussed, which means that one cannot increase gN2 to any value to get the desired quantum correlation.It is because the nonlinearity factors can dramatically change the modified circuit performance (DC and AC characteristics).In other words, by changing the nonlinearity factors such as gm2 and gm3, for example, the DC characteristics of the InP HEMT can be changed.The mentioned variations contribute to the circuit's components changing, such as CN and gN2.To show this point, the circuit illustrated in Fig. 2 is simulated in PySpice (Jupyter Lab) using the data from Table .1.We simulated the internal circuit of InP HEMT, and the related DC characterization is shown in Fig. 5.The left curve as a blue graph shows the InP HEMT drain-source current (Ids) as the function of gN2, whereas the right graph (red) displays the CN as the function of gN2.In the region shown between the dotted line, where the gN2 is changed between 0.8-1.2,Ids is increased to ~30 mA, and CN is reached around 50 pF.It is better to note that Ids is altered in the range of (0.5-2.3) mA for sub-mW cryogenic applications [26].The mentioned trade-off becomes crucially created, and the designers should especially care about this point.
In other words, increasing the current to 30 mA may generate the quantum correlation between states.Still, it forces the system to dissipate energy dramatically, which is a crucial point in cryogenic applications.
As an essential conclusion of this work, it can be mentioned that InP HEMT nonlinearity has the potential to partially generate the quantum correlation between modes of the oscillators coupled through the transistor.Some critical factors play a central role, such as the nonlinear capacitor arising due to the nonlinearity CN and the nonlinearity factor gN2 affected by gm2 and gm3.These parameters give engineers a reasonable degree of freedom to effectively design a cryogenic circuit containing InP HEMT by which the generation of the quantum correlation between modes at 4.2 becomes possible.However, the generation of entangled microwave photons by InP HEMT operating at 4.2 K seems impossible, at least with the recent technology.

Conclusions
This study investigated the quantum correlation of the microwave two-mode squeezed thermal state generated by the nonlinearity of InP HEMT.For this reason, the nonlinear circuit related to InP HEMT was analyzed using quantum theory, and the contributed dynamics equation of the motion of the circuit was theoretically derived.In this work, we intensely concentrated on the nonlinear Hamiltonian by which it was theoretically shown that the two-mode squeezed thermal state was generated.In addition, for the circuit discussed, it was theoretically proved that the circuit just generated the mixed state, and there was no pure state.Because of the facts mentioned, the study focused on quantum discord rather than quantum entanglement.To completely know about the quantum discord, other quantities such as quantum mutual information, classical correlation, and the smaller Symplectic eigenvalue were analyzed.
Some engineering was carried out on the nonlinear circuit to get the desired results, and some critical parameters, such as gm2 and gm3, were manipulated to enhance the quantum correlation between the generated modes.As an important conclusion, the simulation result showed that although it might be possible to improve the quantum correlation between modes, attaining quantum discord greater than unity seems challenging when InP HEMT operated at 4.2 K.
Using the standard approach to solve the equations, the final form of the variables become: where A1 = A1R+iA1j and A2 = A2R+iA2j.Using Eq.A7 the steady state point of the circuit, which is related to the thermal noise and bias point of the circuit, is calculated.

Fig. 2
Fig.2 An approximation of the InP HEMT internal circuit by considering the nonlinearity effects illustration with a scribble-dotted line on the primary circuit.

Finally, the compact
form of quantum discord, classical correlation, and quantum mutual information are given, respectively, by D(ρAB) = h(b) -h (ν-) -h (ν+) + h (τ + η), C(ρAB) = h(a) -h (τ + η), and I(ρAB) = h(a) -h (ν-) -h (ν+), where ν± is the Symplectic eigenvalue of the CM.The Symplectic eigenvalues are defined as ν± = [∆±√(∆ 2 -4D)]/2, where ∆ = det(aI) + det(bI)+2det(do12I) and D = det(VAB); in the recent formulas det{:} stands for the matrix determinant.In the compact form of the equation defined for the quantum discord, the first term stands for the von Neumann entropy of the second oscillator in the system.The second and third terms define the von Neumann conditional entropy of the system.The last term in the equation is the effect of the classical correlation depending on the type of measurement performed on the second oscillator.
as D(ρAB) = h(b) -h (ν-) -h (ν+) + h (τ + η), C(ρAB) = h(a) -h (τ + η), and I(ρAB) = h(a) -h (ν-) -h (ν+)[5-  7].The common point between D(ρAB) and I(ρAB) is -[h (ν-)+h (ν+)], which relates to the von Neumann entropy of the CM; however, the difference between them contributes to the oscillator's von Neumann entropy (h(b) and h(a)), and h (τ + η) which just appears in D(ρAB).In Fact, the term h (τ + η) comes from the classical correlation effect, which minimizes the classical correlation to maximize the quantum correlation.Using the mentioned points and comparing Fig.3awith Fig.3c, one can find that the quantum discord is maximized where the quantum mutual information is minimized.The logarithmic inset graphs in Fig.3c(Fig.3c1andFig.3c2) shows clearly the points.The main difference between the D(ρAB) and C(ρAB) is the last term in the expressions h (τ + η), denoting the measurement effects on the second oscillators.The classical correlation is the maximum information about one subsystem, which depends on the measurement performed on the other subsystem.The result of the classical correlation is demonstrated in Fig.3d.Also, some inset figures in the logarithmic scale are attached to clearly demonstrate the classical correlation as the function of circuit nonlinearity and frequency.This figure also shows that the classical correlation is maximized at the gap where the quantum discord disappeared.This effect contributes to the influence of the first oscillator's von Neumann entropy.

Fig. 4
Fig. 4 quantum discord vs. angular frequency (GHz) and nonlinear factor gN2 (A/V 2 ).a) C2 = 0.5 pF, T = 4.2 K, b) T = 1.2 K; C2 = 0.5 pF.It mentioned that this study calculates the Gaussian quantum A discord, meaning that the measurement is performed on the second oscillator.As a result, if one wants to enhance the quantum A discord, it is necessary to minimize the thermally excited photons in the second oscillator.Thus, if one puts most of the thermal photons on the first oscillator (unmeasured subsystem) and decreases the thermal photons in the second oscillator, this enhances the Gaussian quantum discord.In the latter design, C2 = 1.0 pF, by which the first and second oscillator's mean thermally excited photons for a typical frequency were, respectively, nth1 = 11.74,nth2 = 13.01; also the mean photon number (phase cross-correlation) for two coupled oscillators was around n12 ≡ < δa1δa2> = 2.63e-4.Since C2 is changed to 0.5 pF, the mean number of the thermally excited photons in the same frequency becomes nth1 = 11.73,nth2 = 9.38, and n12 = 9.95e-4.The estimated results show that the quantum A discord should be increased.The simulation result is shown in Fig.4aas a 3D graph.It is shown that the quantum discord amplitude is increased.This is because the second oscillator's thermally excited photons are decreased.Using illustrated 3D figure, one can get more information about the Gaussian quantum discord since the change of the quantifier as the 1D curves in the

Fig. 5
Fig. 5 InP HEMT DC characterization: Ids (mA) vs. gN2 (A/V 2 ), Nonlinear capacitance CN (pF) vs. gN2 (A/V 2 ) expressed in A1 and the related dynamic equation of motion, DC points (constant points) of the circuit can be calculated.Therefore, the steady state equations become: