format: text

text:
# MATHEMATICAL PHYSICS AND LIFE11footnote 1To appear in Mathematical Sciences Series, Vol.4: Computing and Information Sciences: Recent Trends, ed. J.C. Misra, Narosa (2003) 270-293, quant-ph/0202022.

## 1 What is Life?

Often life is characterised on the basis of its fundamental processes: metabolism and reproduction. But perhaps it is easier to characterise what life is not than to define what it is. Life is not a process of thermodynamic equilibrium—all the metabolic activities are driven by a continuous supply of energy which is obtained from the environment as food. Evolution of life is not dictated by precise rules of logic—it is governed by random events where memory plays a role but there is no foresight. Living organisms are made of neither completely regular arrangements as in a solid nor totally random arrangements as in a gas—all the activities take place in a liquid medium with continuous jostling amongst various biomolecules. These properties make the familiar tools of physicists—equilibrium dynamics, axiomatic deductions, periodic structures, perturbation theory—totally useless when it comes to understanding life. Life is termed too nonlinear, too complex, too unwieldy for the simple tools of physicists to handle.

Information theory provides a way out of this apparently hopeless conundrum. Information is the abstract mathematical concept that quantifies the notion of order amongst the building blocks of a message. It extracts just the right features from the apparently chaotic ensemble that is life, and allows the study of its manipulations without going into the nitty-gritty of all the details. Of course, figuring out the appropriate information in a given process requires detailed experiments and modeling. But once the essential features have been extracted, the rules of information processing are precise enough to allow their systematic analysis in a mathematical framework.

It can be argued that information processing is not just a particular characteristic of living organisms, but it is their most important characteristic . The processes of life are aimed towards survival of the organism. Even though we do not quite understand why the organisms want to perpetuate themselves, we have enough evidence to show that they use all available means for this purpose . The damage caused by the disturbances from the environment makes it impossible for a particular individual to survive forever, so the perpetuation is carried out through the process of replication. In this process, the knowledge gathered on how to survive is carried forward from one generation to the next. The pattern of acquiring, interpreting and passing on information, often after using and refining it, has occurred repeatedly during evolution. The hereditary genetic information is the most basic and primitive communication system. As living organisms evolved, the lower levels of information communication have been bypassed in favour of more efficient higher level mechanisms, reducing dependence on direct physical contact at every stage. Electrochemical signaling between cells and nervous systems evolved in multicellular organisms, teaching and imitation arose amongst parents and offsprings, verbal and sign languages originated amongst members of the same species, humans created books and libraries for long term storage of information, we use telecommunication and computers nowadays, and the future (already imagined in science fiction stories) is likely to witness fusion of brains with computers.

## 2 Biological Information

Living organisms absorb free energy from the environment to create order from disorder. All their struggles for raw materials and energy are ultimately for picking up necessary building blocks and organising them in a desired manner. After the organisms die, the building blocks gradually come apart and become random ensembles. Quantitatively the information carried by a set of building blocks is nothing but the entropy (up to an unimportant proportionality constant), and principles of computer science can be used to study its transformations.

Table 1 lists various stages of information processing tasks carried out by living organisms and compares them with similar tasks carried out by our electronic computers. They have been arranged in a hierarchy of levels according to the physical scales involved. Different levels typically have different languages, and the processing translates the information from one level to the next. We are accustomed to looking at our computers from the top level down—from the abstract mathematical equations to the transistors embedded in the silicon chips. On the other hand, living organisms have evolved from the bottom level up—from the biomolecular interactions to multicellular systems.

To understand the working of an information processing system, one must figure out what happens at each level step by step. But where should one begin, when little is known about the system? At this juncture, it is worthwhile to observe that two major simplifications occur as one proceeds from high level information processing to low level:
(1) The languages at the higher levels are abstract software. They get translated from one to another according to established conventions. On the contrary, the languages at the lowest levels are directly connected to the physical responses of the hardware. There is no more translation of the message—the interpretation of the signal is built into the design of the system and not left to an external agency. (e.g. we can programme the mathematical formulae into the computer in a variety of notations, but ultimately they are all converted to pulses of voltages and currents because the transistors in the silicon chips cannot respond to any other language.)
(2) The abstract high levels allow a variety of subjective choices in the languages. That leads to lots of variations, historical adaptations and a large number of possible instructions. At the lowest level, only a handful of instructions related to physical responses of the hardware are possible, and the language tends to become universal. (e.g. for writing the computer programmes, we may choose between Fortran and C, or Unix and Windows, but finally all that is reduced to the operations of Boolean logic.)

These properties make it obvious that it is easiest to decipher an information processing system at the lowest level, where the physical properties of the hardware dictate the language and the instructions . Complex systems are then generated by putting together a large number of simple ingredients. With the progressive miniaturisation of the silicon chips, the lowest level of information processing in our computers nowadays is at the scale of a micron. The lowest level of information processing for the living organisms is at the molecular scale—a scale smaller by a factor of thousand. The molecules involved there are DNA and proteins, and their language is universal—the same all the way from viruses and bacteria to human beings. So our first objective is to identify the information processing tasks carried out by DNA and proteins, and then study the relation between these tasks and physical properties of the molecules. It is important to note that biological systems are not general purpose computers, where the same CPU handles all the tasks. Living organisms have evolved specialised components to perform specific tasks, which makes it easy to identify the task of a component and study its implementation in detail.

DNA is the read-only-memory of the living organisms. It is a linear chain constructed from an alphabet of 4 nucleotide bases. Most of the time, the information carried by it remains nicely protected in the double helical structure. Whenever necessary, the information is “copied” into another DNA molecule during replication and mRNA molecule during transcription. The “copying” is not literal, rather it follows rules of complementary base-pairing. The “copy” is assembled on the master template, by picking up desired nucleotide bases one by one from the surroundings and joining them together in a chain. The available nucleotide bases are floating around randomly in the cell, and the base-pairing rules decide which one is the correct one to insert at a specific location. Thus the only task carried out by DNA is the sequential assembly of a chain of nucleotide bases, selected from a random ensemble.

Proteins carry out almost all the processes required by the living cell, by binding to various molecules. The binding is highly specific, very much like a lock and key mechanism. For this purpose, proteins contain structurally stable features, precisely located on the surface. Each protein has its own unique shape and participates in its own unique process. Proteins are made of one or more tightly folded polypeptide chains (they may contain other ingredients but the role of these other ingredients is essentially chemical and not structural). Polypeptide chains are synthesised from an alphabet of 20 amino acids, and other ingredients get added afterwards. Each polypeptide chain is assembled on the mRNA template during translation when successive segments of 3 nucleotide bases are mapped onto individual amino acids. Every chain subsequently folds up into a three-dimensional structure, uniquely determined by its amino acid sequence. The proteins are thus involved in two information processing tasks. One is similar to that of DNA, i.e. sequential assembly of a chain of amino acids selected from a random ensemble. The other is to create a multitude of three-dimensional structures by folding up linear polypeptide chains.

These tasks are carried out at the molecular scale. We know the physical laws applicable there—classical dynamics is relevant, but atomic structure and quantum dynamics cannot be ignored. We can now pose the question: were we to design a processor to carry out the tasks of DNA and proteins, what design would we come up with knowing all that we do about the physical laws? The rest of the article investigates this question, and compares the results to what is found in nature.

## 3 Survival of the Fittest

It is advantageous to process the information efficiently, and not in any haphazard manner. The first attempt may not provide the best solution to a problem, but the attempts do not stop there. Technological developments are driven as much by new inventions as by untiring attempts to improve and optimise earlier solutions. In general, information processing is optimised following two guidelines: minimisation of physical resources (time, space, energy etc.) and minimisation of errors. These guidelines often impose conflicting demands on the solution, and specific trade-offs are made depending on the details of the problem.

In case of living organisms, the optimisation criteria are paraphrased as “Darwinian evolution”. The environment and competition drive the living organisms to adapt to them by exerting selection pressures, and we can illustrate that by many examples. We now understand the genetic basis behind Darwinian evolution. No information processing system can be perfect, and occasional errors in genetic information processing produce random mutations of living organisms. The error rate is bound to be small in any reasonably stable system, and mutations are small local fluctuations in DNA molecules. The mutated organism is in essentially the same environment as the original one, and both have to compete for the available resources. If the mutation improves the ability of the organism to survive, the mutated organism grows in number, otherwise it fades away. To what extent may all this be quantified?

Let the index $`i`$ label a set of coexisting species in a given environment, and let $`\varphi _i(t)`$ denote their populations at time $`t`$. The simplest evolution scenario is where the future population of each species depends linearly on the present populations of all the species,

$$\varphi _i(t+1)=\underset{j}{}M_{ij}\varphi _j(t),\varphi _i(t)0.$$
(1)
The diagonal terms $`M_{ii}`$ represent the individual rates of growth, while the off-diagonal terms $`M_{ij}`$ represent interactions between species. Every species consumes resources from the environment; when the available resources are limited, the total population cannot exceed a certain value (a normalisation convention can be chosen such that a unit of each species utilises the same amount of resources). Evolution gets really competitive when the total population reaches its maximum value. In this stage,

$$\underset{i}{}\varphi _i(t)=N,\underset{i}{}M_{ij}=1.$$
(2)
Evolutionary models obeying Eqs.(1,2) have often been used together with the constraint $`0M_{ij}1`$. They describe Markovian evolution in the language of classical probability theory. With the constraints on $`M_{ij}`$, one can prove many inequalities and convergence properties. This is not the interesting part of evolution, however.

Fixed amount of available resources correspond to conservation laws. The most general evolution in such circumstances is described by “orthogonal transformations” . The generators for these transformations are antisymmetric matrices. Thus in a general evolutionary setting, it is more appropriate to let $`M_{ij}`$ take positive as well as negative values. As a matter of fact, situations of both positive and negative feedback occur routinely in biological systems (e.g. catalysis and inhibition, symbiosis and parasitic behaviour, defence mechanisms and cancer, etc.). When the resources are limited, one gains only at the expense of someone else—the formal setting is called “zero-sum games”.

We can now look at how evolution changes, when the range of $`M_{ij}`$ is extended to include negative values. Clearly extending the range of $`M_{ij}`$ cannot make the evolutionary algorithm any less efficient. On the contrary, recent developments in quantum computation  offer a hint. Quantum algorithms exploit two features to beat their classical counterparts, superposition of states and cleverly designed destructive interference. Superposition means letting all the states be in the same place at the same time, and it can reduce the spatial resources required for the algorithm exponentially. But superposition is not an option available to coexisting biological species. Even in absence of superposition, destructive interference can be used to eliminate undesired states and to reach the output state more quickly. The execution time is reduced at least by a constant factor if not polynomially. Biological evolution occurs over long time scales, and even a tiny change in the rate of growth is important because it can translate into exponential changes in populations over a long time. Negative values of $`M_{ij}`$ can indeed be interpreted as effects of destructive interference, which help the species reach their asymptotic populations faster. The lessons learnt from quantum computation then tell us that “orthogonal transformations” offer a quicker way of picking a winner amongst the contenders. These arguments are not a rigorous derivation, but they allow us to infer that competition beats altruism and Darwinian selection is a consequence of limited availability of resources .

Exercise: Study evolution algorithms governed by Eqs.(1,2). Estimate the speed-up obtained when $`M_{ij}`$ are allowed to become negative compared to when they are restricted to be positive.

Having justified that optimisation principles do play a role in biological systems, let us analyse what is optimised and how. In designing more and more efficient computers, we have explored the criteria mentioned at the beginning of this section. To optimise spatial resources, we need elementary hardware components that are simple and easily available, and yet versatile enough to be connected together in many different ways. This is the typical choice made at the lowest level of information processing, and complicated systems are then constructed by packing a large number of components in a small volume. It helps to have a small number of elementary components, because that reduces the number of possible instructions and the connections amongst the components. With a small instruction set, individual steps can be implemented quickly and the time required for processing information is cut down. In addition to these, in the implementation of any given task, we have to find a software algorithm which requires smallest number of components and smallest number of execution steps.

For our computers, it is also necessary to minimise the energy consumption during processing. But this feature is surprisingly absent in biological systems. The reason is that the elementary components of biological systems are so simple and cheap, that they carry out their tasks with very little energy—a single biomolecular interaction requires a million times less energy than a Boolean logic operation with modern silicon transistors. As a result, biological systems often exhibit a wasteful feature. Millions of eggs and pollen grains are produced when a few would have sufficed to propagate the species in a secure environment. Such overkills strengthen the competition and enforce survival of the fittest.

Exercise: Observe that biological systems have evolved excellent amplifiers for processing occurring at the molecular scale. Eyes and chlorophyll can detect a few photons, smells can be identified with a few molecules and sound amplitudes with a fraction of atomic size can be heard. Body movements almost always use levers in the configuration of mechanical disadvantage; it is somewhat ironic that mechanical advantage corresponds to gain in power which is the reciprocal of the gain in amplification. Obtain quantitative data for such amplification processes.

Faithful information processing requires that the error rate must be kept in control. One strategy is to shield the system from unwanted external disturbances, but that does not protect against internal fluctuations. Internal errors are minimised by selecting an information processing language based on discrete variables (as opposed to continuous variables). Allowed values of fundamental physical variables are often continuous, in which case a set of non-overlapping intervals of values can be chosen as the discrete variables. This is the common procedure of digitisation. The advantage is that the discrete variables remain unaffected, even when the underlying continuous variables drift, as long as the drifts keep the values within the assigned intervals. As a matter of fact, digitisation allows elimination of bounded errors. When the discrete variables are chosen as far apart from each other as possible in a given range of values, misidentification is minimised and one has the largest protection against errors. Even with digitisation, errors from large fluctuations cannot be eliminated. It is a curious fact that biological systems take advantage of a tiny error rate—with too many errors the organism will not be able to survive, but without mutations there will be no evolution.

## 4 Choice of Language

We are now in a position to look at some examples of information processing systems, and understand how well they implement the optimisation principles. Table 2 lists some of the languages used by living organisms, their basic components and the operations performed on those components. Messages are constructed by linking the basic components—the building blocks of the language—in a variety of arrangements. The information contained in a message depends on the values and locations of the building blocks. Any language that communicates non-trivial information must have the flexibility to arrange its building blocks in different ways to represent different messages.

This notion is quantified by saying that messages are aperiodic chains of the building blocks, and the information contained in a message is its entropy, i.e. a measure of the number of possible forms the message could have taken . This definition tells us that the information content of a message can be increased by eliminating correlations from it and making it more random. It also tells us that local errors in a message can be corrected by building long range correlations into it. But it does not tell us what building blocks are appropriate for a particular message. The choice of building blocks depends on the type of information and not on the amount of information.

A striking feature of all the types of information listed above is that they all have digital form. The simplicity of instructions and the control over errors offered by digitisation are too important to ignore. The fact that digitisation necessarily approximates values of continuous variables does not matter in practical applications. The reason is that no practical result needs infinite precision; as long as the results are obtained within predefined but non-zero tolerance limits they are useful . The outstanding optimisation question then is to figure out the best way of digitising a message, i.e. what should be selected as the building blocks of the aperiodic chain.

When the languages are versatile enough, information can be translated from one language into another by replacing one set of building blocks by another. Nonetheless, physical principles are involved in selecting different building blocks for different information processing tasks. For example, our electronic computers compute using electrical signals but store the results on the disk using magnetic signals; the former encoding is suitable for fast processing while the latter is suitable for long term storage. It is therefore the available hardware and the job to be carried out which dictate the optimal building blocks in any information processing task.

Languages of biological systems have arisen by trial and error exploration, and by accumulation of small modifications, and not by mathematical deduction of the optimal choice. Historical baggage is often carried forward in such situations; one gets stuck in local optima without discovering the globally optimal features. The optimal solution can be reached only when there is a long enough time for exploration . Biologists have often viewed the languages of genes and proteins as “frozen accident”—it is such a vital part of life that any change in it would be highly deleterious . On the other hand, these languages are universal across all living organisms and exhibit efficient packing of information (little correlation or close to maximal entropy), which suggest that the optimal form may have been reached. Keeping all these ideas at the back of our mind, let us step by step search for the optimal way of implementing the tasks of DNA and proteins, and compare that to what is observed in nature.

Exercise: Understand the physical reasons behind various numerical systems we use. Why do we have binary arithmetic for computers, decimal system for counting, and even more complicated system for marking time?

## 5 Structure of Proteins

It is convenient to analyse the language of proteins first, since it is a straightforward exercise in classical geometry . The problem is to find a discrete language than can encode arbitrary shapes in 3-dimensional space, much like designing children’s playing blocks so that they can be linked together to create many different structures. This means discretising the continuous operations of translation and rotation. At the molecular scale, the atomic structure of matter automatically provides such a discretisation of space, i.e. a lattice. Among all possible lattice discretisations, the simplest lattice unit that can faithfully describe the space in any dimension is the simplex—the volume element formed by a set of $`(d+1)`$ points in $`d`$-dimensions .

The changes required to go from the familiar language of 1-dimensional numerical information to the language of 3-dimensional structural information are summarised in Table 3. We can quickly observe the following facts:
$``$ In 3-dimensional space, the simplex is a tetrahedron. One can go from any vertex to any other vertex by a single step. Arbitrary structures can be constructed by joining tetrahedra together.
$``$ A regular tetrahedron describes four equivalent and equidistant states. They correspond to the smallest two representations of the spherical harmonics, $`l=0,1`$. This is the obvious generalisation of the 1-dimensional binary arithmetic to the 3-dimensional space. In quantum theory, these four states are formed by $`sp^3`$-hybridisation of atomic orbitals and provide an orthogonal basis.
$``$ The ideal chemical element for realising tetrahedral geometry in molecular structures is carbon, in the form of a diamond lattice. Carbon has the capability to form aperiodic chains, where different side chemical groups hang onto a backbone, which is a must for encoding information.
$``$ The carbon atom is located at the centre of the tetrahedron, with its four bonds directed towards the four vertices. The tetrahedral symmetry group $`T_d`$ has 24 elements, which can be factored into a group $`T`$ of 12 proper rotations and reflection (or parity). The proper rotations consist of rotations around four 3-fold axes and three 2-fold axes, and describe rotations around bonds of the central carbon atom. The parity transformation flips the chirality of the molecule, and is rare enough to be ignored in the structural analysis of proteins.
$``$ The tetrahedral geometry also plays a crucial role in the properties of water—the essential solvent for all the processes of life. The proteins are often hydrated, with water molecules filling the cavities. The role of water, however, will not be pursued here.

Arbitrary irregular shapes can be constructed by gluing tetrahedra together, similar to how children assemble their playing blocks. But such a process has to rely on an external agency to stop the growing process when the desired structure has been reached. At the lowest level of information processing, there is no luxury of an external agency—the building blocks themselves have to carry the information about the required shape. The solution is to encode the 3-dimensional structural information as a 1-dimensional chain that knows how to bend and fold at every step. (The assembly problem, which cannot do without an external agency, then simplifies—only 1-dimensional chains have to be assembled, and we will look at that later.) This is the solution chosen by nature. The other advantage of constructing proteins as chains which can fold and unfold is that they can cross membranes and cell walls, which they need to do to participate in various cellular processes, without making big holes in them. It is a marvelous piece of engineering that proteins unfold to their chain form, cross the membrane through a small hole that prevents other objects from leaking, and then fold again into their characteristic shapes.

### 5.1 The polypeptide backbone

To have the capability to fold and unfold the backbone of the chain must be strong and its side interactions weak. The simplest backbone formed by a string of carbon atoms is that of polyethylene (i.e. $`(CH_2)_n`$). It is made of rotatable carbon-carbon single bonds, and is too flexible to hold its shape. The rigidity of the backbone can be increased by including in it some double bonds that cannot rotate; the backbone will then have stiff segments alternating with flexible joints similar to a chain of metal rings. The simplest such non-planar backbone is formed by alternating one double bond with two single bonds. That is the backbone of polypeptide chains—the stiff $`CN`$ peptide bond increases its rigidity, while the rotatable $`C_\alpha `$ bonds permit construction of a variety of 3-dimensional structures. The polypeptide chain is synthesised by joining together amino acids, and their chemical structures are shown in Fig.1. Every amino acid contains an acidic $`COOH`$ group and a basic $`NH_2`$ group, and the peptide bond is formed by their acid-base neutralisation. The amino acids differ from each other by their R-groups; these R-groups hang on the sides of the polypeptide backbone and interactions amongst them decide how the backbone would bend and fold.

The next step is to figure out the number of ways a polypeptide chain can be folded on the diamond lattice. Along the backbones of real polypeptide chains, all the bonds are not of the same length, all the bond angles are not exactly equal to the tetrahedral value (in particular the peptide bond is planar), and the nitrogen atoms do not have exactly the same properties as the carbon atoms . The folding of the polypeptide chain on the regular diamond lattice, therefore, is necessarily approximate. But the variations of bond lengths and angles are within 10%, and so we can expect to see features of the diamond lattice in the folded polypeptide chains.

The diamond lattice is a face-centred cubic lattice with a two-point basis. It can accommodate the polypeptide chain reasonably well, with the peptide bonds in the “trans” configuration. The rare “cis” configuration of the peptide bond can be accommodated too, using the well-known local switch between the face-centred cubic lattice and the hexagonal lattice. (Most of the natural peptide bonds occur in “trans” configuration, but a few percent occur in the “cis” configuration.) All possible configurations of the polypeptide chain can be enumerated, by specifying for every peptide bond the location of the next peptide bond in the chain. Each peptide unit of the chain contains two rotatable bonds of the $`C_\alpha `$ atom, labeled by angles $`\varphi `$ and $`\psi `$ in Fig.1. The tetrahedral geometry gives three possible values for each of these angles, separated by $`120^{}`$ steps. There are thus nine possible orientations of each peptide unit in the chain.

The Ramachandran map describes the allowed values of the angles $`\varphi `$ and $`\psi `$ in real polypeptide chains, after taking into account hard core repulsion between the atoms of the polypeptide backbone . As Fig.2 shows, even with all the variations in bond lengths and angles, the discrete diamond lattice approximation is not a bad starting point for describing the preferred orientations of the angles $`\varphi `$ and $`\psi `$. (The Ramachandran map for achiral glycine is different. Unlike the other chiral L-type amino acids, glycine can occupy the region around $`\varphi =60^{},\psi =60^{}`$.) To put it differently, though nature has not been able to perfectly implement the ideal diamond lattice geometry of structural design, evolution has taken her impressively close to it. It is also worthwhile to observe that in general there are many ways a folded chain can cover a 3-dimensional shape. So it is not necessary that all the nine angular orientations are occupied equally—it is sufficient to have all the orientations as possibilities.

The nine angular orientations and the “trans-cis” option exhaust the “elementary logic operations” for embedding the polypeptide backbone on the diamond lattice, i.e. by implementing these operations one can fold the backbone on the lattice in any desired configuration. Real polypeptide chains form other long distance connections, hydrogen and disulfide bonds, when they fold. While these are important for structural stability of the proteins, they do not give rise to new configurational possibilities.

Exercise: In two dimensions, the simplex is a triangle and the simplest connectivity leads to a honeycomb structure. At the molecular scale, such a geometry is also realised by carbon in the form of graphite sheets. Figure out the building blocks required to fold a chain on the graphite lattice in any arbitrary manner.

### 5.2 The R-groups of amino acids

For each polypeptide backbone configuration described above, there are two remaining directions for other chemical groups to attach to the $`C_\alpha `$ atom. One is attached to the R-group, while the other to a hydrogen atom. Two arrangements are possible, and they correspond to opposite chirality. All the amino acids naturally occurring in proteins use only one of these two possibilities, L-type arrangement . Altogether, therefore, there are only nine discrete orientations (and an occasional “trans-cis” switch) for adding a new amino acid to an existing polypeptide chain.

There is no clear association of any individual amino acid with the nine discrete points in the Ramachandran map. Indeed, just the structure of a particular amino acid does not decide its orientation in the polypeptide chain; rather its orientation is fixed by the overall interactions of its R-group with those that precede it and those that follow it. In other words, the structural code of the polypeptide chain is an overlapping one. Even with an overlapping code, every time an amino acid is added to the polypeptide chain, the orientation of one amino acid gets decided. This conservation of information implies that a maximally efficient overlapping code still needs at least nine amino acids (may be 10 to take care of the “trans-cis” option and long distance connections) to construct polypeptide chains that can cover arbitrary shapes. The R-group properties of amino acids have been studied in detail—polar and non-polar, positive and negative charge, straight chains and rings, short and long chains, and so on. Still which sequence of amino acids will lead to which conformation of the polypeptide chain is an exercise, in coding as well as in chemical properties, that has not been solved yet.

At this stage, it is instructive to observe that the 20 amino acids are divided into two classes of 10 each, according to the properties of their aminoacyl-tRNA synthetases . tRNA molecules, with 3 nucleotide bases of anticodon at one end and an amino acid at the other end, translate the information from mRNA molecules to polypeptide chains. Aminoacyl-tRNA synthetases are the bilingual molecules that ensure correct translation of this information, by attaching matching amino acids to tRNA molecules as per their anticodons. The two classes of synthetases totally differ from each other in their active sites and in how they attach amino acids to the tRNA molecules. The lack of any apparent relationship between the two classes of synthetases has led to the conjecture that the two classes evolved independently, and early forms of life could have existed with proteins made up of only 10 amino acids of one type or the other. An inspection of the R-group properties of amino acids in the two classes reveals that each property (polar, non-polar, ring/aromatic, positive and negative charge) is equally divided amongst the two classes, as shown in Table 4. Not only that, the heavier amino acids with each property belong to class I, while the lighter ones belong to class II. This clear binary partition of the amino acids according to the size of their side chains has unambiguous structural significance. The diamond lattice structure is quite loosely packed with many cavities of different sizes. The use of long R-groups to fill up big cavities and short R-groups to fill up small ones can produce a dense compact structure. Proteins indeed have a packing fraction similar to closest packing of identical spheres.

We have thus arrived at a structural explanation for the 20 amino acids as building blocks of proteins, a factor of 10 for folding the polypeptide backbone and a factor of 2 for the size of the R-group . This explanation can be made more concrete by identifying subsequences of amino acids corresponding to different orientations of the polypeptide backbone. We do not have a clear criterion regarding how long the amino acid sequence should be before it assumes a definite orientation, although the short range nature of molecular interactions favours relatively local folding rules. The information accumulated in protein structure databases should help in a detailed analysis.

Exercise: Devise a scheme to discretise the angular distributions of amino acids on the Ramachandran map. Check how the distributions narrow down, when preceding and following amino acids are fixed, step by step.

## 6 Language of DNA

Now we can analyse the other task of synthesising linear chains by joining together the building blocks . This task needs an external agency—called an oracle in the language of computer science—to dictate the required order of building blocks. DNA and polypeptide chains are assembled in presence of preexisting templates; the order of building blocks specified by the template is the blueprint or the master copy.

To understand the optimisation of this process, it is helpful to first look at a familiar game played by school children. In the game, two teams compete to find the names of famous persons. One team selects the name of a famous person and gives it to the referee. The other team has to discover this name by asking a set of questions to the first team. The game is made interesting by the condition that the first team only provides “yes or no” answers to the questions, i.e. the minimal amount of information—one bit—is released in response to every question. The two teams take turns choosing the names and asking the questions. After several rounds, the team which succeeds in discovering the names with a smaller number of questions is the winner.

The players quickly learn that specific questions such as “Is the person Albert Einstein?” are inefficient. They fail most of the time, and when they fail one does not learn much about who the person is. Efficient questions are of the type “Is the person a man or a woman?”, whence there is a substantial reduction in the number of possible names for the next step no matter what the answer to the question is. The best questions are the ones which reduce the number of possible names for the next step by a factor of two. This factor of two reduction is of course a consequence of the fact that the answers provided to the questions are binary.

The computer science paradigm for this game is “database search”. Binary search is the optimal classical algorithm, which finds the desired item in the database using $`\mathrm{log}_2N`$ binary questions. This algorithm can be used when the items in the database are sorted according to some order, so that at every step the items corresponding to “yes” answer can be easily separated from the items corresponding to “no”. Examples are words in a dictionary and names in a telephone directory. On the other hand, if the database is not sorted, e.g. numbers in a telephone directory, then there is no fast search algorithm. The best option is to go through all the items one by one, and that requires $`N/2`$ binary questions on the average.

Now we can return to look at how the DNA and polypeptide chains are assembled by biochemical processes. The questions are the molecular bonds involved in nucleotide base-pairing, the preexisting template is the “other team”, and the answers are binary—either the bonds form or they do not. The available database is unsorted—the building blocks are randomly floating around in the cellular environment. Classically, two nucleotide bases (one complementary pair) would be sufficient to encode genetic information; the number of more complicated building blocks would then follow powers of two. Our computers indeed encode information this way. Moreover, starting from scratch, two nucleotide bases would have evolved before four nucleotide bases. So why did living organisms give up the simple powers of two and evolve more complicated languages?

A surprising answer is provided by the quantum database search algorithm, whose optimisation features are different from their classical analogues.

### 6.1 Quantum database search

Quantum states are unit vectors in a Hilbert space (linear vector space with complex coefficients). They evolve in time by unitary transformations. These properties are totally different from classical states and their evolution, and they form the basis of the subject of quantum computation. Of course, to make contact with our problems defined in classical language, we require a mapping between classical and quantum states. That is achieved by identifying the orthogonal basis vectors of the Hilbert space with the set of distinct classical states. The complex components of a general vector in the Hilbert space can vary continuously, and they are called amplitudes of the quantum state. States with more than one non-zero amplitudes are called superposition states. Quantum algorithms evolve the amplitudes from some initial values to some final values, by a sequence of unitary transformations. The measurement probability of obtaining a particular classical result from a quantum state is given by the absolute value square of the corresponding amplitude.

The quantum database search algorithm works in an $`N`$-dimensional Hilbert space, whose basis vectors are identified with individual items. It takes an initial state whose amplitudes are uniformly distributed over all items, to a final state where all but one amplitudes vanish. Following Dirac’s notation, let $`\{|i\}`$ be the set of basis vectors, $`|s`$ the initial uniform superposition state, and $`|b`$ the final state corresponding to the desired item. Then

$$|i|s|=1/\sqrt{N},i|b=\delta _{ib}.$$
(3)
Since superposition of amplitudes is commutative, the quantum database can be taken to be unsorted without any loss of generality. The optimal solution to the quantum search problem  is based on two properties: (i) the shortest path between two points on the unitary sphere is the geodesic great circle, and (ii) the largest step one can take in a given direction, consistent with unitarity, is the reflection operation. The available reflection operators are,

$$U_b=12|bb|,U_s=2|ss|1.$$
(4)
The former is the binary quantum query for the desired state, while the latter is the sign flip relative to the initial uniform state. The quantum algorithm is the discrete Trotter’s formula with these two operators, which locates the desired item in the database with $`Q`$ queries,

$$(U_sU_b)^Q|s=|b.$$
(5)
This equation is readily solved to give the relation

$$(2Q+1)\mathrm{sin}^1(1/\sqrt{N})=\pi /2.$$
(6)
For a given $`N`$, the solution for $`Q`$ satisfying Eq.(6) may not be an integer. This means that the algorithm will have to stop without the final state being exactly $`|b`$. There will remain a small admixture of other amplitudes in the output, implying an error in the search process. The size of the unwanted admixture is determined by how close one can get to $`\pi /2`$ on the r.h.s. of Eq.(6). Apart from this, the algorithm is fully deterministic.

Exercise: Derive the result in Eq.(6). Show that it is the optimal solution.

For unsorted database search, the result in Eq.(6) is clearly superior to the corresponding classical result (the number of queries is asymptotically $`O(\sqrt{N})`$ vs. $`O(N)`$). The reason behind this improvement is the clever interference amongst superposed amplitudes, which amplifies the desired state amplitude at the expense of the rest. As a matter of fact, this quantum algorithm does not use the full power of quantum dynamics. It can very well be implemented using other systems that obey the superposition principle, e.g. classical waves or coupled pendulums . Fig.3 illustrates the individual steps of the algorithm in the simplest case.

The solutions of Eq.(6) for small values of $`Q`$ are of special significance for the number of building blocks involved in genetic information processing,

$`Q=1`$ $``$ $`N=4,`$
$`Q=2`$ $``$ $`N=10.5,`$
$`Q=3`$ $``$ $`N=20.2,`$ (7)
when the binary quantum query is identified with the nucleotide base-pairing. DNA has an alphabet of 4 nucleotide bases, the well-known triplet genetic code for synthesis of polypeptide chains carries 21 signals (for 20 amino acids plus STOP) , and a primitive doublet genetic code could have labeled the 10 amino acids of a single class (described in the previous section). It is fascinating that these numbers have come out of an algorithm that performs the actual task accomplished by DNA , and that too as optimal solutions . Obviously, the next job is to check whether DNA possesses the hardware components necessary to implement the quantum search algorithm.

### 6.2 A physical scenario

The execution of the quantum search algorithm requires three ingredients: an interaction providing the quantum oracle $`U_b`$, a mechanism to generate the superposition state $`|s`$, and an environment where quantum coherence can survive. Let us look at how these may be realised in the case of DNA.

The formation of a molecular bond is a quantum process. Two reactants come together in an excited state, and the combined system decays to the bound state with the emission of an energy quantum. The total energy is conserved in the process, but the quantum dynamics of the transition between the two molecular energy levels alters the phase of the quantum state by $`\sqrt{1}`$. This phase is a geometric phase, independent of the strength of the transition interaction and the time for the decay.

The DNA molecule is a double helix, and is often schematically represented as a ladder. The sides of the ladder have a periodic structure made of alternating sugar and phosphate groups, while the rungs of the ladder are made of paired nucleotide bases. The atomic structure of paired nucleotide bases is shown in Fig.4. The nucleotide bases are made of planar rings of atoms, and so the base-pairing requires at least two points of contact for structural stability of the helix (otherwise the bases would rotate relative to each other). The contacts between nucleotide bases are hydrogen bonds, which are formed by the quantum process of a proton tunneling between two attractive energy minima. When the nucleotide bases come together with arbitrary initial orientations, the base-pairing is likely to occur as a two step process—the first contact connects the two nucleotide bases and the second one locks their relative orientations. With such a two step base-pairing, the quantum state is transformed precisely by the operator $`U_b`$, i.e. the quantum state changes its sign when the nucleotide bases get paired and it remains unaffected when there is no pairing. From the known binding energy of the hydrogen bond and the uncertainty principle, the time scale for base-pairing can be estimated to be, $`t_b10^{14}`$ sec.

The quantum algorithm starts with the superposition state $`|s`$, so $`|s`$ has to be a stable equilibrium state. Any other initial state has to relax towards $`|s`$, and let $`t_r`$ be the time scale for this relaxation process. In case of DNA, $`|s`$ is a superposition state of chemically distinct nucleotide bases; it can be created only if the environment provides transition matrix elements between its various components. The nucleotide bases differ from each other by 5-10 atoms, and the transition matrix elements between them vanish in free space. Exchange of atomic groups is routine in chemical reactions, however—these reactions take place not with fixed number of atoms but with fixed chemical potential where exchange of atoms with the environment is allowed. (It is important to note that chemical reactions produce mixtures in classical environments while they produce superpositions in quantum environments.) The magnitudes of the transition matrix elements decide the speed at which $`|s`$ cycles through its various components. Let $`t_{osc}`$ be the time scale for these oscillations.

As illustrated in Fig.5, the quantum search algorithm can be executed when the above defined time scales satisfy the hierarchy,

$$t_b<<t_{osc}<<t_r.$$
(8)
Explicitly, (1) the environment relaxes the initial state to $`|s`$, (2) the sudden base-pairing changes the state to $`U_b|s`$ (the energy quantum is not released yet), (3) the relaxation process tries to bring the state back to $`|s`$ causing damped oscillations, and (4) the energy quantum becomes free to wander off at the opposite end of oscillation; its irreversible release confirms the base-pairing and stops the search algorithm, $`U_sU_b|s=|b`$.

This scenario is quite robust, and does not rely on fine-tuned parameters. The all important remaining question is whether the environment of DNA really has the properties ascribed to it, i.e. can it protect quantum coherence for a long enough time and create superposition states?

### 6.3 Decoherence and enzymes

High precision experiments have shown that quantum dynamics is the only valid description of the happenings at the atomic scale. Nevertheless, we hardly observe quantum effects in macroscopic systems; we find that macroscopic systems follow the quite different rules of classical dynamics. It cannot be discounted that the cross-over in the behaviour of physical systems is caused by some new laws of physics, as yet undiscovered. On the other hand, by a careful analysis of quantum dynamics, it is possible to explain the cross-over in a self-consistent manner. This explanation is labeled decoherence. It is a consequence of the fact that no quantum system is perfectly isolated; its inevitable interactions with its environment lead to irreversible loss of information and classical behaviour . Irrespective of any undiscovered laws of physics, therefore, decoherence must be controlled in order to observe quantum effects (e.g. in quantum computers).

It is impossible to give a complete quantitative description of decoherence, because sufficient details of the environment are not available in general. Still we can construct a formulation good enough to provide order-of-magnitude estimates. The dominant interactions of a system with its environment are molecular collisions and long range forces. The cross-sections for these processes can be calculated in terms of global properties of the environment, such as temperature and density. Then the rate at which the system relaxes towards the equilibrium state favoured by the environment is estimated according to Fermi’s golden rule. (This relaxation is often referred to as the “collapse of the quantum wavefunction”.) The relaxation rate is inversely proportional to three factors: the initial flux, the interaction strength and the final density of states. Typical decoherence times are extremely short, beyond direct observation. But we also know specific situations, where quantum states are long-lived due to suppression of one or more of these factors. For instance, low temperatures and shielding reduce the initial flux, the interaction strength is small for lasers and nuclear spins, and the final density of states is suppressed for superconductors and hydrogen bonds due to energy gap.

It must be noted that the relaxation rate cannot become arbitrarily small when these factors become large. It is well-known that classical waves cannot be damped faster than the critical rate of damping, given by their natural period of oscillation. If the damping is increased beyond the critical value, the system stops relaxing instead of relaxing faster. The quantum analogue of this is the Zeno effect, which points out that a continuously monitored system cannot decay. The reason behind this caveat is that Fermi’s golden rule is an approximation, not valid at times smaller than the natural oscillation period. It follows that, in the notation of the previous subsection, $`t_{osc}t_r`$ .

Biochemical reactions take place in a liquid medium at room temperature. In this environment, the decoherence times estimated following conventional thermodynamics and Fermi’s golden rule are miniscule. On the other hand, it is known that highly specific enzymes enhance the rates of biochemical reactions by orders of magnitude (as large as $`10^{12}`$) compared to estimates of kinetic theory and diffusion processes. Enzymes are catalysts, and the standard explanation of their properties is that they speed up the reactions by binding to the transition states and lowering the reaction barriers . What the transition states are and how the reaction barriers are lowered must be explained ultimately in terms of underlying physical laws. A transition state is something in between the reactants and the products, which can be represented by distorted quantum wavefunctions, i.e. a superposition. The intermediate transition state will not be stable on its own, but the enzymes are able to store and supply free energy necessary for stabilising it. The reduction of reaction barrier by stabilisation of superposition states can provide a large speed up, because the tunneling amplitudes depend exponentially on the barrier size. Thus it is entirely possible that the enzymes perform their job by exploiting the quantum superposition phenomenon.

DNA replication takes place only in the presence of polymerase enzymes, so we must focus on the kind of environment they provide. The polymerase enzymes are much larger than the nucleotide bases, and they completely enclose the reaction region during replication. This behaviour automatically shields the reaction from the flux of external disturbances and reduces the final density of states by limiting possible configurations (the allowed phase space is essentially one-dimensional and not three-dimensional). The hydrogen bonds also reduce the final density of states, since their binding energy ($`7kT`$) is considerably larger than the thermal fluctuations. The polymerase enzymes drive the replication process monotonically, using free energy from ATP molecules, which is another indication that the thermal fluctuations are overcome.

Altogether, it is not inconceivable for a cleverly designed polymerase enzyme to create superposition states and protect quantum coherence in accordance with Eq.(8). Moreover, certain aspects of the quantum replication algorithm can be tested experimentally: artificial DNA with varying number of nucleotide bases can show whether 4 is the optimal number of nucleotide bases or not, and isotopic tagging can track whether superposition of atoms takes place or not .

## 7 Outlook

Information theory provides a powerful framework for extracting essential features of complicated processes of life, and then analysing them in a systematic manner. The easiest processes to study are no doubt the ones at the lowest level. We have learnt a lot, both in computer science and in molecular biology, since their early days , and so we can now perform a much more detailed analysis. Physical theories often start out as effective theories, where the predictions of the theories depend on certain parameters. The values of the parameters have to be either assumed or taken from experiments; the effective theory cannot predict them. To understand why the parameters have the values they do, we have to go one level deeper—typically to smaller scales, e.g. viscosity of fluids and crystal shapes of solids can be understood based on atomic forces, structure of atoms can be understood based on properties of electrons and quarks, and so on. When the deeper level reduces the number of unknown parameters, we consider the theory to be more complete and satisfactory. The level below conventional molecular biology is spanned by atomic structure and quantum dynamics, and that is the natural place to look for reasons behind life’s “frozen accident”. It is indeed wonderful that sufficient ingredients exist at this deeper level to explain the frozen accident as the optimal solution.

Counting the number of building blocks in the languages of DNA and proteins is only the first step. The obvious next step to investigate is why the languages use only the observed building blocks and not other similar ones; only a subset of nucleotide bases and amino acids existing in living cells are used as the building blocks. The likely criteria for the selection of particular building blocks are simplicity (for easy availability and quick synthesis) and functional ability (for implementing the desired tasks). In case of proteins, the choices are connected to the solution of the protein folding problem, i.e. which way a particular amino acid sequence will fold. In case of DNA, the choices are possibly connected with the magnitudes of the transition matrix elements. We do not have the answers yet, but we may not have to wait very long. The experimental techniques and the information collected in databases have reached a stage, where it is possible to form hypotheses and then test them. It would certainly be worthwhile to test the arguments presented here in detail, and then build on them to understand more and more complicated processes of life.

The results described here also allow us to speculate about the origin of life. All the complicated processes of life did not come about in one go. The lowest level of information processing not only arose first, but in all likelihood it was much simpler than what exists today. Importance of functionality over memory would put proteins before DNA. Similar codons for similar amino acids and wobble rules found in the present genetic code are possible relics of an earlier simpler system . The solutions of the quantum search algorithm and the amino acid classes suggest that the present triplet genetic code was preceded by a doublet one. (This idea is reinforced by the accidental degeneracy, where 20 amino acids can be coded either by one classical and two quantum queries or by three quantum queries.) A still simpler precursor to the doublet code would be a singlet code of 4 amino acids. In fact such an evolutionary route for the genetic code has been proposed, just based on the chemical properties of the amino acids, GNC $``$ SNS $``$ NNN .

Unraveling the mysteries of life is certainly exciting. It is not merely a theoretical adventure; it has immense potential applications. Medical science will obviously benefit if we learn how to control processes of life that go wrong. Nanotechnology would benefit as well if we learn from biology how to implement specific tasks at the molecular scale.

## Acknowledgements

Many persons, too numerous to name individually, have helped in this study by guiding my thoughts in various directions. I am grateful to the Center for Computational Physics, University of Tsukuba, Japan, and the Optical Physics Research group, Lucent Technologies, USA, for their hospitality during part of this work.


text:
# Scattering in Two Black Hole Moduli Space

## 1 Introduction

Recently the study of black hole moduli space has attracted much attention. Quantum black holes have been studied by means of quantum fields or strings interacting with a single black hole. In the past few years the new quantum mechanics of an arbitrary number $`N`$ of supersymmetric black holes has been focused. Configurations of $`N`$ static black holes parametrize a moduli space. The low-lying quantum states of the system are governed by quantum mechanics on the moduli space. The effective theories of quantum mechanics on a moduli space of Reissner-Nordström multi-black holes were constructed in . Recently, (super)conformal quantum mechanics are constructed on a moduli spaces of five-dimensional multi-black holes in the near-horizon limit in .

The motivation behind these works is the hope that information of these quantum states will lead to the black hole entropy. The quantum states supported in the near-horizon region can be interpreted as internal states of black holes, and the number of such states is related to the black hole entropy. The other motivation is the expectation that investigation of these moduli spaces will lead to the understanding of $`\mathrm{AdS}_2/\mathrm{CFT}_1`$ correspondence and unravel some novel features behind the quantum states in multi-black hole mechanics.

The geometry of black hole moduli spaces was first discussed by Ferrell and Eardley in four dimensions . Further the black hole moduli spaces geometry with dilaton coupling in $`N+1`$ dimensions was discussed by Shiraishi . In these works the structure of the moduli space geometry is different according to dimensions and values of dilaton coupling.

In this work, we discuss quantum mechanics on the moduli space consisting of two maximally charged dilaton black holes. We study the quantum effect resulting from the different structure of the moduli space geometry in the scattering process. In the section 2, we discuss the moduli space structure of two black hole system with dilaton coupling in any dimensions. In the section 3, we study the quantum mechanics on this moduli space and the general view of potential in the (3+1) dimensions. In the section 4, we consider the scattering process on the moduli space. Then we discuss the quantum effects from the different structures of the moduli space geometry. In the section 5, we will give the conclusion and the discussion.

## 2 The Moduli Space Metric for the System Consisting of Maximally Charged Dilaton Black Holes

The Einstein-Maxwell-dilaton system contains a dilaton field $`\varphi `$ coupled to a $`U(1)`$ gauge field $`A_\mu `$ beside the Einstein-Hilbert gravity. In the $`N+1`$ dimensions $`(N3)`$, the action for the fields with particle sources is

$`S`$ $`=`$ $`{\displaystyle d^{N+1}x\frac{\sqrt{g}}{16\pi }\left[R\frac{4}{N1}(\varphi )^2e^{\{4a/(N1)\}\varphi }F^2\right]}`$ (1)
$`{\displaystyle \underset{i=1}{\overset{n}{}}}{\displaystyle 𝑑s_i\left(m_ie^{\{4a/(N1)\}\varphi }+Q_iA_\mu \frac{dx_i^\mu }{ds_i}\right)},`$
where $`R`$ is the scalar curvature and $`F_{\mu \nu }=_\mu A_\nu _\nu A_\mu `$. We set the Newton constant $`G=1`$. The dilaton coupling constant $`a`$ can be assumed to be a positive value.

The metric for the $`N`$-body system of maximally-charged dilaton black holes has been known as

$$ds^2=U^2(𝐱)dt^2+U^{2/(N2)}(𝐱)d𝐱^2,$$
(2)
where

$`U(𝐱)`$ $`=`$ $`(F(𝐱))^{(N2)/(N2+a^2)},`$ (3)
$`F(𝐱)`$ $`=`$ $`1+{\displaystyle \underset{i=1}{\overset{n}{}}}{\displaystyle \frac{\mu _i}{(N1)|𝐱𝐱_i|^{N2}}}.`$ (4)
Using these expressions, the vector one form and dilaton configuration are written as

$`A=\sqrt{{\displaystyle \frac{N1}{2(N2+a^2)}}}\left(F(𝐱)\right)^1dt,`$ (5)
$`e^{4a\varphi /(N1)}=(F(𝐱))^{2a^2/(N2+a^2)}.`$ (6)
In this solution, the asymptotic value of $`\varphi `$ is fixed to be zero.

The electric charge $`Q_i`$ of each black hole are associated with the corresponding mass $`m_i`$ by

$`m_i`$ $`=`$ $`{\displaystyle \frac{A_{N1}(N1)}{8\pi (N2+a^2)}}\mu _i,`$ (7)
$`|Q_i|`$ $`=`$ $`\sqrt{{\displaystyle \frac{N1}{2(N2+a^2)}}}\mu _i,`$ (8)
where $`A_{N1}=2\pi ^{n/2}/\mathrm{\Gamma }(\frac{1}{2}N)`$.

The perturbed metric and potential can be written in the form

$`ds^2=U^2(𝐱)dt^2+2𝐍d𝐱dt+U^{2/(N2)}(𝐱)d𝐱^2,`$ (9)
$`A=\sqrt{{\displaystyle \frac{N1}{2(N2+a^2)}}}(F(𝐱))^1dt+𝐀d𝐱,`$ (10)
where $`U(𝐱)`$ and $`F(𝐱)`$ are defined by (3) and (4). We have only to solve linearized equations with perturbed sources up to $`O(v)`$ for $`N_i`$ and $`A_i`$. (Here $`v`$ represents the velocity of the black hole as a point source.) We should note that each source plays the role of a maximally charged dilaton black hole.

Solving the Einstein-Maxwell equations and substituting the solutions, the perturbed dilaton field and sources to the action (1) with proper boundary terms, we get the effective Lagrangian up to $`O(v^2)`$ for $`N`$-maximally charged dilaton black hole system

$`L`$ $`=`$ $`{\displaystyle \underset{i=1}{\overset{n}{}}}m_i+{\displaystyle \underset{i=1}{\overset{n}{}}}{\displaystyle \frac{1}{2}}m_i(𝐯_i)^2`$ (11)
$`+{\displaystyle \frac{(N1)(Na^2)}{16\pi (N2+a^2)^2}}{\displaystyle d^Nx(F(𝐱))^{2(1a^2)/(N2+a^2)}\underset{i,j}{\overset{n}{}}\frac{(𝐧_i𝐧_j)|𝐯_i𝐯_j|^2\mu _i\mu _j}{2|𝐫_i|^{N1}|𝐫_j|^{N1}}},`$
where $`𝐫_i=𝐱𝐱_i`$ and $`𝐧_i=𝐫_i/|𝐫_i|`$. $`F(𝐱)`$ is defined by (4). In general, a naive integration in equation (11) diverges. Therefore, we regularize that divergent terms proportional to $`d^Nx\delta ^n(x)/|x|^p(p>0)`$ which appear when the integrand is expanded must be regularized . We set them to zero. The prescription is equivalent to carrying out the following replacement in equation (11)

$`(F(x))^{\frac{2(1a^2)}{(N2+a^2)}}1+\left[1+{\displaystyle \frac{8\pi (N2+a^2)}{A_{N1}(N2)(N1)}}{\displaystyle \frac{m_a}{|𝐫_a|^{N2}}}\right]^{\frac{2(1a^2)}{N2+a^2}}`$
$`+\left[1+{\displaystyle \frac{8\pi (N2+a^2)}{A_{N1}(N2)(N1)}}{\displaystyle \frac{m_b}{|𝐫_b|^{N2}}}\right]^{\frac{2(1a^2)}{N2+a^2}}.`$ (12)
After regularization, the effective Lagrangian for two body system (consisting of black hole labeled with $`a`$ and $`b`$) can be rewritten. From this rewriting effective Lagrangian, we obtain the metric of the $`N+1`$ dimensional moduli space for two-body system as

$$g_{ab}=\gamma (r)\delta _{ab},$$
(13)
with

$`\gamma (r)=1`$ $``$ $`{\displaystyle \frac{M}{\mu }}{\displaystyle \frac{8\pi (Na^2)}{A_{N1}(N2)(N1)}}{\displaystyle \frac{M}{r_{ab}^{N2}}}`$ (14)
$`+`$ $`{\displaystyle \frac{M}{m_a}}\left(1+{\displaystyle \frac{8\pi (N2+a^2)}{A_{N1}(N2)(N1)}}{\displaystyle \frac{m_a}{r_{ab}^{N2}}}\right)^{(Na^2)/(N2+a^2)}`$
$`+`$ $`{\displaystyle \frac{M}{m_b}}\left(1+{\displaystyle \frac{8\pi (N2+a^2)}{A_{N1}(N2)(N1)}}{\displaystyle \frac{m_b}{r_{ab}^{N2}}}\right)^{(Na^2)/(N2+a^2)},`$
where $`M=m_a+m_b`$, $`\mu =m_am_b/M`$, and $`r_{ab}=|𝐱_a𝐱_b|`$.

## 3 Quantum mechanics in two-black hole moduli space

We consider the quantum mechanics on moduli space. The quantization of moduli parameters has been discussed in .

Let us introduce a wave function $`\mathrm{\Psi }`$ on the moduli space, which obeys the Schrödinger equation

$$i\mathrm{}\frac{d\mathrm{\Psi }}{dt}=\left(\frac{\mathrm{}^2}{2\mu }^2+\mathrm{}^2\xi R_{(MS)}\right)\mathrm{\Psi },$$
(15)
where $`^2`$ is the covariant Laplacian constructed from the moduli space metric and $`R_{(MS)}`$ is the scalar curvature of the moduli space. We assume $`\xi =0`$ in this paper though this term may be present in most general case.
To simplify, we fix the case of the (3+1) dimensional case. The partial wave in a stationary state is

$$\mathrm{\Psi }=\psi _{ql}(r)Y_{lm}(\theta ,\varphi )\mathrm{exp}(iEt/\mathrm{}),$$
(16)
where $`Y_{lm}(\theta ,\varphi )`$ is the spherical harmonic function and $`E=\mathrm{}^2q^2/(2\mu )`$. We redefine the variables as

$`R`$ $`=`$ $`{\displaystyle \sqrt{\gamma }𝑑r},`$ (17)
$`\psi `$ $`=`$ $`{\displaystyle \frac{\chi }{r\sqrt{\gamma }}}.`$ (18)
The Schrödinger equation (15) is rewritten as

$$\frac{d^2\chi }{dR^2}+(q^2V)\chi =0,$$
(19)
where the potential $`V`$ is

$$V=\frac{\gamma (r\gamma ^{})^{}r(\gamma ^{})^2}{2r\gamma ^3}+\frac{l(l+1)}{r^2\gamma }.$$
(20)
Here  stands for $`\frac{d}{dr}`$.
To study the scattering process, we find the general view of potential on the moduli space. The potential for dilaton coupling $`a^2=0`$, $`1/3`$ and $`1`$ are plotted in Fig.1.

In the case of $`a^2=0`$ for all values of $`l`$ and $`a^2=1/3`$ for $`l=0`$, the potential has a maximal value. Then the incoming particle into two black hole moduli space are scattered or coalesced in the scattering process. On the other hand, in the case of $`a^2=1/3`$ for all values of $`l1`$ and $`a^2=1`$ all values of $`l`$, the incoming particle are always scattered away. We notice that the value of $`a^2=1/3`$ is the critical point of moduli structures . For simplicity, we will study the scattering-away processes in the case of $`a^2=1/3`$ and $`a^2=1`$.

## 4 Scattering on two-black hole moduli space

We consider the case of $`a^2=1/3`$ and $`a^2=1`$. As the first study of quantum effects, we discuss the scattering process in the WKB approximation. For the WKB approximation, the phase shift in the scattering process are

$$\delta _l=qR_0+_{R_0}^{\mathrm{}}𝑑R\left(\sqrt{q^2V(R)\frac{1}{4R^2}}q\right)+\frac{2l+1}{4}\pi ,$$
(21)
where $`R_0`$ is the solution of $`V(R_0)=q^2\frac{1}{4R_0^2}`$.
The partial cross section are

$$\sigma _l=\frac{4\pi }{q^2}(2l+1)\mathrm{sin}^2\delta _l.$$
(22)
The deflection angle are

$$\mathrm{\Theta }=\pi +_{R_0}^{\mathrm{}}𝑑R\frac{}{l}\sqrt{q^2V(R)\frac{1}{4R^2}}.$$
(23)
The phase shift $`\delta _l`$ and the deflection angle $`\mathrm{\Theta }`$ for $`q=0.1`$ and $`0.2`$ are plotted in Fig.2. The solid line in Fig.2 A2, B2 represents the classical deflection angle which are obtained from the impact parameter. For $`a^2=1`$, the effects for the different value of $`q`$ are not found, and the deflection angle are almost correspond to the classical deflection angle. Although, for $`a^2=1/3`$, we can find that the behaviours of the phase shift and the deflection angle depend on the incoming particle energy $`q`$. Then the deflection angle are different from the classical deflection angle at the small value of $`l/q`$. These effects are interpreted as the quantum effects arisen by the difference of the moduli space structures in the scattering process. Although we used the WKB semi-classical approximation in this analysis, the quantum effects can be obtained sufficiently.

## 5 Conclusion

In this work, we studied the quantum effects on the different moduli structure. We obtained the black hole moduli space metrics in any dimensions. In the (3+1) dimensions, we considered the quantum mechanics on the moduli space, and investigated the potential for $`a^2=0`$, $`1/3`$ and $`1`$. For the process in which the incoming particle are always scattered away, $`a^2=1/3`$ and $`1`$, we obtained the phase shift and the deflection angle in the WKB approximation. Then we revealed the quantum effects in scattering process for $`a^2=1/3`$. We considered that these effects due to the difference in the moduli space structure. So, the moduli space structure affect to the quantum scattering process.
In the further study, we will consider the behavior in the different dimensions and for the other values of dilaton coupling $`a^2`$. The moduli space structures are different in the other dimensions and at the other values of dilaton coupling $`a^2`$. In (4+1) dimensions for $`a^2=1`$, the moduli space structure are resemble the structure in (3+1) dimensions for $`a^2=1/3`$ . Then we expect that the similar quantum effects are also arisen. In this work, we consider the quantum effects in the WKB approximation. So, we will more precisely study in quantum mechanically analysis in further study.


text:
# Infrared Space Observatory Polarimetric Imaging of the Egg Nebula (RAFGL 2688)Based on observations with ISO (Kessler et al. 1996), an ESA project with instruments funded by ESA Member States (especially the PI countries: France, Germany, the Netherlands and the United Kingdom) with the participation of ISAS and NASA.

## 1 Introduction

Stellar mass loss during post-main sequence stellar evolution is, more often than not, an axisymmetric, rather than spherically symmetric, process. However, a coherent explanation for the transformation from quasi-spherical red giant to bipolar planetary nebula remains elusive. During pre-main sequence evolution, bipolar outflow is evidently an inevitable consequence of the presence of a circumstellar, protoplanetary disk. But — neglecting the possibility that some evolved stars may retain massive, long-lived, “fossil” protoplanetary disks (Kastner & Weintraub 1995) — there is no obvious natural “collimating agent” during post-main sequence stellar evolution. Hence the shaping of bipolar planetary nebulae (PNs), the endpoints of stellar evolution for intermediate-mass (1-8 $`M_{}`$) stars, is a topic of considerable contemporary interest in astronomy (see reviews in Kastner, Soker, & Rappaport 2000a).

The object RAFGL 2688 (The Egg Nebula) has long been regarded as prototypical of the class of evolved bipolar nebulae. The nebula was first identified as a bright far-IR source associated with axisymmetric optical nebulosity by Ney et al. (1975). Based on its radio molecular emission line profiles, Zuckerman et al. (1976) proposed that the object was representative of the rapid transition from red giant to planetary nebula. Immediately following its discovery, RAFGL 2688 was established as a source of highly polarized optical and near-IR emission (e.g., Jones & Dyck 1978). The large net polarization in these early, large aperture measurements was thought to originate with the scattering of starlight by circumstellar dust, where the dust grain population is distributed in axisymmetric fashion about the central star. The latitudinal dependence of the dust density results in large opacities along the equatorial plane but relatively small opacities toward higher latitudes. Assuming the polar axis is nearly in the plane of the sky, the star itself would then be obscured along our line of sight, whereas starlight can still escape along polar (but not equatorial) directions and be scattered toward Earth.

This qualitative model of the structure of RAFGL 2688 has since been confirmed, first by Monte Carlo scattering models (Yusef-Zadeh et al. 1984) and, most recently, by high spatial resolution Hubble Space Telescope near-infrared polarimetric imaging observations (Sahai et al. 1998a; Fig. 1). The HST imaging demonstrates conclusively that the nebula displays the centrosymmetric patterns and large degree of linear polarization expected for an axisymmetric distribution of dust grains illuminated by a central source. Furthermore, the HST observations have established the position of the central illuminator, relative to the extended nebulosity, to within $`100`$ AU (Weintraub et al. 2000). In addition to the unseen illuminator, RAFGL 2688 contains a prominent, central source of direct (unpolarized) emission that is only detectable at wavelengths $`>1.5\mu `$m (Sahai et al. 1998a; Weintraub et al. 2000). The relationship of this source to the illuminator is uncertain, although the simplest explanation would be that the pair constitutes a widely separated binary.

Although HST has well established the optical and near-infrared morphology of RAFGL 2688 (Sahai et al. 1998a,b), a detailed understanding of the nebula remains elusive. The kinematics of its molecular emission are particularly puzzling (see Cox et al. 2000, Kastner et al. 2001, and references therein). This emission seems to trace multipolar, radially directed outflows, some with no counterparts at optical or infrared wavelengths; alternatively, the emission may be interpreted in terms of a combination of radial and azimuthal velocity fields. Deepening the mystery of the structure of RAFGL 2688 are recent mid-infrared and radio imaging observations of dust emission that appear to trace highly collimated (and perhaps precessing) outflows (Morris & Sahai 2000; Jura et al. 2000).

Polarimetric imaging observations remain central to an understanding of the structure and evolution of RAFGL 2688. While the nebula is well characterized polarimetrically in the optical and near-infrared, little or no spatially resolved polarimetric imaging data exists for wavelengths $`>2\mu `$m. Here, we report on infrared polarimetric imaging of RAFGL 2688 at 4.5 $`\mu `$m with the Infrared Space Observatory (ISO, Kessler et al. 1996). In conjunction with previous HST near-infrared polarimetric imaging and ground-based direct mid-infrared imaging, the ISO polarimetric imaging yields insights into the relative contributions of direct and unpolarized starlight and of the grain size distribution in RAFGL 2688.

## 2 Observations

Data presented here were obtained with the ISO direct imaging instrument, ISOCAM (Cesarsky et al., 1996), using its long wavelength (LW) camera and passband filter LW1. The detector format was 32$`\times `$32 pixels. The response of LW1 is centered at 4.5 $`\mu `$m and is relatively flat, extending from 4 to 5 $`\mu `$m.

ISOCAM was equipped with polarizers oriented at 0, 120, and 240 degrees. Polarimetric imaging was set up according to the standard observing template CAM05 (Siebenmorgen 1996). The corresponding ISO observation numbers, TDT, are: 76800602, 76800703, 76800804, 76800905 76801006,76801107, 76801208, 76801309, 76801410 and 76801511. We used a $`4\times 4`$ raster with a raster step size of $`6^{\prime \prime }`$, the $`1.5^{\prime \prime }`$ lens, a readout time for each exposure of 2.1 s, and a gain of one. Initially, we took 155 exposures on the source to stabilize the detector and then a raster through the free “hole” of the entrance wheel. At each raster position 26 exposures were read out. This procedure was carried out for each of the three polarizers. The polarizer rasters were repeated in three observing cycles.

## 3 Results

### 3.1 Processed ISOCAM LW1 images

The data were reduced with the ISOCAM interactive analysis system (CIA version 4.0, Ott et al. 1996). Only basic reduction steps, such as dark current subtraction, removal of cosmic ray hits and transient correction, are applied to the data. These steps are described together with calibration uncertainties in Blommaert et al. (2001). The coadded images at each raster position were projected on the sky to derive the final mosaics. The mosaics have a $`66^{\prime \prime }\times 68^{\prime \prime }`$ total field of view and a pixel scale of 1.5′′. The polarized signal is found to be consistent between the individual observing cycles. The average of all cycles gives the final mosaic image for each polarizer. The output of this image processing procedure therefore consisted of four images, each of total integration time 1124 s: three polarized images (one for each of the three polarization filters) and one “unpolarized” or total intensity image (obtained through the “hole” position).

As part of the ISO calibration program, the zodiacal light and a set of (unpolarised) photometric standard stars was observed by repeating CAM05 observing templates. From those calibration measurements an instrumental polarization of $`<`$1.5% is derived (Siebenmorgen 1999). This instrumental polarization is small enough to have negligible effect on the results described in this paper.

The total intensity image of RAFGL 2688 is displayed in Fig. 2, alongside a 3420 s ISOCAM LW1 image of the infrared-bright standard star HIC 85317 (spectral type A4V, V$`=6.5`$mag, TDT=35600401). The comparison emphasizes the rather complex, asymmetric point spread function (PSF) of ISOCAM at 4.5 $`\mu `$m. The six-pointed structure of the PSF is due to the support structure for the secondary mirror of the telescope. Although this point spread function dominates the structure of both the RAFGL 2688 and standard star images, it is also apparent that RAFGL 2688 is extended relative to the standard. This is seen most clearly in the comparisons of nonlinear greyscale images (Fig. 2, bottom panels) and plots of intensity vs. radius for RAFGL 2688 and the standard star (Fig. 3), the latter of which reveals twice as much or more light in the RAFGL 2688 images as in the PSF at distances $`10^{\prime \prime }`$ from the bright central PSF core.

### 3.2 ISOCAM 4.5 $`\mu `$m polarimetric map

The polarimetric image set obtained with the three polarizing filters in the LW camera is sufficient to uniquely determine the Stokes parameters for each image pixel $`(x,y)`$. Specifically, we find $`I(x,y),Q(x,y),U(x,y)`$ from the image set $`P_0(x,y),P_{120}(x,y),P_{240}(x,y)`$ as follows:

$`I(x,y)={\displaystyle \frac{2}{3}}[P_0(x,y)+P_{120}(x,y)+P_{240}(x,y)]`$ (1)
$`Q(x,y)={\displaystyle \frac{2}{3}}[2P_0(x,y)P_{120}(x,y)P_{240}(x,y)]`$ (2)
$`U(x,y)={\displaystyle \frac{2}{\sqrt{3}}}P_{120}(x,y)+{\displaystyle \frac{2}{\sqrt{3}}}P_{240}(x,y)`$ (3)
The degree $`p`$ and position angle $`\theta `$ of linear polarization may then be derived from the Stokes parameters in the usual manner, i.e.,

$$p(x,y)=\frac{(Q^2(x,y)+U^2(x,y))^{1/2}}{I(x,y)}$$
(4)
and

$$\theta (x,y)=\frac{1}{2}\mathrm{tan}^1\frac{U(x,y)}{Q(x,y)}.$$
(5)
Polarization maps so obtained from ISOCAM polarimetric imaging have been presented for the Crab nebula (Tuffs et al. 1999), the protostellar system HH 108MMS (Siebenmorgen & Krügel 2001), the starburst NGC 1808 galaxy (Siebenmorgen et al., 2001), and a sample of ultraluminous galaxies (Siebenmorgen & Efstathiou 2001).

The ISO polarization map obtained for RAFGL 2688 (Fig. 4) shows hints of centrosymmetry, particularly in its northern regions. The center of symmetry of this portion of the polarimetric map appears to lie near the bright central region of the image. This spatial polarization structure suggests that much of the northern region of the ISO 4.5 $`\mu `$m image of RAFGL 2688 is, in fact, scattered light, and that the source of illumination lies in the vicinity of the image intensity peak. However, the structure of the ISO polarimetric map is, overall, much less well organized than that of the HST/NICMOS polarization map (Fig. 1). In particular, the relative contributions from nebulosity and PSF are not clear from Fig. 1; nor is it obvious whether any extended 4.5 $`\mu `$m emission is dominated by scattering (as is undoubtedly the case at 2 $`\mu `$m) or by direct, thermal emission. Evidently, either the ISO polarimetric map suffers from systematic distortions, or there are significant differences in the spatial distribution and degree of polarization in RAFGL 2688 between 2 $`\mu `$m and 4.5 $`\mu `$m.

## 4 Reconstruction of Direct and Polarimetric Images

To investigate the extent to which the prominent signature of the PSF in the individual polarimetric images affects the ISO 4.5 $`\mu `$m polarimetric map, we attempted to remove the PSF signature from these images and that of the total intensity image. To this end we have employed two popular deconvolution methods, maximum likelihood (hereafter ML; Richardson 1972; Lucy 1974) and maximum entropy (e.g., Hollis, Dorband, & Yusef-Zadeh 1992; hereafter ME), using programs that are part of the standard IDL Astronomy Library111http://idlastro.gsfc.nasa.gov/.

Both ML and ME techniques require a deconvolution “kernel” consisting of the best known representation of the instrumental PSF. One then makes an initial estimate of the “true” source intensity distribution, i.e., the source intensity distribution that would be detected in the absence of instrumental broadening. This image is convolved with the PSF and is then compared with the observed image intensity distribution. The resulting difference image (or “residual image”) is stored for reference. The “true” image intensity distribution is then modified accordingly (following either the ML or ME algorithm prescription), and the process is repeated. The change in residual image statistics (e.g., residual image maximum and standard deviation) over successive iterations of either the ML or ME schemes can be used to assess the progress of image deconvolution.

### 4.1 Reconstructed ISOCAM 4.5 $`\mu `$m total intensity image

The results of iterative deconvolutions of the ISOCAM LW1 total intensity image using ML and ME methods are displayed in Figs. 5 and 6. For these deconvolutions, we adopted as the instrumental PSF the ISOCAM LW1 image of HIC 85317. It is apparent from both figures that the ML and ME methods produce very similar results and, furthermore, that the two methods asymptotically approach their “best guesses” for the “true” image at similar rates. Specifically, after $`15`$ iterations, the net changes in both the maximum and the standard deviation of the residual image (from their original values) approach $`90`$% of the differences between their original and asymptotic values (Fig. 5). We find that, for either the ML or ME methods, iterating more than about 20 times introduces image artifacts, such as negative values and spurious off-source features.

Notably, both the ML and ME deconvolution methods produce the same fundamental result for the reconstructed total intensity image (Fig. 6). This reconstructed image consists of a compact, bright central source, representing $``$90% of the total source intensity, surrounded by fainter nebulosity. Most of the nebulosity lies to the north-northeast of the bright, compact source. This source structure is consistent with the morphology of the polarimetric map (§3.2), i.e., the bright, central source appears to illuminate a region of reflection nebulosity to its north-northeast.

### 4.2 Reconstructed ISOCAM polarization map

We applied the ML method to the individual LW1 polarimetric images, to obtain a deconvolved polarization image set $`P_0^{}(x,y),P_{120}^{}(x,y),P_{240}^{}(x,y)`$. For these deconvolutions, we again adopted the ISOCAM LW1 image of HIC 85317 as representative of the instrumental PSF, under the assumption that the polarizing filters produce negligible additional point source distortions. To obtain each polarimetric image, we performed 16 iterations of the ML algorithm. The resulting deconvolved LW1 polarimetric images of RAFGL 2688 (not shown) are qualitatively similar to the deconvolved total intensity image described in the preceding section, i.e., a compact central source with nebulosity extending to the north-northeast.

The polarization map constructed from the images $`P_0^{}(x,y),P_{120}^{}(x,y),P_{240}^{}(x,y)`$ (hereafter referred to as the reconstructed polarization map) is presented in Fig. 7. The centrosymmetry evident in the northern regions of the “raw” polarimetric map (Fig. 4) is more pronounced in this reconstructed polarization map. Like the “raw” map, the reconstructed map also displays a minimum of polarization near the intensity peak in the total intensity image, and a somewhat disorganized polarization pattern to the south of the intensity peak.

## 5 Discussion

### 5.1 Can the results of the deconvolutions be trusted?

It is well known that the results of either ML or ME deconvolution are not necessarily unique. Indeed, depending on the quality of both data and PSF model, the results of either technique can be spurious. Furthermore, we are not aware of a rigorous proof that demonstrates that deconvolution of polarimetric images preserves polarization information. Nonetheless, there are several reasons to conclude that the deconvolution results just described are valid.

First, we find strong consistency between the results of the ML and ME methods. The two methods converge to essentially the same end results (Figs. 5,6). Second, the reconstructed ISOCAM 4.5 $`\mu `$m image and polarization map of RAFGL 2688 are qualitatively consistent with the results of previous near-infrared imaging polarimetry and direct mid-infrared imaging. Specifically, in the reconstructed ISOCAM image, RAFGL 2688 is extended along a position angle $`15^{}`$, a result also obtained from high-resolution 2 $`\mu `$m and 9 $`\mu `$m imaging (Weintraub et al. 2000; Morris & Sahai 2000) but that can be only marginally inferred from the raw ISO images (e.g., Fig. 4).

Third, the reconstructed polarization map (Fig. 7) represents an improvement over the polarization map constructed from “raw” images (Fig. 4). In particular, note the stronger centrosymmetry and larger values of polarization in the northern regions of the map in Fig. 7. Both improvements would be expected from a successful deconvolution procedure; such a procedure should correct for instrumental stray light at low spatial frequencies, thereby removing “contaminating” flux that would otherwise suppress the degree of polarization, especially in regions of low source surface brightness.

We conclude that the reconstructed ISOCAM 4.5 $`\mu `$m total intensity image and polarization map of RAFGL 2688 are trustworthy.

### 5.2 The nature of the 4.5 $`\mu `$m emission from RAFGL 2688

The ISOCAM polarimetric map of RAFGL 2688 reveals that the extended nebulosity to the north-northeast is highly polarized, indicative of reflection nebulosity. The centrosymmetry of this polarization pattern about the 4.5 $`\mu `$m intensity peak indicates that the intensity peak is the likely source of illumination of the nebulosity. This source is very likely a combination of direct photospheric emission from the central F-type supergiant (e.g., Cohen & Kuhi 1977) and photospheric emission reprocessed by dust grains at characteristic temperatures of $`\stackrel{>}{}140`$ K (Persi et al. 1999).

The apparent “emergence” of the central source at longer infrared wavelengths is seen in Fig. 8. The source is only seen via reflection at 2 $`\mu `$m (lefthand panel; see also Weintraub et al. 2000) whereas, evidently, it makes a substantial contribution to the nebular morphology (in the form of a strong point-like source and diluted polarization) near the center of the nebula at 4.5 and 8.8 $`\mu `$m (central and righthand panels). Similar behavior is observed in infrared reflection nebulosity surrounding young stellar objects, and is usually interpreted as a decrease in scattering opacity with increasing wavelength (e.g., Fischer 1995). We caution, however, that the registration of the three images in Fig. 8 is not certain to better than $`1^{\prime \prime }`$, the ISOCAM polarization imaging suffers from relatively low spatial resolution, and there is no polarization information in the 8.8 $`\mu `$m image. Furthermore, there are multiple sources of 2 $`\mu `$m photons in RAFGL 2688 (Weintraub et al. 2000). Hence it may be somewhat premature to conclude that the same source that illuminates the nebula at 2 $`\mu `$m also illuminates reflection nebulosity at 4.5 and 8.8 $`\mu `$m.

While it appears we have detected the central illuminating source at 4.5 $`\mu `$m, this source evidently is highly obscured along our line of sight. This conclusion follows from the ratio of nebular surface brightness to point source intensity; specifically, at an angular displacement of $`6^{\prime \prime }`$ along the polar axis of the nebula from the central point source, we find that the surface brightness of scattered light is $``$1% of the integrated intensity of the central source. For a reflection nebula whose dust density is inversely proportional to the square of the angle from the source of illumination ($`\theta `$), the scattering optical depth ($`\tau _s`$) is related to the ratio of scattered to incident starlight, $`I_s(\theta )/I_0`$, via

$$\tau _s(\theta )=8\pi \frac{I_s(\theta )}{I_0}\frac{\theta ^2}{\mathrm{\Delta }\varphi }$$
(6)
where $`\mathrm{\Delta }\varphi `$ is the solid angle subtended by a pixel (Jura & Jacoby 1976; Kastner 1990). Assuming the integrated intensity of the central source as measured in our images accurately represents the available illuminating flux, $`I_0`$, then the measured ratio $`I_s(6^{\prime \prime })/I_00.01`$ would suggest a scattering optical depth $`\tau _s(6^{\prime \prime })4`$. This result is seemingly at odds with our detection of very high levels of polarization ($`60`$%) in this region, however; such a high degree of polarization requires that single scattering is dominant over multiple scattering. We conclude that, along our line of sight, the central object is highly attenuated at 4.5 $`\mu `$m, such that the measured value of $`I_s(6^{\prime \prime })/I_0`$ is an overestimate. A conservative lower limit to the degree of 4.5 $`\mu `$m extinction toward the central source is provided by the requirement that $`\tau _s<1`$, suggesting $`A_{4.5\mu m}>1.5`$.

Fig. 9 shows a comparison of percent polarization vs. radius from the central star, for the 2 $`\mu `$m (HST/NICMOS) and 4.5 $`\mu `$m (ISOCAM) polarimetric imaging observations. The figure demonstrates that percent polarization is more or less independent of wavelength within the nebula, between $`2`$ $`\mu `$m and 5 $`\mu `$m, at radial offsets well displaced from the star. This behavior is similar to that observed in the evolved bipolar reflection nebula OH 231.8+4.2 (Shure et al. 1995), although the OH 231.8+4.2 observations only extended to 3.6 $`\mu `$m. From the large values of polarization in RAFGL 2688 at both 2 $`\mu `$m and 4.5 $`\mu `$m, we may safely conclude that scattering dominates the nebular emission at these wavelengths; i.e., there is little or no contribution from thermal emission from the dust grains in RAFGL 2688, even at wavelengths as long as $`5`$ $`\mu `$m. As no polarization data have been obtained for this nebula at longer wavelengths we are, as yet, unable to ascertain the nature of the lobe emission seen at 8.8 $`\mu `$m (Fig. 8).

### 5.3 Dust grain size distribution

Most detailed modeling of reflection nebulae at $`5`$ $`\mu `$m performed thus far has been geared toward interpeting large-aperture polarization measurements (e.g., Heckert & Smith 1988; Pendleton et al. 1990; Fischer 1995; Shure et al. 1995) rather than spatially resolved polarization imaging. Nevertheless, based on the wavelength independence of polarization of RAFGL 2688 and the typical polarization values of $`60`$% at large offsets from the central star, we conclude that the bulk of the dust grains in the RAFGL 2688 reflection nebula are significantly smaller than the wavelengths employed in the HST/NICMOS and ISOCAM imaging studies. That is, the grains have typical radii $`a<1`$ $`\mu `$m. This result, which is similar to that obtained by Shure et al. (1995) for OH 231.8+4.2, follows from the straightforward analytical treatment of light scattering by small particles (van de Hulst 1957) as well as from the more detailed, but limited, modeling of dust grain populations in infrared reflection nebulae performed to date (Pendleton et al. 1990; Fischer 1995).

## 6 Summary

We have obtained 4.5 $`\mu `$m polarimetric images of the evolved bipolar reflection nebula RAFGL 2688 with the CAM mid-infrared camera aboard the Infrared Space Observatory. The observations are strongly affected by the complex CAM point spread function. The source structure in reconstructed (deconvolved) images is that of a bright central object with faint nebulosity extending to the north-northeast. The position angle of extended nebulosity is the same as that seen in ground-based near-infrared images obtained at other wavelengths. A polarization map constructed from deconvolved, polarized images displays a large degree of polarization and a centrosymmetric pattern of position angle within the faint, extended nebulosity, while the bright central source displays much lower levels of polarization and disorganized patterns of polarization position angle.

The pattern and degree of polarization within the extended nebulosity, as revealed in the reconstructed image and polarization map, confirms that scattering is the dominant source of nebular light. Comparison with previous space-based (HST/NICMOS) polarimetric imaging indicates that the degree of polarization within the extended RAFGL 2688 nebula is independent of wavelength, suggesting that small ($`<1`$ $`\mu `$m) grains do most of the scattering of starlight.

Further, detailed modeling is required to provide additional insight into the circumstellar structure and grain properties of evolved, bipolar nebulae such as RAFGL 2688. It would also be of great value to obtain mid-infrared imaging polarimetry of RAFGL 2688 at higher resolution, to ascertain the precise location of the illuminating source at $`5`$ $`\mu `$m. as well as at longer wavelengths, to determine whether scattering is still important or even dominant beyond 5 $`\mu `$m.

The authors thank Mark Morris and Raghvendra Sahai for the use of their 8.8 $`\mu `$m Keck image of RAFGL 2688. CIA is a joint development by the ESA Astrophysics Division and the ISOCAM Consortium. The ISOCAM Consortium is led by the ISOCAM PI, C. Cesarsky. ISO research by JHK and JL was supported by JPL grants 1230950 and 1211136 to RIT.