Half-explicit Runge-Kutta methods for semi-explicit differential-algebraic equations of index 1

Summary For the numerical solution of non-stiff semi-explicit differentialalgebraic equations (DAEs) of index 1 half-explicit Runge-Kutta methods (HERK) are considered that combine an explicit Runge-Kutta method for the differential part with a simplified Newton method for the (approximate) solution of the algebraic part of the DAE. Two principles for the choice of the initial guesses and the number of Newton steps at each stage are given that allow to construct HERK of the same order as the underlying explicit Runge-Kutta method. Numerical tests illustrate the efficiency of these methods.     

Non-adiabatic calculation of muonic atom-nucleus collisions

In order to calculate the cross sections of the muonic atom-nucleus collisions, we have proposed a precise numerical method, a non-adiabatic coupled-rearrangement-channel method with the use of the Jacobian coordinates for the three-body system in the whole space. The scattering boundary condition is correctly imposed along the coordinates; this method does not suffer from the well known defects of the method of adiabatic representation. We applied our method to the muonic atom-nucleus collisions for an incident c.m. energies of 0.001–100 eV.     

Conversion estimates for gas-solid reactions: Improved estimates for irreversible systems

Solid conversion has been of much interest to research workers and the Cumulative Gas Concentration Method (CGCM) has been developed to provide bounds for the reactant solid in irreversible systems having power-law type kinetics. This CGCM is used in conjunction with the Pseudo-Steady-State (PSS) assumption and these bounds appear to be good under certain conditions and for long times.

In this work, the PSS condition is not assumed when deriving the bounds on the cumulative gas concentration and the effect of the inclusion of the transient part is analysed. A technique, based on an iterative scheme incorporating the CGCM, is developed to provide improved estimates for the reactant solid.     

A high-gain porous silicon metal–semiconductor–metal photodetector through rapid thermal oxidation and rapid thermal annealing

2 tungsten lamp at 1.83 mA excited by a 0.85-mW 675-nm laser diode. We also demonstrate that the dark current could be greatly reduced through rapid thermal oxidation and rapid thermal annealing. Received: 4 August 1997/Accepted: 28 May 1998     

Investigation by EPR and ENDOR spectroscopy of the novel 4Fe ferredoxin fromPyrococcus furiosus

The hyperthermophilic archaeonPyrococcus furiosus contains a four-Fe ferredoxin (Pf- Fd) that differs from most other 4Fe-Fd’s in that its [Fe₄S₄] cluster is anchored to protein by only three cysteinyl residues.Pf- Fd also is of interest because in its reduced form, [Fe₄S₄]⁺, the cluster exhibits bothS = 1/2 andS = 3/2 spin states. Addition of excess cyanide ion converts the cluster exclusively to anS = 1/2 state (g₁ = 2.09, g₂ = 1.95, g₃ = 1.92), however dialysis restores the EPR signal of native reduced protein indicating that the cluster is not irreversibly altered by cyanide. Both the native protein and protein in the presence of excess cyanide ion (Pf- Fd 4Fe-CN) were investigated here using the techniques of electron paramagnetic resonance (EPR) and electron-nuclear double resonance (ENDOR) spectroscopy. In particular,Pf- Fd 4Fe-CN was investigated using¹³CN^? and C¹⁵N^? ligands.¹³C and¹⁵N ENDOR indicated that a single cyanide ion bound directly, with the cluster showing an unusually small contact interaction (a_iso(¹³C)～ ?3 MHz, a_iso(¹⁵N) ～ 0). This is in contrast to cyanide bound to monomeric low-spin Fe(III)-containing proteins such as transferrin and myoglobin, for which the¹³C hyperfine coupling has a large isotropic component (a_iso(¹³C) ≈ ?30 MHz). This small contact interaction is not due to low spin density of Fe, as⁵⁷Fe ENDOR of the singly and triply labeledPf- Fd 4FeCN isotopologs, [⁵⁷FeFe₃S₄]⁺ and [Fe⁵⁷Fe₃S₄]⁺, show hyperfine coupling characteristic for [Fe₄S₄]⁺ clusters, particularly for the Fe to which cyanide binds. Thus, the low spin density on¹³C is not due to low spin density on the Fe ion to which it binds. Further theoretical work is needed to explain the contrast between the strong electronic effect of cyanide ion binding with the low spin density on the ligand.     

CONFLICT AND COOPERATION IN UTILIZING A COMMON PROPERTY RESOURCE

Noncooperative games are used to demonstrate that, while free riding is always an option, the “tragedy of the commons” is not inevitable. When the decision to cooperate or free ride is considered in a dynamic setting, there is no intermediate case where some cooperate and others free ride. The game is only stabilized in either full cooperation or full defection. The important factor in obtaining a cooperative outcome is the critical number of players that decide to cooperate. The concept of commitment is used to demonstrate the necessary and sufficient conditions for full cooperation. Although the analysis is based on a shared water resource, it can be extended to other natural resources with common costs and private benefits, especially in the cases where there are no international authorities or treaties that internalize the externalities involved in privatizing the resource.     

Superconductivity and normal state resistivity of Ba0.6K0.4BiO3: an optical phonon approach

We discuss the nature of the pairing mechanism and the physical properties associated with the normal as well as the superconducting state of cubic perovskites Ba_0.6K_0.4BiO₃using the strong coupling theory. An interaction potential which includes the Coulomb, electron–optical phonon and electron–plasmon interactions is developed to elucidate the superconducting state. A model dielectric function is constructed with these interactions fulfilling thef-sum rule. The screening parameter (μ^* = 0.26) infers the poor screening of charge carriers. The electron–optical phonon strength (λ) estimated as 0.98 is consistent with an attractive electron–electron interaction and supports the moderate to strong coupling theory. The superconducting transition temperature of Ba_0.6K_0.4BiO₃is then estimated as 32 K. Ziman's formula of resistivity is employed to analyse and compare this with the temperature-dependent resistivity of a single crystal. The estimated contribution from the electron–optical phonon together with the residual resistivity clearly infers a difference when a comparison is made with experimental data. The subtracted data infer a quadratic temperature dependence in the temperature domain (30 ≤ T ≤ 200 K). The quadratic temperature dependence of ρ [ = ρ_exp − (ρ₀ + ρ_e–ph)] is understood in terms of 3D electron–electron inelastic scattering. The presence of these el–el and el–ph interactions allows a coherent interpretation of the physical properties. Analysis reveals that a moderate to strong coupling exists in the Ba_0.6K_0.4BiO₃system and the coupling of electrons with the high-energy optical phonons of the oxygen breathing mode will be a reason for superconductivity. The implications of the above analysis are discussed.     

Plasmon resonances in a spherical sodium cluster and in a flat surface with a soft optical edge

Received: 27 August 1997/Revised version: 7 November 1997