Tube bundles in physics

Structures which commonly occur in both man-made applications and nature are bundles of inter-wound tubular shapes. We will all familiar with ropes, used as a device to climb or pull and in addition large steel ropes are a common engineering structure used to stabilise bridges and oil-rigs. Tube bundles also proliferate in the human body. For example, the famed DNA double helixEruption of a coronal flux rope can be thought of as a pair of tubular electron densities, as depicted in the major Figure on my Homepage. More complex bundles from fibred molecules such as Fibronin, Keratin and the Optic nerve (amongst many others). We also see rope-like structures in the Sun, where bundles of magnetic flux are expelled from the Sun's interior to its surface, they then contort in its atmosphere (see Figure 1 to the right) and are often then expelled explosively (a coronal mass ejection). One school of thought holds that this ejection is vital to the so-called dynamo process, which is theorised to generate the Sunís large scale magnetic field. I have and am currently involved with the mathematical modelling of bundle structures in all these areas and have particular interest in their response to applied forces and moments, as well as the analysing the complexity of structure they can exhibit



Mathematical aspects of tube bundles

plectonemeA key consideration in the mathematical modelling of tubes and tube bundles is their changing shape when they are bent, stretched, twisted and folded.  For example, one can attempt to tighten a straightened rope by twisting it (see Figure 2). However, no matter how hard you try it is impossible to tighten a rope beyond a certain limit (try it if you have a rope!). Eventually the rope would have to pass through itself in order to wind any further, a physical iman
        tightening ropepossibility. On the other hand we often have to wrap and knot ropes, once again we use the fact that the rope cannot pass through itself in order to allow the knot to catch. A classic example of this is the formation of a plectoneme, as shown in panels A-E to the left (Figures courtesy of Sebastien Neukirch). The equations used to define the systems I discussed above do not generally ensure this necessary physical constraint is actually enforced.  In order to do so we have to define mathematical measurements of shape which remain unchanged as long as the tubes do not pass through themselves, but change when they do. Ensuring this quantity always stays fixed for our system allows the mathematician to enforce this necessary physical constraint. This is a particular branch of mathematics called Topology (knot and braid theory in this case). The problem is that it is hard to define examples of these measurements which are simple to apply in practice for the mathematical modeller. It is this topic which was the subject of my doctoral work, which concerned the creation of a new topological measurement of the linking of two open ended curves (published in J. Phys. A). This lead to the creation of an open analogue of the celebrated Calugareanu theorem (based on an open curve definition of the writhing). Recently I extended this open topological framework to accounted form knotting and belt trick type deformations, in joint work with Sebastien Neukirch (accepted in J Phys A, a pre-print is here). The code to calculate the writhing of a twisted ribbon/tube can be found in the code section of this website.


Practical Application of Techniques

Biopolymers.
One of the most common methods currently used in identifying the structure of the hundreds of thousands of protein structures which exist in nature are diffraction techniques. Put roughly, electromagnetic radiation (such as x-rays) is fired at a sample of the protein of interest and the radiation is scattered by its atomic structure. The pattern formed after scattering is thus indicative of the protein's shape and the theory gives us the necessary mathematical tool to translate this data, the Fourier transform.  If we can guess a structure whose Fourier transform fits this pattern then we probably have the right structure (in reality itís a little more sophisticated than just guessing!).  Unfortunately proteins can have spectacularly complex structure (just look at haemoglobin). Currently, defining the transform relies on a description of the molecule which depends on defining the hundreds of thousands to millions of orientation angles between the neighbouring molecules composing the structure. For medium resolution diffraction techniques, which are used for the more complex protein structures, this level of detail is unnecessary in order to determine its tertiary and quaternary structure (i.e., its inter-wound shape rather than the specific molecular sequence). With Alain Goriely and Andrew Hausrath I have developed a technique, based on modelling the protein as a continuum tubular (electron) density of arbitrary structural complexity, which can be used to describe protein structures with a vastly reduced number of parameters, but still having sufficient accuracy for the purpose of determining its shape (published in J. Phys. A). This work is now being used to develop efficient methods for identifying new protein structures with Dr Ehmke Pohl.

Elastic tubes
A common method of modelling physical tubes is elastic rod theory. Basically we describe a tube which bends stretches and twists under applied forces (compression and stretching) and moments (twisting and bending forces). The aim of this is to work out what shape such a tube will adopt under the application of a set of forces/moments. I am involved with modelling not single tubes but an ensemble of tubes which are inter-wound, a more complex task which involves some of the mathematical methods mentioned above. The model was used to produce the bent rope-like structure on my homepage (under an applied moment).  This work will in future be used to model both biopolymer structures and wire ropes in specific physical scenarios.

Magnetic flux ropes in the Sunís corona
sigmoidsThe Corona is the Sunís atmosphere directly above its surface. Rope-like structures composed of inter-wound tubes of magnetic field (magnetic flux) emerge though the surface into the corona. The coronal environment is very different from the interior of the Sun and the flux ropes begins to contort adopting coiled structures which appear as S or Z shapes from above (see x-ray images in Figure 4 and model simulations in Figure 5). Often these flux ropes will become unstable and be ejected into space. There is much interest in this process partly because it is twisted tubesvital to understanding what happens in the Sun's interior which, unlike the Corona, we cannot view directly using current observational techniques. In particular there is interest in measuring its helicity, a topological measure of the inter-winding of the flux tubes composing the rope. The results from my thesis  (detailed in J. Phys. A) gave the correct method for measuring the helicity in this case for which the flux ropes are anchored at points on the Sunís surface, rather than forming closed loops. Previously there had only been a measure of helicity of closed loop flux bundles. With Mitchell Berger I showed that this measurement actually indicated certain assumptions about how the helicity could be measured from observational data would not in general be correct (published in Solar Physics).  A future aim will be to use this work to improve the techniques for estimating the amount of helicity of these flux ropes.

Recently, on a more technical note, Anthony Yeates and I have obtained some fundamental results regarding the interpretation of helicity of fields which allow flow through their boundary (Pre-Print). The helicity of such fields seems to be ill-defined as is is not Gauge invariant, however, we show that this is a function of mistake of defining the helicity though its vector potential, rather than basing its definition on the curves composing the field itself.