Numerical integration schemes

A. Ahmad and L. Cohen. A numerical integration scheme for the A -body gravitational problem. J. Comp. Phys., 12 389-402, 1973. 

We focus on so-called symplectic methods [18] for the following reason It has been shown that the preservation of the symplectic structure of phase space under a numerical integration scheme implies a number of very desirable properties. Namely,... 

Multicenter Numerical Integration Scheme for Polyatomic Molecules... 

A numerical integration scheme has produced the following results ... 

Becke, A. D., 1988c, A Multicenter Numerical Integration Scheme for Polyatomic Molecules , J. Chem. Phys.,... 

Use a spreadsheet or equivalent computer program to calculate the concentration of product C as the reaction proceeds with time (/) in a constant-volume batch reactor (try the parameter values supplied below). You may use a simple numerical integration scheme such as Acc =... 

An alternative procedure consists in using a numerical integration scheme to evaluate the exchange-correlation contribution. In this case, no auxiliary basis set is needed for the exchange-correlation terms, and numerically more reliable results can be obtained. 

A. EULER ALGORITHM. The simplest possible numerical-integration scheme (and the most useful) is Euler (pronounced oiler ), illustrated in Fig. 4.7. Assume we wish to solve the ODE... 

Discretizing the population balance in K+1 grid points results in K ordinary differential equations as the population density n(LQ,t) is determined by the algebraic relation Equation 4. The equation for dn(LQ,t)/dt is therefore not required. In the overal model of the crys lllzer, the population density must be integrated in the calculation of c(t) and In the calculation of m (t) in the nucleatlon rate. As the population balance is discretized these integrals have to be replaced by numerical integration schemes. 

Using a Langevin dynamics approach, the stochastic LLG equation [Eq. (3.46)] can be integrated numerically, in the context of the Stratonovich stochastic calculus, by choosing an appropriate numerical integration scheme [51]. This method was first applied to the dynamics of noninteracting particles [51] and later also to interacting particle systems [13] (see Fig. 3.5). 

Treutler, O., Ahlrichs, R., 1995, Efficient Molecular Numerical Integration Schemes , /. Chem. Phys., 102,346. 

If we pick an arbitrary element we can see that it is represented by the xy-coordinates of the four nodal points, as depicted in Fig. 9.16. The figure also shows a -coordinate system embedded within the element. In the r/, or local, coordinate system, we have a perfectly square element of area 2x2, where the element spreads between —1 > < 1 and — 1 > rj < 1. This attribute allows us to easily allows us to use Gauss quadrature as a numerical integration scheme, where the limits vary between -1 and 1. The isoparametric element described in the //-coordinate system is presented in Fig. 9.17. 

Numerical integration schemes allow an opportunity to test the numerical nonempirical pseudopotentials without their fit by analytical functions, which can lead to a considerable reduction in computational efforts. Employment of atomic pseudopotentials only at some selected atoms of a system while treating the rest all-electronically makes an impression of the consistency and reliability of such a combined approach. The results obtained for the MgO clusters embedded in some effective pseudopotential surroundings demonstrate a promise of the approach for compensation of broken bond effects. Specifically, the approach offers a tool for a substantial reduction of the artificially introduced nonequivalence of partial densities and the effective charges for atoms equivalent in the lattice. It is worth to mention that our approach can be modified further in many ways because numerical integration schemes can be easily applied/adapted even in those cases where the analytical methods become too complicated. 

Although H is conserved by the equations of motion, it clearly is not a Hamiltonian for Eqs. [65]. Since non-Hamiltonian systems tend to be more difficult than Hamiltonian systems to integrate stably numerically, the existence of a conserved energy for a non-Hamiltonian system is of vital importance as a check on the stability of the numerical integration scheme employed. 

The electronic state calculation by discrete variational (DV) Xa molecular orbital method is introduced to demonstrate the usefulness for theoretical analysis of electron and x-ray spectroscopies, as well as electron energy loss spectroscopy. For the evaluation of peak energy. Slater s transition state calculation is very efficient to include the orbital relaxation effect. The effects of spin polarization and of relativity are argued and are shown to be important in some cases. For the estimation of peak intensity, the first-principles calculation of dipole transition probability can easily be performed by the use of DV numerical integration scheme, to provide very good correspondence with experiment. The total density of states (DOS) or partial DOS is also useful for a rough estimation of the peak intensity. In addition, it is necessary lo use the realistic model cluster for the quantitative analysis. The... 

A. D. Becke (1988) A multicenter numerical integration scheme for polyatomic molecules. J. Chem. Phys. 88(4), pp. 2547-2553... 

This is the simplest possible numerical integration scheme. It is known as Euler s method. 

All the calculations were carried out using the Amsterdam Density Functional (ADF) code. Version 2.3 (Theoretical Chemistry, Vrije Universiteit, Amsterdam, The Netherlands), developed by Baerends et al. (41), which incorporates the relativistic extensions first proposed by Snijders et al. (42). The code was vectorized by Ravenek (43) and parallelized by Fonseca Guerra et al. (44), and the numerical integration scheme applied for the calculations was developed by te Velde et al. 

This chapter is devoted to numerical integration, and more specifically to the integration of rate expressions encountered in chemical kinetics. For simple cases, integration yields closed-form rate equations, while more complex reaction mechanisms can often be solved only by numerical means. Here we first use some simple reactions to develop and calibrate general numerical integration schemes that are readily applicable to a spreadsheet. We then illustrate several non-trivial applications, including catalytic reactions and the Lotka oscillator. 

The solution of Equations 47, 48, and 49 requires numerical techniques. For such nonlinear equations, it is usually wise to employ a simple numerical integration scheme which is easily understood and pay the price of increased computational time for execution rather than using a complex, efficient, numerical integration scheme where unstable behavior is a distinct possibility. A variety of simple methods are available for integrating a set of ordinary first order differential equations. In particular, the method of Huen, described in Ref. 65, is effective and stable. It is self-starting and consists of a predictor and a corrector step. Let y = f(t,y) be the vector differential equation and let h be the step size. 

Robert, A., A stable numerical integration scheme for the primitive meteorological equations. Atmos Ocean 19, 35, 1981. 

A straightforward method bypasses the introduction of the auxiliary exchange-correlation fitting basis and evaluates those matrix elements directly by three dimensional numerical integration. This not only spares a computational step, it also avoids the limitations of the fitting basis set. There still arises a truncation error due to the numerical integration scheme. 

---