Triangle Evolution-A Hybrid Heuristic for Global Optimization

Abstract This paper presents a hybrid heuristic{triangle evolution (TE) for global optimization. It is a real coded evolutionary algorithm. As in di®erential evolution (DE), TE targets each individual in current population and attempts to replace it by a new better individual. However, the way of generating new individuals is di®erent. TE generates new individuals in a Nelder- Mead way, while the simplices used in TE is 1 or 2 dimensional. The proposed algorithm is very easy to use and e±cient for global optimization problems with continuous variables. Moreover, it requires only one (explicit) control parameter. Numerical results show that the new algorithm is comparable with DE for low dimensional problems but it outperforms DE for high dimensional problems.

Attenuation tomography. Modelling regional love waves: Imperial Valley to Pasadena

Abstract to Part I

The inverse problem of seismic wave attenuation is solved by an iterative back-projection method. The seismic wave quality factor, Q, can be estimated approximately by inverting the S-to-P amplitude ratios. Effects of various uncertain ties in the method are tested and the attenuation tomography is shown to be useful in solving for the spatial variations in attenuation structure and in estimating the effective seismic quality factor of attenuating anomalies.

Back-projection attenuation tomography is applied to two cases in southern California: Imperial Valley and the Coso-Indian Wells region. In the Coso-Indian Wells region, a highly attenuating body (S-wave quality factor (Q_β ≈ 30) coincides with a slow P-wave anomaly mapped by Walck and Clayton (1987). This coincidence suggests the presence of a magmatic or hydrothermal body 3 to 5 km deep in the Indian Wells region. In the Imperial Valley, slow P-wave travel-time anomalies and highly attenuating S-wave anomalies were found in the Brawley seismic zone at a depth of 8 to 12 km. The effective S-wave quality factor is very low (Q_β ≈ 20) and the P-wave velocity is 10% slower than the surrounding areas. These results suggest either magmatic or hydrothermal intrusions, or fractures at depth, possibly related to active shear in the Brawley seismic zone.

No-block inversion is a generalized tomographic method utilizing the continuous form of an inverse problem. The inverse problem of attenuation can be posed in a continuous form , and the no-block inversion technique is applied to the same data set used in the back-projection tomography. A relatively small data set with little redundancy enables us to apply both techniques to a similar degree of resolution. The results obtained by the two methods are very similar. By applying the two methods to the same data set, formal errors and resolution can be directly computed for the final model, and the objectivity of the final result can be enhanced.

Both methods of attenuation tomography are applied to a data set of local earthquakes in Kilauea, Hawaii, to solve for the attenuation structure under Kilauea and the East Rift Zone. The shallow Kilauea magma chamber, East Rift Zone and the Mauna Loa magma chamber are delineated as attenuating anomalies. Detailed inversion reveals shallow secondary magma reservoirs at Mauna Ulu and Puu Oo, the present sites of volcanic eruptions. The Hilina Fault zone is highly attenuating, dominating the attenuating anomalies at shallow depths. The magma conduit system along the summit and the East Rift Zone of Kilauea shows up as a continuous supply channel extending down to a depth of approximately 6 km. The Southwest Rift Zone, on the other hand, is not delineated by attenuating anomalies, except at a depth of 8-12 km, where an attenuating anomaly is imaged west of Puu Kou. The Ylauna Loa chamber is seated at a deeper level (about 6-10 km) than the Kilauea magma chamber. Resolution in the Mauna Loa area is not as good as in the Kilauea area, and there is a trade-off between the depth extent of the magma chamber imaged under Mauna Loa and the error that is due to poor ray coverage. Kilauea magma chamber, on the other hand, is well resolved, according to a resolution test done at the location of the magma chamber.

Abstract to Part II

Long period seismograms recorded at Pasadena of earthquakes occurring along a profile to Imperial Valley are studied in terms of source phenomena (e.g., source mechanisms and depths) versus path effects. Some of the events have known source parameters, determined by teleseismic or near-field studies, and are used as master events in a forward modeling exercise to derive the Green's functions (SH displacements at Pasadena that are due to a pure strike-slip or dip-slip mechanism) that describe the propagation effects along the profile. Both timing and waveforms of records are matched by synthetics calculated from 2-dimensional velocity models. The best 2-dimensional section begins at Imperial Valley with a thin crust containing the basin structure and thickens towards Pasadena. The detailed nature of the transition zone at the base of the crust controls the early arriving shorter periods (strong motions), while the edge of the basin controls the scattered longer period surface waves. From the waveform characteristics alone, shallow events in the basin are easily distinguished from deep events, and the amount of strike-slip versus dip-slip motion is also easily determined. Those events rupturing the sediments, such as the 1979 Imperial Valley earthquake, can be recognized easily by a late-arriving scattered Love wave that has been delayed by the very slow path across the shallow valley structure.

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.

飞秒强激光经氩气和氩团簇的传播

通过数值求解三维电场传播方程，理论模拟了飞秒强激光脉冲（50fs，10^16w／cm^2）在氩气和中等尺寸氩团簇中的传播效应．结果表明，飞秒强激光脉冲经氩气传播将发生频谱蓝移展宽和光束发散；而经中等尺寸氩团簇传播，则存在一定程度的自聚焦效应．

Further acceleration of MeV electrons by a relativistic laser pulse

With the development of photocathode rf electron gun, electrons with high-brightness and mono-energy can be obtained easily. By numerically solving the relativistic equations of motion of an electron generated from this facility in laser fields modelled by a circular polarized Gaussian laser pulse, we find the electron can obtain high energy gain from the laser pulse. The corresponding acceleration distance for this electron driven by the ascending part of the laser pulse is much longer than the Rayleigh length, and the light amplitude experienced on the electron is very weak when the laser pulse overtakes the electron. The electron is accelerated effectively and the deceleration can be neglected. For intensities around 10(19) W(.)mu m(2)/cm(2), an electron's energy gain near 0.1 GeV can be realized when its initial energy is 4.5 MeV, and the final velocity of the energetic electron is parallel with the propagation axis. The energy gain can be up to 1 GeV if the intensity is about 10(21) W(.)mu m(2)/cm(2). The final energy gain of the electron as a function of its initial conditions and the parameters of the laser beam has also been discussed.

Ultraviolet-infrared femtosecond laser-induced damage in fused silica and CaF2 crystals

The damage in fused silica and CaF2 crystals induced by wavelength tunable femtosecond lasers is studied. The threshold fluence is observed to increase rapidly with laser wavelength lambda in the region of 250-800 nm, while it is nearly a constant for 800

Signatures of Topological Superconductors

Topological superconductors are particularly interesting in light of the active ongoing experimental efforts for realizing exotic physics such as Majorana zero modes. These systems have excitations with non-Abelian exchange statistics, which provides a path towards topological quantum information processing. Intrinsic topological superconductors are quite rare in nature. However, one can engineer topological superconductivity by inducing effective p-wave pairing in materials which can be grown in the laboratory. One possibility is to induce the proximity effect in topological insulators; another is to use hybrid structures of superconductors and semiconductors.

The proposal of interfacing s-wave superconductors with quantum spin Hall systems provides a promising route to engineered topological superconductivity. Given the exciting recent progress on the fabrication side, identifying experiments that definitively expose the topological superconducting phase (and clearly distinguish it from a trivial state) raises an increasingly important problem. With this goal in mind, we proposed a detection scheme to get an unambiguous signature of topological superconductivity, even in the presence of ordinarily detrimental effects such as thermal fluctuations and quasiparticle poisoning. We considered a Josephson junction built on top of a quantum spin Hall material. This system allows the proximity effect to turn edge states in effective topological superconductors. Such a setup is promising because experimentalists have demonstrated that supercurrents indeed flow through quantum spin Hall edges. To demonstrate the topological nature of the superconducting quantum spin Hall edges, theorists have proposed examining the periodicity of Josephson currents respect to the phase across a Josephson junction. The periodicity of tunneling currents of ground states in a topological superconductor Josephson junction is double that of a conventional Josephson junction. In practice, this modification of periodicity is extremely difficult to observe because noise sources, such as quasiparticle poisoning, wash out the signature of topological superconductors. For this reason, We propose a new, relatively simple DC measurement that can compellingly reveal topological superconductivity in such quantum spin Hall/superconductor heterostructures. More specifically, We develop a general framework for capturing the junction's current-voltage characteristics as a function of applied magnetic flux. Our analysis reveals sharp signatures of topological superconductivity in the field-dependent critical current. These signatures include the presence of multiple critical currents and a non-vanishing critical current for all magnetic field strengths as a reliable identification scheme for topological superconductivity.

This system becomes more interesting as interactions between electrons are involved. By modeling edge states as a Luttinger liquid, we find conductance provides universal signatures to distinguish between normal and topological superconductors. More specifically, we use renormalization group methods to extract universal transport characteristics of superconductor/quantum spin Hall heterostructures where the native edge states serve as a lead. Interestingly, arbitrarily weak interactions induce qualitative changes in the behavior relative to the free-fermion limit, leading to a sharp dichotomy in conductance for the trivial (narrow superconductor) and topological (wide superconductor) cases. Furthermore, we find that strong interactions can in principle induce parafermion excitations at a superconductor/quantum spin Hall junction.

As we identify the existence of topological superconductor, we can take a step further. One can use topological superconductor for realizing Majorana modes by breaking time reversal symmetry. An advantage of 2D topological insulator is that networks required for braiding Majoranas along the edge channels can be obtained by adjoining 2D topological insulator to form corner junctions. Physically cutting quantum wells for this purpose, however, presents technical challenges. For this reason, I propose a more accessible means of forming networks that rely on dynamically manipulating the location of edge states inside of a single 2D topological insulator sheet. In particular, I show that edge states can effectively be dragged into the system's interior by gating a region near the edge into a metallic regime and then removing the resulting gapless carriers via proximity-induced superconductivity. This method allows one to construct rather general quasi-1D networks along which Majorana modes can be exchanged by electrostatic means.

Apart from 2D topological insulators, Majorana fermions can also be generated in other more accessible materials such as semiconductors. Following up on a suggestion by experimentalist Charlie Marcus, I proposed a novel geometry to create Majorana fermions by placing a 2D electron gas in proximity to an interdigitated superconductor-ferromagnet structure. This architecture evades several manufacturing challenges by allowing single-side fabrication and widening the class of 2D electron gas that may be used, such as the surface states of bulk semiconductors. Furthermore, it naturally allows one to trap and manipulate Majorana fermions through the application of currents. Thus, this structure may lead to the development of a circuit that enables fully electrical manipulation of topologically-protected quantum memory. To reveal these exotic Majorana zero modes, I also proposed an interference scheme to detect Majorana fermions that is broadly applicable to any 2D topological superconductor platform.

Imidazolium-Based Organic Structure Directing Agents for the Synthesis of Microporous Materials

The central theme of this thesis is the use of imidazolium-based organic structure directing agents (OSDAs) in microporous materials synthesis. Imidazoliums are advantageous OSDAs as they are relatively inexpensive and simple to prepare, show robust stability under microporous material synthesis conditions, have led to a wide range of products, and have many permutations in structure that can be explored. The work I present involves the use of mono-, di-, and triquaternary imidazolium-based OSDAs in a wide variety of microporous material syntheses. Much of this work was motivated by successful computational predictions (Chapter 2) that led me to continue to explore these types of OSDAs. Some of the important discoveries with these OSDAs include the following: 1) Experimental evaluation and confirmation of a computational method that predicted a new OSDA for pure-silica STW, a desired framework containing helical pores that was previously very difficult to synthesize. 2) Discovery of a number of new imidazolium OSDAs to synthesize zeolite RTH, a zeolite desired for both the methanol-to-olefins reaction as well as NOX reduction in exhaust gases. This discovery enables the use of RTH for many additional investigations as the previous OSDA used to make this framework was difficult to synthesize, such that no large scale preparations would be practical. 3) The synthesis of pure-silica RTH by topotactic condensation from a layered precursor (denoted CIT-10), that can also be pillared to make a new framework material with an expanded pore system, denoted CIT-11, that can be calcined to form a new microporous material, denoted CIT-12. CIT-10 is also interesting since it is the first layered material to contain 8 membered rings through the layers, making it potentially useful in separations if delamination methods can be developed. 4) The synthesis of a new microporous material, denoted CIT-7 (framework code CSV) that contains a 2-dimensional system of 8 and 10 membered rings with a large cage at channel intersections. This material is especially important since it can be synthesized as a pure-silica framework under low-water, fluoride-mediated synthesis conditions, and as an aluminosilicate material under hydroxide mediated conditions. 5) The synthesis of high-silica heulandite (HEU) by topotactic condensation as well as direct synthesis, demonstrating new, more hydrothermally stable compositions of a previously known framework. 6) The synthesis of germanosilicate and aluminophosphate LTA using a triquaternary OSDA. All of these materials show the diverse range of products that can be formed from OSDAs that can be prepared by straightforward syntheses and have made many of these materials accessible for the first time under facile zeolite synthesis conditions.

Coulomb drag and tunneling studies in quantum Hall bilayers

The bilayer quantum Hall state at total filling factor ν_T=1, where the total electron density matches the degeneracy of the lowest Landau level, is a prominent example of Bose-Einstein condensation of excitons. A macroscopically ordered state is realized where an electron in one layer is tightly bound to a "hole" in the other layer. If exciton transport were the only bulk transportmechanism, a current driven in one layer would spontaneously generate a current of equal magnitude and opposite sign in the other layer. The Corbino Coulomb drag measurements presented in this thesis demonstrate precisely this phenomenon.

Excitonic superfluidity has been long sought in the ν_T=1 state. The tunneling between the two electron gas layers exihibit a dc Josephson-like effect. A simple model of an overdamped voltage biased Josephson junction is in reasonable agreement with the observed tunneling I-V. At small tunneling biases, it exhibits a tunneling "supercurrent". The dissipation is carefully studied in this tunneling "supercurrent" and found to remain small but finite.

Influência da reorientação de spin nas propriedades termomagnéticas dos compostos da série Ho1-yGdyAl2

Neste trabalho foram estudadas as propriedades estruturais e termomagnéticas dos pseudobinários Ho1-yGdyAl2, através de abordagens experimentais e teóricas. A parte experimental envolveu a preparação de cinco amostras, com as concentrações y = 0, 0,25, 0,5, 0,75 e 1, assim como medidas de magnetização, calor especifico e da variação adiabática da temperatura. Na parte teórica usamos um hamiltoniano modelo que leva em consideração a interação dos íons com o campo magnético aplicado, com o campo elétrico cristalino e a troca entre os íons magnéticos. A entropia da rede foi considerada na aproximação de Debye e a entropia eletrônica na aproximação do gás de elétrons livres. A influência das reorientações de spin, espontâneas e induzidas pelo campo magnético, na magnetização e no calor especifico foram investigadas sistematicamente tanto a partir de dados experimentais quanto teoricamente. Também obtemos resultados teóricos para a variação de entropia e variação adiabática da temperatura alterando a intensidade ou a direção do campo magnético.

Quantifying seascape structure: extending terrestrial spatial pattern metrics to the marine realm

Spatial pattern metrics have routinely been applied to characterize and quantify structural features of terrestrial landscapes and have demonstrated great utility in landscape ecology and conservation planning. The important role of spatial structure in ecology and management is now commonly recognized, and recent advances in marine remote sensing technology have facilitated the application of spatial pattern metrics to the marine environment. However, it is not yet clear whether concepts, metrics, and statistical techniques developed for terrestrial ecosystems are relevant for marine species and seascapes. To address this gap in our knowledge, we reviewed, synthesized, and evaluated the utility and application of spatial pattern metrics in the marine science literature over the past 30 yr (1980 to 2010). In total, 23 studies characterized seascape structure, of which 17 quantified spatial patterns using a 2-dimensional patch-mosaic model and 5 used a continuously varying 3-dimensional surface model. Most seascape studies followed terrestrial-based studies in their search for ecological patterns and applied or modified existing metrics. Only 1 truly unique metric was found (hydrodynamic aperture applied to Pacific atolls). While there are still relatively few studies using spatial pattern metrics in the marine environment, they have suffered from similar misuse as reported for terrestrial studies, such as the lack of a priori considerations or the problem of collinearity between metrics. Spatial pattern metrics offer great potential for ecological research and environmental management in marine systems, and future studies should focus on (1) the dynamic boundary between the land and sea; (2) quantifying 3-dimensional spatial patterns; and (3) assessing and monitoring seascape change.

Modeling the effect of snow on seasonal runoff within the Truckee River drainage basin

EXTRACT (SEE PDF FOR FULL ABSTRACT): We estimate monthly runoff for a 2-dimensional solution domain containing those areas tributary to Pyramid Lake, Nevada (the Truckee River drainage basin) at a 1-kilometer grid cell spacing. ... To calculate the effect of snow on the hydrologic system, we perform two experiments. In the first we assume that all precipitation falls as rain; in the second we assume that some precipitation falls as snow, thus available water is a combination of rain and snowmelt. We find that considering the effect of snow results in a more accurate representation of mean monthly flow rates, in particular the peak flow during the melt season in the Sierra Nevada. These preliminary results indicate that a relatively simple snow model can improve the representation of Truckee River basin hydrology, significantly reducing errors in modeled seasonal runoff.

Pipeline uplift response in thawed sandy backfills: Preliminary findings

Offshore and onshore buried pipelines under high operating temperature and pressures may lead to upheaval buckling (UHB) if sufficient soil cover is not present to prevent the upward movement of the pipeline. In regions where seasonal changes involve ground soil undergoing freezing-thawing cycles, the uplift resistance from soil cover may be minimum when the soil is undergoing thawing. This paper presents the results from 2 directly-comparable minidrum centrifuge tests conducted at the Schofield Centre, University of Cambridge, to investigate the difference in uplift resistance responses between fully-saturated and thawed sandy backfill conditions. Both tests were conducted drained at 30g using an 8.6 mm diameter aluminium model pipe, corresponding to a prototype pipe diameter of 258 mm. The soil cover/pipe diameter ratio, H/D, was kept at 1. Fraction E fine silica sand was used as the backfill. Preliminary experimental results indicated that the ultimate uplift resistance of a thawing sand backfill to be lower than that of a fully saturated sand backfill. This suggests that in regions where backfill soil undergoes freeze-thaw cycles, the thawing backfill may be more critical than fully saturated backfill for uplift resistance. The 2-dimensional displacement field during the experiment was accurately measured and analysed using the Particle Image Velocimetry technique. Copyright © 2011 by the International Society of Offshore and Polar Engineers (ISOPE).

Enticing general engineering undergrads into signal processing with Photo-Editing and image search

The University of Cambridge is unusual in that its Department of Engineering is a single department which covers virtually all branches of engineering under one roof. In their first two years of study, our undergrads study the full breadth of engineering topics and then have to choose a specialization area for the final two years of study. Here we describe part of a course, given towards the end of their second year, which is designed to entice these students to specialize in signal processing and information engineering topics for years 3 and 4. The course is based around a photo editor and an image search application, and it requires no prior knowledge of the z-transform or of 2-dimensional signal processing. It does assume some knowledge of 1-D convolution and basic Fourier methods and some prior exposure to Matlab. The subject of this paper, the photo editor, is written in standard Matlab m-files which are fully visible to the students and help them to see how specific algorithms are implemented in detail. © 2011 IEEE.

Eliminating discharge imbalances in predicting shallow water flows with shocks

The shallow water equations are widely used in modelling environmental flows. Being a hyperbolic system of differential equations, they admit shocks that represent hydraulic jumps and bores. Although the water surface can be solved satisfactorily with the modern shock-capturing schemes, the predicted flow rate often suffers from imbalances where shocks occur, eg the mass conservation is violated by failing to maintain a constant discharge rate at every cross-section in a steady open channel flow. A total-variation-diminishing Lax-Wendroff scheme is developed, and used to demonstrate how to achieve an exact flux balance. The performance of the proposed methods is inspected through some test cases, which include 1- and 2-dimensional, flat and irregular bed scenarios. The proposed methods are shown to preserve the mass exactly, and can be easily extended to other shock-capturing models.