String-like propagation of the 5-coordinated defect state in supercooled water: molecular origin of dynamic and thermodynamic anomalies

We find that at low temperature water, large amplitude (similar to 60 degrees) rotational jumps propagate like a string, with the length of propagation increasing with lowering temperature. The strings are formed by mobile 5-coordinated water molecules which move like a Glarum defect (J. Chem. Phys., 1960, 33, 1371), causing water molecules on the path to change from 4-coordinated to 5-coordinated and again back to 4-coordinated water, and in the process cause the tagged water molecule to jump, by following essentially the Laage-Hynes mechanism (Science, 2006, 311, 832-835). The effects on relaxation of the propagating defect causing large amplitude jumps are manifested most dramatically in the mean square displacement (MSD) and also in the rotational time correlation function of the O-H bond of the molecule that is visited by the defect (transient transition to the 5-coordinated state). The MSD and the decay of rotational time correlation function, both remain quenched in the absence of any visit by the defect, as postulated by Glarum long time ago. We establish a direct connection between these propagating events and the known thermodynamic and dynamic anomalies in supercooled water. These strings are found largely in the regions that surround the relatively rigid domains of 4-coordinated water molecules. The propagating strings give rise to a noticeable dynamical heterogeneity, quantified here by a sharp rise in the peak of the four-point density response function, chi(4)(t). This dynamics heterogeneity is also responsible for the breakdown of the Stokes-Einstein relation.

Redundancy resolution using a tractrix and its application to real-time simulations of hyper-redundant manipulators, snakes and tying of knots

To realistically simulate the motion of flexible objects such as ropes, strings, snakes, or human hair,one strategy is to discretise the object into a large number of small rigid links connected by rotary or spherical joints. The discretised system is highly redundant and the rotations at the joints (or the motion of the other links) for a desired Cartesian motion of the end of a link cannot be solved uniquely. In this paper, we propose a novel strategy to resolve the redundancy in such hyper-redundant systems.We make use of the classical tractrix curve and its attractive features. For a desired Cartesian motion of the `head'of a link, the `tail' of the link is moved according to a tractrix,and recursively all links of the discretised objects are moved along different tractrix curves. We show that the use of a tractrix curve leads to a more `natural' motion of the entire object since the motion is distributed uniformly along the entire object with the displacements tending to diminish from the `head' to the `tail'. We also show that the computation of the motion of the links can be done in real time since it involves evaluation of simple algebraic, trigonometric and hyperbolic functions. The strategy is illustrated by simulations of a snake, tying of knots with a rope and a solution of the inverse kinematics of a planar hyper-redundant manipulator.

Potential and electric field distribution in a ceramic disc insulator string with faulty insulators

Ceramic/Porcelain suspension disc insulators are widely used in power systems to provide electrical insulation and mechanically support for high-voltage transmission lines. These insulators are subjected to a variety of stresses, including mechanical, electrical and environmental. These stresses act in unison. The exact nature and magnitude of these stresses vary significantly and depends on insulator design, application and its location. Due to various reasons the insulator disc can lose its electrical insulation properties without any noticeable mechanical failure. Such a condition while difficult to recognize, can enhance the stress on remaining healthy insulator discs in the string further may lead to a flashover. To understand the stress enhancement due to faulty discs in a string, attempt has been made to simulate the potential and electric field profiles for various disc insulators presently used in the country. The results of potential and electric filed stress obtained for normal and strings with faulty insulator discs are presented.

Simulation of potential and electric field across faulty ceramic disc insulator string

Ceramic/Porcelain insulators are widely used in power transmission lines to provide mechanical support for High voltage conductors in addition to withstand electrical stresses. As a result of lightning, switching or temporary over voltages that could initiate flashover under worst weather conditions, and to operate within interference limits. Given that the useful life in service of the individual insulator elements making up the insulator strings is hard to predict, they must be verified periodically to ensure that adequate line reliability is maintained at all times. Over the years utilities have adopted few methods to detect defective discs in a string, subsequently replacement of the faulty discs are being carried out for smooth operation. But, if the insulator is found to be defective in a string at some location that may not create any changes in the field configuration, there is no need to replace to avoid manpower and cost of replacement. Due to deficiency of electric field data for the existing string configuration, utilities are forced to replace the discs which may not be essentially required. Hence, effort is made in the present work to simulate the potential and electric field along the normal and with faults induced discs in a string up to 765 kV system voltages using Surface Charge Simulation Method (SCSM). A comparison is made between simulated results, experimental and field data and it was found that the computed results are quite acceptable and useful.

Correlation between local structural dynamics of proteins inferred from NMR ensembles and evolutionary dynamics of homologues of known structure

Conformational changes in proteins are extremely important for their biochemical functions. Correlation between inherent conformational variations in a protein and conformational differences in its homologues of known structure is still unclear. In this study, we have used a structural alphabet called Protein Blocks (PBs). PBs are used to perform abstraction of protein 3-D structures into a 1-D strings of 16 alphabets (a-p) based on dihedral angles of overlapping pentapeptides. We have analyzed the variations in local conformations in terms of PBs represented in the ensembles of 801 protein structures determined using NMR spectroscopy. In the analysis of concatenated data over all the residues in all the NMR ensembles, we observe that the overall nature of inherent local structural variations in NMR ensembles is similar to the nature of local structural differences in homologous proteins with a high correlation coefficient of .94. High correlation at the alignment positions corresponding to helical and beta-sheet regions is only expected. However, the correlation coefficient by considering only the loop regions is also quite high (.91). Surprisingly, segregated position-wise analysis shows that this high correlation does not hold true to loop regions at the structurally equivalent positions in NMR ensembles and their homologues of known structure. This suggests that the general nature of local structural changes is unique; however most of the local structural variations in loop regions of NMR ensembles do not correlate to their local structural differences at structurally equivalent positions in homologues.

Low tension strings on a cosmological singularity

It has recently been argued that the singularity of the Milne orbifold can be resolved in higher spin theories. In string theory scattering amplitudes, however, the Milne singularity gives rise to ultraviolet divergences that signal uncontrolled backreaction. Since string theory in the low tension limit is expected to be a higher spin theory (although precise proposals only exist in special cases), we investigate what happens to these scattering amplitudes in the low tension limit. We point out that the known problematic ultraviolet divergences disappear in this limit. In addition we systematically identify all divergences of the simplest 2-to-2 tree-level string scattering amplitude on the Milne orbifold, and argue that the divergences that survive in the low tension limit have sensible infrared interpretations.

Strings vs. spins on the null orbifold

We study the null orbifold singularity in 2+1 d flat space higher spin theory as well as string theory. Using the Chern-Simons formulation of 2+1 d Einstein gravity, we first observe that despite the singular nature of this geometry, the eigenvalues of its Chern-Simons holonomy are trivial. Next, we construct a resolution of the singularity in higher spin theory: a Kundt spacetime with vanishing scalar curvature invariants. We also point out that the UV divergences previously observed in the 2-to-2 tachyon tree level string amplitude on the null orbifold do not arise in the at alpha' -> infinity limit. We find all the divergences of the amplitude and demonstrate that the ones remaining in the tensionless limit are physical IR-type divergences. We conclude with a discussion on the meaning and limitations of higher spin (cosmological) singularity resolution and its potential connection to string theory.

Heterotic string on the CHL orbifold of K3

We study N = 2 compactifications of heterotic string theory on the CHL orbifold (K3 x T-2)/Z(N) with N = 2, 3, 5, 7. Z(N) acts as an automorphism on K3 together with a shift of 1/N along one of the circles of T-2. These compactifications generalize the example of the heterotic string on K3 x T-2 studied in the context of dualities in string theories. We evaluate the new supersymmetric index for these theories and show that their expansion can be written in terms of the McKay-Thompson series associated with the Z(N) automorphism embedded in the Mathieu group M-24. We then evaluate the difference in one-loop threshold corrections to the non-Abelian gauge couplings with Wilson lines and show that their moduli dependence is captured by Siegel modular forms related to dyon partition functions of N = 4 string theories.

Heterotic string on the CHL orbifold of K3

We study N = 2 compactifications of heterotic string theory on the CHL orbifold (K3 x T-2)/Z(N) with N = 2, 3, 5, 7. Z(N) acts as an automorphism on K3 together with a shift of 1/N along one of the circles of T-2. These compactifications generalize the example of the heterotic string on K3 x T-2 studied in the context of dualities in string theories. We evaluate the new supersymmetric index for these theories and show that their expansion can be written in terms of the McKay-Thompson series associated with the Z(N) automorphism embedded in the Mathieu group M-24. We then evaluate the difference in one-loop threshold corrections to the non-Abelian gauge couplings with Wilson lines and show that their moduli dependence is captured by Siegel modular forms related to dyon partition functions of N = 4 string theories.

Synthetic molecular machines for active self-assembly : prototype algorithms, designs, and experimental study

Computer science and electrical engineering have been the great success story of the twentieth century. The neat modularity and mapping of a language onto circuits has led to robots on Mars, desktop computers and smartphones. But these devices are not yet able to do some of the things that life takes for granted: repair a scratch, reproduce, regenerate, or grow exponentially fast–all while remaining functional.

This thesis explores and develops algorithms, molecular implementations, and theoretical proofs in the context of “active self-assembly” of molecular systems. The long-term vision of active self-assembly is the theoretical and physical implementation of materials that are composed of reconfigurable units with the programmability and adaptability of biology’s numerous molecular machines. En route to this goal, we must first find a way to overcome the memory limitations of molecular systems, and to discover the limits of complexity that can be achieved with individual molecules.

One of the main thrusts in molecular programming is to use computer science as a tool for figuring out what can be achieved. While molecular systems that are Turing-complete have been demonstrated [Winfree, 1996], these systems still cannot achieve some of the feats biology has achieved.

One might think that because a system is Turing-complete, capable of computing “anything,” that it can do any arbitrary task. But while it can simulate any digital computational problem, there are many behaviors that are not “computations” in a classical sense, and cannot be directly implemented. Examples include exponential growth and molecular motion relative to a surface.

Passive self-assembly systems cannot implement these behaviors because (a) molecular motion relative to a surface requires a source of fuel that is external to the system, and (b) passive systems are too slow to assemble exponentially-fast-growing structures. We call these behaviors “energetically incomplete” programmable behaviors. This class of behaviors includes any behavior where a passive physical system simply does not have enough physical energy to perform the specified tasks in the requisite amount of time.

As we will demonstrate and prove, a sufficiently expressive implementation of an “active” molecular self-assembly approach can achieve these behaviors. Using an external source of fuel solves part of the the problem, so the system is not “energetically incomplete.” But the programmable system also needs to have sufficient expressive power to achieve the specified behaviors. Perhaps surprisingly, some of these systems do not even require Turing completeness to be sufficiently expressive.

Building on a large variety of work by other scientists in the fields of DNA nanotechnology, chemistry and reconfigurable robotics, this thesis introduces several research contributions in the context of active self-assembly.

We show that simple primitives such as insertion and deletion are able to generate complex and interesting results such as the growth of a linear polymer in logarithmic time and the ability of a linear polymer to treadmill. To this end we developed a formal model for active-self assembly that is directly implementable with DNA molecules. We show that this model is computationally equivalent to a machine capable of producing strings that are stronger than regular languages and, at most, as strong as context-free grammars. This is a great advance in the theory of active self- assembly as prior models were either entirely theoretical or only implementable in the context of macro-scale robotics.

We developed a chain reaction method for the autonomous exponential growth of a linear DNA polymer. Our method is based on the insertion of molecules into the assembly, which generates two new insertion sites for every initial one employed. The building of a line in logarithmic time is a first step toward building a shape in logarithmic time. We demonstrate the first construction of a synthetic linear polymer that grows exponentially fast via insertion. We show that monomer molecules are converted into the polymer in logarithmic time via spectrofluorimetry and gel electrophoresis experiments. We also demonstrate the division of these polymers via the addition of a single DNA complex that competes with the insertion mechanism. This shows the growth of a population of polymers in logarithmic time. We characterize the DNA insertion mechanism that we utilize in Chapter 4. We experimentally demonstrate that we can control the kinetics of this re- action over at least seven orders of magnitude, by programming the sequences of DNA that initiate the reaction.

In addition, we review co-authored work on programming molecular robots using prescriptive landscapes of DNA origami; this was the first microscopic demonstration of programming a molec- ular robot to walk on a 2-dimensional surface. We developed a snapshot method for imaging these random walking molecular robots and a CAPTCHA-like analysis method for difficult-to-interpret imaging data.

Evaluation of the capture efficiency and size selectivity of four pot types in the prospective fishery for North Pacific giant octopus (Enteroctopus dofleini)

Over 230 metric tons of octopus is harvested as bycatch annually in Alaskan trawl, long-line, and pot fisheries. An expanding market has fostered interest in the development of a directed fishery for North Pacific giant octopus (Enteroctopus dofleini). To investigate the potential for fishery development we examined the efficacy of four different pot types for capture of this species. During two surveys in Kachemak Bay, Alaska, strings of 16 –20 sablefish, Korean hair crab, shrimp, and Kodiak wooden lair pots were set at depths ranging between 62 and 390 meters. Catch per-unit-of-ef for t estimates were highest for sablefish and lair pots. Sablefish pots caught significantly heavier North Pacific giant octopuses but also produced the highest bycatch of commercially important species, such as halibut (Hippoglossus stenolepis), Pacific cod (Gadus macrocephalus), and Tanner crab (Chionoecetes bairdi).

Dividing population genetic distance data with the software Partitioning Optimization with Restricted Growth Strings (PORGS): an application for Chinook salmon (Oncorhynchus tshawytscha), Vancouver Island, British Columbia

A new method of finding the optimal group membership and number of groupings to partition population genetic distance data is presented. The software program Partitioning Optimization with Restricted Growth Strings (PORGS), visits all possible set partitions and deems acceptable partitions to be those that reduce mean intracluster distance. The optimal number of groups is determined with the gap statistic which compares PORGS results with a reference distribution. The PORGS method was validated by a simulated data set with a known distribution. For efficiency, where values of n were larger, restricted growth strings (RGS) were used to bipartition populations during a nested search (bi-PORGS). Bi-PORGS was applied to a set of genetic data from 18 Chinook salmon (Oncorhynchus tshawytscha) populations from the west coast of Vancouver Island. The optimal grouping of these populations corresponded to four geographic locations: 1) Quatsino Sound, 2) Nootka Sound, 3) Clayoquot +Barkley sounds, and 4) southwest Vancouver Island. However, assignment of populations to groups did not strictly reflect the geographical divisions; fish of Barkley Sound origin that had strayed into the Gold River and close genetic similarity between transferred and donor populations meant groupings crossed geographic boundaries. Overall, stock structure determined by this partitioning method was similar to that determined by the unweighted pair-group method with arithmetic averages (UPGMA), an agglomerative clustering algorithm.

Calibrating Analytical Models for Semilocal Strings

In this work we calibrate two different analytic models of semilocal strings by constraining the values of their free parameters. In order to do so, we use data obtained from the largest and most accurate field theory simulations of semilocal strings to date, and compare several key properties with the predictions of the models. As this is still work in progress, we present some preliminary results together with descriptions of the methodology we are using in the characterisation of semilocal string networks.

Torsional Behaviour of Cello Strings

The behaviour of a bowed string depends, among other things, on the frequency, impedance and internal damping of torsional waves on the string. Very little published information is available about these quantities, especially the torsional damping. Measurements of all relevant torsional properties have been made on cello strings of three different constructions. These show that the torsional modes are harmonically spaced to reasonable accuracy, and that the Q factors are approximately equal for all modes of a given string. These torsional Q factors are roughly an order of magnitude smaller than those of the transverse modes of the same string. The torsional wave speed varies somewhat with the tension in the string, decreasing with higher tension. The damping factors are not significantly influenced by tension. These results have been expressed in terms of a novel "reflection function" [1] suitable for direct incorporation into simulations of the bowing process.