Geometric approach to Fletcher's ideal penalty function

Original article can be found at: www.springerlink.com Copyright Springer. [Originally produced as UH Technical Report 280, 1993]In this note, we derive a geometric formulation of an ideal penalty function for equality constrained problems. This differentiable penalty function requires no parameter estimation or adjustment, has numerical conditioning similar to that of the target function from which it is constructed, and also has the desirable property that the strict second-order constrained minima of the target function are precisely those strict second-order unconstrained minima of the penalty function which satisfy the constraints. Such a penalty function can be used to establish termination properties for algorithms which avoid ill-conditioned steps. Numerical values for the penalty function and its derivatives can be calculated efficiently using automatic differentiation techniques.Peer reviewe

Function reconstruction as a classical moment problem: A maximum entropy approach

We present a systematic study of the reconstruction of a non-negative function via maximum entropy approach utilizing the information contained in a finite number of moments of the function. For testing the efficacy of the approach, we reconstruct a set of functions using an iterative entropy optimization scheme, and study the convergence profile as the number of moments is increased. We consider a wide variety of functions that include a distribution with a sharp discontinuity, a rapidly oscillatory function, a distribution with singularities, and finally a distribution with several spikes and fine structure. The last example is important in the context of the determination of the natural density of the logistic map. The convergence of the method is studied by comparing the moments of the approximated functions with the exact ones. Furthermore, by varying the number of moments and iterations, we examine to what extent the features of the functions, such as the divergence behavior at singular points within the interval, is reproduced. The proximity of the reconstructed maximum entropy solution to the exact solution is examined via Kullback-Leibler divergence and variation measures for different number of moments.Comment: 20 pages, 17 figure

A simplified protocol for detecting two systemic bait markers (Rhodamine B and iophenoxic acid) in small mammals

We developed a method of quantifying levels of fluorescence in the whiskers of wild stoats (Mustela erminea) using fluorescence microscopy and Axiovision 3.0.6.1 software. The method allows for discrimination between natural fluorescence present in or on a whisker, and the fluorescence resulting from the ingestion of the systemic marker Rhodamine B (RB), although some visual judgement is still required. We also developed a new high performance liquid chromatography (HPLC) protocol for detecting the systemic marker iophenoxic acid (IPA) in the blood of laboratory rats (Rattus norvegicus) and wild stoats. With this method, the blood of an animal that has consumed IPA can be tested for the presence of the foreign IPA compound itself. This is a more reliable test than the previous method, which measured the raised level of natural blood protein-bound iodine correlated with IPA absorption. The quantity of blood required from animal subjects is very small (10 μl), so the testing is less intrusive and the method can be extended to smaller species. The extraction technique uses methanol, rather than acids and heavy metal salts, thereby simplifying the procedure. Recovery of IPA is quantitative, giving a highly reliable reading. In experiments on captive rats the IPA method proved successful. Of 12 positively marked carcasses, two that had not been frozen for the 24 h before blood samples were taken showed relatively lower IPA levels. The same IPA detection method, as well as the whisker analysis for RB, was applied successfully to a population of wild stoats to which both Rhodamine B and IPA were made available at bait stations. The presence of both bait markers was detectable in rats for at least 21 days and in stoats for at least 27 days

Tight-binding parameterization for photonic band gap materials

The ideas of the linear combination of atomic orbitals (LCAO) method, well known from the study of electrons, is extended to the classical wave case. The Mie resonances of the isolated scatterer in the classical wave case, are analogous to the localized eigenstates in the electronic case. The matrix elements of the two-dimensional tight-binding (TB) Hamiltonian are obtained by fitting to ab initio results. The transferability of the TB model is tested by reproducing accurately the band structure of different 2D lattices, with and without defects, thus proving that the obtained TB parameters can be used to study other properties of the photonic band gap materials.Comment: 4 pages, 3 postscript figures, sumbitted to Phys. rev. Let

The Physicist's Guide to the Orchestra

An experimental study of strings, woodwinds (organ pipe, flute, clarinet, saxophone and recorder), and the voice was undertaken to illustrate the basic principles of sound production in music instruments. The setup used is simple and consists of common laboratory equipment. Although the canonical examples (standing wave on a string, in an open and closed pipe) are easily reproduced, they fail to explain the majority of the measurements. The reasons for these deviations are outlined and discussed.Comment: 11 pages, 10 figures (jpg files). Submitted to European Journal of Physic

The statistical mechanics of complex signaling networks : nerve growth factor signaling

It is becoming increasingly appreciated that the signal transduction systems used by eukaryotic cells to achieve a variety of essential responses represent highly complex networks rather than simple linear pathways. While significant effort is being made to experimentally measure the rate constants for individual steps in these signaling networks, many of the parameters required to describe the behavior of these systems remain unknown, or at best, estimates. With these goals and caveats in mind, we use methods of statistical mechanics to extract useful predictions for complex cellular signaling networks. To establish the usefulness of our approach, we have applied our methods towards modeling the nerve growth factor (NGF)-induced differentiation of neuronal cells. Using our approach, we are able to extract predictions that are highly specific and accurate, thereby enabling us to predict the influence of specific signaling modules in determining the integrated cellular response to the two growth factors. We show that extracting biologically relevant predictions from complex signaling models appears to be possible even in the absence of measurements of all the individual rate constants. Our methods also raise some interesting insights into the design and possible evolution of cellular systems, highlighting an inherent property of these systems wherein particular ''soft'' combinations of parameters can be varied over wide ranges without impacting the final output and demonstrating that a few ''stiff'' parameter combinations center around the paramount regulatory steps of the network. We refer to this property -- which is distinct from robustness -- as ''sloppiness.''Comment: 24 pages, 10 EPS figures, 1 GIF (makes 5 multi-panel figs + caption for GIF), IOP style; supp. info/figs. included as brown_supp.pd

Impact of Scottish vocational qualifications on residential child care : have they fulfilled the promise?

This article will present findings from a doctoral study exploring the impact of 'SVQ Care: Promoting Independence (level III)' within children's homes. The study focuses on the extent to which SVQs enhance practice and their function within a 'learning society'. A total of 30 staff were selected from seven children's homes in two different local authority social work departments in Scotland. Each member of staff was interviewed on four separate occasions over a period of 9 months. Interviews were structured using a combination of repertory grids and questions. Particular focus was given to the assessment process, the extent to which SVQs enhance practice and the learning experiences of staff. The findings suggest that there are considerable deficiencies both in terms of the SVQ format and the way in which children's homes are structured for the assessment of competence. Rather than address the history of failure within residential care, it appears that SVQs have enabled the status quo to be maintained whilst creating an 'illusion' of change within a learning society

Bound states of edge dislocations: The quantum dipole problem in two dimensions

We investigate bound state solutions of the 2D Schr\"odinger equation with a dipole potential originating from the elastic effects of a single edge dislocation. The knowledge of these states could be useful for understanding a wide variety of physical systems, including superfluid behavior along dislocations in solid $^4$He. We present a review of the results obtained by previous workers together with an improved variational estimate of the ground state energy. We then numerically solve the eigenvalue problem and calculate the energy spectrum. In our dimensionless units, we find a ground state energy of -0.139, which is lower than any previous estimate. We also make successful contact with the behavior of the energy spectrum as derived from semiclassical considerations.Comment: 6 pages, 3 figures, submitted to PR

Ultraviolet Line Profiles of Slowly Rotating Massive Star Winds Using the "Analytic Dynamical Magnetosphere" Formalism

Recent large-scale spectropolarimetric surveys have established that a small but significant percentage of massive stars host stable, surface dipolar magnetic fields with strengths on the order of kG. These fields channel the dense, radiatively driven stellar wind into circumstellar magnetospheres, whose density and velocity structure can be probed using ultraviolet (UV) spectroscopy of wind-sensitive resonance lines. Coupled with appropriate magnetosphere models, UV spectroscopy provides a valuable way to investigate the wind-field interaction, and can yield quantitative estimates of the wind parameters of magnetic massive stars. We report a systematic investigation of the formation of UV resonance lines in slowly rotating magnetic massive stars with dynamical magnetospheres. We pair the Analytic Dynamical Magnetosphere (ADM) formalism with a simplified radiative transfer technique to produce synthetic UV line profiles. Using a grid of models, we examine the effect of magnetosphere size, the line strength parameter, and the cooling parameter on the structure and modulation of the line profile. We find that magnetic massive stars uniquely exhibit redshifted absorption at most viewing angles and magnetosphere sizes, and that significant changes to the shape and variation of the line profile with varying line strengths can be explained by examining the individual wind components described in the ADM formalism. Finally, we show that the cooling parameter has a negligible effect on the line profiles.Comment: 16 pages, 15 figures, accepted to MNRA