The time behavior of DSN (displaced squeezed number state) for a two-dimensional electronic circuit composed of nanoscale elements is investigated using unitary transformation approach. The original Hamiltonian of the system is somewhat complicated. However, through unitary transformation, the Hamiltonian became very simple enough that we can easily treat it. By executing inverse transformation for the wave function obtained in the transformed system, we derived the exact wave function associated to the DSN in the original system. The time evolution of the DSN is described in detail, and its corresponding probability density is illustrated. We confirmed that the probability density oscillates with time like that of a classical state. There are two factors that drive the probability density to oscillate: One is the initial amplitude of complementary functions, and the other is the external power source. The oscillation associated with the initial amplitude gradually disappears with time due to the dissipation raised by resistances of the system. These analyses exactly coincide with those obtained from classical state. The characteristics of quantum fluctuations and uncertainty relations for charges and currents are also addressed.
Keywords:Displaced squeezed number states; Electronic circuits; Nanoscale physics; Unitary transformation; Fluctuations of charges and currents
The technical range of nanoscale is 1 to 999 nm, but people often refer to nanosize when an element is smaller than about 100 nm, where quantum effects are dominant instead of classical ones. Nanophysics and nanoelectronics have been rapidly developed thanks to the advancement of relevant technologies such as crystal growth and lithography, which facilitate sophisticated experiments for nanosystems [1,2]. A recent conspicuous trend in the community of electronic device is that the integrated circuits and components are miniaturized towards atomic-scale dimensions . We can confirm from many experiments and theories associated with nanoscale elements that the quantum effects become prominent when the transport dimension reaches a critical value which is the Fermi wavelength, while at the same situation, the classical theory for the motion of charges and currents is invalid. Not only quantum dot and quantum wire but also the quantum characteristics of electronic circuits involving nanoscale elements are important as a supporting theory for nanometer electronic technology and quantum information technology. For this reason, quantum effects in electronic circuits with nanoscale elements have been widely studied in recent years.
The simple quantum model of a lossless inductor-capacitor (LC) circuit have been suggested firstly by Louisell . Zhang et al. investigated the quantum properties of two-dimensional electronic circuits which have no power source . The quantum behavior of charges and currents for an LC circuit  and a resistor-inductor-capacitor (RLC) linear circuit  driven by a power source have been studied by several researchers. If a circuit contains resistance, the electronic energy of the system dissipates with time. In this case, the system is described by a time-dependent Hamiltonian. Another example of the systems described by time-dependent Hamiltonian is electronic circuits driven by time-varying power sources. The quantum problem of time-dependent Hamiltonian systems attracted great concern in the community of theoretical physics and chemistry for several decades [4,6,7].
The study of electronic characteristics of charge carriers in nanoelectronic circuits is basically pertained to a physical problem. There are plentiful reports associated with the physical properties of miniaturized two-loop (or two-dimensional) circuits [8-12] and more high multi-loop circuits [13-16] including their diverse variants. Various applications which use two-loop circuits include a switch-level resistor-capacitor (RC) model of an n-transistor (see Figure 3 of ), a design of a prototype of current-mode leapfrog ladder filters (Sect. 3 of ), and a port-Hamiltonian system , whereas higher loop circuits can be used as a transmission line model for multiwall carbon nanotube  and a filter circuit for electronic signals (Sect. 5 of ).
In this paper, we derive quantum solutions of a two-dimensional circuit coupled via RL and investigate its displaced squeezed number state (DSN) . We suppose that the system is composed of nanoscale elements and driven by a time-varying power source. The unitary transformation method which is very useful when treating time-dependent Hamiltonian systems in cases like this will be used. We can obtain the wave functions of DSN by first applying the squeezing operator in those of the number state and then applying the unitary displacement operator. Under displaced quantum states of circuit electrodynamics, conducting charges (or currents) exhibit collective classical-like oscillation. The fluctuations and uncertainty relations for charges and currents will be evaluated in the DSN without approximation.
Displaced squeezed number states, which are the main topic in this work, belong to nonclassical states that have been objects of many investigations. The statistical properties of these states exhibit several pure quantum effects which have no classical analogues, including the interference in the phase space , the revival/collapse phenomenon , and sub-Poissonian statistics . The position representation of these states with overall phases is derived by Moller et al. for the simple harmonic oscillator by employing geometric operations in phase space . The effects of quantum interference between two distinct DSNs prepared to be out of phase with respect to each other are investigated by Faisal et. al., discussing various nonclassical properties in connection with quantum number distribution, purity, quadrature squeezing, W-function, etc. .
Methods and results
Simplification via unitary transformation
Let us consider two loops of RLC circuit, whose elements are nanosized, that are coupled with each other via inductance and resistance as shown in Figure 1. Using Kirchhoff’s law, we obtain the classical equations of motion for charges of the system :
Figure 1. Electronic circuit. This is the diagram of a two-dimensional electronic circuit composed of nanoscale elements.
where qj (j=1,2; hereafter, this convention will be used for all j) are charges stored in the capacitances Cj, respectively, and is an arbitrary time-varying voltage source connected in loop 1. If we consider not only the existence of but also the mixed appearance of q1 and q2 in these two equations, it may be not an easy task to treat the system directly. If the scale of resistances are sufficiently large, the system is described by an overdamped harmonic oscillator, whereas the system becomes an underdamped harmonic oscillator in the case of small resistances. In this paper, we consider only the underdamped case.
For convenience, we suppose that R0/L0=R1/L1=R2/L2≡β. Then, the classical Hamiltonian of the system can be written as
where pj are canonical currents of the system, and kj=(1/Lj)(1/L0+1/L1+1/L2)−1/2. From Hamilton’s equations, we can easily see that pj are given by
If we replace classical variables qj and pj in Equation 3 with their corresponding operators, and , the classical Hamiltonian becomes quantum Hamiltonian:
where . Now, we are going to transform into a simple form using the unitary transformation method, developed in  for a two-loop LC circuit, in order to simplify the problem. Let us first introduce a unitary operator
Using Equation 7, we can transform the Hamiltonian such that
A straightforward algebra after inserting Equation 6 into the above equation gives
One can see from Equation 13 that the coupled term involving in the original Hamiltonian is decoupled through this transformation. However, the Hamiltonian still contains linear terms that are expressed in terms of , which are hard to handle when developing a quantum theory of the system. To remove these terms, we introduce another unitary operator of the form
where qjp(t) and pjp(t) are classical particular solutions of the firstly transformed system described by in the charge and the current spaces, respectively. From basic Hamiltonian dynamics with the use of Equation 12, we see that qjp(t) and pjp(t) satisfy the time-dependent classical equations that are given by
Then, the second transformation yields
The finally transformed Hamiltonian, Equation 22, is very simple and no longer involves linear terms that contain . If we neglect , this is exactly the same as that of the two-dimensional simple harmonic oscillator of frequencies ωj. We will use this formula in order to develop DSN, which is a typical nonclassical quantum state.
If we regard that the transformed Hamiltonian is very simple, the quantum dynamics in the transformed system may be easily developed. Let us write the Schrödinger equations for elements of the transformed Hamiltonian as
where represent number state wave functions for each component of the decoupled systems described by .
By means of the usual annihilation operator,
and the creation operator defined as the Hermitian adjoint of , one can identify the initial wave functions of the transformed system in number state such that
This formula of wave functions will be used in the next section in order to derive the DSN of the system.
Displaced squeezed number state
The DSNs are defined by first squeezing the number states and then displacing them. Like squeezed states, DSNs exhibit nonclassical properties of the quantum field in which the fluctuation of a certain observable can be less than that in the vacuum state. This state is a generalized quantum state for dynamical systems and, in fact, equivalent to excited two-photon coherent states in quantum optics. If we consider that DSNs generalize and combine the features of well-known important states such as displaced number states (DNs) , squeezed number states , and two-photon coherent states (non-excited) , the study of DSNs may be very interesting. Different aspects of these states, including quantal statistics, entropy, entanglement, and position space representation with the correct overall phase, have been investigated in [17,23,25].
To obtain the DSN in the original system, we first derive the DSN in the transformed system according to its exact definition. Then, we will transform it inversely into that of the original system. The squeeze operator in the transformed system is given by
Using the Baker-Campbell-Hausdorff relation that is given by 
where , the squeeze operator can be rewritten as
Let us express the DSN in the transformed system in the form
where represent two decoupled states which are drivable from
Here, are displacement operators in the transformed system, which are given by
where αj is an eigenvalue of at initial time. By considering Equation 26, we can confirm that
where qjc(t) and pjc(t) are classical solutions of the equation of motion in charge and current spaces, respectively, for the finally transformed system. If we regard that the complementary functions  of the equation of motion in the firstly transformed system are the same as the classical solutions of the finally transformed system, qjc(t) and pjc(t) can also be complementary functions of the firstly transformed system. The other operators are time-displacement operators:
At first, the action of squeezing operator in wave functions of the initial number state gives
The evaluation of the other actions of the operators in Equation 34 may be easily performed using Equation 31 and the relation 
together with the eighth formula of 7.374 in  (see Appendix Appendix 1), yielding
Here, the time evolution of complementary functions are
The transformed system reduces to a two-dimensional undriven simple harmonic oscillator in the limit . Our result in Equation 44 is exact, and in this limit, we can easily confirm that some errors in Equation 45 in  are corrected (see Appendix Appendix 2).
The wave function associated to the DSN in the transformed system will be transformed inversely to that of the original system in order to facilitate full study in the original system. This is our basic strategy. Thus, we evaluate the DSN in the original system from
Using the unitary operators given in Equations 7 and 16, we derive
This is the full expression of the time evolution of wave functions for the DSN. If we let r→0, the squeezing effects disappear, and consequently, the system becomes DN. Of course the above equation reduces, in this limit, to that of the DN.
To see the time behavior of this state, we take a sinusoidal signal as a power source, which is represented as
Then, the solution of Equations 19 and 20 is given by
The probability densities are plotted in Figures 2 and 3 as a function of q1 and t under this circumstance. As time goes by, the overall probability densities gradually converge to the origin where q1=0 due to the dissipation of energy caused by the existence of resistances in the circuit. If there are no resistances in the circuit, the probability densities no longer converge with time. An electronic system in general loses energy by the resistances, and the lost energy changes to thermal energy. Actually, Figure 2 belongs to DN due to the condition r1=r2=0 supposed in it. The wave function used in Figure 2a is not displaced and is consequently the same as that of the number state. Figure 2b is distorted by the effect of displacement. From Figure 2c,d, you can see that the exertion of a sinusoidal power source gives additional distortion. The frequency of is relatively large for Figure 2c whereas it is small for Figure 2d.
Figure 2. Probability density (A). This represents the probability density as a function of q1 and t. Here, we did not take into account the squeezing effect (i.e., we let r1=r2=0). Various values we have taken are q2=0, n1=n2=2, , R0=R1=R2=0.1, L0=L1=L2=1, C1=1, C2=1.2, p1c(0) = p2c(0) = 0, and δ = 0. The values of are (0,0,0,0) (a), (0.5,0.5,0,0) (b), (0.5,0.5,10,4) (c), and (0.5,0.5,0.5,0.53) (d). All values are taken dimensionlessly for convenience: this convention will be used in all subsequent figures.
Figure 3. Probability density (B). The probability density with squeezing parameters r1 = r2 = 0.7 and ϕ1 = ϕ2 = 1.5 is shown here as a function of q1 and t. Various values we have taken are q2 = 0, n1 = n2 = 2, , R0 = R1 = R2 = 0.1, L0 = L1 = L2 = 1, C1 = 1, C2 = 1.2, p1c(0) = p2c(0) = 0, and δ = 0. The values of are (0,0,0,0) (a), (0.5,0.5,10,4) (b), and (0.5,0.5,0.5,0.53)(c).
You can see the effects of squeezing from Figure 3. The probability densities in the DSN are more significantly distorted than those of the DN. We can see from Figure 3b,c that the time behavior of probability densities is highly affected by external power source. If there is no power source in the circuit, the displacement of charge, specified with an initial condition, may gradually disappear according to its dissipation induced by resistances in the circuit. This is the same as that interpreted from the DN and exactly coincides with classical analysis of the system.
While various means and technologies to generate squeezed and/or displaced light are developed in the context of quantum optics after the seminal work of Slusher et al.  for observing squeezed light in the mid 1980s, (displaced) squeezed number state with sufficient degree of squeezing for charges and currents in a circuit quantum electrodynamics is first realized not long ago by Marthaler et al.  as far as we know. The circuit they designed not only undergoes sufficiently low dissipation but its potential energy also contains a positive quartic term that leads to achieving strong squeezing. Another method to squeeze quantum states of mechanical oscillation of charge carriers in a circuit is to use the technique of back-action evasion [33,34] that is originally devised in order to measure one of two arbitrary conjugate quadratures with high precision beyond the standard quantum limit.
Though it is out of the scope of this work, the superpositions of any two DSNs may also be interesting topics to study, thanks to their nonclassical features that have no classical analogues. The quantum properties such as quadrature squeezing, quantum number distribution, purity, and the Mandel Q parameter for the superposition of two DSNs out of phase with respect to each other are studied in the literatures (see, for example, ).
Now let us see the quantum fluctuations and uncertainty relations for charges and currents in the DSN for the original system. It is well known that quantum energy and any physical observables are temporarily changed due to their quantum fluctuations. The theoretical study for the origin and background physics of quantum fluctuations have been performed in  by introducing stochastic and microcanonical quantizations.
If we consider the method of consecutive unitary transformation, the expectation value for an arbitrary operator in the original system can be evaluated from
Using this relation, the expectation value of charges and currents is derived to be
The expectation value of square of and can also be obtained form the same method, and we listed them in Appendix Appendix 3. In fact, Equations 58 and 59 are the same as the classically predicted amount of charges qcl,1 and qcl,2 in C1 and C2 in the original system, respectively. If we consider that αj are given by Equation 36, qcl,1 and qcl,2 can be rewritten, after a little evaluation, in the form
We illustrated qcl,1 and qcl,2 in Figure 4 as a function of time. To understand the time behavior of these quantities, it may be worth to recall that complementary functions, qjc(t), and particular solutions, qjp(t), are not associated to the original system but to the firstly transformed system. We can also easily confirm from similar evaluation that the time behavior of canonical conjugate currents pcl,j are represented in terms of qjc(t), pjc(t), and pjp(t) (see Appendix Appendix 4).
Figure 4. Classically predicted amount of charges in capacitors. This illustration represents the time behavior of qcl,1 (thick solid line) and qcl,2 (dashed line) where R0 = R1 = R2 = 0.1, L0 = L1 = L2 = 1, C1 = 1, C2 = 1.2, q1c(0) = q2c(0) = 0.5, p1c(0) = p2c(0) = 0, and δ = 0. The values of are (0,0) (a), (10,4) (b), and (0.5,0.53) (c).
The definition of quantum fluctuations for any quantum operator in the DSN is given by
Using this, we obtain the fluctuations of charges and currents as
As we have seen before, the expectation values associated to charges and currents are represented in terms of complementary functions, qjc(t) and pjc(t), and particular solutions qjp(t) and pjp(t). The amplitude of complementary functions is determined from the strength of displacements, whereas the particular solutions are determined by the power source (see Equations 19 and 20). However, all of the fluctuations do not involve such solutions. This means that the displacement and the electric power source do not affect to the fluctuations of charges and currents.
The uncertainty products between charges and their conjugate currents can be easily identified by means of Equations 67 to 70. For the case of the DN that are given from the limit r1=r2→0, we have F1=F2=0 and . Then, the uncertainty products become
These are the same as the uncertainty products in the number states and are always larger than , preserving the uncertainty principle. Thus, we can conclude that the uncertainty products in the DN are the same as those of the ordinary number states. Evidently, the uncertainty principle is inherent in quantum mechanical context described by canonical variables. The results, Equations 71 and 72 with n1=n2=0, are exactly the same as Equations 29 and 30 of , respectively. Moreover, for R1=R2=R3→0 (i.e., β→0), the above two equations reduce to Equations 52 and 53 in , which are evaluated in ordinary number state. Hence, this work includes all the results of both  (no power source) and  (no resistances) as special cases. The fluctuations and uncertainty product in the DN and in the DSN are plotted in Figure 5. We can adjust the uncertainty (or fluctuation) of a quadrature to be small at the expense of broadening that of another quadrature, or vice versa. The uncertainty in the case of this figure is larger than , while is smaller than due to the squeezing effect. Therefore, it is relatively difficult for us to know the precise value of charge q1, while we can find out the conjugate current p1 more precisely. However, the relevant uncertainty product in the DSN is nearly unaltered from that in the DN.
Figure 5. Fluctuations. This inset shows fluctuations (dashed line) and (thick solid line) (a), and (dashed line) and (thick solid line) (b), and uncertainty product (dashed line) and (thick solid line) (c) as a function of t where n1=n2=0, , R0 = R1 = R2 = 0.1, L0 = L1 = L2 = 1, C1 = 1, and C2 = 1.2. The values of squeezing parameters for the DSN are r1 = 0.1, r2 = 0.3, ϕ1 = 1.2, and ϕ2 = 0.6.
In summary, the time evolution of the DSN for the two-dimensional electronic circuit composed of nanoscale elements and driven by a power source is investigated using unitary transformation method. Two steps of the unitary transformation are executed: We removed the cross term involving in the original Hamiltonian from the first step, and the linear terms represented in terms of in the firstly transformed Hamiltonian are eliminated by second unitary transformation.
We can see from Equation 6 that the original Hamiltonian is time-dependent. When treating a time-dependent Hamiltonian system dynamically, one usually employs classical solutions of the equation of motion for a given system (or for a system similar to a given system) [6,7]. We also introduced such classical solutions in Equations 19 to 20 and in Equations 47 to 48. Among them, particular solutions qjp and pjp are important in developing quantum theory of the system involving external power source since they are crucial factors that lead the transformed Hamiltonian to be simple so that we can easily treat it.
Since the transformed system is just the same as the one that consists of two independent simple harmonic oscillators, provided that we can neglect the trivial terms in the transformed Hamiltonian, we easily identified the complete quantum solutions in the DSN in the transformed system. We also obtained the wave functions of the DSN in the original system via the technique of inverse transformation, as shown in Equation 50. If we regard the fact that the probability does not reflect the phase of a wave function, the overall phase of these states is relatively unimportant for many cases. However, in some applications such as the computation of expectation values using generating or characteristic functions given in , the exact knowledge of overall phase is crucial. For r1=r2=0, the wave function in the DSN exactly reduces to that of the DN.
We analyzed the probability densities in the DN and in the DSN from Figures 2 and 3, respectively, with the choice of sinusoidal signal source. The probability densities in the DN given in Figure 2b,c,d oscillate with time. Moreover, their time behaviors are more or less distorted. The probability density, however, does not oscillate when there are no displacement and no signal of power source (see Figure 2a). The probability densities in the DSN are distorted much more significantly than those of the DN.
The time behavior of probability densities of quantum states, both the DN and the DSN, is highly affected by external driving power source. When there is no external power source( =0), the displacement of charges, specified with a certain initial condition, gradually disappears as time goes by like a classical state.
The fluctuations and uncertainty products of charges and currents are derived in the DSN, and it is shown that their value is independent of the size of the particular solutions qjp(t) and pjp(t). From this, together with the fact that qjp(t) and pjp(t) are determined by the characteristics of , it is clear that the electric power source does not affect on the fluctuation of canonical variables. If we ignore the time dependence of Fj(t) and , decrease exponentially with time, whereas increase exponentially.
From Equations 64 and 65, we can see that the time behavior of qj is determined by two factors: One is displacement and the other is the signal of power source. For better understanding of this, recall that the amplitude of complementary functions gives displacement of the system, and the particular solutions are closely related to external driving force (i.e., in this case, the power source).
In this paper, we did not consider thermal effects for the system. The thermal effects, as well as dissipation, may be worth to be considered in the studies of quantum fluctuations of electronic circuits with nanosize elements because the practical circuits are always working in thermal states with the presence of damping. It may therefore be a good theme to investigate DSNs with thermalization as a next task, and we plan to investigate it in the near future.
The eighth formula of 7.374 in 
Correction of Equation 45 of 
The second line of Equation 45 of  needs to be corrected as
Besides, among various functions that appeared in Equation 45 of , (Equation 23) and A (Equation 46) should be altered as
For the convenience of comparison, we provide a list of correspondences between our notations and the notations used in :
Expectation value of and
According to the rule, Equation 57, for evaluating expectation values, we also have the expectation value of square of charges and currents as
Through the same vein as that of the calculation of qcl,1 and qcl,2 given in Equations 64 and 65, we can evaluate classical currents pcl,1 and pcl,2 from their quantum expectation value given in Equations 60 and 61. Thus, we have
DN: displaced number state.
The authors declare that they have no competing interests.
Physica A 1993, 197:364-370. Publisher Full Text
Phys Scr 2006, 73:587-595. Publisher Full Text
Integration-VLSI J 1996, 21:1-94. Publisher Full Text
Int J Electron Commun 2010, 64:503-511. Publisher Full Text
Automatica 2003, 39:969-979. Publisher Full Text
Mechatronics 2010, 20:201-212. Publisher Full Text
Microelectron Reliab 2012, 52:3026-3034. Publisher Full Text
Microelectron J 2011, 42:63-73. Publisher Full Text
Microelectron J 2010, 41:834-839. Publisher Full Text
Physica A 2001, 291:449-466. Publisher Full Text
Phys Rev A 1981, 23:236-247. Publisher Full Text
J Mod Opt 1993, 40:1331-1350. Publisher Full Text
J Opt B: Quant Semiclass Opt 2003, 5:60-72. Publisher Full Text
Phys Rev D 1985, 32:400-404. Publisher Full Text
J Mod Opt 1987, 34:709-759. Publisher Full Text
Chaos, Solitons Fractals 2001, 12:2455-2470. Publisher Full Text
J Phys A: Math Gen 1998, 31:4329-4336. Publisher Full Text
Quantum Semiclass Opt 1996, 8:1061-1066. Publisher Full Text
Phys Lett A 1997, 229:135-143.
[Equation 45]Publisher Full Text
New J Phys 2008, 10:095010. Publisher Full Text
Phys Lett A 2004, 324:145-151. Publisher Full Text