Energy conserving schemes for the simulation of musical instrument contact dynamics

Collisions are an innate part of the function of many musical instruments. Due to the nonlinear nature of contact forces, special care has to be taken in the construction of numerical schemes for simulation and sound synthesis. Finite difference schemes and other time-stepping algorithms used for musical instrument modelling purposes are normally arrived at by discretising a Newtonian description of the system. However because impact forces are non-analytic functions of the phase space variables, algorithm stability can rarely be established this way. This paper presents a systematic approach to deriving energy conserving schemes for frictionless impact modelling. The proposed numerical formulations follow from discretising Hamilton׳s equations of motion, generally leading to an implicit system of nonlinear equations that can be solved with Newton׳s method. The approach is first outlined for point mass collisions and then extended to distributed settings, such as vibrating strings and beams colliding with rigid obstacles. Stability and other relevant properties of the proposed approach are discussed and further demonstrated with simulation examples. The methodology is exemplified through a case study on tanpura string vibration, with the results confirming the main findings of previous studies on the role of the bridge in sound generation with this type of string instrument.

Discrete-Time Conserved Quantities for Damped Oscillators

Numerical sound synthesis is often carried out using the finite difference time domain method. In order to analyse the stability of the derived models, energy methods can be used for both linear and nonlinear settings. For Hamiltonian systems the existence of a conserved numerical energy-like quantity can be used to guarantee the stability of the simulations. In this paper it is shown how to derive similar discrete conservation laws in cases where energy is dissipated due to friction or in the presence of an energy source due to an external force. A damped harmonic oscillator (for which an analytic solution is available) is used to present the proposed methodology. After showing how to arrive at a conserved quantity, the simulation of a nonlinear single reed shows an example of an application in the context of musical acoustics.

Polarisation analysis of VUV synchrotron radiation emitted from a bending magnet source in the energy range 20-50 eV: A comparison between measurements and theoretical predictions

The results of a study to characterise the polarisation properties of the photon beam emerging from beamline 5D, mounted on a bending magnet source at the Synchrotron Radiation Source, Daresbury Laboratory, are presented. The expectation values for the Stokes parameters corresponding to the light transmitted by the beamline have been calculated by combining ray-tracing and optical methods. The polarisation of the light at the source is modified both by the beamline geometry and by the reflections at the optical components. Although it is often assumed that the polarising properties of grazing incidence optics are negligible, this assumption leads to rather inaccurate results in the VUV region. A study of the reflectivity shows that even at incidence angles (theta(i) = 80-85degrees) which are far from the Brewster angle (theta(B) similar to 45degrees for VUV and soft X-ray radiation) the residual changes in the amplitudes of the reflected light can result in non-negligible polarisation effects. Furthermore, reflection at grazing incidence gives rise to a substantial change in the phase, and this has the effect of rotating the elliptically polarised state. Theoretical Stokes parameters have been compared with full polarisation measurements obtained using a reflection polarimeter in the energy range 20-40 eV. (C) 2003 Elsevier B.V. All rights reserved.

The transfer of energy between electrons and ions in solids

In this review we consider those processes in condensed matter that involve the irreversible flow of energy between electrons and nuclei that follows from a system being taken out of equilibrium. We survey some of the more important experimental phenomena associated with these processes, followed by a number of theoretical techniques for studying them. The techniques considered are those that can be applied to systems containing many nonequivalent atoms. They include both perturbative approaches (Fermi's Golden Rule and non-equilibrium Green's functions) and molecular dynamics based (the Ehrenfest approximation, surface hopping, semi-classical Gaussian wavefunction methods and correlated electron-ion dynamics). These methods are described and characterized, with indications of their relative merits.

A new lead for nonpeptidic active-site-directed inhibitors of the severe acute respiratory syndrome coronavirus main protease discovered by a combination of screening and docking methods.

The coronavirus main protease, Mpro, is considered to be a major target for drugs suitable for combating coronavirus infections including severe acute respiratory syndrome (SARS). An HPLC-based screening of electrophilic compounds that was performed to identify potential Mpro inhibitors revealed etacrynic acid tert-butylamide (6a) as an effective nonpeptidic inhibitor. Docking studies suggested a binding mode in which the phenyl ring acts as a spacer bridging the inhibitor's activated double bond and its hydrophobic tert-butyl moiety. The latter is supposed to fit into the S4 pocket of the target protease. Furthermore, these studies revealed etacrynic acid amide (6b) as a promising lead for nonpeptidic active-site-directed Mpro inhibitors. In a fluorimetric enzyme assay using a novel fluorescence resonance energy transfer (FRET) pair labeled substrate, compound 6b showed a Ki value of 35.3 M. Since the novel lead compound does not target the S1', S1, and S2 subsites of the enzyme's substrate-binding pockets, there is room for improvement that underlines the lead character of compound 6b.

Energy levels, radiative rates, and excitation rates for transitions in OIV

Aims. In this paper we report calculations for energy levels, radiative rates, and excitation rates for transitions in O IV. Methods. The grasp (general-purpose relativistic atomic structure package) and FAC (flexible atomic code) were adopted for calculating energy levels and radiative rates, and the Dirac atomic R-matrix code (DARC) used to determine the excitation rates. Results. Oscillator strengths and radiative rates are reported for all E1, E2, M1, and M2 transitions among the lowest 75 levels of O IV. Additionally, lifetimes are reported for all levels and comparisons made with those available in the literature. Finally, effective collision strengths are reported for all transitions over a wide temperature range below 106 K. Comparisons are made with earlier results and the accuracy of the data is assessed.

Energy levels, radiative rates and excitation rates for transitions in Ni XI

Aims. In this paper we report calculations for energy levels, radiative rates and excitation rates for transitions in Ni xi.
Methods. The grasp (General-purpose Relativistic Atomic Structure Package) and fac (Flexible Atomic Code) have been adopted
for calculating energy levels and radiative rates, and the Dirac Atomic R-matrix Code (darc) has been used to determine the excitation
rates.
Results. Oscillator strengths, radiative rates and line strengths are reported for all E1, E2, M1 and M2 transitions among the lowest
250 levels of Ni xi. Additionally, lifetimes are also reported for all levels. However, results for excitation rates are presented only for
transitions among the lowest 17 levels.

Measurement of Quenching Coefficients and Developement of Calibration Methods for Quantitative Spectroscopy of Plasma at Elevated Pressures

Measurements of collisional de-excitation (quenching) coefficients required for the interpretation of emission and fluorescence spectroscopic measurements are reported. Particular attention is turned on argon transitions which are of interest for actinometric determinations of atomic ground state populations and on fluorescence lines originating from excited atoms and noble gases in connection with two-photon excitation (TALIF) of atomic radicals. A novel method is described which allows to infer quenching coefficients for collisions with molecular hydrogen of noble gas states in the energy range up to 24 eV. The excitation is performed in these experiments by collisions of energetic electrons in the sheath of an RF excited hydrogen plasma during the field reversal phase which lasts about 10 ns. We describe in addition a calibration method - including quenching effects - for the determination by TALIF of absolute atomic radical densities of hydrogen, nitrogen and oxygen using two-photon resonances in noble gases close by the resonances of the species mentioned. The paper closes with first ideas on a novel technique to bypass quenching effects in TALIF by introducing an additional, controllable loss by photoionization that will allow quenching-free determination of absolute atomic densities with prevalent nanosecond laser systems in situations where collisional de-excitation dominates over spontaneous emission.

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.

Energy levels, radiative rates, and electron impact excitation rates for transitions in OVII

Aims.
In this paper we report calculations for energy levels, radiative rates, and electron impact excitation rates for transitions in O vii.
Methods.
The grasp (general-purpose relativistic atomic structure package) is adopted for calculating energy levels and radiative
rates. For determining the collision strengths and subsequently the excitation rates, the Dirac atomic R-matrix code (darc) and the
flexible atomic code (fac) are used.
Results.
Oscillator strengths, radiative rates, and line strengths are reported for all E1, E2, M1, and M2 transitions among the lowest
49 levels of O vii. Collision strengths have been averaged over a Maxwellian velocity distribution, and the resulting effective collision
strengths are reported over a wide temperature range below 2 × 106 K. Additionally, lifetimes are also listed for all levels.
Key words.

The control of wave energy converters using active bipolar damping

A novel method for controlling wave energy converters using active bipolar damping is described and compared with current control methods. The performance of active bipolar damping is modelled numerically for two distinct types of wave energy converter and it is found that in both cases the power capture can be significantly increased relative to optimal linear damping. It is shown that this is because active bipolar damping has the potential for providing a quasi-spring or quasi-inertia, which improves the wave energy converter's tuning and amplitude of motion, resulting in the increase in power capture observed. The practical implementation of active bipolar damping is also discussed. It is noted that active bipolar damping does not require a reactive energy store and thereby reduces the cost and eliminates losses due to the cycling of reactive energy. It is also noted that active bipolar damping could be implemented using a single constant pressure double-acting hydraulic cylinder and so potentially represents a simple, efficient, robust and economic solution to the control of wave energy converters.

Resting energy expenditure in non-ventilated, non-sedated patients recovering from serious traumatic brain injury: Comparison of prediction equations with indirect calorimetry values

Background & aims: Little is known about energy requirements in brain injured (TBI) patients, despite evidence suggesting adequate nutritional support can improve clinical outcomes. The study aim was to compare predicted energy requirements with measured resting energy expenditure (REE) values, in patients recovering from TBI.

Methods: Indirect calorimetry (IC) was used to measure REE in 45 patients with TBI. Predicted energy requirements were determined using FAO/WHO/UNU and Harris–Benedict (HB) equations. Bland– Altman and regression analysis were used for analysis.

Results: One-hundred and sixty-seven successful measurements were recorded in patients with TBI. At an individual level, both equations predicted REE poorly. The mean of the differences of standardised areas of measured REE and FAO/WHO/UNU was near zero (9 kcal) but the variation in both directions was substantial (range 591 to þ573 kcal). Similarly, the differences of areas of measured REE and HB demonstrated a mean of 1.9 kcal and range 568 to þ571 kcal. Glasgow coma score, patient status, weight and body temperature were signi?cant predictors of measured REE (p < 0.001; R²⁼ 0.47).

Conclusions: Clinical equations are poor predictors of measured REE in patients with TBI. The variability in REE is substantial. Clinicians should be aware of the limitations of prediction equations when estimating energy requirements in TBI patients.

Lowest energy potentials of Li2 + F accessible at ultracold temperatures

Ab initio calculations for the strongly exoergic Li2+F harpoon reaction are presented using density-functional theory, complete active space self-consistent field, and multireference configuration interaction methods to argue that this reaction would be an ideal candidate for investigation with ultracold molecules. The lowest six states are calculated with the aug-correlation-consistent polarized valence triple-zeta basis set and at least two can be accessed by a ground rovibronic Li2 molecule with zero collision energy at all reaction geometries. The large reactive cross section (characteristic of harpoon reactions) and chemiluminescent products are additional attractive features of these reactions.

Collective excitations of trapped Bose-Einstein condensates in the energy and time domains

A time-dependent method for calculating the collective excitation frequencies and densities of a trapped, inhomogeneous Bose-Einstein condensate with circulation is presented. The results are compared with time-independent solutions of the Bogoliubov-de Gennes equations. The method is based on time-dependent linear-response theory combined with spectral analysis of moments of the excitation modes of interest. The technique is straightforward to apply, extremely efficient in our implementation with parallel fast Fourier transform methods, and produces highly accurate results. For high dimensionality or low symmetry the time-dependent approach is a more practical computational scheme and produces accurate and reliable data. The method is suitable for general trap geometries, condensate flows and condensates permeated with defects and vortex structures.

An understanding of chemoselective hydrogenation on crotonaldehyde over Pt(111) in the free energy landscape: The microkinetics study based on first-principles calculations

Microkinetic model is developed in the free energy landscape based on density functional theory (DFT) to quantitatively investigate the reaction mechanism of chemoselective partial hydrogenation of crotonaldehyde to crotyl alcohol over Pt(1 1 1) at the temperature of 353 K. Three different methods (mobile, immobile and collision theory models) were carried out to obtain free energy barrier of adsorption/desorption processes. The results from mobile and collision theory models are similar. The calculated TOFs from both models are close to the experiment value. However, for the immobile model, in which the free energy barrier of desorption approaches the energy barrier, the calculated TOF is 2 orders of magnitude lower than the other models. The difficulty of adsorption/ desorption may be overestimated in the immobile model. In addition, detailed analyses show that for the surface hydrogenation elementary steps, the entropy and internal energy effects are small under the reaction condition, while the zero-point-energy (ZPE) correction is significant, especially for the multi-step hydrogenation reaction. The total energy with the ZPE correction approaches to the full free energy calculation for the surface reaction under the reaction condition. (c) 2011 Elsevier B.V. All rights reserved.