Efficient adaptive integration of functions with sharp gradients and cusps in n-dimensional parallelepipeds

In this paper, we study the efficient numerical integration of functions with sharp gradients and cusps. An adaptive integration algorithm is presented that systematically improves the accuracy of the integration of a set of functions. The algorithm is based on a divide and conquer strategy and is independent of the location of the sharp gradient or cusp. The error analysis reveals that for a $C^0$ function (derivative-discontinuity at a point), a rate of convergence of $n+1$ is obtained in $R^n$. Two applications of the adaptive integration scheme are studied. First, we use the adaptive quadratures for the integration of the regularized Heaviside function---a strongly localized function that is used for modeling sharp gradients. Then, the adaptive quadratures are employed in the enriched finite element solution of the all-electron Coulomb problem in crystalline diamond. The source term and enrichment functions of this problem have sharp gradients and cusps at the nuclei. We show that the optimal rate of convergence is obtained with only a marginal increase in the number of integration points with respect to the pure finite element solution with the same number of elements. The adaptive integration scheme is simple, robust, and directly applicable to any generalized finite element method employing enrichments with sharp local variations or cusps in $n$-dimensional parallelepiped elements.Comment: 22 page

CONTINUITY OF CONDITION SPECTRUM AND ITS LEVEL SET IN BANACH ALGEBRA

For 0 < � < 1 and a Banach algebra element a, this thesis aims to establish the results related to continuity of condition spectrum and its level set correspondence at (�; a). Here we propose a method of study to achieve the continuity. We first identify the Banach algebras at which the interior of the level set of condition spectrum is empty and then we obtain the continuity results. This thesis consists of four chapters. Chapter 1 contains all the prerequisites which are crucial for the development of the thesis. In particular, this chapter has a quick review of the basic properties of spectrum, condition spectrum, upper and lower hemicontiuous correspondences. We also concentrate on analytic vector valued maps and generalized maximum modulus theorem for them. For an element a in A, Chapter 2 has the results related to interior of the level of set of the condition spectrum of a. At first, we focus on

Multi-Channel Inverse Scattering Problem on the Line: Thresholds and Bound States

We consider the multi-channel inverse scattering problem in one-dimension in the presence of thresholds and bound states for a potential of finite support. Utilizing the Levin representation, we derive the general Marchenko integral equation for N-coupled channels and show that, unlike to the case of the radial inverse scattering problem, the information on the bound state energies and asymptotic normalization constants can be inferred from the reflection coefficient matrix alone. Thus, given this matrix, the Marchenko inverse scattering procedure can provide us with a unique multi-channel potential. The relationship to supersymmetric partner potentials as well as possible applications are discussed. The integral equation has been implemented numerically and applied to several schematic examples showing the characteristic features of multi-channel systems. A possible application of the formalism to technological problems is briefly discussed.Comment: 19 pages, 5 figure

PT-symmetric square well and the associated SUSY hierarchies

The PT-symmetric square well problem is considered in a SUSY framework. When the coupling strength $Z$ lies below the critical value $Z_0^{\rm (crit)}$ where PT symmetry becomes spontaneously broken, we find a hierarchy of SUSY partner potentials, depicting an unbroken SUSY situation and reducing to the family of $\sec^2$-like potentials in the $Z \to 0$ limit. For $Z$ above $Z_0^{\rm (crit)}$, there is a rich diversity of SUSY hierarchies, including some with PT-symmetry breaking and some with partial PT-symmetry restoration.Comment: LaTeX, 18 pages, no figure; broken PT-symmetry case added (Sec. 6

Generalized Duffy transformation for integrating vertex singularities

For an integrand with a 1/r vertex singularity, the Duffy transformation from a triangle (pyramid) to a square (cube) provides an accurate and efficient technique to evaluate the integral. In this paper, we generalize the Duffy transformation to power singularities of the form p(x)/r α , where p is a trivariate polynomial and α &gt; 0 is the strength of the singularity. We use the map (u, v, w) → (x, y, z) : x = u β , y = x v, z = x w, and judiciously choose β to accurately estimate the integral. For α = 1, the Duffy transformation (β = 1) is optimal, whereas if α ≠ 1, we show that there are other values of β that prove to be substantially better. Numerical tests in two and three dimensions are presented that reveal the improved accuracy of the new transformation. Higher-order partition of unity finite element solutions for the Laplace equation with a derivative singularity at a re-entrant corner are presented to demonstrate the benefits of using the generalized Duffy transformation

Controls of soil spatial variability in a dry tropical forest

We examined the roles of lithology, topography, vegetation and fire in generating local-scale (= 1 cm diameter at breast height (DBH), and spatial variation in fire frequency (times burnt during the 17 years preceding soil sampling) in a permanent 50-ha plot. Unlike classic catenas, lower elevation soils had lesser moisture, plant-available Ca, Cu, Mn, Mg, Zn, B, clay and total C. The distribution of plant-available Ca, Cu, Mn and Mg appeared to largely be determined by the whole-rock chemical composition differences between amphibolites and hornblende-biotite gneisses. Amphibolites were associated with summit positions, while gneisses dominated lower elevations, an observation that concurs with other studies in the region which suggest that hillslope-scale topography has been shaped by differential weathering of lithologies. Neither NO3--N nor NH4+-N was explained by the basal area of trees belonging to Fabaceae, a family associated with N-fixing species, and no long-term effects of fire on soil parameters were detected. Local-scale lithological variation is an important first-order control over soil variability at the hillslope scale in this SDTF, by both direct influence on nutrient stocks and indirect influence via control of local relief

Pseudospin, Spin, and Coulomb Dirac-Symmetries: Doublet Structure and Supersymmetric Patterns

Relativistic symmetries of the Dirac Hamiltonian with a mixture of spherically symmetric Lorentz scalar and vector potentials, are examined from the point of view of supersymmetric quantum mechanics. The cases considered include the Coulomb, pseudospin and spin limits relevant, respectively, to atoms, nuclei and hadrons.Comment: 8 pages, 1 figure, Proc. Int. Workshop on "Blueprints for the Nucleus: From First Principles to Collective Motion", May 17-23, 2004, Feza Gursey Institute, Istanbul, Turke

Phase shift effective range expansion from supersymmetric quantum mechanics

Supersymmetric or Darboux transformations are used to construct local phase equivalent deep and shallow potentials for $\ell \neq 0$ partial waves. We associate the value of the orbital angular momentum with the asymptotic form of the potential at infinity which allows us to introduce adequate long-distance transformations. The approach is shown to be effective in getting the correct phase shift effective range expansion. Applications are considered for the $^1P_1$ and $^1D_2$ partial waves of the neutron-proton scattering.Comment: 6 pages, 3 figures, Revtex4, version to be publised in Physical Review