An application of quasi-Newton methods for the numerical solution of interface problems

Quasi-Newton methods are applied to solve interface problems which arise from domain decomposition methods. These interface problems are usually sparse systems of linear or nonlinear equations. We are interested in applying these methods to systems of linear equations where we are not able or willing to calculate the Jacobian matrices as well as to systems of nonlinear equations resulting from nonlinear elliptic problems in the context of domain decomposition. Suitability for parallel implementation of these algorithms on coarse-grained parallel computers is discussed.

Neutron quasi-elastic scattering in disordered solids: a Monte Carlo study of metal-hydrogen systems

The dynamic structure factor of neutron quasi-elastic scattering has been calculated by Monte Carlo methods for atoms diffusing on a disordered lattice. The disorder includes not only variation in the distances between neighbouring atomic sites but also variation in the hopping rate associated with each site. The presence of the disorder, particularly the hopping rate disorder, causes changes in the time-dependent intermediate scattering function which translate into a significant increase in the intensity in the wings of the quasi-elastic spectrum as compared with the Lorentzian form. The effect is particularly marked at high values of the momentum transfer and at site occupancies of the order of unity. The MC calculations demonstrate how the degree of disorder may be derived from experimental measurements of the quasi-elastic scattering. The model structure factors are compared with the experimental quasi-elastic spectrum of an amorphous metal-hydrogen alloy.

Monitoring of a quasi-stationary eddy in the Bay of Biscay by means of satellite, in situ and model results

The presence of a quasi-stationary anticyclonic eddy within the southeastern Bay of Biscay (centred around 44°30′N-4°W) has been reported on various occasions in the bibliography. The analysis made in this study for the period 2003–2010, by using in situ and remote sensing measurements and model results shows that this mesoscale coherent structure is present almost every year from the end of winter-beginning of spring, to the beginning of fall. During this period it remains in an area limited to the east by the Landes Plateau, to the west by Le Danois Bank and Torrelavega canyon and to the northwest by the Jovellanos seamount. All the observations and analysis made in this contribution, suggest that this structure is generated between Capbreton and Torrelavega canyons. Detailed monitoring from in situ and remote sensing data of an anticyclonic quasi-stationary eddy, in 2008, shows the origin of this structure from a warm water current located around 43°42′N-3°30′W in mid-January. This coherent structure is monitored until August around the same area, where it has a marked influence on the Sea Level Anomaly, Sea Surface Temperature and surface Chlorophyll-a concentration. An eddy tracking method, applied to the outputs of a numerical model, shows that the model is able to reproduce this type of eddy, with similar 2D characteristics and lifetimes to that suggested by the observations and previous works. This is the case, for instance, of the simulated MAY04 eddy, which was generated in May 2004 around Torrelavega canyon and remained quasi-stationary in the area for 4 months. The diameter of this eddy ranged from 40 to 60 km, its azimuthal velocity was less than 20 cm s−1, its vertical extension reached 3000–3500 m depth during April and May and it was observed to interact with other coherent structures.

Monitoring of a quasi-stationary eddy in the Bay of Biscay by means of satellite, in situ and model results

The presence of a quasi-stationary anticyclonic eddy within the southeastern Bay of Biscay (centred around 44°30′N-4°W) has been reported on various occasions in the bibliography. The analysis made in this study for the period 2003–2010, by using in situ and remote sensing measurements and model results shows that this mesoscale coherent structure is present almost every year from the end of winter-beginning of spring, to the beginning of fall. During this period it remains in an area limited to the east by the Landes Plateau, to the west by Le Danois Bank and Torrelavega canyon and to the northwest by the Jovellanos seamount. All the observations and analysis made in this contribution, suggest that this structure is generated between Capbreton and Torrelavega canyons. Detailed monitoring from in situ and remote sensing data of an anticyclonic quasi-stationary eddy, in 2008, shows the origin of this structure from a warm water current located around 43°42′N-3°30′W in mid-January. This coherent structure is monitored until August around the same area, where it has a marked influence on the Sea Level Anomaly, Sea Surface Temperature and surface Chlorophyll-a concentration. An eddy tracking method, applied to the outputs of a numerical model, shows that the model is able to reproduce this type of eddy, with similar 2D characteristics and lifetimes to that suggested by the observations and previous works. This is the case, for instance, of the simulated MAY04 eddy, which was generated in May 2004 around Torrelavega canyon and remained quasi-stationary in the area for 4 months. The diameter of this eddy ranged from 40 to 60 km, its azimuthal velocity was less than 20 cm s−1, its vertical extension reached 3000–3500 m depth during April and May and it was observed to interact with other coherent structures.

The impact of nurse-directed protocolised-weaning from mechanical ventilation on nursing practice: A quasi-experimental study

Background:
Internationally, nurse-directed protocolised-weaning has been evaluated by measuring its impact on patient outcomes. The impact on nurses’ views and perceptions has been largely ignored.

Aim:
To determine the change in intensive care nurses’ perceptions, satisfaction, knowledge and attitudes following the introduction of nurse-directed weaning. Additionally, views were obtained on how useful protocolised-weaning was to practice.

Methods:
The sample comprised nurses working in general intensive care units in three university-affiliated hospitals. Nurse-directed protocolised-weaning was implemented in one unit (intervention group); two ICUs continued with usual doctor-led practice (control group). Nurses’ perceptions, satisfaction, knowledge and attitudes were measured by self-completed questionnaires before (Phase I) and after the implementation of nurse-directed weaning (Phase II) in all units.

Results:
Response rates were 79% (n=140n=140) for Phase 1 and 62% (n=132n=132) for Phase II. Regression-based analyses showed that changes from Phase I to Phase II were not significantly different between the intervention and control groups. Sixty-nine nurses responded to both Phase I and II questionnaires. In the intervention group, these nurses scored their mean perceived level of knowledge higher in Phase II (6.39 vs 7.17, p=0.01p=0.01). In the control group, role perception (4.41 vs 4.22, p=0.01p=0.01) was lower and, perceived knowledge (6.03 vs 6.63, p=0.04p=0.04), awareness of weaning plans (6.09 vs 7.06, p=0.01p=0.01) and satisfaction with communication (5.28 vs 6.19, p=0.01p=0.01) were higher in Phase II. The intervention group found protocolised weaning useful in their practice (75%): this was scored significantly higher by junior and senior nurses than middle grade nurses (p=0.02p=0.02).

Conclusion

We conclude that nurse-directed protocolised-weaning had no effect on nurses’ views and perceptions due to the high level of satisfaction which encouraged nurses’ participation in weaning throughout. Control group changes are attributed to a ‘reactive effect’ from being study participants. Weaning protocols provide a uniform method of weaning practice and are particularly beneficial in providing safe guidance for junior staff.

Ice-binding proteins adsorb to their ligand via anchored clathrate waters

The main success of my thesis has been to establish the mechanism by which antifreeze proteins (AFPs) bind irreversibly to ice crystals, and hence prevent their growth. AFPs organize ice-like water on their ice-binding site, which then merges and freezes with the quasi-liquid layer of ice. This was revealed from studying the exceptionally large (ca. 1.5-MDa) Ca 2+-dependent AFP from the Antarctic bacterium Marinomonas primoryensis (MpAFP). The 34-kDa antifreeze- active region of MpAFP was predicted to fold as a novel Ca 2+-binding β-helix. Site-directed mutagenesis confirmed the model and demonstrated that its ice-binding site (IBS) consisted of solvent-exposed Thr and Asx parallel arrays on the Ca 2+-binding turns. The X-ray crystal structure of the antifreeze region was solved to a resolution of 1.7 Å. Two of the four molecules within the unit cell of the crystal had portions of their IBSs freely exposed to solvent. Identical clathrate-like cages of water molecules were present on each IBS. These waters were organized by the hydrophobic effect and anchored to the protein via hydrogen bonds. They matched the spacing of water molecules in an ice lattice, demonstrating that anchored clathrate waters bind AFPs to ice. This mechanism was extended to other AFPs including the globular type III AFP from fishes. Site-directed mutagenesis and a modified ice-etching technique demonstrated this protein uses a compound ice-binding site, comprised of two flat and relatively hydrophobic surfaces, to bind at least two planes of ice. Reinvestigation of several crystal structures of type III AFP identified anchored clathrate waters on the solvent-exposed portion of its compound IBS that matched the spacing of waters on the primary prism plane of ice. Ice nucleation proteins (INPs), which can raise the temperature at which ice forms in solution to just slightly below 0oC, have the opposite effect to AFPs. A novel dimeric β-helical model was proposed for the INP produced by the bacterium Pseudomonas borealis. Molecular dynamics simulations showed that INPs are also capable of ordering water molecules into an ice- like lattice. However, their multimerization brings together sufficient ordered waters to form an ice nucleus and initiate freezing.

Three-dimensional photonic crystals based on macroporous silicon with modulated pore diameter

We report on the fabrication and optical characterization of a three-dimensional (3D) photonic crystal on the basis of macroporous silicon. The structure consists of a 2D array of air pores in silicon whose diameter is varied (modulated) periodically with depth. The bandstructure of the resulting 3D hexagonal photonic crystal is calculated and compared with transmission measurements. The described structure allows to adjust the dispersion relation along the pore axis almost independently from the dispersion relation in the plane perpendicular to the pore axis.

Optical bistability in nonlinear surface-plasmon polaritonic crystals

Nonlinear optical transmission through periodically nanostructured metal films (surface-plasmon polaritonic crystals) has been studied. The surface polaritonic crystals have been coated with a nonlinear polymer. The optical transmission of such nanostructures has been shown to depend on the control-light illumination conditions. The resonant transmission exhibits bistable behavior with the control-light intensity. The bistability is different at different resonant signal wavelengths and for different wavelengths of the control light. The effect is explained by the strong sensitivity of the surface-plasmon mode resonances at the signal wavelength to the surrounding dielectric environment and the electromagnetic field enhancement due to plasmonic excitations at the controlled light wavelengths.

Ferroelectricity and isotope effects in hydrogen-bonded KDP crystals

Based on an accurate first principles description of the energetics in H-bonded potassium-dihydrogen-phosphate crystals, we conduct a first study of nuclear quantum effects and of the changes brought about by deuteration. Tunneling is allowed only for clusters involving correlated protons and heavy ion displacements, the main effect of deuteration being a depletion of the proton probability density at the O-H-O bridge center, which in turn weakens its proton-mediated covalent bonding. The ensuing lattice expansion couples self-consistently with the proton off-centering, thus explaining both the giant isotope effect and its close connection with geometrical effects.

Interplay between proton ordering and ferroelectric polarization in H-bonded KDP-type crystals

The origin of ferroelectricity in KH2PO4 (KDP) is studied by first-principles electronic structure calculations. In the low-temperature phase, the collective off-centre ordering of the protons is accompanied by an electronic charge delocalization from the near and localization at the far oxygen within the O-H...O bonds. Electrostatic forces. then, push the K+ ions towards off-centre positions, and induce a macroscopic polarization. The analysis of the correlation between different geometrical and electronic quantities, in connection with experimental data. supports the idea that the role of tunnelling in isotopic effects is irrelevant. Instead, geometrical quantum effects appear to play a central role. (C) 2001 Elsevier Science B.V. All rights reserved.

First-principles study of ferroelectricity and isotope effects in H-bonded KH2PO4 crystals

By means of extensive first-principles calculations we studied the ferroelectric phase transition and the associated isotope effect in KH2PO4 (KDP). Our calculations revealed that the spontaneous polarization of the ferroelectric phase is due to electronic charge redistributions and ionic displacements which are a consequence of proton ordering, and not vice versa. The experimentally observed double-peaked proton distribution in the paraelectric phase cannot be explained by a dynamics of only protons. This requires, instead, collective displacements within clusters that include also the heavier ions. These tunneling clusters can explain the recent evidence of tunneling obtained from Compton scattering measurements. The sole effect of mass change upon deuteration is not sufficient to explain the huge isotope effect. Instead, we find that structural modifications deeply connected with the chemistry of the H bonds produce a feedback effect on tunneling that strongly enhances the phenomenon. The resulting influence of the geometric changes on the isotope effect agrees with experimental data from neutron scattering. Calculations under pressure allowed us to analyze the issue of universality in the disappearance of ferroelectricity upon compression. Compressing DKDP so that the distance between the two peaks in the deuteron distribution is the same as for protons in KDP, corresponds to a modification of the underlying double-well potential, which becomes 23 meV shallower. This energy difference is what is required to modify the O-O distance in such a way as to have the same distribution for protons and deuterons. At the high pressures required experimentally, the above feedback mechanism is crucial to explain the magnitude of the geometrical effect.

Luminous blue variables as the progenitors of supernovae with quasi-periodic radio modulations

The interaction between supernova ejecta and circumstellar matter, arising from previous episodes of mass loss, provides us with a means of constraining the progenitors of supernovae. Radio observations of a number of supernovae show quasi-periodic deviations from a strict power-law decline at late times. Although several possibilities have been put forward to explain these modulations, no single explanation has proven to be entirely satisfactory. Here we suggest that Luminous blue variables undergoing S-Doradus type variations give rise to enhanced phases of mass loss that are imprinted on the immediate environment of the exploding star as a series of density enhancements. The variations in mass loss arise from changes in the ionization balance of Fe, the dominant ion that drives the wind. With this idea, we find that both the recurrence timescale of the variability and the amplitude of the modulations are in line with the observations. Our scenario thus provides a natural, single-star explanation for the observed behaviour that is, in fact, expected on theoretical grounds.