Eigenvalue bounds for a class of singular potentials in N dimensions

The eigenvalue bounds obtained earlier [J. Phys. A: Math. Gen. 31 (1998) 963] for smooth transformations of the form V(x) = g(x^2) + f(1/x^2) are extended to N-dimensions. In particular a simple formula is derived which bounds the eigenvalues for the spiked harmonic oscillator potential V(x) = x^2 + lambda/x^alpha, alpha > 0, lambda > 0, and is valid for all discrete eigenvalues, arbitrary angular momentum ell, and spatial dimension N.Comment: 10 pages (plain tex with 2 ps figures). J.Phys.A:Math.Gen.(In Press

Effective Elastic Moduli in Solids with High Crack Density

We investigate the weakening of elastic materials through randomly distributed circles and cracks numerically and compare the results to predictions from homogenization theories. We find a good agreement for the case of randomly oriented cracks of equal length in an isotropic plane-strain medium for lower crack densities; for higher densities the material is weaker than predicted due to precursors of percolation. For a parallel alignment of cracks, where percolation does not occur, we analytically predict a power law decay of the effective elastic constants for high crack densities, and confirm this result numerically.Comment: 8 page

Branching Structures in Elastic Shape Optimization

Fine scale elastic structures are widespread in nature, for instances in plants or bones, whenever stiffness and low weight are required. These patterns frequently refine towards a Dirichlet boundary to ensure an effective load transfer. The paper discusses the optimization of such supporting structures in a specific class of domain patterns in 2D, which composes of periodic and branching period transitions on subdomain facets. These investigations can be considered as a case study to display examples of optimal branching domain patterns. In explicit, a rectangular domain is decomposed into rectangular subdomains, which share facets with neighbouring subdomains or with facets which split on one side into equally sized facets of two different subdomains. On each subdomain one considers an elastic material phase with stiff elasticity coefficients and an approximate void phase with orders of magnitude softer material. For given load on the outer domain boundary, which is distributed on a prescribed fine scale pattern representing the contact area of the shape, the interior elastic phase is optimized with respect to the compliance cost. The elastic stress is supposed to be continuous on the domain and a stress based finite volume discretization is used for the optimization. If in one direction equally sized subdomains with equal adjacent subdomain topology line up, these subdomains are consider as equal copies including the enforced boundary conditions for the stress and form a locally periodic substructure. An alternating descent algorithm is employed for a discrete characteristic function describing the stiff elastic subset on the subdomains and the solution of the elastic state equation. Numerical experiments are shown for compression and shear load on the boundary of a quadratic domain.Comment: 13 pages, 6 figure

Spectra generated by a confined softcore Coulomb potential

Analytic and approximate solutions for the energy eigenvalues generated by a confined softcore Coulomb potentials of the form a/(r+\beta) in d>1 dimensions are constructed. The confinement is effected by linear and harmonic-oscillator potential terms, and also through `hard confinement' by means of an impenetrable spherical box. A byproduct of this work is the construction of polynomial solutions for a number of linear differential equations with polynomial coefficients, along with the necessary and sufficient conditions for the existence of such solutions. Very accurate approximate solutions for the general problem with arbitrary potential parameters are found by use of the asymptotic iteration method.Comment: 17 pages, 2 figure

Solutions for certain classes of Riccati differential equation

We derive some analytic closed-form solutions for a class of Riccati equation y'(x)-\lambda_0(x)y(x)\pm y^2(x)=\pm s_0(x), where \lambda_0(x), s_0(x) are C^{\infty}-functions. We show that if \delta_n=\lambda_n s_{n-1}-\lambda_{n-1}s_n=0, where \lambda_{n}= \lambda_{n-1}^\prime+s_{n-1}+\lambda_0\lambda_{n-1} and s_{n}=s_{n-1}^\prime+s_0\lambda_{k-1}, n=1,2,..., then The Riccati equation has a solution given by y(x)=\mp s_{n-1}(x)/\lambda_{n-1}(x). Extension to the generalized Riccati equation y'(x)+P(x)y(x)+Q(x)y^2(x)=R(x) is also investigated.Comment: 10 page

Energies and wave functions for a soft-core Coulomb potential

For the family of model soft Coulomb potentials represented by V(r) = -\frac{Z}{(r^q+\beta^q)^{\frac{1}{q}}}, with the parameters Z>0, \beta>0, q \ge 1, it is shown analytically that the potentials and eigenvalues, E_{\nu\ell}, are monotonic in each parameter. The potential envelope method is applied to obtain approximate analytic estimates in terms of the known exact spectra for pure power potentials. For the case q =1, the Asymptotic Iteration Method is used to find exact analytic results for the eigenvalues E_{\nu\ell} and corresponding wave functions, expressed in terms of Z and \beta. A proof is presented establishing the general concavity of the scaled electron density near the nucleus resulting from the truncated potentials for all q. Based on an analysis of extensive numerical calculations, it is conjectured that the crossing between the pair of states [(\nu,\ell),(\nu',\ell')], is given by the condition \nu'\geq (\nu+1) and \ell' \geq (\ell+3). The significance of these results for the interaction of an intense laser field with an atom is pointed out. Differences in the observed level-crossing effects between the soft potentials and the hydrogen atom confined inside an impenetrable sphere are discussed.Comment: 13 pages, 5 figures, title change, minor revision

Memory of the Unjamming Transition during Cyclic Tiltings of a Granular Pile

Discrete numerical simulations are performed to study the evolution of the micro-structure and the response of a granular packing during successive loading-unloading cycles, consisting of quasi-static rotations in the gravity field between opposite inclination angles. We show that internal variables, e.g., stress and fabric of the pile, exhibit hysteresis during these cycles due to the exploration of different metastable configurations. Interestingly, the hysteretic behaviour of the pile strongly depends on the maximal inclination of the cycles, giving evidence of the irreversible modifications of the pile state occurring close to the unjamming transition. More specifically, we show that for cycles with maximal inclination larger than the repose angle, the weak contact network carries the memory of the unjamming transition. These results demonstrate the relevance of a two-phases description -strong and weak contact networks- for a granular system, as soon as it has approached the unjamming transition.Comment: 13 pages, 15 figures, soumis \`{a} Phys. Rev.

Spectral characteristics for a spherically confined -1/r + br^2 potential

We consider the analytical properties of the eigenspectrum generated by a class of central potentials given by V(r) = -a/r + br^2, b>0. In particular, scaling, monotonicity, and energy bounds are discussed. The potential $V(r)$ is considered both in all space, and under the condition of spherical confinement inside an impenetrable spherical boundary of radius R. With the aid of the asymptotic iteration method, several exact analytic results are obtained which exhibit the parametric dependence of energy on a, b, and R, under certain constraints. More general spectral characteristics are identified by use of a combination of analytical properties and accurate numerical calculations of the energies, obtained by both the generalized pseudo-spectral method, and the asymptotic iteration method. The experimental significance of the results for both the free and confined potential V(r) cases are discussed.Comment: 16 pages, 4 figure