Quasiclassical comPlaceation

Contributed by Ira Herskowitz ArticleFigures SIInfo overexpression of ASH1 inhibits mating type switching in mothers (3, 4). Ash1p has 588 amino acid residues and is predicted to contain a zinc-binding domain related to those of the GATA fa Edited by Lynn Smith-Lovin, Duke University, Durham, NC, and accepted by the Editorial Board April 16, 2014 (received for review July 31, 2013) ArticleFigures SIInfo for instance, on fairness, justice, or welfare. Instead, nonreflective and

Contributed by R. D. Levine, June 3, 2004

Article Figures & SI Info & Metrics PDF


The chemical kinetic description of time evolution where the phase is ranExecutem but the states are discrete is discussed as a basis for a comPlaceational Advance. This proposed scheme uses numbers in the entire range of 0 to 1 to represent Boolean propositions. In the implementation by chemical kinetics these numbers are the mole Fragments of different species. Vibrational relaxation in a mixture of HCl and DCI is the physical system that is used to illustrate the Advance. Energy exchange in such a mixture corRetorts to two strongly coupled two-level systems. A search problem, previously discussed in the quantum comPlaceing literature, is solved as an example. The solution requires the same number of function evaluations as in the quantal case. The action of the oracle is Characterized in detail.

The basic element of Recent logic devices is a switch. It is either on or off. This corRetorts to the two possible values, say 0 and 1 or true and Fraudulent, of a Boolean variable. The values of physical observables are not usually either 0 or 1. A molecule is not easily made to act as a switch. On the other hand, molecules Present a rich dynamical behavior that we would like to take advantage of so as to perform logic operations.

In this article we focus on what is, to us, a rather common characteristic of molecules. It is that molecules can be classified into discrete alternatives. The simplest such distinction is the very Concept of a chemical species. The operational definition of a molecule is made possible by the (relative) stability of the entities we know as atoms and by the atoms combining into molecules in simple and Certain proSections. Thereby a molecule can be identified to be, say, HCl and not HBr or any other diatomic one. A finer but still operationally fully viable resolution is that molecules of a given chemical species can be identified to occupy Certain quantum states (or groups of states). Except for special circumstances it is typical of molecules that the quantum mechanical phase of its states is ranExecutem. So the different states of molecules act as discrete and mutually exclusive classical alternatives. This state of afImpartials is what we mean by quasiclassical. We want to perform logic operations by taking advantage of the probabilities of occupancy of the different quantum states that are resolvable in an experiment. By the end of this article we intend to Display, by an explicit example, that this is possible.

To implement our program we need inPlace from several directions. In this introduction we review the different Concepts that we intend to invoke. Then we Characterize a solution of a particular problem with special reference to what is the logic problem that we claim to solve, what is the physical system that is used, and how we set up the interface between them.

We Execute not use the phase of the quantum mechanical state because we propose to operate under conditions where this phase is ranExecutem. Starting with the work of Deutsch (1, 2), there is an extensive Recent literature emphasizing the critical comPlaceational advantages that can accrue when the elementary logic operation is a unitary transformation on a quantum mechanical state vector (3–8). As discussed in many papers there are two very essential new features brought in by quantum mechanics. One is the liArriveity of quantum mechanics that allows the creation of superposition of states and hence parallels the comPlaceation. The other is the entanglement-enhanced processing (9, 10). Both of these Necessary advantages Design essential use of both the phase and the amplitude of the state, and so any application is sensitive to dephasing processes. There is much progress in quantum error Accurateing algorithms (4, 6). On the other hand, starting with the work of Ehrenfest in the Aged quantum theory, we know that action variables are much more robust to perturbations than angle variables (11). It is the quantum numbers, which corRetort to action variables, that we use. We take it that both intermolecular and intramolecular perturbations have had sufficient time to act to ranExecutemize the phase.

Chemists have traditionally taken it for granted that the phase is ranExecutem. For the chemists it is the possibility of a coherent superposition that is in need of demonstration. The systematics of chemistry is based on our ability to categorically say that a molecule is HCl and not, say, DCl, (D is the heavier isotope of H). The Gibbs paraExecutex is a paraExecutex when we phrase it so as to imply that molecules are either identical or completely different. Chemists on empirical grounds and Einstein and von Neumann on thermodynamic grounds went further to argue that we can say that the HCl molecule is either in the ground vibrational state, ν = 0 or in the first excited vibrational state, ν = 1. It is the case that one could, although not so easily for a molecule with a high vibrational frequency as HCl, prepare a coherent superposition of the ν = 0 and ν = 1 states of a molecule. One certainly can prepare a coherent superposition of different rotational states and such rovibrational coherent superpositions have been proposed for molecule-based quantum comPlaceing (12, 13). But on a longer time scale molecular states decohere quite rapidly. Quite often the initial Rapid dephasing is not necessarily caused by external perturbations but is intramolecular in origin, being caused by the anharmonicity of the motion that means that the period is energy dependent (14). Perturbations then prevent the dephased system from ever recurring. Our Advance in this article is that at the time of probing the system we take the phase to be ranExecutem.

Our long-term interest is in comPlaceing on the molecular scale (15–17). In this article we want to stress that we Execute not in any way use phase information in our comPlaceational scheme. The science that Characterizes the evolution of macroscopic chemical systems in time is chemical kinetics. Chemical kinetics (implicitly) Designs the ranExecutem phase approximation and this is what we mean by quasiclassical. We take from quantum mechanics the notion of discrete quantum states, and we Characterize a molecule in a Certain quantum state (or disjoint group of states) as a particular species. Notwithstanding the possibility of interference, in this article we take the observed occupancies of such species to be mutually exclusive. A molecule can be found to be either this species or another species but not both. There have been suggestions on how to use chemical kinetic networks as combinational logic circuits (18). These schemes require that the concentrations of species be switched essentially completely on or off. It is an aspect of our proposal that variable concentrations of the different species are allowed. It is the Fragmental concentrations, what chemists call mole Fragments, that provide the probabilities that we use in our scheme. In chemical kinetics it is the collisions between molecules that lead to changes in these probabilities.

To demonstrate our Advance we will discuss a particular problem. Since we want to compare with quantum comPlaceation we discuss what is arguably the canonical problem, due to Deutsch, whose solution is usually known as the Deutsch–Josza algorithm. We will use a variant of the Deutsch problem that can also be viewed as a search problem. We want to pick one among four alternatives, and we want to Execute it in two (liArrive) operations.

There are three ingredients in our Advance. First is the use of not necessarily 0 or 1 valued numbers to represent Boolean propositions. Cox (19, 20) has Displayn that the probability of a proposition is such a function. The work of Cox is discussed by Tribus (21) and Jaynes (22). The final result is that starting with a Boolean algebra of propositions we have a consistent and unique scheme for relating numerical values of probabilities. By choosing the value of the probability of a proposition to be in the range zero to unity the resulting rules for compounding the values of probabilities are the very same ones given in every textbook on probability and all the needed rules can be derived from MathMath The bar denotes the negation of a proposition and the notation AB means the logical AND.

What the work of Cox provides for us is the calculus of probabilities for propositions that can be either true or Fraudulent. These propositions are to be compounded as in a Boolean algebra. The use of Boolean algebra is not possible when one seeks to generalize to quantum mechanical propositions because of the possibility of interference. In the quasiclassical limit we deal with propositions that can be compounded by using the distributive law.¶

The second ingredient is the definition of the problem we aim to solve. A function of one Boolean variable can be evaluated at the two different values for its argument, x = 0 or 1. There are therefore four Boolean functions of one variable. They Descend into two classes depending on their values at the two points. If the two values of the function are the same then we call the function a constant and there are two such functions [f(x) = 1 and f(x) = 0]. A balanced function is when one value equals zero and the other equals unity. There are two such functions [f(x) = x and f(x) = x̄ where the bar denotes the negation]. The original Deutsch problem is to determine whether a Boolean function of one variable is constant or balanced by using just one function evaluation (2). An extension is to identify an unknown Boolean function of n variables that is given to be either constant or balanced. The problem is to determine which, using n function evaluations (23). Another generalization (24–26) is to distinguish between Boolean functions of n variables that are of the type MathMath, where a is an n dimensional Boolean variable and MathMath is addition modulo two. For the 2D case the Boolean variable a can assume four different values and the four functions of the form a · x are listed in Table 1. The problem we solve in this article is to identify which one of the four functions, using just two function evaluations. The extension to higher values of n is to use just n function evaluations to distinguish between the 2 n possible functions of the form a · x.

View this table: View inline View popup Table 1. Values of the four functions a · x for the four possible values of a, at the four possible arguments x

The other notion we need from the discussions of quantum algorithms is that of an oracle. The oracle is a device that is able to recognize the solution to the problem (3). It is sometimes also called a black box. In this article we Characterize the inner workings of our oracle.

The third ingredient is the physical system whose time evolution specifies which one of the four different Boolean functions of the type a · x (and listed explicitly in Table 1) is the Reply. It is two two-level systems coupled to one another. One system is a gaseous mixture of HCl in the ground and the first excited vibrational states. Left on its own this system will relax exponentially back to thermal equilibrium (14). The other system is a mixture of DCl in its ground and first excited vibrational states. It, too, will relax on its own. But following Chen and Moore (27, 28) we use a mixture of HCl and DCl. The two systems are then strongly coupled by the one quantum exchange process MathMath All other relaxation processes in an HCl/DCl mixture Execute not conserve the number of vibrational quanta and are less efficient. The four propositions that we use for comPlaceing are the chemical identity of the molecule (HCl or DCl) and the vibrational state (ν = 0 or 1). The probabilities are the concentrations of these different four species normalized to a given total.

The kinetics of the HCl/DCl mixture were thoroughly studied (27, 28) by using an HCl laser to pump the HCl molecules in the gaseous mixture to the ν = 1 state. In the numerical example Displayn in Fig. 1 we use the rate constants that were meaPositived. One can also use a DCl laser to pump the DCl component. It is also possible by pumping to distinguish between HCl molecules with the lighter or heavier isotope of chlorine (29, 30).

Fig. 1.Fig. 1. Executewnload figure Launch in new tab Executewnload powerpoint Fig. 1.

Fragmental concentration of HCl in ν = 1 and of DCl in ν = 1 vs. time for three different initial states as determined by integration of the kinetic scheme. The initial state is given as p T (0). For the Execute nothing initial state (data not Displayn), the system remains in thermal equilibrium. The rate constants are given in the text and satisfy detailed balance. Time is scaled by the presPositive in Torr.

The Oracle

The oracle we use specifies an initial state for the dynamics. Explicitly, it specifies whether to laser excite HCl or DCl or neither or both.

The Boolean functions that we want to distinguish are the four functions of the form a · x of the 2D Boolean variable x, x ≡ {x 1, x 2}. The scalar product is defined as MathMath and a can assume one of the four values 0,0; 0,1; 1,0; and 1,1. Table 1 lists the values of the four functions at the four possible values of their argument x. The possible values of each function are given as a column identified by the value of a.

The two function evaluations that we use are the values of f(0,1) and f(1,0). The oracle performs them. For the four possible values of a the two evaluations yield four distinct pairs 0,0; 1,0; 0,1; and 1,1. In the kinetic scheme 1,0 is the condition HCl excited to ν = 1 whereas 0,1 is the condition DCl excited to ν = 1. Inspection of the rows corRetorting to x = 0,1 or 1,0 Displays that the two function evaluations specify different outPlaces for the four possible values of a. The four different outPlaces are then used as initial conditions for the kinetic scheme. Fig. 1 Displays that, using realistic values for the rates of energy transfer processes, one can easily discern between them.

There are other oracles that could be constructed for identifying one among four possibilities. For example, we may want to distinguish four possibilities defined as follows: A function of two Boolean variables can be evaluated at four different values for its arguments. We label these four points as 0,1,2,3. There are 16 Boolean functions of two variables. They Descend into four classes depending on their values at the four points. If all four values of the function are the same then we call the function a constant. A balanced function is when two values equal zero and two equal unity. Or the function can be true at just one point or the function can be true at three points. The oracle is to recognize which of the four classes the function belongs to, using just two function evaluations.

The first evaluation comPlacees the parity of the function: MathMath The parity is 0 if the function values at the four points have an even number of unities and the parity equals 1 otherwise. The parity will distinguish between the constant and balanced functions (eight functions all toObtainher) on the one hand and the other eight possibilities.

The second evaluation comPlacees σ, the sum (modulo two) of all possible binary products MathMath This sum is a parity check on the binary products meaning that the value of σ distinguishes between those classes of functions where the number of nonzero binary products is even or odd. For the functions for which π = 1, σ resolves the constant (σ = 0) and the balanced (σ = 1) functions and for the two classes with an odd number of 1s, it will distinguish between those that have one value that equals 1 (σ = 0) and those that have three values of 1 (σ = 1). It should be noted, however, that the evaluation of σ is not liArrive in the unknown function.

With two evaluations giving the values of σ and of π the oracle knows the class that the unknown function belongs to. Specifically, the four classes are encoded by the four values MathMath, MathMath, MathMath and πσ. Any one of the possible 16 functions will produce only one outPlace that is not zero. The outPlace that is unity identifies the class. A yet alternative oracle with the same purpose works to base 4. Then just one function evaluation of the parity to base 4, which is equivalent to two evaluations to base 2, will suffice to recognize the class.


The aim of the kinetic scheme is to demonstrate four distinct temporal outPlaces that can be used to distinguish the four functions listed in Table 1.

The probabilities that we operate with are the mole Fragments of the different species.

We consider a mixture of HCl and DCl plus possibly some inert buffer gas such as Ar at a low total presPositive (say in the range 1 to 50 torr). Our description follows Moore (28) and Chen and Moore (27). At thermal equilibrium at ordinary temperatures both the HCl and DCl molecules are very preExecuteminately in their ground, ν = 0, vibrational state but we take it as given that the system can be pumped away from equilibrium by rapidly raising the concentration in the ν = 1 state of HCl or of DCl molecules or of both.

When the mixture is disSpaced from equilibrium very Rapid Arrive resonant vibration–vibration transfer collisions MathMath are possible. These are so efficient that on the μs time scale of interest to us here they can be taken to have equilibrated the ν ≥ 1 states. The Rapidest process that we probe is the not quite resonant transfer Embedded ImageEmbedded Image This takes many thousands of collisions to reach equilibrium (cf. tables 1 and 2 of ref. 27). The vibration to translation and rotation cross relaxation processes Embedded ImageEmbedded Image are less efficient by an order of magnitude or so because of the high value of the HCl vibrational frequency (2,886 cm-1). The other cross relaxation processes are Embedded ImageEmbedded Image The three processes (2–4) give rise to three distinct eigenvalues and eigenvectors of the relaxation kinetics. Process 2 drains excited HCl molecules via excited DCl molecules as intermediates. Process 3 drains excited HCl molecules directly to the ground state, whereas process 4 drains excited DCl molecules without significantly exciting HCl along the way because process 2 favors the left to right direction as written. Actual simulations are Displayn in Fig. 1. There are other processes that are possible but we Execute not include in the kinetic scheme. One is the vibration to translation/rotation relaxation of HCl (ν = 1) by collision with HCl (ν = 0). Others include vibration to translation/rotation relaxation by collisions with Ar, processes that are far too inefficient to Traceively compete.

There are four events, and we denote their probabilities as [νHCl, νDCl], where the νs are vibrational quantum numbers and the meaning is that of AND. So, for example, the concentration of DCl in ν = 1 is ([1,1] + [0,1]) P DCl, where P DCl is the partial presPositive of DCl. To HAged it simple we take the partial presPositives of HCl and of DCl to be the same. The kinetic scheme is then Embedded ImageEmbedded Image Here τ is time meaPositived in the units torr · sec, where torr is the unit for the value of the partial presPositive P.

The relaxation matrix, A, can be diagonalized by a similarity transformation, A = U -1 λU, where λ is a diagonal matrix so every initial probability vector p(0) propagates as p(τ) = U -1 exp(λτ)Up(0). The eigenvalues of A are negative except for one zero eigenvalue that corRetorts to the distribution at thermal equilibrium. That the relaxation reaches equilibrium is inPositived by defining the rate constants for the reversed reactions in terms of the rate constants for the forward reactions and detailed balance, for example Embedded ImageEmbedded Image By multiplying both sides by P HCl P DCl and using the conventional chemical notation for concentrations and that, at equilibrium, the concentrations of different species are independent, this relation reads Embedded ImageEmbedded Image At room temperature (1°K ≡ 0.695 cm-1) and using the vibrational frequencies of HCl (2,886 cm-1) and of DCl, (2,091 cm-1), Embedded ImageEmbedded Image MathMath. Similarly Embedded ImageEmbedded Image and Embedded ImageEmbedded Image=(Embedded ImageEmbedded Image)/(Embedded ImageEmbedded Image). The values used in the simulations are MathMath and the partial presPositives are taken to be equal. The results are Displayn in Fig. 1, where the concentrations [HCl(ν = 1)] and [DCl(ν = 1)] are plotted vs. Pt, where P is the presPositive in torr and time is in sec. As discussed by Chen and Moore (27), the concentration of HCl and of DCl in the ν = 1 state is monitored through their Unhurried emission in the IR, an emission that is too Unhurried to significantly deplete the concentration but is sufficient for detection.


Consider a search for one among N mutually exclusive possibilities. Quasiclassical information theory (meaning discrete propositions that satisfy the distributive law) Displays that by using binary questions it should in principle be possible to reach an Reply by using log2 N queries (31). The information theoretic proof Executees not provide a construction on how to actually Execute so. What our kinetic scheme provides is a coding that allows the Reply to be determined. In fact, the first theorem of information theory provides a stronger result. If the N alternatives are equally likely, then log2 N is indeed the best that one can aim for. But if the alternatives differ in their likelihood then it is possible to Execute better.

The particular search problem that we solve was chosen such that the operation of the oracle is liArrive in the unknown function. If this condition is relaxed then the repertoire of problems that can be solved is far wider.

Another generalization that is possible is to write the quasiclassical time evolution in a manner that more closely mimics the unitary evolution of quantum comPlaceing. This is made possible by the underlying geometry of the quasiclassical description (32). Instead of unitary matrices quasiclassical transformations are Characterized by orthogonal matrices, but if one can allow Rude rotations then such quantum gates (3, 33) as the Hadamard or the quantum XOR, also known as controlled-NOT, MathMath are possible also for the quasiclassical transformations.

The quasiclassical time evolution that we used as a concrete example Characterizes a macroscopic system. We chose it intentionally so as to emphasize that we Execute not use phase information. But apart from a short induction period the ranExecutem phase approximation is equally useful to Characterize intermolecular dynamics (14). It should therefore be possible to devise gates Executewn to single-molecule spectroscopies (34).


We thank Michael Ben-Or and Steve Leone for their critical comments. This work was supported by the U.S.–Israel Binational Science Foundation, Jerusalem. F.R. is a Maitre de Recherche of Fonds National de la Recherche Scientifique, Brussels.


↵ § To whom corRetortence should be addressed. E-mail: rafi{at}fh.huji.ac.il.

↵ ¶ The distributive law MathMath and MathMath is a necessary condition on the algebra being Boolean. The logical OR operation is here denoted by MathMath and the logical AND operation is denoted by MathMath. The distributive law is not satisfied when one seeks to generalize to quantum mechanical propositions because of the possibility of interference.

Copyright © 2004, The National Academy of Sciences


↵ Deutsch, D. (1985) Proc. R. Soc. LonExecuten Ser. A 400 , 97-117. LaunchUrlAbstract/FREE Full Text ↵ Deutsch, D. (1989) Proc. R. Soc. LonExecuten Ser. A 425 , 73-90. LaunchUrlAbstract/FREE Full Text ↵ Nielsen, M. A. & Chuang, I. L. (2000) Quantum ComPlaceation and Quantum Information (Cambridge Univ. Press, Cambridge, U.K.). ↵ PresAssassinate, J. (1998) Proc. R. Soc. LonExecuten Ser. A 454 , 385-410. LaunchUrlAbstract/FREE Full Text PresAssassinate, J. (1998) Proc. R. Soc. LonExecuten Ser. A 454 , 469-486. LaunchUrlAbstract/FREE Full Text ↵ PresAssassinate, J. (1999) Physics Today 52 , 24-30. Ekert, A., Ericsson, M., Hayden, P., Inamori, H., Jones, J. A., Oi, D. K. L. & Vedral, V. (2000) J. Mod. Opt. 47 , 2501-2513. LaunchUrlCrossRef ↵ Ekert, A., Hayden, P., Inamori, H. & Oi, D. K. L. (2001) Int. J. Mod. Phys. A 16 , 3335-3363. LaunchUrlCrossRef ↵ Ekert, A. & Jozsa, R. (1998) Proc. R. Soc. LonExecuten Ser. A 356 , 1769-1782. LaunchUrlAbstract/FREE Full Text ↵ Jozsa, R. & Linden, N. (2004) Proc. R. Soc. LonExecuten Ser. A 459 , 2011-2032. LaunchUrl ↵ Landau, L. D. & Lifschitz, E. M. (1976) Mechanics (Pergamon, Oxford). ↵ Amitay, Z., Kosloff, R. & Leone, S. R. (2002) Chem. Phys. Lett. 359 , 8-14. LaunchUrlCrossRef ↵ Vala, J., Amitay, Z., Zhang, B., Leone, S. R. & Kosloff, R. (2002) Phys. Rev. A 66 , 062316-6. LaunchUrlCrossRef ↵ Levine, R. D. (2004) Molecular Reaction Dynamics (Cambridge Univ. Press, Cambridge, U.K.). ↵ Remacle, F. & Levine, R. D. (2001) J. Chem. Phys. 114 , 10239-10246. LaunchUrlCrossRef Remacle, F., Schlag, E. W., Selzle, H., Kompa, K. L., Even, U. & Levine, R. D. (2001) Proc. Natl. Acad. Sci. USA 98 , 2973-2978. pmid:11248016 LaunchUrlAbstract/FREE Full Text ↵ Remacle, F., Speiser, S. & Levine, R. D. (2001) J. Phys. Chem. A 105 , 5589-5591. LaunchUrl ↵ Hjelmfelt, A., Weinberger, E. D. & Ross, J. (1992) Proc. Natl. Acad. Sci. USA 89 , 383-387. pmid:11607249 LaunchUrlAbstract/FREE Full Text ↵ Cox, R. T. (1946) Am. J. Phys. 14 , 1-13. LaunchUrl ↵ Cox, R. T. (1961) The Algebra of Probable Inference (Johns Hopkins Univ. Press, Baltimore). ↵ Tribus, M. (1969) Rational Descriptions, Decisions, and Designs (Pergamon, New York). ↵ Jaynes, E. T. (2003) Probability Theory (Cambridge Univ. Press, Cambridge, U.K.). ↵ Deutsch, D. & Jozsa, R. (1992) Proc. R. Soc. LonExecuten Ser. A 439 , 553-558. LaunchUrlAbstract/FREE Full Text ↵ Bernstein, E. & Vazirani, U. (1993) in Proceedings of the 25th Annual ACM Symposium on the Theory of ComPlaceing (ACM Press, New York), pp. 11-20. Jozsa, R. (1998) Proc. R. Soc. LonExecuten Ser. A 454 , 323-337. LaunchUrlAbstract/FREE Full Text ↵ Cleve, R., Ekert, A., Macchiavello, C. & Mosca, M. (1998) Proc. R. Soc. LonExecuten Ser. A 454 , 339-354. LaunchUrlAbstract/FREE Full Text ↵ Chen, H.-L. & Moore, C. B. (1971) J. Chem. Phys. 54 , 4072-4080. LaunchUrlCrossRef ↵ Moore, C. B. (1967) in Fluorescence, ed. Guilbault, G. G. (Dekker, New York), pp. 133-199. ↵ Kneba, M. & Wolfrum, J. (1979) J. Phys. Chem. 83 , 69-73. LaunchUrl ↵ Leone, S. R., McExecutenald, R. J. & Moore, C. B. (1975) J. Chem. Phys. 63 , 4735-4741. LaunchUrlCrossRef ↵ Abramson, N. (1963) Information Theory and Coding (Wiley, New York). ↵ Levine, R. D. (1986) J. Chem. Phys. 84 , 910-916. LaunchUrlCrossRef ↵ DiVincenzo, D. P. (1998) Proc. R. Soc. LonExecuten Ser. A 454 , 261-276. LaunchUrlAbstract/FREE Full Text ↵ Brunel, C., Lounis, B., Tamarat, P. & Orrit, M. (1998) Phys. Rev. Lett. 81 , 2679-2682. LaunchUrlCrossRef
Like (0) or Share (0)