Product Integration Rules at Clenshaw-Curtis and Related Points: A Robust Implementation

Product integration rules generalizing the Fej?r, Clenshaw-Curtis,and Filippi quadrature rules respectively are derived for integralswith trigonometric and hyperbolic weight factors. The Chebyshevmoments of the weight functions are found to be given by well-conditionedexpressions, in terms of hypergeometric functions ₀F₁. An a priori error estimator is discussed which is shown bothto avoid wasteful invocation of the integration rule and toincrease significantly the robustness of the automatic quadratureprocedure. Then, specializing to extended Clenshaw-Curtis (ECC) rules,three types of a posteriori error estimates are considered andthe existence of a great risk of their failure is demonstratedby large scale validation tests. An empirical error estimator,superseding them for slowly varying integrands, is found toresult in a spectacular increase in the output reliability. Finally, enhancements in the control of the interval subdivisionstrategy aiming at increasing code robustness is discussed.Comparison with the code DQAWO of QUADPACK, with about a hundredthousand solved integrals, is illustrative of the increasedrobustness and error estimate reliability of our computer codeimplementation of the ECC rules.     

A First-Passage Kinetic Monte Carlo algorithm for complex diffusion–reaction systems

We develop an asynchronous event-driven First-Passage Kinetic Monte Carlo (FPKMC) algorithm for continuous time and space systems involving multiple diffusing and reacting species of spherical particles in two and three dimensions. The FPKMC algorithm presented here is based on the method introduced in Oppelstrup et al. [10] and is implemented in a robust and flexible framework. Unlike standard KMC algorithms such as the n-fold algorithm, FPKMC is most efficient at low densities where it replaces the many small hops needed for reactants to find each other with large first-passage hops sampled from exact time-dependent Green’s functions, without sacrificing accuracy. We describe in detail the key components of the algorithm, including the event-loop and the sampling of first-passage probability distributions, and demonstrate the accuracy of the new method. We apply the FPKMC algorithm to the challenging problem of simulation of long-term irradiation of metals, relevant to the performance and aging of nuclear materials in current and future nuclear power plants. The problem of radiation damage spans many decades of time-scales, from picosecond spikes caused by primary cascades, to years of slow damage annealing and microstructure evolution. Our implementation of the FPKMC algorithm has been able to simulate the irradiation of a metal sample for durations that are orders of magnitude longer than any previous simulations using the standard Object KMC or more recent asynchronous algorithms.     

Finite generation of powers of ideals

Suppose is a maximal ideal of a commutative integral domain and that some power of is finitely generated. We show that is finitely generated in each of the following cases: (i) is of height one, (ii) is integrally closed and , (iii) is a monoid domain over a field , where is a cancellative torsion-free monoid such that , and is the maximal ideal . We extend the above results to ideals of a reduced ring such that is Noetherian. We prove that a reduced ring is Noetherian if each prime ideal of has a power that is finitely generated. For each with , we establish existence of a -dimensional integral domain having a nonfinitely generated maximal ideal of height such that is -generated.

     

Thermodynamic Complexation of Dopamine with Zinc(II) in Media with Different Dielectric Constants

The complexation of zinc(II) with dopamine has been investigated by spectrophotometric measurements in mixed solvent system at an ionic strength of 0.2 mol•dm－3 sodium chloride, employed at at (15±0.1), (25±0.1), (35±0.1) ℃ at inin a pH range of ca. 6 to ca. 7 with a high ratio of ligand to metal. The effect of solvent systems on protonation and complexation are was was discussed. Linear relationships are werewere observed by plotting lg K versus 1/D, where K and D show stability and dielectric constants, respectively.