Correction Term (correction + term)

Distribution by Scientific Domains


Selected Abstracts


Estimating diesel degradation rates from N2, O2 and CO2 concentration versus depth data in a loamy sand

EUROPEAN JOURNAL OF SOIL SCIENCE, Issue 1 2007
J. Van De Steene
Summary The degradation rate of the pollutant is often an important parameter for designing and maintaining an active treatment system or for determining the rate of natural attenuation. A quasi-steady-state gas transport model based on Fick's law with a correction term for advective flux, for estimating diesel degradation rates from N2, O2 and CO2 concentration versus depth data, was evaluated in a laboratory column study. A loamy sand was spiked with diesel fuel at 0, 1000, 5000 and 10 000 mg kg,1 soil (dry weight basis) and incubated for 15 weeks. Soil gas was sampled weekly at 6 selected depths in the columns and analysed for O2, CO2 and N2 concentrations. The agreement between the measured and the modelled concentrations was good for the untreated soil (R2= 0.60) and very good for the soil spiked with 1000 mg kg,1 (R2= 0.96) and 5000 mg kg,1 (R2= 0.97). Oxygen consumption ranged from ,0.15 to ,2.25 mol O2 m,3 soil day,1 and CO2 production ranged from 0.20 to 2.07 mol CO2 m,3 soil day,1. A significantly greater mean O2 consumption (P < 0.001) and CO2 production (P < 0.005) over time was observed for the soils spiked with diesel compared with the untreated soil, which suggests biodegradation of the diesel substrate. Diesel degradation rates calculated from respiration data were 1.5,2.1 times less than the change in total petroleum hydrocarbon content. The inability of this study to correlate respiration data to actual changes in diesel concentration could be explained by volatilization, long-term sorption of diesel hydrocarbons to organic matter and incorporation of diesel hydrocarbons into microbial biomass, aspects of which require further investigation. [source]


Long-life torsion fatigue with normal mean stresses

FATIGUE & FRACTURE OF ENGINEERING MATERIALS AND STRUCTURES, Issue 4 2000
Marquis
Relatively simple fatigue tests have been performed on two common engineering materials, cast ductile iron and low-carbon steel, using two stress states, cyclic torsion and cyclic torsion with static axial and hoop stresses. Tests were designed to discriminate between normal stress and hydrostatic stress as the most suitable mean stress correction term for high cycle fatigue analysis. Microscopy shows that cracks in low-carbon steel nucleate and grow on maximum shear planes, while for cast iron pre-existing flaws grow on maximum normal stress planes. The data illustrate that tensile normal stress acting on a shear plane significantly reduced fatigue life and is an appropriate input for fatigue analysis of ductile materials. Static normal stresses did not significantly affect the fatigue life for the cast iron because the net mean stress on the maximum normal stress plane was zero. Mean torsion significantly reduced the fatigue strength of the cast iron. A critical plane long-life parameter for nodular iron which accounts for both stress state and mean stress is proposed, and is found to accurately correlate experimental data. [source]


Uniform Head in Horizontal and Vertical Wells

GROUND WATER, Issue 1 2006
David R. Steward
The steady-state head within a fully penetrating well may be estimated by evaluating the Thiem equation at the radius of the well. A method is presented here to extend results from the Thiem equation to horizontal wells and to partially penetrating wells. The particular model used in this investigation is based upon the analytic element method; it accurately reproduces a boundary condition of uniform head along the cylindrical surface at the perforated face of the well. This model is exercised over a representative range of parameters including the well's length, radius, and pumping rate, and the aquifer's hydraulic conductivity and thickness. Results are presented in a set of figures and tables that compare the well's drawdown to the drawdown that would have been obtained using the Thiem solution with the same pumping rate and radius. A methodology is presented to estimate the head within a horizontal or partially penetrating well by adding a correction term to results that can be readily obtained from computer models of vertical fully penetrating wells. This approach may also be used to contrast the differences in head between horizontal and vertical wells of various lengths, radii, and placement elevations. [source]


A hybrid Padé ADI scheme of higher-order for convection,diffusion problems

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 5 2010
Samir KaraaArticle first published online: 8 SEP 200
Abstract A high-order Padé alternating direction implicit (ADI) scheme is proposed for solving unsteady convection,diffusion problems. The scheme employs standard high-order Padé approximations for spatial first and second derivatives in the convection-diffusion equation. Linear multistep (LM) methods combined with the approximate factorization introduced by Beam and Warming (J. Comput. Phys. 1976; 22: 87,110) are applied for the time integration. The approximate factorization imposes a second-order temporal accuracy limitation on the ADI scheme independent of the accuracy of the LM method chosen for the time integration. To achieve a higher-order temporal accuracy, we introduce a correction term that reduces the splitting error. The resulting scheme is carried out by repeatedly solving a series of pentadiagonal linear systems producing a computationally cost effective solver. The effects of the approximate factorization and the correction term on the stability of the scheme are examined. A modified wave number analysis is performed to examine the dispersive and dissipative properties of the scheme. In contrast to the HOC-based schemes in which the phase and amplitude characteristics of a solution are altered by the variation of cell Reynolds number, the present scheme retains the characteristics of the modified wave numbers for spatial derivatives regardless of the magnitude of cell Reynolds number. The superiority of the proposed scheme compared with other high-order ADI schemes for solving unsteady convection-diffusion problems is discussed. A comparison of different time discretizations based on LM methods is given. Copyright © 2009 John Wiley & Sons, Ltd. [source]


Turbulence modelling of problem aerospace flows

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2006
Paul G. Tucker
Abstract Unsteady Reynolds averaged Navier,Stokes (URANS) and detached eddy simulation (DES) related approaches are considered for high angle of attack NACA0012 airfoil, wing,flap, generic tilt-rotor airfoil and double-delta geometry flows. These are all found to be problem flows for URANS models. For DES fifth-order upwinding is found too dissipative and the use of, for high speed flows, instability prone centred differencing essential. An existing hybrid ILES,RANS modelling approach, intended for flexible geometry, relatively high numerical dissipation codes is tested along with differential wall distance algorithms. The former gives promising results. The standard turbulence modelling approaches are found to give perhaps a surprising results variation. Results suggest that for the problem flows, the explicit algebraic stress and Menter shear stress transport (SST) URANS models are more accurate than the economical Spalart,Allmaras (SA). However, the explicit algebraic stress model (EASM) in its k,, form is impractically expensive to converge. Here, SA predictions lack a rotation correction term and this is likely to improve these results. Copyright © 2005 John Wiley & Sons, Ltd. [source]


An eddy viscosity model with near-wall modifications

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 9 2005
M. M. Rahman
Abstract An extended version of the isotropic k,, model is proposed that accounts for the distinct effects of low-Reynolds number (LRN) and wall proximity. It incorporates a near-wall correction term to amplify the level of dissipation in nonequilibrium flow regions, thus reducing the kinetic energy and length scale magnitudes to improve prediction of adverse pressure gradient flows, involving flow separation and reattachment. The eddy viscosity formulation maintains the positivity of normal Reynolds stresses and the Schwarz' inequality for turbulent shear stresses. The model coefficients/functions preserve the anisotropic characteristics of turbulence. The model is validated against a few flow cases, yielding predictions in good agreement with the direct numerical simulation (DNS) and experimental data. Comparisons indicate that the present model is a significant improvement over the standard eddy viscosity formulation. Copyright © 2005 John Wiley & Sons, Ltd. [source]


One-particle density matrix functional for correlation in molecular systems

INTERNATIONAL JOURNAL OF QUANTUM CHEMISTRY, Issue 6 2003
Mario Piris
Abstract Based on the analysis of the general properties for the one- and two-particle reduced density matrices, a new natural orbital functional is obtained. It is shown that by partitioning the two-particle reduced density matrix in an antisymmeterized product of one-particle reduced density matrices and a correction ,c we can derive a corrected Hartree,Fock theory. The spin structure of the correction term from the improved Bardeen,Cooper,Schrieffer theory is considered to take into account the correlation between pairs of electrons with antiparallel spins. The analysis affords a nonidempotent condition for the one-particle reduced density matrix. Test calculations of the correlation energy and the dipole moment of several molecules in the ground state demonstrate the reliability of the formalism. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 317,323, 2003 [source]


Flight control of a rotary wing UAV using backstepping

INTERNATIONAL JOURNAL OF ROBUST AND NONLINEAR CONTROL, Issue 6 2010
Bilal Ahmed
Abstract This paper presents a novel application of backstepping controller for autonomous landing of a rotary wing UAV (RUAV). This application, which holds good for the full flight envelope control, is an extension of a backstepping algorithm for general rigid body velocity control. The nonlinear RUAV model used in this paper includes the flapping and servo dynamics. The backstepping-based controller takes advantage of the ,decoupling' of the translation and rotation dynamics of the rigid body, resulting in a two-step procedure to obtain the RUAV control inputs. The first step is to compute desired thrusts and flapping angles to achieve the commanded position and the second step is to compute control inputs, which achieve the desired thrusts and flapping angles. This paper presents a detailed analysis of the inclusion of a flapping angle correction term in control. The performance of the proposed algorithm is tested using a high-fidelity RUAV simulation model. The RUAV simulation model is based on miniature rotorcraft parameters. The closed-loop response of the rotorcraft indicates that the desired position is achieved after a short transient. The Eagle RUAV control inputs, obtained using high-fidelity simulation results, clearly demonstrate that this algorithm can be implemented on practical RUAVs. Copyright © 2009 John Wiley & Sons, Ltd. [source]


Ab initio structure solution by iterative phase-retrieval methods: performance tests on charge flipping and low-density elimination

JOURNAL OF APPLIED CRYSTALLOGRAPHY, Issue 1 2010
Frank Fleischer
Comprehensive tests on the density-modification methods charge flipping [Oszlányi & Süt, (2004). Acta Cryst. A60, 134,141] and low-density elimination [Shiono & Woolfson (1992). Acta Cryst. A48, 451,456] for solving crystal structures are performed on simulated diffraction data of periodic structures and quasicrystals. A novel model-independent figure of merit, which characterizes the reliability of the retrieved phase of each reflection, is introduced and tested. The results of the performance tests show that the quality of the phase retrieval highly depends on the presence or absence of an inversion center and on the algorithm used for solving the structure. Charge flipping has a higher success rate for solving structures, while low-density elimination leads to a higher accuracy in phase retrieval. The best results can be obtained by combining the two methods, i.e. by solving a structure with charge flipping followed by a few cycles of low-density elimination. It is shown that these additional cycles dramatically improve the phases not only of the weak reflections but also of the strong ones. The results can be improved further by averaging the results of several runs and by applying a correction term that compensates for a reduction of the structure-factor amplitudes by averaging of inconsistently observed reflections. It is further shown that in most cases the retrieved phases converge to the best solution obtainable with a given method. [source]


On singular behaviors of impedance-based repeatable control for redundant robots

JOURNAL OF FIELD ROBOTICS (FORMERLY JOURNAL OF ROBOTIC SYSTEMS), Issue 4 2001
Chau-Chang Wang
This article addresses the association between the unstiffening phenomena in structural mechanics and the algorithmic singularities encountered in the impedance-based repeatable control algorithms used to command redundant manipulators. It is well known that velocity control schemes such as the pseudoinverse control schemes do not guarantee repeatability for redundant manipulators. In other words, for a closed end-effector trajectory, the joints do not, in general, exhibit a closed trajectory. One way to overcome this problem is to model each joint with compliance and incorporate a second-order correction term for the pseudoinverse. With this model, the joint configuration adopted by the manipulator at a given point in task space is one which minimizes the artificial potential energy of the system and is locally unique. In terms of statics, this is equivalent to saying that the elastic structure reaches its static equilibrium under external load. Keep this analogy in mind. We know that the impedance control commands the manipulator to mimic the behavior of an elastic articulated chain. For any phenomena observable on a real elastic structure, we should be able to find its counterpart embedded in the impedance control. In this article, we analyze the performance of such repeatable control algorithms from the point of view of structure mechanics. Singularities in the algorithm are examined and their significance in mechanics are also discussed. © 2001 John Wiley & Sons, Inc. [source]


On the cotangent cohomology of rational surface singularities with almost reduced fundamental cycle

MATHEMATISCHE NACHRICHTEN, Issue 11 2006
Trond Střlen Gustavsen
Abstract We prove dimension formulas for the cotangent spaces T1 and T2 for a class of rational surface singularities by calculating a correction term in the general dimension formulas. We get that it is zero if the dual graph of the rational surface singularity X does not contain a particular type of configurations, and this generalizes a result of Theo de Jong stating that the correction term c (X ) is zero for rational determinantal surface singularities. In particular our result implies that c (X ) is zero for Riemenschneiders quasi-determinantal rational surface singularities, and this also generalizes results for quotient singularities. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim) [source]


Effective series-resistance model of spiral inductors

MICROWAVE AND OPTICAL TECHNOLOGY LETTERS, Issue 2 2005
Heng-Ming Hsu
Abstract A new resistance model, which considers both the metal-skin effect and eddy-current losses of rectangular spiral-inductor structures, is proposed. Specifically, a correction term is added to the in-series resistance model so as to more accurately model the high-frequency behavior of metal losses. The accuracy of this analytical model is validated by comparing its simulated results with the measured ones. © 2005 Wiley Periodicals, Inc. Microwave Opt Technol Lett 46: 107,109, 2005; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/mop.20915 [source]


A comparative method to evaluate and validate stochastic parametrizations

THE QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY, Issue 642 2009
Leon Hermanson
Abstract There is a growing interest in using stochastic parametrizations in numerical weather and climate prediction models. Previously, Palmer (2001) outlined the issues that give rise to the need for a stochastic parametrization and the forms such a parametrization could take. In this article a method is presented that uses a comparison between a standard-resolution version and a high-resolution version of the same model to gain information relevant for a stochastic parametrization in that model. A correction term that could be used in a stochastic parametrization is derived from the thermodynamic equations of both models. The origin of the components of this term is discussed. It is found that the component related to unresolved wave,wave interactions is important and can act to compensate for large parametrized tendencies. The correction term is not proportional to the parametrized tendency. Finally, it is explained how the correction term could be used to give information about the shape of the random distribution to be used in a stochastic parametrization. Copyright © 2009 Royal Meteorological Society [source]


A Semiparametric Estimate of Treatment Effects with Censored Data

BIOMETRICS, Issue 3 2001
Ronghui Xu
Summary. A semiparametric estimate of an average regression effect with right-censored failure time data has recently been proposed under the Cox-type model where the regression effect ,(t) is allowed to vary with time. In this article, we derive a simple algebraic relationship between this average regression effect and a measurement of group differences in K -sample transformation models when the random error belongs to the Gp family of Harrington and Fleming (1982, Biometrika69, 553,566), the latter being equivalent to the conditional regression effect in a gamma frailty model. The models considered here are suitable for the attenuating hazard ratios that often arise in practice. The results reveal an interesting connection among the above three classes of models as alternatives to the proportional hazards assumption and add to our understanding of the behavior of the partial likelihood estimate under nonproportional hazards. The algebraic relationship provides a simple estimator under the transformation model. We develop a variance estimator based on the empirical influence function that is much easier to compute than the previously suggested resampling methods. When there is truncation in the right tail of the failure times, we propose a method of bias correction to improve the coverage properties of the confidence intervals. The estimate, its estimated variance, and the bias correction term can all be calculated with minor modifications to standard software for proportional hazards regression. [source]


Modeling of Kinetic Expressions for the Reduction of NOx by Hydrogen in Oxygen-Rich Exhausts Using a Gradient-Free Loop Reactor

CHEMICAL ENGINEERING & TECHNOLOGY (CET), Issue 6 2003
E. Frank
Abstract The reduction of NOx by hydrogen under lean conditions is investigated in a gradient-free loop reactor. Using this computer-controlled reactor, the reaction rates can be measured under exact isothermal conditions. Systematic variation of the input concentrations of hydrogen, nitric oxide, oxygen as well as reaction temperature provides a complete data set of reaction rates for the given reaction system. A number of kinetic rate expressions were evaluated for their ability to fit the experimental data by using toolboxes of MATLAB. The temperature influence on reaction rate constants and adsorption equilibrium constants were correlated simultaneously using Arrhenius and van't Hoff equations, respectively. The kinetic rate expression based on a Langmuir-Hinshelwood-type model describes the data and the model can be improved by introducing a correction term in square root of hydrogen partial pressure over the range of conditions investigated. [source]


Ab Initio Group Contribution Method for Activation Energies of Hydrogen Abstraction Reactions

CHEMPHYSCHEM, Issue 1 2006
Mark Saeys Prof.
Abstract The group contribution method for activation energies is applied to hydrogen abstraction reactions. To this end an ab initio database was constructed, which consisted of activation energies calculated with the ab initio CBS-QB3 method for a limited set of well-chosen homologous reactions. CBS-QB3 is shown to predict reaction rate coefficients within a factor of 2,4 and Arrhenius activation energies within 3,5 kJ,mol,1of experimental data. Activation energies in the set of homologous reactions vary over 156 kJ,mol,1with the structure of the abstracting radical and over 94 kJ,mol,1with the structure of the abstracted hydrocarbon. The parameters required for the group contribution method, the so-called standard activation group additivity values, were determined from this database. To test the accuracy of the group contribution method, a large set of 88 additional activation energies were calculated from first principles and compared with the predictions from the group contribution method. It was found that the group contribution method yields accurate activation energies for hydrogen-transfer reactions between hydrogen molecules, alkylic hydrocarbons, and vinylic hydrocarbons, with the largest deviations being less than 6 kJ,mol,1. For reactions between allylic and propargylic hydrocarbons, the transition state is believed to be stabilized by resonance effects, thus requiring the introduction of an appropriate correction term to obtain a reliable prediction of the activation energy for this subclass of hydrogen abstraction reactions. [source]


Universality at the edge of the spectrum for unitary, orthogonal, and symplectic ensembles of random matrices

COMMUNICATIONS ON PURE & APPLIED MATHEMATICS, Issue 6 2007
Percy Deift
We prove universality at the edge of the spectrum for unitary (, = 2), orthogonal (, = 1), and symplectic (, = 4) ensembles of random matrices in the scaling limit for a class of weights w(x) = e,V(x) where V is a polynomial, V(x) = ,2mx2m + · · ·, ,2m > 0. The precise statement of our results is given in Theorem 1.1 and Corollaries 1.2 and 1.4 below. For the same class of weights, a proof of universality in the bulk of the spectrum is given in [12] for the unitary ensembles and in [9] for the orthogonal and symplectic ensembles. Our starting point in the unitary case is [12], and for the orthogonal and symplectic cases we rely on our recent work [9], which in turn depends on the earlier work of Widom [46] and Tracy and Widom [42]. As in [9], the uniform Plancherel-Rotach-type asymptotics for the orthogonal polynomials found in [12] plays a central role. The formulae in [46] express the correlation kernels for , = 1, 4 as a sum of a Christoffel-Darboux (CD) term, as in the case , = 2, together with a correction term. In the bulk scaling limit [9], the correction term is of lower order and does not contribute to the limiting form of the correlation kernel. By contrast, in the edge scaling limit considered here, the CD term and the correction term contribute to the same order: this leads to additional technical difficulties over and above [49]. © 2006 Wiley Periodicals, Inc. [source]


Flux compactification of M-theory on compact manifolds with Spin(7) holonomy,

FORTSCHRITTE DER PHYSIK/PROGRESS OF PHYSICS, Issue 11-12 2005
D. Constantin
At the leading order, M-theory admits minimal supersymmetric compactifications if the internal manifold has exceptional holonomy. The inclusion of non-vanishing fluxes in M-theory and string theory compactifications induce a superpotential in the lower dimensional theory, which depends on the fluxes. In this work, we check the conjectured form of this superpotential in the case of warped M-theory compactifications on Spin(7) holonomy manifolds. We perform a Kaluza-Klein reduction of the eleven-dimensional supersymmetry transformation for the gravitino and we find by direct comparison the superpotential expression. We check the conjecture for the heterotic string compactified on a Calabi-Yau three-fold as well. The conjecture can be checked indirectly by inspecting the scalar potential obtained after the compactification of M-theory on Spin(7) holonomy manifolds with non-vanishing fluxes. The scalar potential can be written in terms of the superpotential and we show that this potential stabilizes all the moduli fields describing deformations of the metric except for the radial modulus. All the above analyses require the knowledge of the minimal supergravity action in three dimensions. Therefore we calculate the most general causal ,, =1 three-dimensional, gauge invariant action coupled to matter in superspace and derive its component form using Ectoplasmic integration theory. We also show that the three-dimensional theory which results from the compactification is in agreement with the more general supergravity construction. The compactification procedure takes into account higher order quantum correction terms in the low energy effective action. We analyze the properties of these terms on a Spin(7) background. We derive a perturbative set of solutions which emerges from a warped compactification on a Spin(7) holonomy manifold with non-vanishing flux for the M-theory field strength and we show that in general the Ricci flatness of the internal manifold is lost, which means that the supergravity vacua are deformed away from the exceptional holonomy. Using the superpotential form we identify the supersymmetric vacua out of this general set of solutions. [source]


Learning to recognize vulnerable patterns due to undesirable Zone-3 relay operations

IEEJ TRANSACTIONS ON ELECTRICAL AND ELECTRONIC ENGINEERING, Issue 3 2009
Koji Yamashita Member
Abstract Undesirable zone 3 relay operations caused by unexpected loading conditions can contribute to the cascaded events, leading to catastrophic outages. Identifying the basic patterns of zone 3 relay operations in advance is an effective way to help prevent cascaded events. The postcontingency impedances seen by zone 3 relays can be calculated on line in a steady state security assessment framework. However, their accuracy is inadequate compared with the postcontingency apparent impedance obtained from off-line time domain dynamic simulations. This paper proposes a fuzzy inference system (FIS) to correct discrepancies between the postcontingency apparent impedances obtained from the results of steady state security assessment and the corresponding values obtained by time-domain simulations. The postcontingency apparent impedances obtained from the results of steady state security assessment can be corrected on line using correction terms provided by the FIS. The dynamic model of a 200-bus system is used to validate the performance of the proposed FIS. Copyright © 2009 Institute of Electrical Engineers of Japan. Published by John Wiley & Sons, Inc. [source]


Coupled-cluster corrected MR-CISD method with noniterative evaluation of connected triples,

INTERNATIONAL JOURNAL OF QUANTUM CHEMISTRY, Issue 12 2008
Leszek Meissner
Abstract The single-reference coupled-cluster (CC) analysis of the multi-reference configuration interaction singles and doubles (MR-CISD) wave function shows how the MR-CISD results can be used to mimic more effective exponential cluster expansion. The CC corrections have been constructed relying on this analysis and used to improve the accuracy of the MR-CISD energy and reduce the inextensivity error of the MR-CISD method. Disconnected contributions that approximate those components which are not present in the MR-CISD wave function generate the correction terms. Although the linear CI expansion is usually improved through the inclusion of disconnected contributions the coupled-cluster methodology introduces a variety of corrections that, contrary to the CI case, account for higher excitation rank connected components. In this article, we present a combination of both approaches and illustrate its performance on some numerical examples. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008 [source]


Logarithmic perturbation theory for a spiked oscillator and sum rules

INTERNATIONAL JOURNAL OF QUANTUM CHEMISTRY, Issue 1 2005
S. K. Bandyopadhyay
Abstract We show that logarithmic perturbation theory nicely yields the wavefunction correction terms in closed forms for the spiked perturbation ,/x2 on the first excited state of the harmonic oscillator, where the conventional Rayleigh-Schrödinger scheme is known to encounter serious problems. The study also provides a direct route to obtain several sum rules, some of which appear to be new. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005 [source]


A simple approach for improving the hybrid MMVB force field: Application to the photoisomerization of s - cis butadiene

JOURNAL OF COMPUTATIONAL CHEMISTRY, Issue 11 2003
Marco Garavelli
Abstract MMVB is a QM/MM hybrid method, consisting of a molecular mechanics force field coupled to a valence bond Heisenberg Hamiltonian parametrized from ab initio CASSCF calculations on several prototype molecules. The Heisenberg Hamiltonian matrix elements Qij and Kij, whose expressions are partitioned here into a primary contribution and second-order correction terms, are calculated analytically in MMVB. When the original MMVB force field fails to produce potential energy surfaces accurate enough for dynamics calculations, we show that significant improvements can be made by refitting the second-order correction terms for the particular molecule(s) being studied. This "local" reparametrization is based on values of Kij extracted (using effective Hamiltonian techniques) from CASSCF calculations on the same molecule(s). The method is demonstrated for the photoisomerization of s - cis butadiene, and we explain how the correction terms that enabled a successful MMVB dynamics study [Garavelli, M.; Bernardi, F.; Olivucci, M.; Bearpark, M. J.; Klein, S.; Robb, M. A. J Phys Chem A 2001, 105, 11496] were refitted. © 2003 Wiley Periodicals, Inc. J Comput Chem 24: 1357,1363, 2003 [source]


Multiple spin-echo spectroscopic imaging for rapid quantitative assessment of N-acetylaspartate and lactate in acute stroke

MAGNETIC RESONANCE IN MEDICINE, Issue 2 2004
Astrid Stengel
Abstract Monitoring the signal levels of lactate (Lac) and N-acetylaspartate (NAA) by chemical shift imaging can provide additional knowledge about tissue damage in acute stroke. Despite the need for this metabolic information, spectroscopic imaging (SI) has not been used routinely for acute stroke patients, mainly due to the long acquisition time required. The presented data demonstrate that the application of a fast multiple spin-echo (MSE) SI sequence can reduce the measurement time to 6 min (four spin echoes per echo train, 32 × 32 matrix). Quantification of Lac and NAA in terms of absolute concentrations (i.e., mmol/l) can be achieved by means of the phantom replacement approach, with correction terms for the longitudinal and transversal relaxation adapted to the multiple spin-echo sequence. In this pilot study of 10 stroke patients (symptom onset < 24 hr), metabolite concentrations obtained from MSE-SI add important information regarding tissue viability that is not provided by other sequences (e.g., diffusion-weighted imaging (DWI) and perfusion-weighted imaging (PWI)). Metabolic changes extended beyond the borders of the apparent diffusion coefficient (ADC) lesion in nine of the 10 patients, showing a rise in Lac concentrations up to 18 mmol/l, while NAA levels sometimes dropped below the detection level. Considerable differences among the patients in terms of the Lac concentrations and the size of the SI-ADC mismatch were observed. Magn Reson Med 52:228,238, 2004. © 2004 Wiley-Liss, Inc. [source]


First-Order Schemes in the Numerical Quantization Method

MATHEMATICAL FINANCE, Issue 1 2003
V. Bally
The numerical quantization method is a grid method that relies on the approximation of the solution to a nonlinear problem by piecewise constant functions. Its purpose is to compute a large number of conditional expectations along the path of the associated diffusion process. We give here an improvement of this method by describing a first-order scheme based on piecewise linear approximations. Main ingredients are correction terms in the transition probability weights. We emphasize the fact that in the case of optimal quantization, many of these correcting terms vanish. We think that this is a strong argument to use it. The problem of pricing and hedging American options is investigated and a priori estimates of the errors are proposed. [source]


Activation energy of thermally grown silicon dioxide layers on silicon substrates

PHYSICA STATUS SOLIDI (B) BASIC SOLID STATE PHYSICS, Issue 10 2009
G. Gerlach
Abstract A detailed numerical consideration is used as basic approach for calculating profiles of activation energy versus oxide thickness for various temperatures between 780 and 930,°C. Results presented here are intentionally not based on models of diffusion and reaction kinetics to avoid introducing correction terms due to the expansion of theory still under discussion. The statistical calculation gives the mean activation energy of 2.01,eV with standard deviation of 0.10,eV, very close to the overall activation energy of 2.05,eV [M. A. Rabie, Y. M. Haddara, and J. Carette, J. Appl. Phys. 98, 074904 (2005)]. More instructive features of the thermal oxidation of silicon have been disclosed directly from measurements of oxide thickness with time [M. A. Hopper, R. A. Clarke, and L. Young, J. Electrochem. Soc. 122, 1216 (1975) and J. Blanc, Philos. Mag. B 55, 685 (1987)]. Graphs of the natural logarithm of the growth rate versus oxide thickness, in the range between 2 and 65,nm, show that the oxide thickness influences the activation energy EA between 1.4 and 2.7,eV. [source]