Newton`s iterates can converge to non-stationary points

In this note we discuss the convergence of Newton`s method for minimization. We present examples in which the Newton iterates satisfy the Wolfe conditions and the Hessian is positive definite at each step and yet the iterates converge to a non-stationary point. These examples answer a question posed by Fletcher in his 1987 book Practical methods of optimization.

Noether and topological currents equivalence and soliton/particle correspondence in affine sl(2)((1)) Toda theory coupled to matter

A submodel of the so-called conformal affine Toda model coupled to the matter field (CATM) is defined such that its real Lagrangian has a positive-definite kinetic term for the Toda field and a usual kinetic term for the (Dirac) spinor field. After spontaneously broken the conformal symmetry by means of BRST analysis, we end up with an effective theory, the off-critical affine Toda model coupled to the matter (ATM). It is shown that the ATM model inherits the remarkable properties of the general CATM model such as the soliton solutions, the particle/soliton correspondence and the equivalence between the Noether and topological currents. The classical solitonic spectrum of the ATM model is also discussed. (C) 2001 Elsevier B.V. B.V. All rights reserved.

Viabilidade econômica de quatro intensidades de pastejo do capim Tanzânia na pecuária de corte1

Were evaluated the economic effects of four grazing heights (20, 40, 60 and 80 cm) of Tanzânia grass for beef cattle. The total area was 12 acres divided into paddocks of one hectare, with three replications. We used three animals, males, Nelore per paddock, as more animals need additional adjustment to the height you want. The reduction in sward height allowed higher stocking rate, which, even with a reduction in individual performance, there was more beef production per area. The interest rate on capital given the property represented the largest share in the final cost of production. The participation of fixed variables showed a positive linear function of the increase in height from grazing, showing significant reduction of the production scale. A reduction in cost of production per hectare per year with increasing grazing height. No differences were found in certain economic indicators, and the four systems remunerate the capital invested.

Zero sets of bivariate Hermite polynomials

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Reproducing kernel Hilbert spaces associated with kernels on topological spaces

We analyze reproducing kernel Hilbert spaces of positive definite kernels on a topological space X being either first countable or locally compact. The results include versions of Mercer's theorem and theorems on the embedding of these spaces into spaces of continuous and square integrable functions.

Weighted Fourier–Laplace transforms in reproducing kernel Hilbert spaces on the sphere

We study the action of a weighted Fourier–Laplace transform on the functions in the reproducing kernel Hilbert space (RKHS) associated with a positive definite kernel on the sphere. After defining a notion of smoothness implied by the transform, we show that smoothness of the kernel implies the same smoothness for the generating elements (spherical harmonics) in the Mercer expansion of the kernel. We prove a reproducing property for the weighted Fourier–Laplace transform of the functions in the RKHS and embed the RKHS into spaces of smooth functions. Some relevant properties of the embedding are considered, including compactness and boundedness. The approach taken in the paper includes two important notions of differentiability characterized by weighted Fourier–Laplace transforms: fractional derivatives and Laplace–Beltrami derivatives.

Novel algorithms for electronic structure based molecular dynamics with linear system-size scaling

Die vorliegende Arbeit behandelt die Entwicklung und Verbesserung von linear skalierenden Algorithmen für Elektronenstruktur basierte Molekulardynamik. Molekulardynamik ist eine Methode zur Computersimulation des komplexen Zusammenspiels zwischen Atomen und Molekülen bei endlicher Temperatur. Ein entscheidender Vorteil dieser Methode ist ihre hohe Genauigkeit und Vorhersagekraft. Allerdings verhindert der Rechenaufwand, welcher grundsätzlich kubisch mit der Anzahl der Atome skaliert, die Anwendung auf große Systeme und lange Zeitskalen. Ausgehend von einem neuen Formalismus, basierend auf dem großkanonischen Potential und einer Faktorisierung der Dichtematrix, wird die Diagonalisierung der entsprechenden Hamiltonmatrix vermieden. Dieser nutzt aus, dass die Hamilton- und die Dichtematrix aufgrund von Lokalisierung dünn besetzt sind. Das reduziert den Rechenaufwand so, dass er linear mit der Systemgröße skaliert. Um seine Effizienz zu demonstrieren, wird der daraus entstehende Algorithmus auf ein System mit flüssigem Methan angewandt, das extremem Druck (etwa 100 GPa) und extremer Temperatur (2000 - 8000 K) ausgesetzt ist. In der Simulation dissoziiert Methan bei Temperaturen oberhalb von 4000 K. Die Bildung von sp²-gebundenem polymerischen Kohlenstoff wird beobachtet. Die Simulationen liefern keinen Hinweis auf die Entstehung von Diamant und wirken sich daher auf die bisherigen Planetenmodelle von Neptun und Uranus aus. Da das Umgehen der Diagonalisierung der Hamiltonmatrix die Inversion von Matrizen mit sich bringt, wird zusätzlich das Problem behandelt, eine (inverse) p-te Wurzel einer gegebenen Matrix zu berechnen. Dies resultiert in einer neuen Formel für symmetrisch positiv definite Matrizen. Sie verallgemeinert die Newton-Schulz Iteration, Altmans Formel für beschränkte und nicht singuläre Operatoren und Newtons Methode zur Berechnung von Nullstellen von Funktionen. Der Nachweis wird erbracht, dass die Konvergenzordnung immer mindestens quadratisch ist und adaptives Anpassen eines Parameters q in allen Fällen zu besseren Ergebnissen führt.

Analysis and numerical solution of the Peterlin viscoelastic model

Liquids and gasses form a vital part of nature. Many of these are complex fluids with non-Newtonian behaviour. We introduce a mathematical model describing the unsteady motion of an incompressible polymeric fluid. Each polymer molecule is treated as two beads connected by a spring. For the nonlinear spring force it is not possible to obtain a closed system of equations, unless we approximate the force law. The Peterlin approximation replaces the length of the spring by the length of the average spring. Consequently, the macroscopic dumbbell-based model for dilute polymer solutions is obtained. The model consists of the conservation of mass and momentum and time evolution of the symmetric positive definite conformation tensor, where the diffusive effects are taken into account. In two space dimensions we prove global in time existence of weak solutions. Assuming more regular data we show higher regularity and consequently uniqueness of the weak solution. For the Oseen-type Peterlin model we propose a linear pressure-stabilized characteristics finite element scheme. We derive the corresponding error estimates and we prove, for linear finite elements, the optimal first order accuracy. Theoretical error of the pressure-stabilized characteristic finite element scheme is confirmed by a series of numerical experiments.

Hyalotekite, (Ba,Pb,K)(4)(Ca,Y)(2)Si-8(B,Be)(2)(Si,B)(2)O28F, a Tectosilicate Related to Scapolite: New Structure Refinement, Phase Transitions and a Short-Range Ordered 3B Superstructure

Hyalotekite, a framework silicate of composition (Ba,Pb,K)(4)(Ca,Y)(2)Si-8(B,Be)(2) (Si,B)(2)O28F, is found in relatively high-temperature(greater than or equal to 500 degrees C) Mn skarns at Langban, Sweden, and peralkaline pegmatites at Dara-i-Pioz, Tajikistan. A new paragenesis at Dara-i-Pioz is pegmatite consisting of the Ba borosilicates leucosphenite and tienshanite, as well as caesium kupletskite, aegirine, pyrochlore, microcline and quartz. Hyalotekite has been partially replaced by barylite and danburite. This hyalotekite contains 1.29-1.78 wt.% Y2O3, equivalent to 0.172-0.238 Y pfu or 8-11% Y on the Ca site; its Pb/(Pb+Ba) ratio ranges 0.36-0.44. Electron microprobe F contents of Langban and Dara-i-Pioz hyalotekite range 1.04-1.45 wt.%, consistent with full occupancy of the F site. A new refinement of the structure factor data used in the original structural determination of a Langban hyalotekite resulted in a structural formula, (Pb1.96Ba1.86K0.18)Ca-2(B1.76Be0.24)(Si1.56B0.44)Si8O28F, consistent with chemical data and all cations with positive-definite thermal parameters, although with a slight excess of positive charge (+57.14 as opposed to the ideal +57.00). An unusual feature of the hyalotekite framework is that 4 of 28 oxygens are non-bridging; by merging these 4 oxygens into two, the framework topology of scapolite is obtained. The triclinic symmetry of hyalotekite observed at room temperature is obtained from a hypothetical tetragonal parent structure via a sequence of displacive phase transitions. Some of these transitions are associated with cation ordering, either Pb-Ba ordering in the large cation sites, or B-Be and Si-B ordering on tetrahedral sites. Others are largely displacive but affect the coordination of the large cations (Pb, Ba, K, Ca). High-resolution electron microscopy suggests that the undulatory extinction characteristic of hyalotekite is due to a fine mosaic microstructure. This suggests that at least one of these transitions occurs in nature during cooling, and that it is first order with a large volume change. A diffuse superstructure observed by electron diffraction implies the existence of a further stage of short-range cation ordering which probably involves both (Pb,K)-Ba and (BeSi,BB)-BSi.

Regulation of Ferulate-5-Hydroxylase Expression in Arabidopsis in the Context of Sinapate Ester Biosynthesis1

Sinapic acid is an intermediate in syringyl lignin biosynthesis in angiosperms, and in some taxa serves as a precursor for soluble secondary metabolites. The biosynthesis and accumulation of the sinapate esters sinapoylglucose, sinapoylmalate, and sinapoylcholine are developmentally regulated in Arabidopsis and other members of the Brassicaceae. The FAH1 locus of Arabidopsis encodes the enzyme ferulate-5-hydroxylase (F5H), which catalyzes the rate-limiting step in syringyl lignin biosynthesis and is required for the production of sinapate esters. Here we show that F5H expression parallels sinapate ester accumulation in developing siliques and seedlings, but is not rate limiting for their biosynthesis. RNA gel-blot analysis indicated that the tissue-specific and developmentally regulated expression of F5H mRNA is distinct from that of other phenylpropanoid genes. Efforts to identify constructs capable of complementing the sinapate ester-deficient phenotype of fah1 mutants demonstrated that F5H expression in leaves is dependent on sequences 3′ of the F5H coding region. In contrast, the positive regulatory function of the downstream region is not required for F5H transcript or sinapoylcholine accumulation in embryos.

Sequential and parallel synchronous alternating iterative methods

The so-called parallel multisplitting nonstationary iterative Model A was introduced by Bru, Elsner, and Neumann [Linear Algebra and its Applications 103:175-192 (1988)] for solving a nonsingular linear system Ax = b using a weak nonnegative multisplitting of the first type. In this paper new results are introduced when A is a monotone matrix using a weak nonnegative multisplitting of the second type and when A is a symmetric positive definite matrix using a P -regular multisplitting. Also, nonstationary alternating iterative methods are studied. Finally, combining Model A and alternating iterative methods, two new models of parallel multisplitting nonstationary iterations are introduced. When matrix A is monotone and the multisplittings are weak nonnegative of the first or of the second type, both models lead to convergent schemes. Also, when matrix A is symmetric positive definite and the multisplittings are P -regular, the schemes are also convergent.

Quantum dynamics with stochastic gauge simulations

The general idea of a stochastic gauge representation is introduced and compared with more traditional phase-space expansions, like the Wigner expansion. Stochastic gauges can be used to obtain an infinite class of positive-definite stochastic time-evolution equations, equivalent to master equations, for many systems including quantum time evolution. The method is illustrated with a variety of simple examples ranging from astrophysical molecular hydrogen production, through to the topical problem of Bose-Einstein condensation in an optical trap and the resulting quantum dynamics.

On the possibility of elastic strain localisation in a fault

The phenomenon of strain localisation is often observed in shear deformation of particulate materials, e.g., fault gouge. This phenomenon is usually attributed to special types of plastic behaviour of the material (e.g., strain softening or mismatch between dilatancy and pressure sensitivity or both). Observations of strain localisation in situ or in experiments are usually based on displacement measurements and subsequent computation of the displacement gradient. While in conventional continua the symmetric part of the displacement gradient is equal to the strain, it is no longer the case in the more realistic descriptions within the framework of generalised continua. In such models the rotations of the gouge particles are considered as independent degrees of freedom the values of which usually differ from the rotation of an infinitesimal volume element of the continuum, the latter being described for infinitesimal deformations by the non-symmetric part of the displacement gradient. As a model for gouge material we propose a continuum description for an assembly of spherical particles of equal radius in which the particle rotation is treated as an independent degree of freedom. Based on this model we consider simple shear deformations of the fault gouge. We show that there exist values of the model parameters for which the displacement gradient exhibits a pronounced localisation at the mid-layers of the fault, even in the absence of inelasticity. Inelastic effects are neglected in order to highlight the role of the independent rotations and the associated additional parameters. The localisation-like behaviour occurs if (a) the particle rotations on the boundary of the shear layer are constrained (this type of boundary condition does not exist in a standard continuum) and (b) the contact moment-or bending stiffness is much smaller than the product of the effective shear modulus of the granulate and the square of the width of the gouge layer. It should be noted however that the virtual work functional is positive definite over the range of physically meaningful parameters (here: contact stiffnesses, solid volume fraction and coordination number) so that strictly speaking we are not dealing with a material instability.

Phase-space analysis of bosonic spontaneous emission

We present phase-space techniques for the modelling of spontaneous emission in two-level bosonic atoms. The positive-P representation is shown to give a full and complete description within the limits of our model. The Wigner representation, even when truncated at second order, is shown to need a doubling of the phase-space to allow for a positive-definite diffusion matrix in the appropriate Fokker-Planck equation and still fails to agree with the full quantum results of the positive-P representation. We show that quantum statistics and correlations between the ground and excited states affect the dynamics of the emission process, so that it is in general non-exponential. (c) 2005 Elsevier B.V. All rights reserved.

Nonclassicality of a photon-subtracted Gaussian field

We investigate the nonclassicality of a photon-subtracted Gaussian field, which was produced in a recent experiment, using negativity of the Wigner function and the nonexistence of well-behaved positive P function. We obtain the condition to see negativity of the Wigner function for the case including the mixed Gaussian incoming field, the threshold photodetection and the inefficient homodyne measurement. We show how similar the photon-subtracted state is to a superposition of coherent states.