Kinetics and equilibrium of small metallic clusters: Ab initio confinement molecular dynamics study of 4

The ab initio molecular dynamics (AIMD) [1] is combined with the heuristic, successive confinement method of surveying a potential energy surface (PES) [2], thereby offering a framework for the simulation study of kinetics and equilibrium properties of metallic clusters. This approach is applied to the study of Au₄, a cluster possessing a simple but specific PES, which consists of very shallow and deep basins and due to this presents a challenge to the conventional AIMD methods. Among other things, the probabilities of the transitions between isomers have been found, and on this basis, both the time-dependent and equilibrium populations of the isomers have been calculated for the conditions typical of the NeNePo experiments [3] in the femtosecond pump-probe spectroscopy.     

An Extension of the Coordinate Transformation Method for Open-Loop Nash Equilibria

In Leitmann (Ref. 1), a coordinate transformation method was introduced to obtain global solutions for free problems in the calculus of variations. This direct method was extended and broadened in Carlson (Ref. 2) and later in Leitmann (Ref. 3). The applicability of the original work of Leitmann (Ref. 1) was further developed in Dockner and Leitmann (Ref. 4) to include the class of open-loop dynamic games. In the present work, we improve the results of Ref. 4 in two directions. First, we enlarge the class of open-loop dynamic games to permit coupling among the dynamic equations via the states of the players; second, we incorporate the modifications given in Refs. 2 and 3. Our results greatly increase the applicability of this method. An example arising from the harvesting of a renewable resource is presented to illustrate the utility of our results.     

Numerical analysis of reversible A B C reaction-diffusion systems

We develop an effective numerical method of studying large-time properties of reversible reaction-diffusion systems of type A + B ↔ C with initially separated reactants. Using it we find that there are three types of asymptotic reaction zones. In particular we show that the reaction rate can be locally negative and concentrations of species A and B can be nonmonotonic functions of the space coordinate x, locally significantly exceeding their initial values. Received 6 June 2002 / Received in final form 20 January 2003 Published online 7 May 2003     

Permeability,permselectivity, and penetrant‐induced plasticization in fluorinated polyimides studied by positron lifetime measurements

The ortho‐positronium (o‐Ps) lifetime τ₃ and its intensity I₃ in various fluorinated polyimides were determined by the positron annihilation technique and were studied with the spin–lattice relaxation time T₁ and the propylene permeability, solubility, diffusivity, and permselectivity for propylene/propane in them. τ₃, I₃, and the distribution of τ₃ changed when the bulky moieties in the polyimides were changed. The polyimides, having both large τ₃ and I₃ values, exhibited a short T₁ and a high permeability with a low permselectivity. The propylene permeability and diffusivity were exponentially correlated with the product of I₃ and the average free‐volume hole size estimated from τ₃. In highly plasticized states induced by the sorption of propylene, the permeability increased with the propylene pressure in excellent agreement with the change in the free‐volume hole properties probed by o‐Ps. The large and broad distribution of the free‐volume holes and increased local chain mobility for the 2,2‐bis(3,4‐decarboxyphenyl) hexafluoropropane dianhydride‐based polyimides are thought to be important physical properties for promoting penetrant‐induced plasticization. These results suggest that o‐Ps is a powerful probe of not only the free‐volume holes but also the corresponding permeation mechanism and penetrant‐induced plasticization phenomenon. © 2003 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 41: 308–318, 2003     

The modified characteristic finite difference fractional steps method for the coupled system of fluid dynamics in porous media and its analysis

For the coupled system of multilayer fluid dynamics in porous media, the modified characteristic finite difference fractional steps method applicable to parallel arithmetic is put forward and two‐dimensional and three‐dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, energy method, piecewise biquadratic interpolation, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in L² norm are derived to determine the error in the approximate solution. This method has already been applied to the numerical simulation of multilayer fluid dynamics in porous media. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 665–681, 2003.     

The dynamics of a two-component fluid in the presence of capillary forces

In the present paper we study the qualitative behavior ast→∞ of the solution of the Cauchy problem for a system of equations describing a dynamics of a two-component viscous fluid. The model under consideration takes into account the mutual diffusion of the fluid components as well as their capillary interaction. We describe the ω-limit set of trajectories of the dynamical system generated by the problem. It is proved that the stationary solution of the problem, is a homogeneous stationary distribution of one of the components, is asymptotically stable. Any other stationary solution is not asymptotically stable and is even unstable if there are no close stationary solutions corresponding to a smaller energy level. Translated fromMatematicheskie Zametki, Vol. 62, No. 2, pp. 293–305, August, 1997. Translated by A. M. Chebotarev     

AN ENCOMPASSING APPROACH TO MODELING FISHERY DYNAMICS: MODELING DYNAMIC NONLINEAR SYSTEMS

Uncertainty about the nature and significance of nonlinearities and the manner in which dynamics affect future realizations makes model specification the most difficult aspect of modeling dynamic systems. By interpreting several popular fishery models as subcases of a nesting dynamic Taylor series approximation, we isolate the specification differences between these models in a way that accounts for commonalities. On the argument that the differences due to alternative nonlinear forms are likely to be small compared to more mundane considerations such as delay difference and general dynamic lag specification, we propose an alternative model that uses the terms from the first order approximation common to all models combined with a data-based determination of the appropriate lags using the methods of state space time series analysis. Finally, the success of the alternative models is judged in an application to Pacific halibut data.     

Low-energy ion-assisted control of interfacial structures in metallic multilayers

A molecular dynamics method has been used to simulate the argon ion-assisted deposition of Cu/Co/Cu multilayers and to explore ion beam assistance strategies that can be used during or after the growth of each layer to control interfacial structures. A low-argon ion energy of 5–10 eV was found to minimize a combination of interfacial roughness and interlayer mixing (alloying) during the ion-assisted deposition of multilayers. However, complete flattening with simultaneous ion assistance could not be achieved without some mixing between the layers when a constant ion energy approach was used. It was found that multilayers with lower interfacial roughness and intermixing could be grown either by modulating the ion energy during the growth of each metal layer or by utilizing ion assistance only after the completion of each layers deposition. In these latter approaches, relatively high-energy ions could be used since the interface is buried and less susceptible to intermixing. The interlayer mixing dependence upon the thickness of the over layer has been determined as a function of ion energy.     

Particle multiplicities in the fission-like reactions of 340 MeV Si on Th

Pre-scission and post-scission multiplicities of neutrons and alpha particles have been simultaneously measured for the fission-like reactions of 340 MeV ²⁸Si on ²³²Th. Dynamical model calculations using HICOL code predict that about 90% of the observed events are of quasi-fission type while the remaining 10% are from compound nucleus fission decay. Moving source fits were carried out to the observed neutron and alpha particle spectra, measured at different angles with respect to the fragment directions. The pre-scission and post-scission neutron multiplicities are deduced to be 8.7±2.0 and 9.4±2.0, respectively. The corresponding multiplicity values for alpha particles are found to be 0.22±0.08 and 0.1±0.03. From the measured post-scission neutron multiplicity, it is inferred that about 65±20 MeV of the initial excitation energy remains at scission. This may be compared to the value of 85±30 MeV estimated from PACE2 statistical model calculations, adjusted to reproduce the measured pre-scission neutron multiplicity. From a comparison of the Statistical Model predictions with the measured pre-scission neutron multiplicity, the fission delay is estimated to be of 5⁺⁷₋₃×10⁻²⁰ s which overlaps with the average duration of fission-like process from the contact to the scission point (2×10⁻²⁰ s) as determined from HICOL-based dynamical calculations. For the delay time deduced as above, the pre-scission alpha particle multiplicity calculated by the PACE2 code is about a factor two larger than the experimental one, demonstrating the difficulties in modelling the alpha particle emission from highly elongated shapes that characterize the fissioning system from the contact point to scission.