Approximating the Set of Local Minima in Partial RNA Folding Landscapes

Motivation: We study a stochastic method for approximating the set of local minima in partial RNA folding landscapes associated with a bounded-distance neighbourhood of folding conformations. The conformations are limited to RNA secondary structures without pseudoknots. The method aims at exploring partial energy landscapes pL induced by folding simulations and their underlying neighbourhood relations. It combines an approximation of the number of local optima devised by Garnier and Kallel (2002) with a run-time estimation for identifying sets of local optima established by Reeves and Eremeev (2004).

Results: The method is tested on nine sequences of length between 50 nt and 400 nt, which allows us to compare the results with data generated by RNAsubopt and subsequent barrier tree calculations. On the nine sequences, the method captures on average 92% of local minima with settings designed for a target of 95%. The run-time of the heuristic can be estimated by O(n2D?ln?), where n is the sequence length, ? is the number of local minima in the partial landscape pL under consideration and D is the maximum number of steepest descent steps in attraction basins associated with pL.

Density functional and Monte Carlo studies of sulfur. I. Structure and bonding in S-n rings and chains (n=2-18) (vol 118, pg 9257, 2003)

Density functional calculations have been performed for ring isomers of sulfur with up to 18 atoms, and for chains with up to ten atoms. There are many isomers of both types, and the calculations predict the existence of new forms. Larger rings and chains are very flexible, with numerous local energy minima. Apart from a small, but consistent overestimate in the bond lengths, the results reproduce experimental structures where known. Calculations are also performed on the energy surfaces of S8 rings, on the interaction between a pair of such rings, and the reaction between one S8 ring and the triplet diradical S8 chain. The results for potential energies, vibrational frequencies, and reaction mechanisms in sulfur rings and chains provide essential ingredients for Monte Carlo simulations of the liquid–liquid phase transition. The results of these simulations will be presented in Part II.

Path-integral study of a two-dimensional Lennard-Jones glass

The glass transition in a quantum Lennard-Jones mixture is investigated by constant-volume path-integral simulations. Particles are assumed to be distinguishable, and the strength of quantum effects is varied by changing h from zero (the classical case) to one (corresponding to a highly quantum-mechanical regime). Quantum delocalization and zero point energy drastically reduce the sensitivity of structural and thermodynamic properties to the glass transition. Nevertheless, the glass transition temperature T-g can be determined by analyzing the phase space mobility of path-integral centroids. At constant volume, the T-g of the simulated model increases monotonically with increasing h. Low temperature tunneling centers are identified, and the quantum versus thermal character of each center is analyzed. The relation between these centers and soft quasilocalized harmonic vibrations is investigated. Periodic minimizations of the potential energy with respect to the positions of the particles are performed to determine the inherent structure of classical and quantum glassy samples. The geometries corresponding to these energy minima are found to be qualitatively similar in all cases. Systematic comparisons for ordered and disordered structures, harmonic and anharmonic dynamics, classical and quantum systems show that disorder, anharmonicity, and quantum effects are closely interlinked.

HOS-Based Semi-Blind Spatial Equalization for MIMO Rayleigh Fading Channels

In this paper we concentrate on the direct semi-blind spatial equalizer design for MIMO systems with Rayleigh fading channels. Our aim is to develop an algorithm which can outperform the classical training based method with the same training information used, and avoid the problems of low convergence speed and local minima due to pure blind methods. A general semi-blind cost function is first constructed which incorporates both the training information from the known data and some kind of higher order statistics (HOS) from the unknown sequence. Then, based on the developed cost function, we propose two semi-blind iterative and adaptive algorithms to find the desired spatial equalizer. To further improve the performance and convergence speed of the proposed adaptive method, we propose a technique to find the optimal choice of step size. Simulation results demonstrate the performance of the proposed algorithms and comparable schemes.

Evidence for two-phase regions in BST thin films from capacitance-voltage data

The Curie-Weiss plots of reciprocal dielectric constant versus temperature, in Ba0.5Sr0.5TiO3 films grown onto SrRuO3 lower electrodes by pulsed-laser deposition, show two minima below film thicknesses of 280 nm. This double minima implies possible mixed phases in the thin films. A graphical plot of capacitance for decreasing dc voltage versus that of increasing dc voltage shows a well-defined triangular shape for both Pb(Zr0.4Ti0.6)O-3 and SrBi2Ta2O9 thin films. However, for a 175-nm-thick Ba0.5Sr0.5TiO3 thin film, the plot shows an overlapping of two triangles, suggesting mixed phases. This graphical method appears to be effective in detecting structural subtleties in ferroelectric capacitors.

Stochastic protein folding simulation in the three-dimensional HP-model

We present results from three-dimensional protein folding simulations in the HP-model on ten benchmark problems. The simulations are executed by a simulated annealing-based algorithm with a time-dependent cooling schedule. The neighbourhood relation is determined by the pull-move set. The results provide experimental evidence that the maximum depth D of local minima of the underlying energy landscape can be upper bounded by D < n(2/3). The local search procedure employs the stopping criterion (In/delta)(D/gamma) where m is an estimation of the average number of neighbouring conformations, gamma relates to the mean of non-zero differences of the objective function for neighbouring conformations, and 1-delta is the confidence that a minimum conformation has been found. The bound complies with the results obtained for the ten benchmark problems. (c) 2008 Elsevier Ltd. All rights reserved.

Robust fuzzy Gustafson-Kessel clustering for nonlinear system identification

This paper deals with Takagi-Sugeno (TS) fuzzy model identification of nonlinear systems using fuzzy clustering. In particular, an extended fuzzy Gustafson-Kessel (EGK) clustering algorithm, using robust competitive agglomeration (RCA), is developed for automatically constructing a TS fuzzy model from system input-output data. The EGK algorithm can automatically determine the 'optimal' number of clusters from the training data set. It is shown that the EGK approach is relatively insensitive to initialization and is less susceptible to local minima, a benefit derived from its agglomerate property. This issue is often overlooked in the current literature on nonlinear identification using conventional fuzzy clustering. Furthermore, the robust statistical concepts underlying the EGK algorithm help to alleviate the difficulty of cluster identification in the construction of a TS fuzzy model from noisy training data. A new hybrid identification strategy is then formulated, which combines the EGK algorithm with a locally weighted, least-squares method for the estimation of local sub-model parameters. The efficacy of this new approach is demonstrated through function approximation examples and also by application to the identification of an automatic voltage regulation (AVR) loop for a simulated 3 kVA laboratory micro-machine system.

The invasive genus Asparagopsis (Bonnemaisoniaceae, Rhodophyta): Molecular systematics, morphology, and ecophysiology of Falkenbergia isolates

The genus Asparagopsis was studied using 25 Falkenbergia tetrasporophyte strains collected worldwide. Plastid (cp) DNA RFLP revealed three groups of isolates, which differed in their small subunit rRNA gene sequences, temperature responses, and tetrasporophytic morphology (cell sizes). Strains from Australia, Chile, San Diego, and Atlantic and Mediterranean Europe were identifiable as A. armata Harvey, the gametophyte of which has distinctive barbed spines. This species is believed to be endemic to cold-temperate waters of Australia and New Zealand and was introduced into Europe in the 1920s. All isolates showed identical cpDNA RFLPs, consistent with a recent introduction from Australia. Asparagopsis taxiformis (Delile) Trevisan, the type and only other recognized species, which lacks spines, is cosmopolitan in warm-temperate to tropical waters. Two clades differed morphologically and ecophysiologically and in the future could be recognized as sibling species or subspecies. A Pacific/Italian clade had 4-8degrees C lower survival minima and included a genetically distinct apomictic isolate from Western Australia that corresponded to the form of A. taxiformis originally described as A. sanfordiana Harvey. The second clade, from the Caribbean and the Canaries, is stenothermal (subtropical to tropical) with some ecotypic variation. The genus Asparagopsis consists of two or possibly three species, but a definitive taxonomic treatment of the two A. taxiformis clades requires study of field-collected gametophytes.

Population-based local search for protein folding simulation in the MJ energy model and cubic lattices

We present experimental results on benchmark problems in 3D cubic lattice structures with the Miyazawa-Jernigan energy function for two local search procedures that utilise the pull-move set: (i) population-based local search (PLS) that traverses the energy landscape with greedy steps towards (potential) local minima followed by upward steps up to a certain level of the objective function; (ii) simulated annealing with a logarithmic cooling schedule (LSA). The parameter settings for PLS are derived from short LSA-runs executed in pre-processing and the procedure utilises tabu lists generated for each member of the population. In terms of the total number of energy function evaluations both methods perform equally well, however. PLS has the potential of being parallelised with an expected speed-up in the region of the population size. Furthermore, both methods require a significant smaller number of function evaluations when compared to Monte Carlo simulations with kink-jump moves. (C) 2009 Elsevier Ltd. All rights reserved.

Design of asymmetrical RF and microwave bandpass filters by computer optimization

This paper presents an optimization-based approach to the design of asymmetrical filter structures having the maximum number of return- or insertion-loss ripples in the passband such as those based upon Chebyshev function prototypes. The proposed approach. has the following advantages over the general purpose optimization techniques adopted previously such as: less frequency sampling is required, optimization is carried out with respect to the Chebyshev (or minimax) criterion, the problem of local minima does not arise, and optimization is usually only required for the passband. When implemented around an accurate circuit simulation, the method can be used to include all the effects of discontinuities, junctions, fringing, etc. to reduce the amount of tuning required in the final filter. The design of asymmetrical ridged-waveguide bandpass filters is considered as an example. Measurements on a fabricated filter confirm the accuracy of the design procedure.

Geometric bonding effects in the X2A1, A2∑+u, and B2IIg states of Li2F

Published ab-initio and pseudopotential calculations for the dialkali halide systems suggest that the preferred co-linear geometry is for the metal to approach the metal end of the alkali halide. Here, ab-initio calculations on the Li2F system reveal that the well depth on the halide side in this radical is much deeper and is a local saddle-point associated with the ionic non-linear global minima. Although many features of the pseudopotential surfaces are confirmed, significant differences are apparent including the existence of a linear excited state instead of a triangular one, a considerably deeper global minimum some 50% lower in energy and a close approach between the X2A1 and the states, with the minimum 87 kJ mol-1 below the ground state asymptote. All the results can be rationalised as the avoided crossings between a long range, covalent potential dominant within the LiLiF geometry and an ionic state that forms the global minimum. Calculations on the 3rd 2A' potential indicate that even for Li + LiF collisions at ultracold temperatures the collision dynamics could involve as many as three electronic states.

Changes in solar activity and Holocene climatic shifts derived from C-14 wiggle-match dated peat deposits

Closely spaced sequences of accelerator mass spectrometer (AMS) C-14 dates of peat deposits display century-scale wiggles which can be fitted to the radiocarbon calibration curve. By wiggle-matching such sequences, high-precision calendar age chronologies can be generated which show that changes in mire surface wetness during the Bronze Age/Iron Age transition (c. 850 cal. BC) and the 'Little Ice Age' (Wolf, Sporer, Maunder and Dalton Minima) occurred during periods of suddenly increasing atmospheric concentration of C-14. Replicate evidence from peat-based proxy climate indicators in northwest Europe suggest these changes in climate may have been driven by temporary declines of solar activity. Carbon-accumulation rates of two raised peat bogs in the UK and Denmark record low values during the 'Little Ice Age' which reflects reduced primary productivity of the peat-forming vegetation during these periods of climatic deterioration.

Simple method for the rapid simultaneous screening of photocatalytic activity over multiple positions of self-cleaning films

An intelligent ink, previously shown to be capable of rapidly assessing photocatalytic activity, was simply applied via a felt-pen onto a commercially available piece of Activ (TM) self-cleaning glass. The ink, comprising of redox dye resazurin and the sacrificial electron donor glycerol within an aqueous hydroxy ethyl cellulose (HEC) polymer media, was photocatalytically degraded in a two-step process. The key initial stage was the photo-reductive conversion of resazurin to resorufin, whereby a colour change from blue to pink occurred. The latter stage was the subsequent photo-reduction of the resorufin, where a slower change from pink to colourless was seen. Red and green components of red-green-blue colour extracted from flat-bed scanner digital images of resazurin ink coated photocatalytic films at intervals during the photocatalysis reaction were inversely proportional to the changes seen via UV-visible absorption spectroscopy and indicative of reaction kinetics. A 3 x 3 grid of intelligent ink was drawn onto a piece of Activ (TM) and a glass blank. The photocatalysis reaction was monitored solely by flat-bed digital scanning. Red-green-blue values of respective positions on the grid were extracted using a custom-built program entitled RGB Extractor (c). The program was capable of extracting a number of 5 x 5 pixel averages of red-green-blue components simultaneously. Allocation of merely three coordinates allowed for the automatic generation of a grid, with scroll-bars controlling the number of positions to be extracted on the grid formed. No significant change in red and green components for any position on the glass blank was observed; however, the Activ (TM) film displayed a homogenous photo-reduction of the dye, reaching maxima in red and minima in green components in 23 +/- 3 and 14 +/- 2 min, respectively. A compositionally graded N-doped titania film synthesised in house via a combinatorial APCVD reaction was also photocatalytically tested by this method where 247 positions on a 13 x 19 grid were simultaneously analysed. The dramatic variation in photocatalysis observed was rapidly quantified for all positions (2-3 hours) allowing for correlations to be made between thicknesses and N : Ti% compositions attained from Swanepoel and WDX analysis, respectively. N incorporation within this system was found to be detrimental to film activity for the photocatalysis reaction of intelligent ink under 365 nm light.

Partial and total multiphoton detachment rates for H- in a perturbative B-spline-based approach

An ab initio approach has been applied to study multiphoton detachment rates for the negative hydrogen ion in the lowest nonvanishing order of perturbation theory. The approach is based on the use of B splines allowing an accurate treatment of the electronic repulsion. Total detachment rates have been determined for two- to six-photon processes as well as partial rates for detachment into the different final symmetries. It is shown that B-spline expansions can yield accurate continuum and bound-state wave functions in a very simple manner. The calculated total rates for two- and three-photon detachment are in good agreement with other perturbative calculations. For more than three-photon detachment little information has been available before now. While the total cross sections show little structure, a fair amount of structure is predicted in the partial cross sections. In the two-photon process, it is shown that the detached electrons mainly have s character. For four- and six-photon processes, the contribution from the d channel is the most important. For three- and five-photon processes p electrons dominate the electron emission spectrum. Detachment rates for s and p electrons show minima as a function of photon energy. © 1994 The American Physical Society.

Structure and interactions of ultracold Yb ions and Rb atoms

In order to study ultracold charge-transfer processes in hybrid atom-ion traps, we have mapped out the potential-energy curves and molecular parameters for several low-lying states of the Rb, Yb+ system. We employ both a multireference configuration interaction and a full configuration interaction (FCI) approach. Turning points, crossing points, potential minima, and spectroscopic molecular constants are obtained for the lowest five molecular states. Long-range parameters, including the dispersion coefficients, are estimated from our ab initio data. The separated-atom ionization potentials and atomic polarizability of the ytterbium atom (ad=128.4 atomic units) are in good agreement with experiment and previous calculations. We present some dynamical calculations for (adiabatic) scattering lengths for the two lowest (Yb, Rb+) channels that were carried out in our work. However, we find that the pseudopotential approximation is rather limited in validity and only applies to nK temperatures. The adiabatic scattering lengths for both the triplet and singlet channels indicate that both are large and negative in the FCI approximation.