Towards slime mould electrical logic gates with optical coupling

Abstract: Physarum polycephalum is a macroscopic single celled plasmodial slime mould. We employ plasmodial phototactic responses to construct laboratory prototypes of NOT and NAND logical gates with electrical inputs/outputs and optical coupling; the slime mould plays dual roles of computing device and electrical conductor. Slime mould logical gates are fault tolerant and resettable. The results presented here advance our understanding of how biological computing substrates may be manipulated to implement logical operations and demonstrate the feasibility of integrating living substrates into silicon hardware.

\section{Introduction}

	Biological computing devices are currently a source of intense interest in both research and industry: massively-parallel, amorphous computing substrates such as slime moulds or mammalian brains are far superior to conventional solid-state computing architectures in a trade-off between energy consuming, self-growth and  processing ability, and our rapid approach to the physical limitations of inorganic material, the comparative rarity of earth metals used in electronic components and the toxic by-products produced in solid-state computer manufacture are all prime reasons to investigate unconventional computing devices. Furthermore, if sufficiently developed, such devices may allow us insight into other novel bio-computing technologies, e.g. neurally-integrated medical or augmentative prosthetic implants. \emph{P. polycephalum} is an ideal unconventional computing substrate as it is easy to culture, tolerant to abuse and has no ethical issues surrounding its use~\cite{Keller_Everhart_2010}.

The computing potential of the myxomycete slime mould \emph{Physarum polycephalum's} plasmodium (vegetative life-cycle stage, pl. plasmodia) has become a hot topic in computer science, biology and biophysics. Behaving as an excitable or a reaction-diffusion system encapsulated in an elastic membrane~\cite{Adamatzky_2007}, the slime mould may approximate shortest paths between two points~\cite{Nakagaki_2001} and solve maze problems \cite{Adamatzky_2012a,Nakagaki_2001}. Slime mould can be regarded as a proliferating massively parallel processor (a 'Physarum machine' \cite{Adamatzky_2010a}), capable of simultaneous processing of inputs (light, chemical gradients, temperature etc.) \cite{Adamatzky_2007}, concurrent decision making and distributed actuation. In experimental situations, slime mould based computing devices employ chemotaxis --- positive, towards chemoattractants, or negative, away from repellents --- or photoavoidance; the plasmoidum is motile via the contraction of muscle-like proteins invoking rhythmic shuttle streaming of its hydrodynamic core \cite{Adamatzky_2012b,Ishikawa_etal_1991}.

Previously we have demonstrated how \emph{P. polycephalum} may be hybridised with nano-scale metallic particles in order to manually alter the electrical properties of the plasmodium, e.g. facilitate the conduction of electricity and generate permanent, solid-state wires which persist following the plasmodium's migration and/or death \cite{Mayne_etal_2013,Mayne_Adamatzky_2013}.

In this study, we explore the use of light as a programmable input into a Physarum machine, and use the results to implement a basic yet functional device. Significant progress has been made programming a Physarum machine with spatially distributed nutrient sources~\cite{Nakagaki_2001,Adamatzky_2012b,Adamatzky_2010a}, but these are undynamic and problematic to `reprogram' (i.e. alter the input mid-experiment). Indeed, viable logic gates have been designed and produced by a number of authors using different variations of chemical attractant and repellent models~\cite{Tsuda_Aono_Gunji_2004,Adamatzky_2010b,Jones_Adamatzky_2010}, but all suffer from the aforementioned restrictions. These limitations may potentially be overcome by utilising light as an input, or rather a coupling medium between the sub-circuits of logical gates.

Two logic gates were designed: a NOT gate (more commonly known as an inverter), chosen for being the simplest logic operation, and a NAND gate, for being a universal gate, i.e. any logical operation can be carried out by multiple NAND gates.

\section{Selecting inputs: Phototaxis experiments}
\label{phototaxis}

\begin{figure}[!tbp]
\centering
\subfigure[]{\includegraphics[width=0.6\textwidth]{figs/Fig_1a}}
\subfigure[]{\includegraphics[width=0.9\textwidth]{figs/Fig_1bii_arrowed}}
\caption{Phototaxis experiments.
(a)~A scheme of experimental set-up for phototaxis experiments. Gray squares represent agar cubes: the slime mould is inoculated on the central cube (pentagon) and is left to proliferate. Black rectangles represent card barriers with a 0.5~mm gap underneath, which the plasmodium is able to migrate underneath to access the peripheral agar cubes, which are loaded with chemoattractants and illuminated by 2$\times$ LEDs (circles).
(b)~Photograph of a completed phototaxis experiment: note how the plasmodium has migrated towards the left 'blue' side (arrowed).}
\label{fig_taxis}
\end{figure}

Plasmodia were inoculated on to 2.5~cm\textsuperscript{2} squares of 2\% non-nutrient agar in the centre of 9~mm Petri dishes. Similar squares of agar loaded with chemoattractants (oat flakes) were situated at the 3 and 9 o'clock poles of the dishes  and were each illuminated by 2 lid-mounted LEDs (Fig.~\ref{fig_taxis}a). 4 colours of LED were used --- blue (466~nm), green (568~nm), yellow (585~nm) and red (626~nm) --- and each pole had a different colour assigned to it in each experiment. The peripheral regions were separated from the centre by a strip of thick card glued to the dish lid, with a 0.5~mm gap underneath: this prevented most of the light reaching the inoculation point, but allowed the plasmodium to migrate through the gap underneath.

Test plasmodia were therefore given three basic choices in each experiment: migrate towards colour A, migrate towards colour B or do neither (and hence stay in the centre/sclerotinize/attempt to colonise the card/escape the dish etc.). The choices each plasmodium made were compiled into an order of preference for each colour: for each experiment, the colour which was avoided was given 1 'phobia point'; the colour with the most points after all experiments had been completed was determined to be the strongest repellent.

	A completed phototaxis experiment is shown in Fig.~\ref{fig_taxis}b. This pattern of avoidance was identical for every repetition. On no occasion did the plasmodium chose to stay on the central agar square and sclerotinize, although on two occasions it migrated up the card barrier and attempted to consume the glue securing it to the lid of the Petri dish. Phototaxis experiments also took considerably more time to complete than anticipated, with each lasting 5--6 days (whereas similar experiments with only chemoattractants/repellents typically last less time): the plasmodium would typically spend several days exploring the central agar piece and occasionally extend pseudopodia several millimetres along the plastic base of the Petri dish in random directions. This would seem to imply that the plasmodium's actions were a 'last resort' --- a protistic equivalent of reluctance.

\emph{P. polycephalum's} preference to each colour of LED is summarised below, arranged in order of most avoided $\rightarrow$ least avoided:
\begin{center}Green $\rightarrow$ Red $\rightarrow$ Yellow $\rightarrow$ Blue \\\end{center}

Thus green LEDs were used in further experiments.

\section{NOT Gate}

\begin{figure}[!tbp]
\centering
\subfigure[]{\includegraphics[width=0.9\textwidth]{figs/Fig_2a}}
\subfigure[]{\includegraphics[width=0.49\textwidth]{figs/Fig_2b}}
\subfigure[]{\includegraphics[width=0.49\textwidth]{figs/Fig_2c}}
\caption{Physarum NOT logic gate.
(a)~Schematic diagram of a PNOT logic gate. A plasmodium (pentagon) is placed on agar (grey squares) overlying electrode X and chemoattractants are placed on Y. If the input [a] to the LEDs (circles)  is at logic 0, the LED does not switch on and a plasmodial tube is formed between electrodes X and Y; \emph{vice versa} also occurs. The electrodes are wired into an independently-powered output circuit [o].
(bc)~Photographs of a functional PNOT gate.
(b)~Experimental setup with input at logic 1; the plasmodium did not migrate and the output was left at logic 0.
(c)~With input at logic 0, the plasmodium has migrated across the gap and created a `wire' which was electrically conductive, leaving the output at logic 1.
}
\label{fig_pnot}
\end{figure}

	The design for the Physarum-NOT gate (PNOT) is shown in Fig.~\ref{fig_pnot}a. Two aluminium tape electrodes (X and Y) were stuck to the base of a 9~mm plastic Petri dish, separated by a central 10~mm gap. 2~ml agar blobs were placed on top of each electrode, one of which was loaded with plasmodium (X) and the other a chemoattractanct (an oat flake) (Y). The gate's 'input' consisted of a 9~V power supply (Lascar Electronics, Salisbury, UK) connected to a lid-mounted LED overlying electrode Y. At input equals 1 (logic levels of 0 and 1 are represented here as 0~V and 9~V, respectively), the LED was active and prevented the plasmodium from migrating from X to Y, preventing the completion of the circuit. If the input equals 0, the LED was not active and therefore the plasmodium was free to migrate to electrode Y, forming a continuous tube of plasmodium between both electrodes.

	Both electrodes were connected to a separate 'output' circuit with a constant potential difference applied across it. The plasmodium, when spanning two electrodes, acted as a physical conductor of electricity and hence completed the output circuit. As such, the device's functionality is the same as a conventional NOT gate, i.e. if the input is live, the output is not and \emph{vice versa}.

	Experimental photographs of a PNOT gate are shown in Fig.~\ref{fig_pnot}. The time it took for the plasmodium to complete its logic function was extremely variable, ranging from 1--4 days; approximately 25\% entirely failed to propagate. With reference to the results of the phototaxis experiments, this was reasoned to be due to the plasmodium's tendency to stay still for as long as possible before migrating, during which time the agar blobs may have desiccated to the point that they could no longer support its demand for moisture.

\section{NAND Gate}

\begin{figure}[!tbp]
\centering
\subfigure[]{\includegraphics[width=0.9\textwidth]{figs/Fig_3aii}}
\subfigure[]{\includegraphics[width=0.49\textwidth]{figs/Fig_3b}}
\subfigure[]{\includegraphics[width=0.49\textwidth]{figs/Fig_3c}}
\subfigure[]{\includegraphics[width=0.49\textwidth]{figs/Fig_3d}}
\caption{Physarum NAND gate.
(a)~Schematic diagram of a PNAND logic gate. Slime mould (pentagon) is placed on agar (grey squares) overlying electrode X and attractants are placed on Y and Z. If inputs A or B are at logic 1, the corresponding LEDs (circles) turn on, preventing slime mould growth to the corresponding electrode and hence stop a circuit from forming. If both A and B are on, no circuit forms, resulting in no output and \emph{vice versa}. The electrodes are wired into an independently-powered output circuit [o]. A card barrier separates electrodes Y and Z to prevent light contamination.
(bcd)~Photographs of a functional PNAND gate.
(b)~T=6~h; the inoculation point is fully colonised and the left electrode is illuminated.
(c)~ T=12~h; the plasmodium has migrated to the left of the inoculation point.
(d)~T=24~h; the plasmodium has migrated to the unilluminated electrode.
}
\label{Fig_nand}
\end{figure}

Physarum-NAND (PNAND) gates were constructed based upon the same principles of the PNOT gate, but with a second electrode and LED added to the device (Fig.~\ref{Fig_nand}a). LEDs A and B were linked to separate power sources (i.e. inputs) so that their corresponding LEDs worked independently of each other. A plasmodium was once again loaded onto electrode X, but was given the added option of migrating towards electrode Z as well as Y. If both inputs equal 0, the plasmodium could migrate towards either Y or Z at random, forming a circuit between the initial and new electrode, causing the output  to equal 1. If input A is 1 but B is 0, the plasmodium at electrode X is still free to complete the circuit by migrating to electrode Z, and \emph{vice versa}. If both A and B equal 1, the plasmodium would not migrate towards either electrode, preventing a circuit from forming and hence leaving the output at 0. Electrodes X and Y were physically separated by a barrier similar to those described in Sec.~\ref{phototaxis} to prevent light contamination from opposing LEDs. Experimental photos of PNAND gate are shown in Fig.~\ref{Fig_nand}b--d.

\section{Reusability}

	Both PNOT and PNAND gates were tested for reusability. Once a logical operation had been carried out, the inputs were changed and the consequent responses of experimental plasmodia were recorded.

Both PNOT and PNAND gates are reusable, but only a finite number of times. Resetting gates --- i.e. allowing the plasmodium to propagate between two unilluminated points then illuminating one of them in expectation that the plasmodium would withdraw --- was successful in that it caused the plasmodium to withdraw from the target electrode in all instances. This process usually took 2-6 hours and involved the contents of the plasmodium flowing back down the tubes spanning both electrodes. Occasionally, the plasmodium split into two separate fragments, both of which sclerotinized. With only two electrodes present, the PNOT gate was therefore effectively reset as the conductive trail between both electrodes was disrupted. PNOT gates were not successfully re-reset, i.e. the slime mould would not migrate back to an area it had already colonised and withdrawn from.

PNAND gates were successfully reset and reprogrammed, as shown in Fig.~\ref{fig_reset}: once a logical operation had been completed, the inputs were altered and the slime mould migrated accordingly to complete the new operation. PNAND gates that had been reset had a higher failure rate than `fresh' gates, but their propagation delay was significantly shorter (typically $<$~24 hours).

\begin{figure}[!tbp]
\centering
\subfigure[]{\includegraphics[width=0.47\textwidth]{figs/Fig_4a}}
\subfigure[]{\includegraphics[width=0.51\textwidth]{figs/Fig_4b}}
\subfigure[]{\includegraphics[width=0.51\textwidth]{figs/Fig_4c}}
\caption{Resetting and reprogramming a PNAND gate.
(a)~T=36~h. A logical operation has been completed by the slime mould: the right LEDs are on and so it has migrated to the left electrode.
(b)~T=38~h. The inputs are switched; the right LEDs are now off and the left are on. The slime mould begins to withdraw towards the inoculation point and sends pseudopodia in multiple directions looking for a new place to migrate to.
(c)~T=42~h. The slime mould migrates towards the unilluminated right electrode and begins to withdraw its protoplasm from the left electrode and the tube between left and centre electrodes gradually ceases to be conductive.
}
\label{fig_reset}
\end{figure}

\section{Fault-tolerance}

The fault tolerance of Physarum gates was also investigated by altering physical variables in the test environment, including:
\begin{enumerate}
\item Increasing the luminosity of the LEDs by replacing them with brighter variants (each was increased by approximately 50~mcd).
\item Increasing the distance between electrodes to 20~mm
\item Increasing the potential difference between electrodes to 24~V
\end{enumerate}

None of these variables were `decreased', as the initial variables were originally designed as minimal operational parameters; dimmer light sources were anticipated to be ignored by the plasmodium, decreasing the electrode-electrode gap to below 10~mm could result in the plasmodium easily spanning all three (resulting in gate mis-operation) and applying a lower voltage through the electrodes would likely have resulted in an output so low as to be undetectable.

	Whilst adjusting these variables is of little practical use in the design of these specific prototype devices, they were chosen to (a) investigate whether fine-tuning any of these variables would have a beneficial effect on the operation of the device, and (b) demonstrate how a Physarum machine may adapt to unfavourable conditions in a way that traditional computing architectures cannot --- e.g. it is doubtful that any solid state architectures would function properly following exposure to prolonged overvoltage. By introducing extra elements of environmental instability, \emph{P. polycephalum} was given the opportunity to demonstrate how resilient it is, and therefore how robust a Physarum machine could potentially be.

Increasing the luminosity of the LEDs appeared to have no beneficial or detrimental effect on their operation; failure rates were not significantly different and neither was the operation time/propagation delay. Whether \emph{P. polycephalum} is generally phobic to any intensity of green light, or whether the increase in incandescence was too low to provoke any stronger responses was not ascertained.

Doubling the distance between the electrodes markedly increased the failure rate of both varieties of Physarum logic gate. This was likely to result from the combined effects of \emph{P. polycephalum's} reluctance to cross the arid plastic regions of Petri dishes and a lack of any detectable nutrient gradients due to the distant proximity of the supplied chemoattractants.

Experimental plasmodia appeared to be extremely tolerant of high voltages. Whilst the current that would have run through the plasmodium would have been comparatively low due to the resistance of the agar blobs and the innate resistance of the plasmodium its self, it was nevertheless tolerant to applied potential differences of up to and including 24~V. The morphology of plasmodia exposed to high voltages was slightly different to the wild type; rather than spanning the gap between electrodes with a single tube of plasmodium, plasmodia exposed to higher voltages tended to occupy a significant amount of space between electrodes with multiple interlinking tubules (Fig.`\ref{Fig_5}b).

\section{Conclusions}

The principle finding of this study was that it is indeed possible to `program' an unaltered \emph{P. polycephalum} plasmodium to carry out basic logic functions with optical inputs. Whilst the prototype devices presented here have significant limitations, they demonstrate a genuine implementation of unconventional computation in which a living computing substrate functions as both a sensing, decision making computing device and a physical conductor of electricity.

\emph{P. polycephalum's} photoavoidance is widely mentioned in literature, and indeed the plasmodium is best cultivated in the dark. Sauer~\cite{Sauer_1982}  reports that the plasmodium is phobic to blue light (which may also induce sporulation whilst the plasmodium is also starved of nutrients and/or moisture), but may sometimes migrate towards red light. Green light has been reported to inhibit sporulation, and seemingly random wavelengths of every colour (and some UV wavelengths) induce a range of metabolic changes \cite{Hato_etal_1976,Sauer_1982}. Whilst \emph{P. polycephalum} is not photosynthetic, it contains pigments similar to phytochromes which are thought to act as photoreceptors \cite{Kakiuchi_etal_2001}. A migrating plasmodium may be diverted by placing illuminated barriers in its path \cite{Adamatzky_2009}.

There were clear discrepancies between the results of the phototaxis experiments conducted for this investigation and those in published literature. There were, however, stark differences in experimental setup between the experiments described here and those in literature; LED diffusion, different wavelengths of light used, different intensities of light and different experimental conditions could all have altered plasmodial responses drastically. It should also be noted that if, for example, the plasmodium migrated towards blue light, this does not necessarily imply that \emph{P. polycephalum} is preferentially moving towards the blue light, but is more likely moving towards a nutrient source out of necessity and can tolerate blue light better than the alternative. Functionally, however, these distinctions are unimportant as the most-avoided colour of light is the most appropriate for use as a repellent. LEDs were chosen for this investigation for being cheap, utilitarian, small and power efficient: these are all highly desirable characteristics in emergent computing devices.

The most significant limitation of this model was the vast propagation delays which prevent the construction of combinational logic circuits, e.g. adders. The prototype gates produced here therefore hold no specific practical value other than demonstrating an implementation of optically-programmable logic. Aside from optimising the experimental environment (reducing the effects of desiccation, changing the size of the gaps between electrodes etc.), the most viable research direction for producing gates of practical value would appear to be minimising the role of physical movement. In Ref. \cite{Adamatzky_Jones_2011}, Adamatzky and Jones demonstrated that \emph{P. polycephalum} responds to a range of environmental stimuli with predictable patterns of electrical activity (when measured at the outer membrane); if the electrical responses of the plasmodium to light are similarly predictable, culturing the plasmodium on a large array of electrically conductive surfaces would therefore allow the researcher to track plasmodial responses electrically. This could potentially vastly reduce propagation delays and hence allow for gate cascading.

Another limitation of this model was physical size. The test environments used (Petri dishes) were 9~cm in diameter, meaning that a single gate was comparatively large, especially when including the breadboard and power supply. If the aforementioned propagation delays could be negated, combinational logic circuits would still be large and ungainly --- e.g. a one-bit half adder comprised of 7 PNAND gates would use approximately 0.5~m\textsuperscript{2} on a bench top --- and therefore of limited practical value. If future prototypes still relied on slime mould taxis, each test environments could theoretically be scaled down to approximately 3~cm in diameter. Multi-gate testing environments with in-built mechanisms for cascading gates and reducing solid-state components could further scale the devices down.

	The electrical output of these gates was extremely low, even when the output circuit voltage was increased to over 24 Volts. This was likely in part due to the high resistance created by the agar blobs (which were measured in Ref. \cite{Mayne_Adamatzky_2013} and found to have an average resistance of 18K$\Omega$ each); any next-generation device would therefore benefit from using a more electrically-conductive substrate. This could present significant barriers to research as conductive gels, such as those used for skin-electrode contacts in medical science, typically contain high amounts of salt to facilitate electrical conductivity which would likely be harmful to the plasmodium.

		After discussing the limitations and failings of these devices, it is pertinent to mention their redeeming features. The devices produced here were all cheap and comparatively easy to fabricate; they are reusable and produce no hazardous waste. The power consumption was also relatively low, although the power to processing power ratio was far inferior to conventional solid-state hardware. Furthermore, the apparent robustness of slime mould devices to adverse conditions such as overvoltage highlights how Physarum-based computing may be adapted to roles to which standard computing architectures are not suited, e.g. in harsh conditions applications and research. The experiments designed here may also be used to model and approximate the growth and behaviour of other similarly structured (i.e. branching) cells or organisms, such as neurons, melanocytes, plant roots and fungi.

	It is relevant to note that although the devices presented here principally relied on light as a controlling input, other inputs are entirely feasible to explore. For example, \emph{P. polycephalum} also exhibits magnetotaxis \cite{Pazur_Schimek_Galland_2007}, thermotaxis (to an ideal temperature of 29$^\circ$C) \cite{Tso_Mansour_1975}, negative chemotaxis away from repellents \cite{Costello_Adamatzky_2013} and negative thigmotaxis \cite{Adamatzky_2013b}, to name a few of the other senses the plasmodium is able to perceive and coordinate a response to. Any of these may be used as viable attractants/repellents in slime mould devices. Indeed, a putative next-generation slime mould logic device could employ multiple inputs.

	Although revisiting moving-part computing elements is essentially retrogressive, it must be remembered that a Physarum machine is not a linear device, i.e. its actions are not comparable to, for example, a vacuum tube or other moving component whose functionality is extremely limited. The Physarum machine's genetically-encoded behaviour - analagous to firmware - bestow it with a form of judgement which is expressed in its highly efficient foraging behaviour. If a computer is given a command with one figure wrong, the program will terminate; the Physarum machine will adapt and dynamically react to any range of inputs it is given. This ability to reason (automated reasoning) is not a property which has been fully achieved with solid-state hardware.\\

\section{Introduction}

	Biological computing devices are currently a source of intense interest in both research and industry: massively-parallel, amorphous computing substrates such as slime moulds or mammalian brains are far superior to conventional solid-state computing architectures in a trade-off between energy consuming, self-growth and  processing ability, and our rapid approach to the physical limitations of inorganic material, the comparative rarity of earth metals used in electronic components and the toxic by-products produced in solid-state computer manufacture are all prime reasons to investigate unconventional computing devices. Furthermore, if sufficiently developed, such devices may allow us insight into other novel bio-computing technologies, e.g. neurally-integrated medical or augmentative prosthetic implants. \emph{P. polycephalum} is an ideal unconventional computing substrate as it is easy to culture, tolerant to abuse and has no ethical issues surrounding its use~\cite{Keller_Everhart_2010}.

The computing potential of the myxomycete slime mould \emph{Physarum polycephalum's} plasmodium (vegetative life-cycle stage, pl. plasmodia) has become a hot topic in computer science, biology and biophysics. Behaving as an excitable or a reaction-diffusion system encapsulated in an elastic membrane~\cite{Adamatzky_2007}, the slime mould may approximate shortest paths between two points~\cite{Nakagaki_2001} and solve maze problems \cite{Adamatzky_2012a,Nakagaki_2001}. Slime mould can be regarded as a proliferating massively parallel processor (a 'Physarum machine' \cite{Adamatzky_2010a}), capable of simultaneous processing of inputs (light, chemical gradients, temperature etc.) \cite{Adamatzky_2007}, concurrent decision making and distributed actuation. In experimental situations, slime mould based computing devices employ chemotaxis --- positive, towards chemoattractants, or negative, away from repellents --- or photoavoidance; the plasmoidum is motile via the contraction of muscle-like proteins invoking rhythmic shuttle streaming of its hydrodynamic core \cite{Adamatzky_2012b,Ishikawa_etal_1991}.

Previously we have demonstrated how \emph{P. polycephalum} may be hybridised with nano-scale metallic particles in order to manually alter the electrical properties of the plasmodium, e.g. facilitate the conduction of electricity and generate permanent, solid-state wires which persist following the plasmodium's migration and/or death \cite{Mayne_etal_2013,Mayne_Adamatzky_2013}.

In this study, we explore the use of light as a programmable input into a Physarum machine, and use the results to implement a basic yet functional device. Significant progress has been made programming a Physarum machine with spatially distributed nutrient sources~\cite{Nakagaki_2001,Adamatzky_2012b,Adamatzky_2010a}, but these are undynamic and problematic to `reprogram' (i.e. alter the input mid-experiment). Indeed, viable logic gates have been designed and produced by a number of authors using different variations of chemical attractant and repellent models~\cite{Tsuda_Aono_Gunji_2004,Adamatzky_2010b,Jones_Adamatzky_2010}, but all suffer from the aforementioned restrictions. These limitations may potentially be overcome by utilising light as an input, or rather a coupling medium between the sub-circuits of logical gates.

Two logic gates were designed: a NOT gate (more commonly known as an inverter), chosen for being the simplest logic operation, and a NAND gate, for being a universal gate, i.e. any logical operation can be carried out by multiple NAND gates.

\section{Selecting inputs: Phototaxis experiments}
\label{phototaxis}

\begin{figure}[!tbp]
\centering
\subfigure[]{\includegraphics[width=0.6\textwidth]{figs/Fig_1a}}
\subfigure[]{\includegraphics[width=0.9\textwidth]{figs/Fig_1bii_arrowed}}
\caption{Phototaxis experiments.
(a)~A scheme of experimental set-up for phototaxis experiments. Gray squares represent agar cubes: the slime mould is inoculated on the central cube (pentagon) and is left to proliferate. Black rectangles represent card barriers with a 0.5~mm gap underneath, which the plasmodium is able to migrate underneath to access the peripheral agar cubes, which are loaded with chemoattractants and illuminated by 2$\times$ LEDs (circles).
(b)~Photograph of a completed phototaxis experiment: note how the plasmodium has migrated towards the left 'blue' side (arrowed).}
\label{fig_taxis}
\end{figure}

Plasmodia were inoculated on to 2.5~cm\textsuperscript{2} squares of 2\% non-nutrient agar in the centre of 9~mm Petri dishes. Similar squares of agar loaded with chemoattractants (oat flakes) were situated at the 3 and 9 o'clock poles of the dishes  and were each illuminated by 2 lid-mounted LEDs (Fig.~\ref{fig_taxis}a). 4 colours of LED were used --- blue (466~nm), green (568~nm), yellow (585~nm) and red (626~nm) --- and each pole had a different colour assigned to it in each experiment. The peripheral regions were separated from the centre by a strip of thick card glued to the dish lid, with a 0.5~mm gap underneath: this prevented most of the light reaching the inoculation point, but allowed the plasmodium to migrate through the gap underneath.

Test plasmodia were therefore given three basic choices in each experiment: migrate towards colour A, migrate towards colour B or do neither (and hence stay in the centre/sclerotinize/attempt to colonise the card/escape the dish etc.). The choices each plasmodium made were compiled into an order of preference for each colour: for each experiment, the colour which was avoided was given 1 'phobia point'; the colour with the most points after all experiments had been completed was determined to be the strongest repellent.

	A completed phototaxis experiment is shown in Fig.~\ref{fig_taxis}b. This pattern of avoidance was identical for every repetition. On no occasion did the plasmodium chose to stay on the central agar square and sclerotinize, although on two occasions it migrated up the card barrier and attempted to consume the glue securing it to the lid of the Petri dish. Phototaxis experiments also took considerably more time to complete than anticipated, with each lasting 5--6 days (whereas similar experiments with only chemoattractants/repellents typically last less time): the plasmodium would typically spend several days exploring the central agar piece and occasionally extend pseudopodia several millimetres along the plastic base of the Petri dish in random directions. This would seem to imply that the plasmodium's actions were a 'last resort' --- a protistic equivalent of reluctance.

\emph{P. polycephalum's} preference to each colour of LED is summarised below, arranged in order of most avoided $\rightarrow$ least avoided:
\begin{center}Green $\rightarrow$ Red $\rightarrow$ Yellow $\rightarrow$ Blue \\\end{center}

Thus green LEDs were used in further experiments.

\section{NOT Gate}

\begin{figure}[!tbp]
\centering
\subfigure[]{\includegraphics[width=0.9\textwidth]{figs/Fig_2a}}
\subfigure[]{\includegraphics[width=0.49\textwidth]{figs/Fig_2b}}
\subfigure[]{\includegraphics[width=0.49\textwidth]{figs/Fig_2c}}
\caption{Physarum NOT logic gate.
(a)~Schematic diagram of a PNOT logic gate. A plasmodium (pentagon) is placed on agar (grey squares) overlying electrode X and chemoattractants are placed on Y. If the input [a] to the LEDs (circles)  is at logic 0, the LED does not switch on and a plasmodial tube is formed between electrodes X and Y; \emph{vice versa} also occurs. The electrodes are wired into an independently-powered output circuit [o].
(bc)~Photographs of a functional PNOT gate.
(b)~Experimental setup with input at logic 1; the plasmodium did not migrate and the output was left at logic 0.
(c)~With input at logic 0, the plasmodium has migrated across the gap and created a `wire' which was electrically conductive, leaving the output at logic 1.
}
\label{fig_pnot}
\end{figure}

	The design for the Physarum-NOT gate (PNOT) is shown in Fig.~\ref{fig_pnot}a. Two aluminium tape electrodes (X and Y) were stuck to the base of a 9~mm plastic Petri dish, separated by a central 10~mm gap. 2~ml agar blobs were placed on top of each electrode, one of which was loaded with plasmodium (X) and the other a chemoattractanct (an oat flake) (Y). The gate's 'input' consisted of a 9~V power supply (Lascar Electronics, Salisbury, UK) connected to a lid-mounted LED overlying electrode Y. At input equals 1 (logic levels of 0 and 1 are represented here as 0~V and 9~V, respectively), the LED was active and prevented the plasmodium from migrating from X to Y, preventing the completion of the circuit. If the input equals 0, the LED was not active and therefore the plasmodium was free to migrate to electrode Y, forming a continuous tube of plasmodium between both electrodes.

	Both electrodes were connected to a separate 'output' circuit with a constant potential difference applied across it. The plasmodium, when spanning two electrodes, acted as a physical conductor of electricity and hence completed the output circuit. As such, the device's functionality is the same as a conventional NOT gate, i.e. if the input is live, the output is not and \emph{vice versa}.

	Experimental photographs of a PNOT gate are shown in Fig.~\ref{fig_pnot}. The time it took for the plasmodium to complete its logic function was extremely variable, ranging from 1--4 days; approximately 25\% entirely failed to propagate. With reference to the results of the phototaxis experiments, this was reasoned to be due to the plasmodium's tendency to stay still for as long as possible before migrating, during which time the agar blobs may have desiccated to the point that they could no longer support its demand for moisture.

\section{NAND Gate}

\begin{figure}[!tbp]
\centering
\subfigure[]{\includegraphics[width=0.9\textwidth]{figs/Fig_3aii}}
\subfigure[]{\includegraphics[width=0.49\textwidth]{figs/Fig_3b}}
\subfigure[]{\includegraphics[width=0.49\textwidth]{figs/Fig_3c}}
\subfigure[]{\includegraphics[width=0.49\textwidth]{figs/Fig_3d}}
\caption{Physarum NAND gate.
(a)~Schematic diagram of a PNAND logic gate. Slime mould (pentagon) is placed on agar (grey squares) overlying electrode X and attractants are placed on Y and Z. If inputs A or B are at logic 1, the corresponding LEDs (circles) turn on, preventing slime mould growth to the corresponding electrode and hence stop a circuit from forming. If both A and B are on, no circuit forms, resulting in no output and \emph{vice versa}. The electrodes are wired into an independently-powered output circuit [o]. A card barrier separates electrodes Y and Z to prevent light contamination.
(bcd)~Photographs of a functional PNAND gate.
(b)~T=6~h; the inoculation point is fully colonised and the left electrode is illuminated.
(c)~ T=12~h; the plasmodium has migrated to the left of the inoculation point.
(d)~T=24~h; the plasmodium has migrated to the unilluminated electrode.
}
\label{Fig_nand}
\end{figure}

Physarum-NAND (PNAND) gates were constructed based upon the same principles of the PNOT gate, but with a second electrode and LED added to the device (Fig.~\ref{Fig_nand}a). LEDs A and B were linked to separate power sources (i.e. inputs) so that their corresponding LEDs worked independently of each other. A plasmodium was once again loaded onto electrode X, but was given the added option of migrating towards electrode Z as well as Y. If both inputs equal 0, the plasmodium could migrate towards either Y or Z at random, forming a circuit between the initial and new electrode, causing the output  to equal 1. If input A is 1 but B is 0, the plasmodium at electrode X is still free to complete the circuit by migrating to electrode Z, and \emph{vice versa}. If both A and B equal 1, the plasmodium would not migrate towards either electrode, preventing a circuit from forming and hence leaving the output at 0. Electrodes X and Y were physically separated by a barrier similar to those described in Sec.~\ref{phototaxis} to prevent light contamination from opposing LEDs. Experimental photos of PNAND gate are shown in Fig.~\ref{Fig_nand}b--d.

\section{Reusability}

	Both PNOT and PNAND gates were tested for reusability. Once a logical operation had been carried out, the inputs were changed and the consequent responses of experimental plasmodia were recorded.

Both PNOT and PNAND gates are reusable, but only a finite number of times. Resetting gates --- i.e. allowing the plasmodium to propagate between two unilluminated points then illuminating one of them in expectation that the plasmodium would withdraw --- was successful in that it caused the plasmodium to withdraw from the target electrode in all instances. This process usually took 2-6 hours and involved the contents of the plasmodium flowing back down the tubes spanning both electrodes. Occasionally, the plasmodium split into two separate fragments, both of which sclerotinized. With only two electrodes present, the PNOT gate was therefore effectively reset as the conductive trail between both electrodes was disrupted. PNOT gates were not successfully re-reset, i.e. the slime mould would not migrate back to an area it had already colonised and withdrawn from.

PNAND gates were successfully reset and reprogrammed, as shown in Fig.~\ref{fig_reset}: once a logical operation had been completed, the inputs were altered and the slime mould migrated accordingly to complete the new operation. PNAND gates that had been reset had a higher failure rate than `fresh' gates, but their propagation delay was significantly shorter (typically $<$~24 hours).

\begin{figure}[!tbp]
\centering
\subfigure[]{\includegraphics[width=0.47\textwidth]{figs/Fig_4a}}
\subfigure[]{\includegraphics[width=0.51\textwidth]{figs/Fig_4b}}
\subfigure[]{\includegraphics[width=0.51\textwidth]{figs/Fig_4c}}
\caption{Resetting and reprogramming a PNAND gate.
(a)~T=36~h. A logical operation has been completed by the slime mould: the right LEDs are on and so it has migrated to the left electrode.
(b)~T=38~h. The inputs are switched; the right LEDs are now off and the left are on. The slime mould begins to withdraw towards the inoculation point and sends pseudopodia in multiple directions looking for a new place to migrate to.
(c)~T=42~h. The slime mould migrates towards the unilluminated right electrode and begins to withdraw its protoplasm from the left electrode and the tube between left and centre electrodes gradually ceases to be conductive.
}
\label{fig_reset}
\end{figure}

\section{Fault-tolerance}

The fault tolerance of Physarum gates was also investigated by altering physical variables in the test environment, including:
\begin{enumerate}
\item Increasing the luminosity of the LEDs by replacing them with brighter variants (each was increased by approximately 50~mcd).
\item Increasing the distance between electrodes to 20~mm
\item Increasing the potential difference between electrodes to 24~V
\end{enumerate}

None of these variables were `decreased', as the initial variables were originally designed as minimal operational parameters; dimmer light sources were anticipated to be ignored by the plasmodium, decreasing the electrode-electrode gap to below 10~mm could result in the plasmodium easily spanning all three (resulting in gate mis-operation) and applying a lower voltage through the electrodes would likely have resulted in an output so low as to be undetectable.

	Whilst adjusting these variables is of little practical use in the design of these specific prototype devices, they were chosen to (a) investigate whether fine-tuning any of these variables would have a beneficial effect on the operation of the device, and (b) demonstrate how a Physarum machine may adapt to unfavourable conditions in a way that traditional computing architectures cannot --- e.g. it is doubtful that any solid state architectures would function properly following exposure to prolonged overvoltage. By introducing extra elements of environmental instability, \emph{P. polycephalum} was given the opportunity to demonstrate how resilient it is, and therefore how robust a Physarum machine could potentially be.

Increasing the luminosity of the LEDs appeared to have no beneficial or detrimental effect on their operation; failure rates were not significantly different and neither was the operation time/propagation delay. Whether \emph{P. polycephalum} is generally phobic to any intensity of green light, or whether the increase in incandescence was too low to provoke any stronger responses was not ascertained.

Doubling the distance between the electrodes markedly increased the failure rate of both varieties of Physarum logic gate. This was likely to result from the combined effects of \emph{P. polycephalum's} reluctance to cross the arid plastic regions of Petri dishes and a lack of any detectable nutrient gradients due to the distant proximity of the supplied chemoattractants.

Experimental plasmodia appeared to be extremely tolerant of high voltages. Whilst the current that would have run through the plasmodium would have been comparatively low due to the resistance of the agar blobs and the innate resistance of the plasmodium its self, it was nevertheless tolerant to applied potential differences of up to and including 24~V. The morphology of plasmodia exposed to high voltages was slightly different to the wild type; rather than spanning the gap between electrodes with a single tube of plasmodium, plasmodia exposed to higher voltages tended to occupy a significant amount of space between electrodes with multiple interlinking tubules (Fig.`\ref{Fig_5}b).

\section{Conclusions}

The principle finding of this study was that it is indeed possible to `program' an unaltered \emph{P. polycephalum} plasmodium to carry out basic logic functions with optical inputs. Whilst the prototype devices presented here have significant limitations, they demonstrate a genuine implementation of unconventional computation in which a living computing substrate functions as both a sensing, decision making computing device and a physical conductor of electricity.

\emph{P. polycephalum's} photoavoidance is widely mentioned in literature, and indeed the plasmodium is best cultivated in the dark. Sauer~\cite{Sauer_1982}  reports that the plasmodium is phobic to blue light (which may also induce sporulation whilst the plasmodium is also starved of nutrients and/or moisture), but may sometimes migrate towards red light. Green light has been reported to inhibit sporulation, and seemingly random wavelengths of every colour (and some UV wavelengths) induce a range of metabolic changes \cite{Hato_etal_1976,Sauer_1982}. Whilst \emph{P. polycephalum} is not photosynthetic, it contains pigments similar to phytochromes which are thought to act as photoreceptors \cite{Kakiuchi_etal_2001}. A migrating plasmodium may be diverted by placing illuminated barriers in its path \cite{Adamatzky_2009}.

There were clear discrepancies between the results of the phototaxis experiments conducted for this investigation and those in published literature. There were, however, stark differences in experimental setup between the experiments described here and those in literature; LED diffusion, different wavelengths of light used, different intensities of light and different experimental conditions could all have altered plasmodial responses drastically. It should also be noted that if, for example, the plasmodium migrated towards blue light, this does not necessarily imply that \emph{P. polycephalum} is preferentially moving towards the blue light, but is more likely moving towards a nutrient source out of necessity and can tolerate blue light better than the alternative. Functionally, however, these distinctions are unimportant as the most-avoided colour of light is the most appropriate for use as a repellent. LEDs were chosen for this investigation for being cheap, utilitarian, small and power efficient: these are all highly desirable characteristics in emergent computing devices.

The most significant limitation of this model was the vast propagation delays which prevent the construction of combinational logic circuits, e.g. adders. The prototype gates produced here therefore hold no specific practical value other than demonstrating an implementation of optically-programmable logic. Aside from optimising the experimental environment (reducing the effects of desiccation, changing the size of the gaps between electrodes etc.), the most viable research direction for producing gates of practical value would appear to be minimising the role of physical movement. In Ref. \cite{Adamatzky_Jones_2011}, Adamatzky and Jones demonstrated that \emph{P. polycephalum} responds to a range of environmental stimuli with predictable patterns of electrical activity (when measured at the outer membrane); if the electrical responses of the plasmodium to light are similarly predictable, culturing the plasmodium on a large array of electrically conductive surfaces would therefore allow the researcher to track plasmodial responses electrically. This could potentially vastly reduce propagation delays and hence allow for gate cascading.

Another limitation of this model was physical size. The test environments used (Petri dishes) were 9~cm in diameter, meaning that a single gate was comparatively large, especially when including the breadboard and power supply. If the aforementioned propagation delays could be negated, combinational logic circuits would still be large and ungainly --- e.g. a one-bit half adder comprised of 7 PNAND gates would use approximately 0.5~m\textsuperscript{2} on a bench top --- and therefore of limited practical value. If future prototypes still relied on slime mould taxis, each test environments could theoretically be scaled down to approximately 3~cm in diameter. Multi-gate testing environments with in-built mechanisms for cascading gates and reducing solid-state components could further scale the devices down.

	The electrical output of these gates was extremely low, even when the output circuit voltage was increased to over 24 Volts. This was likely in part due to the high resistance created by the agar blobs (which were measured in Ref. \cite{Mayne_Adamatzky_2013} and found to have an average resistance of 18K$\Omega$ each); any next-generation device would therefore benefit from using a more electrically-conductive substrate. This could present significant barriers to research as conductive gels, such as those used for skin-electrode contacts in medical science, typically contain high amounts of salt to facilitate electrical conductivity which would likely be harmful to the plasmodium.

		After discussing the limitations and failings of these devices, it is pertinent to mention their redeeming features. The devices produced here were all cheap and comparatively easy to fabricate; they are reusable and produce no hazardous waste. The power consumption was also relatively low, although the power to processing power ratio was far inferior to conventional solid-state hardware. Furthermore, the apparent robustness of slime mould devices to adverse conditions such as overvoltage highlights how Physarum-based computing may be adapted to roles to which standard computing architectures are not suited, e.g. in harsh conditions applications and research. The experiments designed here may also be used to model and approximate the growth and behaviour of other similarly structured (i.e. branching) cells or organisms, such as neurons, melanocytes, plant roots and fungi.

	It is relevant to note that although the devices presented here principally relied on light as a controlling input, other inputs are entirely feasible to explore. For example, \emph{P. polycephalum} also exhibits magnetotaxis \cite{Pazur_Schimek_Galland_2007}, thermotaxis (to an ideal temperature of 29$^\circ$C) \cite{Tso_Mansour_1975}, negative chemotaxis away from repellents \cite{Costello_Adamatzky_2013} and negative thigmotaxis \cite{Adamatzky_2013b}, to name a few of the other senses the plasmodium is able to perceive and coordinate a response to. Any of these may be used as viable attractants/repellents in slime mould devices. Indeed, a putative next-generation slime mould logic device could employ multiple inputs.

	Although revisiting moving-part computing elements is essentially retrogressive, it must be remembered that a Physarum machine is not a linear device, i.e. its actions are not comparable to, for example, a vacuum tube or other moving component whose functionality is extremely limited. The Physarum machine's genetically-encoded behaviour - analagous to firmware - bestow it with a form of judgement which is expressed in its highly efficient foraging behaviour. If a computer is given a command with one figure wrong, the program will terminate; the Physarum machine will adapt and dynamically react to any range of inputs it is given. This ability to reason (automated reasoning) is not a property which has been fully achieved with solid-state hardware.\\


The Hercules-Aquila Cloud

Abstract: We present evidence for a substantial overdensity of stars in the direction of the constellations of Hercules and Aquila. The Cloud is centered at a Galactic longitude of about 40 degrees and extends above and below the Galactic plane by at least 50 degrees. Given its off-centeredness and height, it is unlikely that the Hercules-Aquila Cloud is related to the bulge or thick disk. More likely, this is a new structural component of the Galaxy that passes through the disk. The Cloud stretches about 80 degrees in longitude. Its heliocentric distance lies between 10 and 20 kpc so that the extent of the Cloud in projection is roughly 20 kpc by 15 kpc. It has an absolute magnitude of -13 and its stellar population appears to be comparable to, but somewhat more metal-rich than, M92.

\section{Introduction}

The construction of the Milky Way via the hierarchical aggregation of
small galaxies has left behind debris in the form of streams,
satellites and substructure in the Galactic halo. There are many
examples now known of such detritus, including the the disrupting
Sagittarius galaxy and its stream~\citep[e.g.,][]{Ib01,Ma03,Be06b},
the ``Monoceros Ring''~\citep[e.g.,][]{Ya03} and the Orphan
Stream~\citep[e.g.,][]{Be07}.  In addition to the easily recognizable
streams, the depth and wide area coverage of the Sloan Digital Sky
Survey (SDSS) enables identification of other less well-defined
features.  A substantial overdensity in the direction of
Virgo~\citep{Ju06} and a possible underdensity in the direction of
Ursa Major~\citep{Xu06} have been found. Here, we identify and
characterize another overdensity of stars, called the {\it Hercules-Aquila
Cloud}.

\section{The Hercules-Aquila Cloud}

\begin{figure*}[t]
\begin{center}
\includegraphics[height=8.5cm]{f1.eps}
\caption{\label{fig:bigmap} The areal density of SDSS stars with $0.1
< g-i < 0.7$ and $20 < i < 22.5$ in Galactic coordinates. The color
plot is an RGB composite with blue representing $0.1 \le g-i < 0.3$,
green $0.3 \le g-i <0.5$ and red $0.5\le g-i < 0.7$. Structures
visible on the figure include the arch of the Sagittarius stream and
the Orphan Stream in blue and the multiple wraps of the Monoceros ring
in red. Behind, and offset from, the Galactic bulge, there is a
yellow-colored structure visible in SDSS data in both hemispheres,
centered at $\ell \approx 40^\circ$. This is the Hercules-Aquila
Cloud.  Overlaid on the plot are contours of constant extinction
(black is 0.1 mag and grey is 0.25 mag in the $i$ band). The two red
vertical lines show the SEGUE stripes at $\ell = 31^\circ$ (runs 5378
and 5421) and $50^\circ$ (run 4828). The inset is a color-magnitude
diagram of the SDSS data, with blue, green and red bands encompass the
stars used to create the large figure.}
\end{center}
\end{figure*}
\begin{figure*}
\begin{center}
\includegraphics[height=8.5cm]{f2a.eps}
\includegraphics[height=8.5cm]{f2b.eps}
\caption{\label{fig:slicea} The upper panels show greyscale contours
of the spatial density of SDSS stars with $0.25 < g -i < 0.4$ in two
slices at longitude $\ell = 31^\circ$ (left) and $\ell = 50^\circ$
(right). The spatial density is measured in ${\rm mag}^{-1} {\rm
deg}^{-2}$.  In each case, the leftmost (rightmost) colored lines
shows the variation of extinction (density) with latitude in each
hemisphere. Green denotes the current profile, whilst red refers to
the opposite hemisphere and blue the difference. There are a number of
prominent enhancements of density along the line of sight, as
manifested by green-colored profiles, which are not smoothly decaying
curves but show bumps.  The southernmost overdensities are associated
with the Sagittarius stream ($\ell = 31^\circ, b \approx -35^\circ, i
\approx 20.5$ mag and $\ell = 50^\circ, b \approx -40^\circ, i \approx
20.5$ mag). The other enhancement is part of the Hercules-Aquila
Cloud, which is revealed in the difference plots shown in the lower
panels. These are created by subtracting the density in the two
hemispheres so that black (white) indicates relative overdensity
(underdensity). The Hercules-Aquila Cloud is at $\ell = 31^\circ, b
\approx 20^\circ, i \approx 19.$ mag and $\ell = 50^\circ, b \approx
-15^\circ, i \approx 20.$ mag.}
\end{center}
\end{figure*}
\begin{figure}[t]
\begin{center}
\includegraphics[width=8.cm]{f3a.eps}
\includegraphics[width=8.cm]{f3b.eps}
\includegraphics[width=8.cm]{f3c.eps}
\caption{\label{fig:lonmap} Three density slices at $b = 70^\circ$
(top), $55^\circ$ (middle) and $40^\circ$ (bottom), each of which is
$15^\circ$ wide. The key gives the scale in number of stars ${\rm
mag}^{-1} {\rm deg}^{-2}$ satisfying the color cut $0.25 <g-i
<0.4$. The variation of the extinction along the strips and the
fraction (between 0 and 1) covered in DR5 are shown as a function of
longitude at the bottom of the panels. Three structures are
discernible in the panels and they are clearly separated. In the top
panel, the Sagittarius stream and Virgo Overdensity are seen to cross
each other obliquely, whilst the tip of the Hercules-Aquila Cloud is
just visible; in the middle panel, the Sagittarius stream is
bifurcated and cleanly differentiated from the Virgo Overdensity,
whilst the Hercules-Aquila Cloud is becoming more substantial; in the
lowest panel only the Hercules-Aquila Cloud is seen, due to SDSS DR5
coverage.}
\end{center}
\end{figure}
\begin{figure}[t]
\begin{center}
\includegraphics[width=9.cm]{f4.eps}
\caption{\label{fig:cmd} Left: Hess diagrams for the on-Cloud field at
$\ell = 30^\circ, b=40^\circ$ (dotted lines), and the off-Cloud field
at $\ell = 15^\circ, b=40^\circ$ (solid).  The units are $10^4$ stars
per square magnitude.  Right: Difference of the Hess diagrams divided
by the square root of their sum.  This shows essentially the
signal-to-noise ratio in the difference. There is one obvious
overdensity corresponding to upper main sequence and turn-off stars in
the Hercules-Aquila Cloud. The white lines show M92 ridgelines shifted
to the distance of 10 and 20 kpc, which bracket the main sequence
location.}
\end{center}
\end{figure}

SDSS imaging data are produced in five photometric bands, namely $u$,
$g$, $r$, $i$, and $z$~\citep[see e.g.,][]{Fu96,Ho01,Sm02,Gu06} and
are automatically processed through pipelines to measure photometric
and astrometric properties and to select targets for spectroscopic
follow-up~\citep{Lu99,St02,Pi03,Iv04}.  In this paper, all the
photometric data are de-reddened using the maps of~\citet*{Sc98}.  We
use primarily Data Release 5 (DR5), which covers $\sim 8000$ square
degrees around the Galactic North Pole, and 3 stripes in the Galactic
southern hemisphere~\citep{AM07}.  Whilst the motivation for the
original SDSS survey was largely extragalactic, a substantial portion
of its successor survey SDSS-II is devoted to Galactic structure (the
Sloan Extension for Galactic Understanding and Exploration or
SEGUE). There are a number of SEGUE imaging stripes that run from
positive to negative latitudes through the disk at constant longitude.

Fig.~\ref{fig:bigmap} shows an RGB composite density map of stars with
$20 < i < 22.5$ and $0.1 < g-i < 0.7$ in Galactic coordinates centered
at $(\ell =0^\circ, b= 0^\circ$).  As illustrated in the inset, the
blue color represents the density of stars with $0.1 \le g-i < 0.3$,
green those with $0.3 \le g-i <0.5$ and red those with $0.5\le g-i <
0.7$. A number of the prominent substructures that dominate the halo
are marked, including the Sagittarius Stream, the Monoceros Ring, the
Virgo Overdensity and the Orphan Stream.  The subject of this {\it
Letter} is the gold-colored overdensity that is most prominent in DR5
at longitudes of $\sim 40^\circ$ and is visible in both hemispheres.
This is the Hercules-Aquila Cloud.  Overlaid on the plot are contours
of constant total extinction from~\citet{Sc98}.

To investigate this structure, we exploit the SEGUE imaging at $\ell =
31^\circ$ (runs 5378 and 5421) and $50^\circ$ (run 4828), shown as
the two red dashed lines in Fig.~\ref{fig:bigmap}. Each SDSS and
SDSS-II stripe comprises two interlocking runs, each of which have a
width of $2.34^\circ$ but covers half the area of the stripe.  In
Fig.~\ref{fig:slicea}, we show the density of stars in the runs with
$0.25 < g-i < 0.4$ in units of ${\rm mag}^{-1}{\rm deg}^{-2}$. The
color cut is designed to identify primarily halo stars. We also show
profiles of density and extinction as well as their differences
(in red, green and blue lines respectively). We see that, comparing
the Northern and Southern hemispheres, there is dramatic asymmetry at
$\ell = 50^\circ$ with a substantial overdensity at negative
latitudes. The asymmetry is still present at $\ell = 31^\circ$, but
the overdensity now occurs at positive latitudes.

The structure is seen most cleanly in the difference maps in the lower
two panels of Fig.~\ref{fig:slicea}. Here, overdensities relative to
the other hemisphere are dark, while underdensities are white.  The
Hercules-Aquila Cloud is seen clearly at $\ell = 31^\circ, b \approx
20^\circ, i \approx 19$ mag and $\ell = 50^\circ, b \approx
-15^\circ, i \approx 20$ mag. It is therefore closer to the Sun at
$\ell =31^\circ$ than at $\ell = 50^\circ$.  There are two further
overdensities visible in Fig~\ref{fig:slicea} at $\ell = 31^\circ, b
\approx -35^\circ, i \approx 20.5$ mag and $\ell = 50^\circ, b \approx
-40^\circ, i \approx 20.5$ mag. They are most likely associated with
the Sagittarius stream.  The orbital plane of the Sagittarius dSph as
computed in \citet{Fe06} is shown as a blue dotted line in
Fig~\ref{fig:bigmap}. Its tidal debris has a thickness of $\sim
15^\circ$ on the sky and so should intersect the SEGUE imaging runs at
their southernmost latitudinal points.

In the slice at $\ell = 50^\circ$, the total number of excess stars
with the color cut $0.25 < g-i < 0.4$, the latitude cut $-35^\circ <
b<0^\circ$, and within the contour corresponding to $n_{\rm stars} =
50$ in the lower right panel of Fig~\ref{fig:slicea} is $10^4$. This
ignores the contribution above the plane, and so is a conservative
estimate.  Scaling to the overall longitudinal extent of the Cloud, we
estimate that there are $5 \times 10^5$ stars with $0.25 < g-i <
0.4$. Comparing with the CMD of the metal-poor globular cluster M92,
we etsimate an absolute magnitude of $M_v = -13$ for the overdensity

Fig~\ref{fig:lonmap} shows longitudinal cuts through the
Hercules-Aquila Cloud at three latitudes $b= 70^\circ, 55^\circ$ and
$40^\circ$, each $15^\circ$ wide.  The SDSS data sometimes only partly
cover the slice. To take this incompleteness into account, we
normalize the density in each longitude bin by the fraction covered.
There are three structures visible in the panels.  The bifurcated
Sagittarius stream and the Virgo Overdensity dominate the top two
panels. However, in both panels, a new structure is becoming visible;
this is the Hercules-Aquila Cloud.  In the bottom panel, the Cloud has
come clearly into view, extending over $90^\circ$ in longitude, whilst
the Sagittarius and the Virgo Overdensity are largely missing in the
area covered by the DR5 footprint.

Fig.~\ref{fig:cmd} shows Hess diagrams for a $8^\circ \times 8^\circ$
field centered on $\ell = 30^\circ, b=40^\circ$, and a $16^\circ
\times 16^\circ$ field centered on $\ell = 15^\circ, b=40^\circ$
(left), together with their difference divided by the square root of
their sum (right). The two fields are chosen to lie interior to and
exterior to the Cloud at the same latitude. In the difference panel,
there is a clear upper main sequence and turn-off, corresponding to
stars in the Cloud.  Also shown in white are two M92 ridgelines, using
data from \citet{Cl05}, that are used to bracket the distance range in
the Hercules-Aquila Cloud. These suggest that the Cloud in SDSS DR5 is
primarily composed of upper main sequence and turn-off stars at
heliocentric distances of between 10 and 20 kpc. The turn-off lies
redward of the M92 ridgeline, suggesting that the Cloud is composed of
more metal-rich stars than M92 ([Fe/H] = -2.2).

\begin{figure}[t]
\begin{center}
\includegraphics[width=7cm]{f5.eps}
\caption{\label{fig:vel} Histograms of velocities (converted to the
Galactic standard of rest) for stars with $g-r <1$ and $r < 19$. Thick
grey lines refers to all stars in the field centered on the
Aquila-Hercules Cloud. We select red giant branch candidates from the
CMD in three fields, two off-Cloud (red and blue) and one on-Cloud
(black). The red and blue histograms peak at $\sim 80$ kms$^{-1}$, as
expected for thick disk stars. The black histogram peaks at $\sim 180$
kms$^{-1}$, which we interpret as the velocity signature of the
Cloud. The total numbers of stars in the fields after the CMD cut, as
well as those picked up by the CMD mask, are marked in color above the
histograms.}
\end{center}
\end{figure}

Fig.~\ref{fig:vel} shows the kinematic evidence for the
Hercules-Aquila Cloud using SDSS radial velocities for stars with $g-r
< 1$ and $r < 19$. All the velocities are converted to the Galactic
standard of rest and have a typical uncertainty of $\sim 10$
kms$^{-1}$. The first field covers the Cloud ($20^\circ \le b \le
55^\circ$ and $20^\circ \le \ell \le 75^\circ$). There are two
comparison fields at larger longitude and higher latitude, namely
field 2 at $40^\circ \le b \le 55^\circ$ and $110^\circ \le \ell \le
180^\circ$, and field 3 at $65^\circ \le b \le 80^\circ$ and $20^\circ
\le \ell \le 80^\circ$.  Most of the stars are at high latitude and
selected to be bluer than the thin disk dwarfs. The signal is
therefore dominated by halo and thick disk stars. This is clear from
the thick grey histogram in Fig.~\ref{fig:vel}, which reports
velocities for all stars in Field 1. The distribution is built from
one component with almost zero mean velocity and another with mean at
$\sim 80$ kms$^{-1}$. The all-star velocity distribution is the same
for Fields 2 and 3, and not shown. To pick out the stars in the
structure, we use a mask based on the M92 ridgeline, shifted to 10 and
20 kpc. The distribution for field 1 centered on the Hercules-Aquila
Cloud is shown as a black solid line and can be compared to the
histograms for the off-Cloud fields 2 and 3 (blue and red
respectively). Nearby thick disk stars in these locations have a line
of sight velocity of $\sim 150 \sin \ell \cos b \sim 80$ kms$^{-1}$
with respect to the Galactic standard of rest~\citep[e.g.,][]{Gi89},
which is consistent with the peaks in the red and blue
histograms. There is a clear kinematic signal from stars moving
coherently with a radial velocity of $\sim 180$ kms$^{-1}$ in the
on-Cloud field. This cannot be due to either the halo or the thick
disk.

Could the Hercules-Aquila Cloud be related to the Galactic bulge or
bar? There is an asymmetry in the integrated light of the bulge,
famously seen by the DIRBE instrument on the COBE satellite
~\citep{Dw95}, but it is only of the order of a few degrees at
most~\citep[see e.g.,][]{We94}. This asymmetry would be larger if
confined to stars in a range of heliocentric distances between 4 to 8
kpc, as we would then be primarily picking up stars on the near-side
of the bar. However, our estimate of the distance to the
Hercules-Aquila Cloud is 10 to 20 kpc.  At $\ell = 40^\circ$, the
closest the Cloud can be to the Galactic Center is in the range $6$ to
$15$ kpc. Although the Cloud passes through the disk, it seemingly
lies well beyond the Galactic bar. Also, the Cloud extends to Galactic
latitudes of $\pm 50^\circ$, again much larger than the vertical
scale-height of the bulge or of the thin and thick disks.
\citet{Pa03,Pa04} also discovered some asymmetries in the distribution
of bright, probably thick disc, stars using POSS 1 data, but these
stars are much closer at $\sim 2$ kpc from the Sun.

Could the Hercules-Aquila Cloud be a smooth component that is centered
on $\ell=0^\circ$ but made to appear asymmetric by extinction?  Most
of the extinction is caused by dust with a scalelength of $\sim 100$
pc. However, there is a prominent North Polar Spur of dust that
extends to high Galactic latitudes~\citep[see e.g.,][]{Fr81}. This can
be seen in Fig.~\ref{fig:bigmap} as the extension to the black contour
over latitudes $30^\circ$ to $60^\circ$, where it overlaps part of the
Sagittarius stream. However, a similar amount of extinction appears to
have no effect on the underlying density of detected stars in the
Sagittarius stream and the Virgo Overdensity at $\ell = -100^\circ$,
and so variable extinction is not a plausible explanation of the
offset structure. {\it We conclude that the Hercules-Aquila Cloud is a
real, hitherto unknown, structural component of the inner halo of the
Galaxy.}

\acknowledgments Funding for the SDSS and SDSS-II has been provided by
the Alfred P.  Sloan Foundation, the Participating Institutions, the
National Science Foundation, the U.S. Department of Energy, the
National Aeronautics and Space Administration, the Japanese
Monbukagakusho, the Max Planck Society, and the Higher Education
Funding Council for England. The SDSS Web Site is
http://www.sdss.org/.

The SDSS is managed by the Astrophysical Research Consortium for the
Participating Institutions. The Participating Institutions are the
American Museum of Natural History, Astrophysical Institute Potsdam,
University of Basel, Cambridge University, Case Western Reserve
University, University of Chicago, Drexel University, Fermilab, the
Institute for Advanced Study, the Japan Participation Group, Johns
Hopkins University, the Joint Institute for Nuclear Astrophysics, the
Kavli Institute for Particle Astrophysics and Cosmology, the Korean
Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos
National Laboratory, the Max-Planck-Institute for Astronomy (MPIA),
the Max-Planck-Institute for Astrophysics (MPA), New Mexico State
University, Ohio State University, University of Pittsburgh,
University of Portsmouth, Princeton University, the United States
Naval Observatory, and the University of Washington.

The Hercules-Aquila Cloud