Source author record

Paolo Amore

Paolo Amore appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

45works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

45 published item(s)

preprint2022arXiv

Efficient algorithms for the dense packing of congruent circles inside a square

We study dense packings of a large number of congruent non-overlapping circles inside a square by looking for configurations which maximize the packing density, defined as the ratio between the area occupied by the disks and the area of the square container. The search for these configurations is carried out with the help of two algorithms that we have devised: a first algorithm is in charge of obtaining sufficiently dense configurations starting from a random guess, while a second algorithm improves the configurations obtained in the first stage. The algorithms can be used sequentially or independently. The performance of these algorithms is assessed by carrying out numerical tests for configurations with a large number of circles.

preprint2022arXiv

Quantum particles in a suddenly accelerating potential

We study the behavior of a quantum particle trapped in a confining potential in one dimension under multiple sudden changes of velocity and/or acceleration. We develop the appropriate formalism to deal with such situation and we use it to calculate the probability of transition for simple problems such as the particle in an infinite box and the simple harmonic oscillator. For the infinite box of length $L$ under two and three sudden changes of velocity, where the initial and final velocity vanish, we find that the system undergoes quantum revivals for $Δt = τ_0 \equiv \frac{4mL^2}{π\hbar}$, regardless of other parameters ($Δt$ is the time elapsed between the first and last change of velocity). For the simple harmonic oscillator we find that the states obtained by suddenly changing (one change) the velocity and/or the acceleration of the potential, for a particle initially in an eigenstate of the static potential, are {\sl coherent} states. For multiple changes of acceleration or velocity we find that the quantum expectation value of the Hamiltonian is remarkably close (possibly identical) to the corresponding classical expectation values. Finally, the probability of transition for a particle in an accelerating harmonic oscillator (no sudden changes) calculated with our formalism agrees with the formula derived long time ago by Ludwig and recently modified by Dodonov~\cite{Dodonov21}, but with a different expression for the dimensionless parameter $γ$. Our probability agrees with the one of ref.~\cite{Dodonov21} for $γ\ll 1$ but is not periodic in time (it decays monotonously), contrary to the result derived in ref.~\cite{Dodonov21}.

preprint2021arXiv

Comment on: "Quantum aspects of the Lorentz symmetry violation on an electron in a nonuniform electric field'' Eur. Phys. J. Plus (2020) 135:623

We analyze recent results concerning the hypothesis of a privileged direction in the space-time that is made by considering a background of the Lorentz symmetry violation determined by a fixed spacelike vector field and the analysis of quantum effects of this background on the interaction of a nonrelativistic electron with a nonuniform electric field produced by a uniform electric charge distribution. We show that the conclusions derived by the authors are an artifact of the truncation of the Frobenius series by means of the tree-term recurrence relation for the expansion coefficients. Thus, the existence of allowed angular frequencies stemming from this procedure is meaningless and unphysical.

preprint2021arXiv

Computing the solutions of the van der Pol equation to arbitrary precision

We describe an extension of the Taylor method for the numerical solution of ODEs that uses Padé approximants to obtain extremely precise numerical results. The accuracy of the results is essentially limited only by the computer time and memory, provided that one works in arbitrary precision. In this method the stepsize is adjusted to achieve the desired accuracy (variable stepsize), while the order of the Taylor expansion can be either fixed or changed at each iteration (variable order). As an application, we have calculated the periodic solutions (limit cycle) of the van der Pol equation with an unprecedented accuracy for a large set of couplings (well beyond the values currently found in the literature) and we have used these numerical results to validate the asymptotic behavior of the period, of the amplitude and of the Lyapunov exponent reported in the literature. We have also used the numerical results to infer the formulas for the asymptotic behavior of the fast component of the period and of the maximum velocity, which have never been calculated before.

preprint2020arXiv

Comment on "Bidimensional bound states for charged polar nanoparticles''

In a recent paper Castellanos-Jaramillo and Castellanos-Moreno proposed a simple quantum-mechanical model for an electron in the vicinity of an ionized nanostructure with a permanent electric dipole. They chose the interaction of the electron with the charge and the dipole in such a way that the resulting Schrödinger equation is separable into radial and angular parts. In this comment we show that those authors did not solve the angular eigenvalue equation with proper periodic boundary conditions and that they also made a mistake in the elimination of the first derivative in the radial equation. Such errors invalidate their results of the Einstein coefficients for the $\left( GaAs\right)_{3}$ system considered.

preprint2020arXiv

On some conditionally solvable quantum-mechanical problems

We analyze two conditionally solvable quantum-mechanical models: a one-dimensional sextic oscillator and a perturbed Coulomb problem. Both lead to a three-term recurrence relation for the expansion coefficients. We show diagrams of the distribution of their exact eigenvalues with the addition of accurate ones from variational calculations. We discuss the symmetry of such distributions. We also comment on the wrong interpretation of the exact eigenvalues and eigenfunctions by some researchers that has led to the prediction of allowed cyclotron frequencies and field intensities.

preprint2020arXiv

Spectral sum rules for the Schrödinger equation

We study the sum rules of the form $Z(s) = \sum_n E_n^{-s}$, where $E_n$ are the eigenvalues of the time--independent Schrödinger equation (in one or more dimensions) and $s$ is a rational number for which the series converges. We have used perturbation theory to obtain an explicit formula for the sum rules up to second order in the perturbation and we have extended it non--perturbatively by means of a Padé--approximant. For the special case of a box decorated with one impurity in one dimension we have calculated the first few sum rules of integer order exactly; the sum rule of order one has also been calculated exactly for the problem of a box with two impurities. In two dimensions we have considered the case of an impurity distributed on a circle of arbitrary radius and we have calculated the exact sum rules of order two. Finally we show that exact sum rules can be obtained, in one dimension, by transforming the Schrödinger equation into the Helmholtz equation with a suitable density.

preprint2019arXiv

Exact sum rules for heterogeneous spherical drums

We have obtained explicit integral expressions for the sums of inverse powers of the eigenvalues of the Laplacian on a unit sphere, in presence of an arbitrary variable density. The exact expressions for the sum rules are obtained by properly "renormalizing" the series, excluding the divergent contribution of the vanishing lowest eigenvalue. For a non--trivial example of a variable density we have applied our formulas to calculate the exact sum rules of order two and three, and we have verified these results calculating the sum rules numerically using the eigenvalues obtained with the Rayleigh-Ritz method.

preprint2016arXiv

Bound states in weakly deformed waveguides: numerical vs analytical results

We have studied the emergence of bound states in weakly deformed and/or heterogeneous waveguides, comparing the analytical predictions obtained using a recently developed perturbative method, with precise numerical results, for different configurations (a homogeneous asymmetric waveguide, a heterogenous asymmetric waveguide and a homogeneous broken-strip). In all the examples considered in this paper we have found excellent agreement between analytical and numerical results, thus providing a numerical verification of the analytical approach.

preprint2016arXiv

On two perturbation series for short-range attractive potentials

We compare two alternative expansions for finite attractive wells. One of them is known from long ago and is given in terms of powers of the strength parameter. The other one is based on the solution of the equations of the Rayleigh-Schrödinger perturbation theory in a basis set of functions of period $L$. The analysis of exactly solvable models shows that although the exact solution of the problem with periodic boundary conditions yields the correct result when $L\rightarrow \infty$ the coefficients of the series for this same problem blow up and fail to produce the correct asymptotic expansion.

preprint2015arXiv

High order eigenvalues for the Helmholtz equation in complicated non-tensor domains through Richardson Extrapolation of second order finite differences

We apply second order finite difference to calculate the lowest eigenvalues of the Helmholtz equation, for complicated non-tensor domains in the plane, using different grids which sample exactly the border of the domain. We show that the results obtained applying Richardson and Padé-Richardson extrapolation to a set of finite difference eigenvalues corresponding to different grids allows to obtain extremely precise values. When possible we have assessed the precision of our extrapolations comparing them with the highly precise results obtained using the method of particular solutions. Our empirical findings suggest an asymptotic nature of the FD series. In all the cases studied, we are able to report numerical results which are more precise than those available in the literature.

preprint2015arXiv

On the symmetry of three identical interacting particles in a one-dimensional box

We study a quantum-mechanical system of three particles in a one-dimensional box with two-particle harmonic interactions. The symmetry of the system is described by the point group $D_{3d}$. Group theory greatly facilitates the application of perturbation theory and the Rayleigh-Ritz variational method. A great advantage is that every irreducible representation can be treated separately. Group theory enables us to predict the connection between the states for the small box length and large box length regimes of the system. We discuss the crossings and avoided crossings of the energy levels as well as other interesting features of the spectrum of the system.

preprint2015arXiv

Weakly bound states in heterogeneous waveguides

We study the spectrum of the Helmholtz equation in a two-dimensional infinite waveguide, containing a weak heterogeneity localized at an internal point, and obeying Dirichlet boundary conditions at its border. We prove that, when the heterogeneity corresponds to a locally denser material, the lowest eigenvalue of the spectrum falls below the continuum threshold and a bound state appears, localized at the heterogeneity. We devise a rigorous perturbation scheme and derive the exact expression for the energy to third order in the heterogeneity.

preprint2014arXiv

Accurate calculation of the solutions to the Thomas-Fermi equations

We obtain highly accurate solutions to the Thomas-Fermi equations for atoms and atoms in very strong magnetic fields. We apply the Padé-Hankel method, numerical integration, power series with Padé and Hermite-Padé approximants and Chebyshev polynomials. Both the slope at origin and the location of the right boundary in the magnetic-field case are given with unprecedented accuracy.

preprint2014arXiv

Heterogeneous systems in $d$ dimensions: lower spectrum

We show that the properties of the lower part of the spectrum of the Helmholtz equation for an heterogeneous system in a finite region in $d$ dimensions, where the solutions to the homogeneous problems are known, can be systematically approximated by means of iterative methods. These methods only require the specification of an arbitrary ansatz and necessarily converge to the desired solution, regardless of the strength of the inhomogeneities, provided that the ansatz has a finite overlap with it. Different boundary conditions at the borders of the domain can be assumed. Applications in one and two dimensions are used to illustrate the methods.

preprint2014arXiv

Is space-time symmetry a suitable generalization of parity-time symmetry?

We discuss space-time symmetric Hamiltonian operators of the form $% H=H_{0}+igH^{\prime}$, where $H_{0}$ is Hermitian and $g$ real. $H_{0}$ is invariant under the unitary operations of a point group $G$ while $H^{\prime}$ is invariant under transformation by elements of a subgroup $G^{\prime}$ of $G$. If $G$ exhibits irreducible representations of dimension greater than unity, then it is possible that $H$ has complex eigenvalues for sufficiently small nonzero values of $g$. In the particular case that $H$ is parity-time symmetric then it appears to exhibit real eigenvalues for all $% 0<g<g_{c}$, where $g_{c}$ is the exceptional point closest to the origin. Point-group symmetry and perturbation theory enable one to predict whether $% H $ may exhibit real or complex eigenvalues for $g>0$. We illustrate the main theoretical results and conclusions of this paper by means of two- and three-dimensional Hamiltonians exhibiting a variety of different point-group symmetries.

preprint2014arXiv

Non-Hermitian oscillators with $T_{d}$ symmetry

We analyse some PT-symmetric oscillators with $T_{d}$ symmetry that depend on a potential parameter $g$. We calculate the eigenvalues and eigenfunctions for each irreducible representation and for a range of values of $g$. Pairs of eigenvalues coalesce at exceptional points $g_c$; their magnitude roughly decreasing with the magnitude of the eigenvalues. It is difficult to estimate whether there is a phase transition at a nonzero value of $g$ as conjectured in earlier papers. Group theory and perturbation theory enable one to predict whether a given space-time symmetry leads to real eigenvalues for sufficiently small nonzero values of $g$.

preprint2014arXiv

On an approximation to the Schrödinger equation with the Hellmann potential

We tried to determine the range of validity of a recently proposed modification of the Hellmann potential that leads to analytical eigenvalues and eigenfunctions. We discuss the difficulties that we found in the analysis of the main equations and results. We conclude that the eigenvalues reported by the authors do not exhibit the same order as those of the Hellmann potential thus leading to a different underlying physics. What is more: the spectrum of the modified model is qualitatively different from the one supported by the Hellmann potential.

preprint2014arXiv

On the eigenvalues of some non-Hermitian Hamiltonians with space-time symmetry

We calculate the eigenvalues of some two-dimensional non-Hermitian Hamiltonians by means of a pseudospectral method and straightforward diagonalization of the Hamiltonian matrix in a suitable basis set. Both sets of results agree remarkably well but differ considerably from the eigenvalues obtained some time ago by other authors. In particular, we do not observe the multiple phase transitions claimed to occur in one of the anharmonic oscillators.

preprint2014arXiv

Small-energy series for one-dimensional quantum-mechanical models with non-symmetric potentials

We generalize a recently proposed small-energy expansion for one-dimensional quantum-mechanical models. The original approach was devised to treat symmetric potentials and here we show how to extend it to non-symmetric ones. Present approach is based on matching the logarithmic derivatives for the left and right solutions to the Schrödinger equation at the origin (or any other point chosen conveniently) . As in the original method, each logarithmic derivative can be expanded in a small-energy series by straightforward perturbation theory. We test the new approach on four simple models, one of which is not exactly solvable. The perturbation expansion converges in all the illustrative examples so that one obtains the ground-state energy with an accuracy determined by the number of available perturbation corrections.

preprint2013arXiv

$\mathcal{PT}$-symmetric strings

We study both analytically and numerically the spectrum of inhomogeneous strings with $\mathcal{PT}$-symmetric density. We discuss an exactly solvable model of $\mathcal{PT}$-symmetric string which is isospectral to the uniform string; for more general strings, we calculate exactly the sum rules $Z(p) \equiv \sum_{n=1}^\infty 1/E_n^p$, with $p=1,2,\dots$ and find explicit expressions which can be used to obtain bounds on the lowest eigenvalue. A detailed numerical calculation is carried out for two non-solvable models depending on a parameter, obtaining precise estimates of the critical values where pair of real eigenvalues become complex.

preprint2013arXiv

Exact sum rules for inhomogeneous drums

We derive general expressions for the sum rules of the eigenvalues of drums of arbitrary shape and arbitrary density, obeying different boundary conditions. The formulas that we present are a generalization of the analogous formulas for one dimensional inhomogeneous systems that we have obtained in a previous paper. We also discuss the extension of these formulas to higher dimensions. We show that in the special case of a density depending only on one variable the sum rules of any integer order can be expressed in terms of a single series. As an application of our result we derive exact sum rules for the a homogeneous circular annulus with different boundary conditions, for a homogeneous circular sector and for a radially inhomogeneous circular annulus with Dirichlet boundary conditions.

preprint2013arXiv

Exact sum rules for inhomogeneous strings

We derive explicit expressions for the sum rules of the eigenvalues of inhomogeneous strings with arbitrary density and with different boundary conditions. We show that the sum rule of order $N$ may be obtained in terms of a diagrammatic expansion, with $(N-1)!/2$ independent diagrams. These sum rules are used to derive upper and lower bounds to the energy of the fundamental mode of an inhomogeneous string; we also show that it is possible to improve these approximations taking into account the asymptotic behaviour of the spectrum and applying the Shanks transformation to the sequence of approximations obtained to the different orders. We discuss three applications of these results.

preprint2013arXiv

Exact sum rules for inhomogeneous systems containing a zero mode

We show that the formulas for the sum rules for the eigenvalues of inhomogeneous systems that we have obtained in two recent papers are incomplete when the system contains a zero mode. We prove that there are finite contributions of the zero mode to the sum rules and we explicitly calculate the expressions for the sum rules of order one and two. The previous results for systems that do not contain a zero mode are unaffected.

preprint2013arXiv

One cannot hear the density of a drum (and further aspects of isospectrality)

It is well known that certain pairs of planar domains have the same spectra of the Laplacian operator. We prove that these domains are still isospectral for a wider class of physical problems, including the cases of heterogeneous drums and of quantum billiards in an external field. In particular we show that the isospectrality is preserved when the density or the potential are symmetric under reflections along the folding lines of the domain. These results are also confirmed numerically using the finite difference method: we find that the pairs of numerical matrices obtained in the discretization are exactly isospectral up to machine precision.

preprint2012arXiv

A perturbative approach to the spectral zeta functions of strings, drums and quantum billiards

We have obtained an explicit expression for the spectral zeta functions and for the heat kernel of strings, drums and quantum billiards working to third order in perturbation theory, using a generalization of the binomial theorem to operators. The perturbative parameter used in the expansion is either the small deformation of a reference domain (for instance a square), or a small variation of the density around a constant value (in two dimensions both cases can apply). This expansion is well defined even in presence of degenerations of the unperturbed spectrum. We have discussed several examples in one, two and three dimensions, obtaining in some cases the analytic continuation of the series, which we have then used to evaluate the corresponding Casimir energy. For the case of a string with piecewise constant density, subject to different boundary conditions, and of two concentric cylinders of very close radii, we have reproduced results previously published, thus obtaining a useful check of our method.

preprint2012arXiv

Bound states for the quantum dipole moment in two dimensions

We calculate accurate eigenvalues and eigenfunctions of the Schrödinger equation for a two-dimensional quantum dipole. This model proved useful for the study of elastic effects of a single edge dislocation. We show that the Rayleigh-Ritz variational method with a basis set of Slater-type functions is considerably more efficient than the same approach with the basis set of point-spectrum eigenfunctions of the two-dimensional hydrogen atom used in earlier calculations.

preprint2012arXiv

Comment on: "Analytical approximations for the collapse of an empty spherical bubble"

We analyze the Rayleigh equation for the collapse of an empty bubble and provide an explanation for some recent analytical approximations to the model. We derive the form of the singularity at the second boundary point and discuss the convergence of the approximants. We also give a rigorous proof of the asymptotic behavior of the coefficients of the power series that are the basis for the approximate expressions.

preprint2012arXiv

Comment on: `Numerical estimates of the spectrum for anharmonic PT symmetric potentials' [Phys. Scr. \textbf{85} (2012) 065005]

We show that the authors of the commented paper draw their conclusions from the eigenvalues of truncated Hamiltonian matrices that do not converge as the matrix dimension increases. In one of the studied examples the authors missed the real positive eigenvalues that already converge towards the exact eigenvalues of the non-Hermitian operator and focused their attention on the complex ones that do not. We also show that the authors misread Bender's argument about the eigenvalues of the harmonic oscillator with boundary conditions in the complex-$x$ plane (Rep. Prog. Phys. {\bf 70} (2007) 947).

preprint2012arXiv

Particle correlation from uncorrelated non Born-Oppenheimer SCF wavefunctions

We analyse a nonadiabatic self-consistent field method by means of an exactly-solvable model. The method is based on nuclear and electronic orbitals that are functions of the cartesian coordinates in the laboratory-fixed frame. The kinetic energy of the center of mass is subtracted from the molecular Hamiltonian operator in the variational process. The results for the simple model are remarkably accurate and show that the integration over the redundant cartesian coordinates leads to couplings among the internal ones.

preprint2011arXiv

Bound states in open coupled asymmetrical waveguides and quantum wires

The behavior of bound states in asymmetric cross, T and L shaped configurations is considered. Because of the symmetries of the wavefunctions, the analysis can be reduced to the case of an electron localized at the intersection of two orthogonal crossed wires of different width. Numerical calculations show that the fundamental mode of this system remains bound for the widths that we have been able to study directly; moreover, the extrapolation of the results obtained for finite widths suggests that this state remains bound even when the width of one arm becomes infinitesimal. We provide a qualitative argument which explains this behavior and that can be generalized to the lowest energy states in each symmetry class. In the case of odd-odd states of the cross we find that the lowest mode is bounded when the width of the two arms is the same and stays bound up to a critical value of the ratio between the widths; in the case of the even-odd states we find that the lowest mode is unbound up to a critical value of the ratio between the widths. Our qualitative arguments suggest that the bound state survives as the width of the vertical arm becomes infinitesimal.

preprint2011arXiv

Comment on "Bound states of edge dislocations: The quantum dipole problem in two dimensions"

We show that the numerical results contained in a recent paper are affected by a non optimal implementation of the methods which have been used to obtain these results. A careful analysis done using the Rayleigh-Ritz method provides a rigorous upper bound for the energy of the ground state of an electron in a two dimensional potential generated by the edge dislocation, as well as precise values for the excited states. The extrapolation of the results corresponding to different subspaces is used to obtain a precise estimate of the fundamental energy of the model. The energies of the first 500 states that we have calculated are in perfect agreement with the expected asymptotic behavior.

preprint2011arXiv

Solution to the Equations of the Moment Expansions

We develop a formula for matching a Taylor series about the origin and an asymptotic exponential expansion for large values of the coordinate. We test it on the expansion of the generating functions for the moments and connected moments of the Hamiltonian operator. In the former case the formula produces the energies and overlaps for the Rayleigh-Ritz method in the Krylov space. We choose the harmonic oscillator and a strongly anharmonic oscillator as illustrative examples for numerical test. Our results reveal some features of the connected-moments expansion that were overlooked in earlier studies and applications of the approach.

preprint2010arXiv

Spectroscopy of annular drums and quantum rings: perturbative and nonperturbative results

We obtain systematic approximations to the states (energies and wave functions) of quantum rings (annular drums) of arbitrary shape by conformally mapping the annular domain to a simply connected domain. Extending the general results of Ref.\cite{Amore09} we obtain an analytical formula for the spectrum of quantum ring of arbirtrary shape: for the cases of a circular annulus and of a Robnik ring considered here this formula is remarkably simple and precise. We also obtain precise variational bounds for the ground state of different quantum rings. Finally we extend the Conformal Collocation Method of \cite{Amore08,Amore09} to the class of problems considered here and calculate precise numerical solutions for a large number of states ($\approx 2000$).

preprint2010arXiv

The string of variable density: further results

We analyze the problem of calculating the solutions and the spectrum of a string with arbitrary density and fixed ends. We build a perturbative scheme which uses a basis of WKB-type functions and obtain explicit expressions for the eigenvalues and eigenfunctions of the string. Using this approach we show that it is possible to derive the asymptotic (high energy) behavior of the string, obtaining explicit expressions for the first three coefficients (the first two can also be obtained with the WKB method). Finally using an iterative approach we also obtain analytical expressions for the low energy behavior of the eigenvalues and eigenfunctions of a string with rapidly oscillating density, recovering (in a simpler way) results in the literature.

preprint2010arXiv

The string of variable density: perturbative and non-perturbative results

We obtain systematic approximations for the modes of vibration of a string of variable density, which is held fixed at its ends. These approximations are obtained iteratively applying three theorems which are proved in the paper and which hold regardless of the inhomogeneity of the string. Working on specific examples we obtain very accurate approximations which are compared both with the results of WKB method and with the numerical results obtained with a collocation approach. Finally, we show that the asymptotic behaviour of the energies of the string obtained with perturbation theory, worked to second order in the inhomogeinities, agrees with that obtained with the WKB method and implies a different functional dependence on the density that in two and higher dimensions.

preprint2009arXiv

Can one hear the density of a drum? Weyl's law for inhomogeneous media

We generalize Weyl's law to inhomogeneous bodies in $d$ dimensions. Using a perturbation scheme recently obtained by us in Ref. \cite{Amore09}, we have derived an explicit formula, which describes the asymptotic behavior of the eigenvalues of the negative laplacian on a closed $d$-dimensional cubic domain, either with Dirichlet or Neumann boundary conditions. For homogeneous bodies, the leading term in our formula reduces to the standard expression for Weyl's law. We have also used Weyl's conjecture to obtain a non-perturbative extension of our formula and we have compared our analytical results with the precise numerical results obtained using the Conformal Collocation Method of Refs. \cite{Amore08,Amore09}.

preprint2009arXiv

Collocation method for fractional quantum mechanics

We show that it is possible to obtain numerical solutions to quantum mechanical problems involving a fractional Laplacian, using a collocation approach based on Little Sinc Functions (LSF), which discretizes the Schrödinger equation on a uniform grid. The different boundary conditions are naturally implemented using sets of functions with the appropriate behavior. Good convergence properties are observed. A comparison with results based on a WKB analysis is performed.

preprint2009arXiv

Spectroscopy of drums and quantum billiards: perturbative and non-perturbative results

We develop powerful numerical and analytical techniques for the solution of the Helmholtz equation on general domains. We prove two theorems: the first theorem provides an exact formula for the ground state of an arbirtrary membrane, while the second theorem generalizes this result to any excited state of the membrane. We also develop a systematic perturbative scheme which can be used to study the small deformations of a membrane of circular or square shapes. We discuss several applications, obtaining numerical and analytical results.

preprint2009arXiv

Variational collocation for systems of coupled anharmonic oscillators

We have applied a collocation approach to obtain the numerical solution to the stationary Schrödinger equation for systems of coupled oscillators. The dependence of the discretized Hamiltonian on scale and angle parameters is exploited to obtain optimal convergence to the exact results. A careful comparison with results taken from the literature is performed, showing the advantages of the present approach.

preprint2003arXiv

Improved Lindstedt-Poincare method for the solution of nonlinear problems

We apply the Linear Delta Expansion (LDE) to the Lindstedt-Poincare (``distorted time'') method to find improved approximate solutions to nonlinear problems. We find that our method works very well for a wide range of parameters in the case of the anharmonic oscillator (Duffing equation) and of the non-linear pendulum. The approximate solutions found with this method are better behaved and converge more rapidly to the exact ones than in the simple Lindstedt-Poincaré method.