Original Research ARTICLE
Molecular quantum robotics: particle and wave solutions, illustrated by “leg-over-leg” walking along microtubules
- Forschungszentrum Informatik (Centre of Computer Science), Intelligent System and Production Engineering, Interactive Diagnosis and Service Systems, Karlsruhe, Germany
Remarkable biological examples of molecular robots are the proteins kinesin-1 and dynein, which move and transport cargo down microtubule “highways,” e.g., of the axon, to final nerve nodes or along dendrites. They convert the energy of ATP hydrolysis into mechanical forces and can thereby push them forwards or backwards step by step. Such mechano-chemical cycles that generate conformal changes are essential for transport on all different types of substrate lanes. The step length of an individual molecular robot is a matter of nanometers but the dynamics of each individual step cannot be predicted with certainty (as it is a random process). Hence, our proposal is to involve the methods of quantum field theory (QFT) to describe an overall reliable, multi–robot system that is composed of a huge set of unreliable, local elements. The methods of QFT deliver techniques that are also computationally demanding to synchronize the motion of these molecular robots on one substrate lane as well as across lanes. Three different challenging types of solutions are elaborated. The impact solution reflects the particle point of view; the two remaining solutions are wave based. The second solution outlines coherent robot motions on different lanes. The third solution describes running waves. Experimental investigations are needed to clarify under which biological conditions such different solutions occur. Moreover, such a nano-chemical system can be stimulated by external signals, and this opens a new, hybrid approach to analyze and control the combined system of robots and microtubules externally. Such a method offers the chance to detect mal-functions of the biological system.
Molecular robotics, which operates on a nano scale, has in the last decade witnessed impressive growth, (e.g., Murata et al., 2013). The current topics in this field are molecular machines (Balzani et al., 2008; Roux, 2011; Fukuda et al., 2012; Seeman, 2014) and, even more important in the context of this contribution, biped DNA walkers (Sherman and Seeman, 2004; Shin and Pierce, 2004; Omabegho et al., 2009; Lund et al., 2010). In a previous contribution, the biological process of muscle contraction by the “one-legged” motion of myosin II along actin-filaments has been described using QFT methods (Haken and Levi, 2012). The results of different types of synchronization of multi-molecular systems have been conclusive in the sense that beside the more classical impact solution, two quantum mechanical solutions exist that can describe the synchronization of billions of unreliable molecules. This result can usually not be modeled by pure classical particle solutions. Despite the theoretical existence of these innovative solutions, knowledge of a great set of experimentally confirmed data is unfortunately still now very limited. Nevertheless, these theoretical outcomes encourage us to continue with our quantum theoretical approach, since we are convinced that in the near future confirming experimental data will be available.
This paper focuses on a description of the “leg-over-leg” motion of molecular robots along tubulin strands. These biological nano-robots (motor proteins) walk along such lanes in a four-step process, and are fueled by ATP consumption (in physics, energy exchange with a heat bath) after the first and third steps. From a system perspective, this denotes that we are modeling an open system that is not in a thermal equilibrium. The C-terminal domain, e.g., of kinesin-1, act as a gripper and is attached to cargo, e.g., vesicles that are transported along axons and dendrites, which are both connected to neuron cell bodies. Very similar processes occur for the much larger dynein, with two legs (in biological terms two “heads”), which also transports vesicles along axons and dendrites. Due to this similarity, this paper focuses on a general description of the motion processes of both molecular robots. We abandon the modeling of cargo transport in this paper for the sake of understandability.
In this paper, a set of such molecules is regarded as a swarm of molecular robots that must be configured, synchronized, and in real time supported by energy to perform the next step (Levi and Haken, 2010). Due to the utilization of a quantum theoretical approach, such molecules can be described as particles or as waves (according to matter/field dualism). The following three different solutions are presented:
(1) Impact (stroke) solution of one molecular robot on a single substrate lane,
(2) Coherent motion of many robots on and between parallel substrate lanes,
(3) Running wave that synchronizes the motions of many robots on one substrate lane.
In all three solutions, a local B-field is activated within a lane and across different lanes. Such a field not only controls the process of energy consumption but also, even more importantly, acts as a signal field that synchronizes the steps of the molecular robots. This field is generated by the molecular robots themselves and is not externally injected.
The impact solution pertains to the motion of a single molecular robot on one lane. Damping effects transform the typical wave characteristics of a QFT approach to particle behavior. The appearance of a sequence of impacts pushes the molecule (particle) forwards or backwards step by step.
The next two solutions plainly reveal the wave features of our approach. The coherent motion solution distinguishes between the individual robots walking on a lane and across the lanes. In contrast to the second solution, the running wave solution presents a result that can be obtained in special restrictive conditions.
The “leg-over-leg” walking of a kinesin molecule (dynein) is modeled as a bipedal molecular robot r, where the two heads (light chains) are considered as two legs, and the “coiled-coil tail” as an effector (Alberts et al., 2008). The track is a microtubule surface (substrate lane), along which r moves in discrete steps. During such a walk, a robot r can take one of four leg states and one of two walking states:
Leg state a: leg 1 or leg 2 points backwards.
Leg state b: leg 1 or leg 2 points forwards.
Leg-over-leg state c1: leg 1 is loosely bound to the top of leg 2,
which is tightly attached to the substrate.
Leg-over-leg state c2: leg 2 is loosely bound to the top of leg 1,
which is tightly attached to the substrate.
Walking state 1: r is moving forward with leg 1.
Walking state 2: r is moving forward with leg 2.
Figure 1 shows the four different leg states.
Figure 1. Representation of the two leg states a and b and the two leg-over-leg states c1 and c2. The walking direction is from left to right. For simplicity, we present in this diagram each leg as tightly connected to the substrate (crossed circle). The loosely bound states connected to the substrate are omitted (see Figure 2).
The two legs are attached to the substrate in two different connection states:
Connection state 1: leg 1 is tightly attached to position pk of the tubulin substrate and leg 2 is weakly attached to site pk + 1.
Connection state 2: leg 2 is tightly attached to site pk and leg 1 is weakly attached to position pk + 1 of the substrate.
The level of attachment to the substrate in these two connection modes is defined by the two possible states of the substrate:
Ground state g: represents a weak attachment of a leg to the substrate molecule.
Excited state e: represents a strong attachment of a leg to the substrate molecule.
The discrete periodic movement of the molecular robot r is characterized by a four-step cycle (Figure 2) that starts with the two leg states (1a, 2b), continues with the transference of these states in step 1 into the leg-over-leg state (1c1, 2c1), and is maintained by step 2, which produces the leg states (2a, 1b). In step 3, the leg-over-leg state (1c2, 2c2) is achieved, and finally, in step 4, the states (1a, 2b) are again established.
Figure 2. Demonstration of one complete “leg-over-leg” cycle of a molecular robot r which walks along a substrate lane. The loosely bound states are marked by crossed circles . The tightly bound states are marked by striped circles . The first position is fixed by k = 1; the initial state is (1a, 2b).
The following notation is employed for the creation (or annihilation) operators:
Molecular robot r: r†leg1 statej positionk; leg2 statej′ positionk′; e.g., r†l1ap1; l2bp2.
Substrate s: s†statejpositionk; e.g., s†gp1.
Table 1 subsumes the different individual motion patterns of the two legs of r.
Table 1. Representation of the motion patterns of a “leg-over-leg” walking process of a molecular robot r.
The energy transfer process (s → r) occurs in steps 1 and 3 and is combined with the excited status of the substrate. In addition, the heat-bath operators (“fueling” and synchronization operators) B†pk and B†pk are applied in this walking phase. The operator transfers during the four steps are summarized as follows (readout from Table 1):
This transfer list of operators delivers the definition of the interaction Hamiltonian Hint and finally the elaboration of the resulting equations of motion.
The local interaction Hamiltonian Hint at positions to p1 to p4 is defined by expression 3.1, where the first step starts at the initial site k marked in Figure 2. Further, two real coupling constants, g1 and g2, are introduced, where the parameter g1 describes the uneven steps and g2 denotes the even steps.
By the replacements p1 → pk, p2 → pk + 1, etc., Hint is transferred into Hint, k and has to be summed up by k = 1, 2, …. All operators that are quoted here are Bose operators. Generally, we should involve both Fermi operators (anti-commutation rule) and Bose operators (commutation rule) for the definition of Hint and the resulting calculations of the Heisenberg equations of motions. Here, we use only Bose operators because in this case aggregations of molecules are allowed since the Pauli Exclusion Principle does not have to be applied. We neglect Fermi operators that depict supplementary interactions of Fermions with Bosons.
Endowed with the Hamiltonian Hint, the corresponding full equations of motion have been calculated and are presented in the next sub-chapter.
Full Heisenberg Equations of Motion
The full set of Heisenberg equations of motion is completed by 7 damping constants, γab, etc., 7 fluctuating forces, F†ab, etc., and two coupling constants, g1, g2, and is defined by the following expressions, representing coupled nonlinear, delayed, and complex operator equations.
To solve the Equations (4.1–4.7), we utilize the semi classical approach and replace the operators by their expectation values in a coherent state representation. In doing so, the operators become complex numbers and the expectation values of the fluctuating forces are zero.
We start the description of the three a fore-mentioned solution with the impact solution and are firstly searching for solutions for each step separately. Later, one overall solution will be outlined by combining all four separate solutions. We are hereby concentrating on one robot walking on one lane. The walking processes of several robots on different lanes will be explained in the succeeding sections.
Impact Solution: First Step
As mentioned above, the procedure begins with the first step (Figure 2), at position k. The initial state is defined as (starting with position k = 1):
where | ϕ0 〉 delineates the vacuum state. The initial conditions of all relevant states are:
We set = 0 in Equation (4.1) because this state is not present. The resulting modified Equation (5.1) now reads:
In a similar way, going through all remaining equations the following reduced set of equations is obtained:
Before presenting the complete set of numerical solutions for all Equations (5.4–5.10), we prefer to assure ourselves that the right solution can also be generated by analytical methods and not only by numerical calculations. This should convince us that we are on the right track.
The following approach (setting equally the robot-relevant damping constants γ = γab = γba = γc1 = γc2, and neglecting for the moment the damping of the B operator: γB = 0), delivers a consistent, periodic solution for a walking biped molecular robot:
By inserting these expressions into the Equations (5.4–5.10), four coupled equations for the four variables f, g, h, b are finally obtained:
The last equation can be directly solved by the technique of separation of variables:
The calculation is started with the assumption that γ = 0 is valid in order to start with “perfect” symmetry. Afterwards, we will switch to the damping process in order to observe how this symmetry will be broken.
Incidentally, it is not surprising if all solutions of the Equations (5.14–5.17) are periodic (wave-like) if the damping process is switched off. Figure 3 demonstrates this prediction in a phase portrait (Holmes et al., 2012) that also includes the expected periodicity of the variable b (B-field):
Figure 3. Three-dimensional limit cycle defined by the three variables b, g, and h. The initial point is marked by a cross.
Even if the general complex solutions of the Equations (5.4–5.10) are calculated, symmetrical patterns are visible. Figure 4 shows the real parts of all operators with the expected symmetry. The imaginary parts show a very similar behavior; therefore they are not presented here or in the next sections.
Figure 4. Temporal dependence of the real parts of all variables of step 1, with γ = 0. The algebraic symbols are: . The coupling constants are g1 = g2 = 0.1.
It is well known that processes without any damping are artifacts. Therefore, we will now include damping processes. Even if “perfect symmetry” is not realistic, it is worth taking it into consideration since it serves as a kind of “roadmap” to the new individual trajectories if γ deviates from zero and becomes positive.
Even if the damping constant is assigned the small value γ = 0.0005 the symmetry (periodicity) will be broken. The limit cycle shown in Figure 3 is dissolved and the trajectory goes directly to zero (fixed point). The same happens with the general solution—it goes directly to zero. Thus, in both cases a large γ value should be chosen.
If the value of γ is continuously decreased, the ringing effect increases along the trajectory to γ = 0. Therefore, we can achieve all possible trajectory patterns, from a direct path to zero to complete oscillation (no damping).
Even if the value of the damping constant slightly diminishes, to γ = 0.0001, the first effect can be observed; after the first peak the trajectory converges to zero (Figure 5).
Figure 5. Damped temporal curves representing the real parts of the following operators of the first step: . The value of the damping constant is set to γ = 0.0001. The coupling constants are g1 = g2 = 0.1.
The results of all remaining steps, 2–4, will not be presented in the same detail as in step 1. Only the actual initial conditions and the corresponding equations of motion are given. At the end of this subchapter, we present all four steps in sequence.
Impact Solution: Second Step
The description continues with a presentation of the second step. For reasons of clarity, we set k + 1 = 2; the initial state is defined by:
The initial conditions are:
The Heisenberg equations of motion of the expectation values of step 2 are:
Here, the damping process is turned on again at the beginning, with γ = 0.0005. Figure 6 shows the expected “ringing” effects now represented in phase portraits.
Figure 6. Phase portraits of the real parts of the second step of the following operators: . The value of the damping constant is γ = 0.0005. All trajectories vary after their starting points (marked by a cross in the first figure) until they end up at the fixed point 0. The coupling constants are g1 = g2 = 0.1.
Impact Solution: Third Step
The presentation continues with a description of the third step. For reasons of comprehensibility better readability, we set k + 2 = 3. The initial state is fixed by:
In this case, the initial conditions of all states are:
The Heisenberg equations of motion of the expectation values of step 3 are:
Attention is drawn to the symmetry that exists between all equations of motion of the first and third steps. Further, the initial condition for s†gp2 in step 3 is s†gp2(t2) = 0, while in step 1 s†gp2(t0) = 1 was valid. Similarly, the value s†ep2(t2) = 1 is assumed in step 3, whereas s†ep2(t0) = 0 was correct in step 1.
Impact Solution: Fourth Step
The presentation continues with a description of the fourth step. Here, we set +2 = 4, and the initial state is given by:
where | ϕ0 〉 defines the vacuum state. In this case, the initial conditions of the two robot states are:
The initial states of the substrate at the three positions p2, p3 and p4 are:
The heat-bath operator satisfies the initial condition B†p3(t3) = 1.
The relevant Heisenberg equations of motion of step 4 are:
Impact Solution: All Four Steps in Succession
Now all four partial solutions are collected and put together into one sequence. The first consecutive view describes the changes in the four robot states during the four separate walking steps. Figure 7 illustrates the corresponding real parts (imaginary parts are again very similar, therefore their representation is skipped) of the transitions of these walking states on a reduced scale of 5000 steps for better visibility. Clearly, two typical effects are observable. The similar steps 1 and 3 are regularly performed and go directly from their initial values to zero. Steps 2 and 4 also conform to this pattern, but one on a broader scale, because step 4 shows for slightly more oscillations than for in step 2.
Figure 7. Damped temporal curves representing the real parts of the four molecular robot operators that are consecutively performed by the steps 1–4. The following algebraic equivalents for the selected operators are introduced: . The value of the damping constant is set to γ = 0.0005. The coupling constants are g1 = g2 = 0.1.
Figure 8 summarizes the sequence of the real parts of the ground states of the substrate for all four steps. Here another effect comes to light. The behavior of s†gp3 is different in step 2, and during step 2 (see Figure 2), the ground state s†gp3 does not exist and has to be generated, whereas during step 3 this ground state already exists and must be annihilated and transformed into an exited state.
Figure 8. Damped temporal curves representing the real parts of three selected substrate operators (ground states) that are consecutively performed by steps 1–4. The following algebraic equivalents for the selected operators are introduced: gp2 = s†gp2, gp3 = s†gp3, gp4 = s†gp4. The value of the damping constant is set to γ = 0.0005. The coupling constants are g1 = g2 = 0.1.
From a purely mathematical point of view, we have to compare the two Equations (5.25) and (5.34) for s†gp3 in the second and third steps. In the second step, there is a strong mutual dependence between s†gp3 and s†gp1. In Equation (5.34), we mainly have to consider the dependence of s†gp3 and s†ep3. The operator product sep2s†gp2 is time-independent and self-consistent, while substrate states g and e produced themselves reciprocally (see also Figure 10). The three remaining bordering operators (2 r-operators, one B-operator) are equivalent in both expressions.
The next focus of attention is the behavior of the synchronizing heat-bath operators B†p2 and B†p3. Here, we consider at first the effect of the damping value γB = 0.001 and fix all other damping constants to zero. Figure 9 reveals the behavior of the overlay of the real parts (gray) and imaginary parts (red) of these two operators. Concerning the real parts, both operators are “strong” during steps 2 and 4 (the fueling process) and “weak” in the two other steps. During step 1, B†p2 starts with B†p2(t0) = 0 and ends with B†p2(t1) = 1. In step 2, it starts with this initial value, B†p2(t2) = 1, and ends permanently with B†p2(t2) = 0. Such interplay is valid in the same manner for B†p3(t3). The trajectories of the imaginary parts of these two operators are distinct in that way that the real parts of these two operators in steps 2 and 4 are no longer dominant.
Figure 9. Damped temporal curves showing the overlay of the real (gray) and imaginary (red) parts of the operators B†p2 and B†p3, which are consecutively outlined through steps 1–4. The algebraic equivalents of the selected operators are: Bp2 = B†p2, Bp3 = B†p3. The damping constant is set to γB = 0.001. The coupling constants are g1 = g2 = 0.1.
The situation greatly changes if we decrease the corresponding damping constant. Figure 10 portrays the temporal diagram of the overlay of the real and imaginary parts of B†p2 and B†p3. The contours of each step can be considered as an envelope around the oscillating solutions.
Figure 10. Damped temporal curves showing the overlay of the real (gray) and imaginary (red) parts of the operators B†p2 and B†p3, which are consecutively outlined through steps 1–4. The algebraic equivalents of the selected operators are Bp2 = B†p2 and Bp3 = B†p3. The value of the damping constant is γB = 0.0001. The coupling constants are g1 = g2 = 0.1.
The next two sections focus on wave solutions. Firstly, we combine modes that synchronize within a lane and between lanes (similar to light modes or laser modes, e.g., Haken, 1970). Secondly, we present a solution that constitutes a running wave.
Synchronized Wave-based Motion
Coherent Motion on Parallel Substrate Lanes
It is assumed that there are M different lanes, and on each lane up to L positions (lane length) are available. The different walking lanes of robots are distinguished by the index m = 1, …, M, and the position of a molecular robot on lane m is denoted by pm = 1, …, L. The B fields are also coupled in a lane and across different lanes. Therefore, we introduce a new two-dimensional vector l by combining both parameters l = (m, pm). This combination modifies the Equations (4.1–4.7) by an additional index, as e.g., demonstrated by the reformulation of Equation (4.1):
Further, we define another two-dimensional vector which defines a two-dimensional wave number vector. This second vector has been introduced since the operator B†l will be developed by a plane wave approximation:
The synchronization activities of the two heat-bath operators B†l, Bl on the molecular robots are described as external signal field modes that are coupled by the constant Jll′.
We insert this hypothesis in the “parallel” extension of Equation (4.7), indicating that our approach now applies to parallel M lanes. This therefore implies that l is fixed and over all paths we sum up m at the various positions with respect to the fixed position pm of l. The use of the coupling constant Jll′ for requests a summation along a lane and between lanes:
We multiply both sides of (5.48) with and use the orthogonal relation :
The last, rather complicated term can be simplified if we introduce the center of gravity x = l + l′/2 and the distance d = l − l′. Further, it is supposed that the coupling constant J depends only on the distance Jll′ = Jd and the lane length L is defined by = LxLd. It follows that:
where ωk has the dimension of a frequency. We continue with the simplification process and set
if all damping constants are set equal to
Finally, after all these simplifications the resulting formula is obtained:
The basic solution to Equation (3.67) can be achieved if only the term k′ = 0 that corresponds to the mode with an infinite wavelength is kept. We will now solve this “reduced” equation.
The last term, , can be deleted in the interaction representation and the expectation value of can be ignored. These two modifications lead to the following equation:
Comparing the result of Equation (5.56) with the original expression (4.7), two main differences are immediately observable. At first, the two coupling constants, g1 and g2, have to be replaced by and . Second, there is a summation over all lanes m = 1, …, M. The elaborate calculations of this subchapter show that the synchronization of motion can easily be accomplished (at least for by replacing the coupling constants and performing an addition of all lanes at the same positions.
Running Wave Solutions
We again concentrate on a single robot that moves on a lane and restrict ourselves for simplicity to step 1 since the remaining steps can be handled very similarly. Furthermore, the damping constants with respect to the molecular robot are equaled: γab = γba = γ, the damping constants of the substrate are inoperative: γg = γe = 0, and the damping constant of the B-field is assumed to be effectively γB ≠ 0. The positions on the lane are pk = k, pk + 1 = k + 1, … Under these assumptions, the approach that guides us to a running wave solution is formulated as follows:
where K is a real variable that will be later defined.
Insertion of these expressions into the Equations (5.4) and (5.5) provides the following formulas:
In doing so, we introduced the following abbreviations:
The corresponding bosonic operators (and Hermitean adjoint operators) commutate, thus they can be interchanged.
Figure 11 shows the -like peaks of these operator products at the same time slots. At the same time, e.g., the ground state g will be annihilated and the exited state e created. The periodicity of this self-consistent process is due to its calculation without damping (revealing again the original mathematical symmetry). In a more realistic case, damping effects are activated and the curves converge after the first peak to zero. Therefore, for a fixed time t, both the real and imaginary parts of D and/or D† are constant expressions.
Figure 11. Undamped temporal curves, showing the real parts of the operator products D1, D†1 (abridged to D1k), D2, D†2 (abridged to D2k), D, D† (abridged to Dk) and B†p2Bp2, consecutively outlined for step 1. The damping constant is γ = 0, the coupling constants are constant: g1 = g2 = 0.1.
In the next step, we reformulate the Equations (5.6–5.10) by including the D–terms:
The calculations are continued by substituting the three operators R†ab, R†c1 and B†b by
The derivatives of the first two expressions read as follows:
Insertion of in (5.64) provides the expression
where = 0 and both φ and G are real.
By inserting in (5.65), the following expression is obtained:
It follows, by insertion of B†b, given by Equation (5.71) into expression (5.70), that
To reproduce the r.h.s. of these three Equations (5.73–5.75), the following assumption is pursued:
With this ansatz the derivatives of the three operators ρ†ab, ρ†c1, ρ†b can be cast in the following form:
From expression (5.83), the value of G is concluded as:
Both sides of this expression should be real, therefore the exponents must be zero: φ = − K/(k + 1) and ρ0Re(D) = G. A second solution is formed by the same φ but with . This last result comes from the two Equations (5.74) and (5.75).
Due to the determination of the two variables φ and K, running wave solutions of the molecular “leg-over-leg” walking can be expressed for the three operators , B†p2 in the conventional form (k = 1, 2,..; G = ρ0Re(D) = const.):
Such a solution is not surprising because it is well known that a composition of different modes (more than only the ground mode k = 0) can generate running waves (see Section Coherent Motion on Parallel Substrate Lanes).
The central goal of this contribution was to support the hypothesis that quantum mechanical effects in molecular biology—especially in human brains—cannot be neglected. This assumption is mainly justified by the size of the interacting objects (nano size) and by the experimental results showing that even greater molecules up to about 100 atoms can demonstrate quantum behavior, e.g., in double-slit experiments (Haken and Levi, 2012). The nano size argument is further supported by the established experimental technology to construct molecular machines and walkers from DNA.
The most striking results to emerge from the produced solutions is that neural processes like the motion of molecular robots (kinesin and dynein) in axons and dendrites can be modeled by an approach that is particle-based (impact solution) or purely wave-oriented (synchronized modes solution and running wave solution). In addition, wave-based solutions also seem to be very relevant for the interactions between different neural layers.
We are aware that our predictions suffer from limitations and have to be considered as a first, basic approach due to the following three reasons. Firstly, the number of the many parameters (in total 16) was reduced to four (two coupling constants and two damping constants). Secondly, there are still not enough experimental data to fix the values of the last mentioned four parameters or even the great set of all relevant available biological data. Thirdly, a greater number of experimentally approved parameters could lead to a higher generalization of the achieved results.
A first extension of our approach could be the integration of tunneling effects, the inclusion of fluctuation forces, and the consideration of the cargo transport. All these additional considerations would greatly increase the predictive power of the outlined model. In this way, it would be feasible to describe e.g., the walk of a molecular robot along the DNA and the transcription of a nucleotide into RNA, whereby for each step the effective potential barrier must be tunneled.
Conflict of Interest Statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The research leading to these results has received funding from the European Union Seventh Frame Program (FP7/2007-2013) under grant agreement no. 604102 (Human Brain Project, Neurorobotics Platform, SP10). We express also our gratitude to Dr. B. Schenke for his numerical calculation support of this contribution.
Holmes, P. H., Lumley, J. L., Berkooz, G., and Rowley, C. W. (2012). Turbulence, Coherent Structures, Dynamical Systems and Symmetry, Cambridge Monograph on Mechanics. Cambridge: Cambridge University Press.
Levi, P., and Haken, H. (2010). “Towards a synergetic quantum field theory for evolutionary symbiotic multi-robotics,” in Symbiotic Multi–robot Organisms, eds P. Levi and S. Kernbach (Berlin: Springer), 25–54.
Keywords: molecular robotics, biped walking on microtubules, motion of kinesin-1 and dynein alog axons and dendrites, particles and wave solutions, neuro-robotics
Citation: Levi P (2015) Molecular quantum robotics: particle and wave solutions, illustrated by “leg-over-leg” walking along microtubules. Front. Neurorobot. 9:2. doi: 10.3389/fnbot.2015.00002
Received: 19 September 2014; Accepted: 23 April 2015;
Published: 08 May 2015.
Edited by:Marc-Oliver Gewaltig, Ecole Polytechnique Federale de Lausanne, Switzerland
Reviewed by:Hiroaki Wagatsuma, Kyushu Institute of Technology, Japan
Valery E. Karpov, National Research University Higher School of Economics, Russia
Copyright © 2015 Levi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Paul Levi, Forschungszentrum Informatik (Centre of Computer Science), Intelligent System and Production Engineering, Interactive Diagnosis and Service Systems, Haid-und-Neu-Street 10-14, 76131 Karlsruhe, Baden-Württemberg, Germany, firstname.lastname@example.org