Source author record

Amlan K. Roy

Amlan K. Roy 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

35works
7topics
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

35 published item(s)

preprint2026arXiv

Quantum information and statistical complexity of hydrogen-like ions in Dunkl-Schrödinger system

In this work, we present analytical solutions of Schrödinger equation for Coulomb potential in presence of a Dunkl reflection operator. Expressions are offered for eigenvalues, eigenfunctions and radial densities for H-isoelectronic series (Z=1-3). The degeneracy in energy in absence and presence of the reflection has been discussed. The standard deviation, Shannon entropy, Rényi entropy in position space have been derived for arbitrary quantum states. Then several important complexity measures like López-Ruiz-Mancini-Calbet (LMC), Shape-Rényi complexity (SRC), Generalized Rényi complexity (GRC), Rényi complexity ratio (RCR) are considered in the analytical framework. Representative results are given for three one-electron atomic ions in tabular and graphical format. Changes in these measures with respect to parity and Dunkl parameter have been given in detail. Most of these results are offered here for the first time.

preprint2024arXiv

Information entropy in excited states in confined quantum systems

The present contribution constitutes a brief account of information theoretical analysis in several representative model as well as real quantum mechanical systems. There has been an overwhelming interest to study such measures in various quantum systems, as evidenced by a vast amount of publications in the literature that has taken place in recent years. However, while such works are numerous in so-called \emph{free} systems, there is a genuine lack of these in their constrained counterparts. With this in mind, this chapter will focus on some of the recent exciting progresses that has been witnessed in our laboratory \cite{sen06,roy14mpla,roy14mpla_manning,roy15ijqc, roy16ijqc, mukherjee15,mukherjee16,majumdar17,mukherjee18a,mukherjee18b,mukherjee18c,mukherjee18d,majumdar20,mukherjee21,majumdar21a, majumdar21b}, and elsewhere, with special emphasis on following prototypical systems, namely, (i) double well (DW) potential (symmetric and asymmetric) (ii) \emph{free}, as well as a \emph{confined hydrogen atom} (CHA) enclosed in a spherical impenetrable cavity (iii) a many-electron atom under similar enclosed environment.

preprint2022arXiv

A real-time TDDFT scheme for strong-field interaction in Cartesian coordinate grid

In this communication, we present a new approach towards RT-TDDFT through time-dependent KS equations based on an \emph{adiabatic eigenstate subspace} (AES) procedure. It introduces a second-order split operator technique in energy representation to implement the approximate TD propagator in AES. Most of the elements in TDKS matrix are directly computed in Cartesian coordinate grid (CCG). To demonstrate the internal consistency of our proposed scheme, we computed the TD dipole moment and high harmonic generation spectra using an adiabatic local density approximation. The comparison with available theoretical results ensures the feasibility of this proposed route.

preprint2022arXiv

A self-consistent systematic optimization of range-separated hybrid functionals from first principles

In this communication, we represent a self-consistent systematic optimization procedure for the development of optimally tuned (OT) range-separated hybrid (RSH) functionals from \emph{first principles}. This is an offshoot of our recent work, which employed a purely numerical approach for efficient computation of exact exchange contribution in the conventional global hybrid functionals through a range-separated (RS) technique. We make use of the size-dependency based ansatz i.e., RS parameter, $γ$, is a functional of density, $ρ(\mathbf{r})$, of which not much is known. To be consistent with this ansatz, a novel procedure is presented that relates the characteristic length of a given system (where $ρ(\mathbf{r})$ exponentially decays to zero) with $γ$ self-consistently via a simple mathematical constraint. In practice, $γ_{\mathrm{OT}}$ is obtained through an optimization of total energy as follows: $γ_{\mathrm{OT}} \equiv \underset{γ}{\mathrm{opt}} \ E_{\mathrm{tot},γ}$. It is found that the parameter $γ_{\mathrm{OT}}$, estimated as above can show better performance in predicting properties (especially from frontier orbital energies) than conventional respective RSH functionals, of a given system. We have examined the nature of highest fractionally occupied orbital from exact piece-wise linearity behavior, which reveals that this approach is sufficient to maintain this condition. A careful statistical analysis then illustrates the viability and suitability of the current approach. All the calculations are done in a Cartesian-grid based pseudopotential (G)KS-DFT framework.

preprint2022arXiv

Analytical solution of $D$ dimensional Schrödinger equation for Eckart potential with a new improved approximation in centrifugal term

Analytical solutions are presented for eigenvalues, eigenfunctions of {\color{red} D-dimensional Schrodinger equation having Eckart potential} within Nikiforov-Uvarov method. This uses a new, improved approximation for centrifugal term, from a combination of Greene-Aldrich and Pekeris approximations. Solutions are obtained in terms of hypergeometric functions. It facilitates an accurate representation in entire domain. Its validity is illustrated for energies in an arbitrary $\ell \neq 0$ quantum state. Results are compared for a chosen set of potential parameters in different dimensions. In short, a simple accurate approximation is offered for Eckart and other potentials in quantum mechanics, in higher dimension.

preprint2022arXiv

Confined H$^-$ ion within a density functional framework

Ground and excited states of a confined negative Hydrogen ion has been pursued under Kohn-Sham density functional approach by invoking a physically motivated work-function-based exchange potential. The exchange-only results are of near Hartree-Fock quality. Local parameterised Wigner-type, and gradient- and Laplacian-dependent non-local Lee-Yang-Parr functionals are chosen to investigate the electron correlation effects. Eigenfunctions and eigenvalues are extracted by using a generalized pseudospectral method obeying Dirichlet boundary condition. Energy values are reported for 1s$^{2}$ ($^{1}$S), 1s2s ($^{3,1}$S) and 1s2p ($^{3,1}$P) states. Performance of the correlation functionals in the context of confinement is examined critically. The present results are in excellent agreement with available literature. Additionally, Shannon entropy and Onicescu energy are offered for ground and low lying singly excited 1s2s ($^{3}$S) and 1s2p ($^{3}$P) states. The influence of electron correlation is more predominant in the weaker confinement limit and it decays with an increase in confinement strength. In essence, energy and some information measures are estimated using a newly formulated density functional strategy.

preprint2022arXiv

Density functional study of atoms spatially confined inside a hard sphere

An atom placed inside a cavity of finite dimension offers many interesting features, and thus has been a topic of great current activity. This work proposes a density functional approach to pursue both ground and excited states of a multi-electron atom under a spherically impenetrable enclosure. The radial Kohn-Sham (KS) equation has been solved by invoking a physically motivated work-function-based exchange potential, which offers near-Hartree-Fock-quality results. Accurate numerical eigenfunctions and eigenvalues are obtained through a generalized pseudospectral method (GPS) fulfilling the Dirichlet boundary condition. Two correlation functionals, \emph{viz.,} (i) simple, parametrized local Wigner-type, and (ii) gradient- and Laplacian-dependent non-local Lee-Yang-Parr (LYP) functionals are adopted to analyze the electron correlation effects. Preliminary exploratory results are offered for ground states of He-isoelectronic series ($Z=2-4$), as well as Li and Be atom. Several low-lying singly excited states of He atom are also reported. These are compared with available literature results -- which offers excellent agreement. Radial densities as well as expectation values are also provided. The performance of correlation energy functionals are discussed critically. In essence, this presents a simple, accurate scheme for studying atomic systems inside a \emph{hard} spherical box within the rubric of KS density functional theory.

preprint2022arXiv

Hydrogen-like ions in plasma environment

The behavior of H-like ions embedded in astrophysical plasmas in the form of \emph{dense, strongly and weakly coupled} plasmas are investigated. In these, the increase and decrease in temperature is impacted with a change in confinement radius $(r_{c})$. Two independent and generalized scaling ideas have been applied to modulate the effect of plasma screening constant ($λ$) and charge of ion ($Z$) on such systems. Several new relations are derived to interconnect the original Hamiltonian and two scaled Hamiltonians. In exponential cosine screened Coulomb potential (ECSCP) (dense) and weakly coupled plasma (WCP) these scaling relations have provided a linear equation connecting the critical screening constant $(λ^{(c)})$ and $Z$. Their ratio offers a state-dependent constant, beyond which, a particular state vanishes. Shannon entropy has been employed to understand the plasma effect on the ion. With increase in $λ$, the accumulation of opposite charge surrounding the ion increases leading to a reduction in number of bound states. However, with rise in ionic charge $Z$, this effect can be delayed. The competing effect of plasma charge density ($n_e$) and temperature in WCP and ECSCP is investigated. A recently proposed simple virial-like theorem has been established for these systems. Multipole ($k=1-4$) oscillator strength (OS) and polarizabilities for these are studied considering $1s, 2s$ states. As a bonus, analytical closed-form expressions are derived for $f^{(k)}$ and $α^{(k)} (k=1-4)$ involving $1s$ and $2s$ state, for \emph{free H-like ion}.

preprint2022arXiv

Ro-vibrational energy analysis of Manning-Rosen and Pöschl-Teller potentials with a new improved approximation in the centrifugal term

Two physically important potentials (Manning-Rosen and Pöschl-Teller) are considered for the ro-vibrational energy in diatomic molecules. An improved new approximation is invoked for the centrifugal term, which is then used for their solution within the Nikiforov-Uvarov framework. This employs a recently proposed scheme, which combines the two widely used Greene-Aldrich and Pekeris-type approximations. Thus, approximate analytical expressions are derived for eigenvalues and eigenfunctions. The energies are examined with respect to two approximation parameters, $λ$ and $ν$. The original approximations are recovered for certain specials values of these two parameters. This offers a simple effective scheme for these and other relevant potentials in quantum mechanics.

preprint2022arXiv

Ro-vibrational energy and thermodynamic properties of molecules subjected to Deng-Fan potential through an improved approximation

Accurate solution of the Schrödinger equation with Deng-Fan potential is presented by means of Nikiforov-Uvarov method. A modified Pekeris-type approximation is proposed for the centrifugal term, from a linear combination of the $r \to 0$ and $r \to r_e$ limits. It can potentially offer a series of approximations (depending on an adjustable parameter $λ$). The existing approximations in the literature can then be recovered in certain special cases. Its efficiency and feasibility is demonstrated by a critical comparison of eigenvalues produced at various $λ$'s for four molecules, \emph{viz.}, H$_2$, LiH, HCl and CO. Analytical expressions are derived for energies, eigenfunctions and the thermodynamic properties such as vibrational mean free energy, vibrational free energy, vibrational entropy and vibrational specific heat. The effect of quantum correction on partition function and thermodynamic properties is discussed by including the correction up to 10th-order, for H$_2$ and LiH. The effect of $λ$ parameter on these properties is also studied.

preprint2022arXiv

Shell-confined atom and plasma: incidental degeneracy, metallic character and information entropy

Shell confined atom can serve as a generalized model to explain both \emph{free} and \emph{confined} condition. In this scenario, an atom is trapped inside two concentric spheres of inner $(R_{a})$ and outer $(R_{b})$ radius. The choice of $R_{a}, R_{b}$ renders four different quantum mechanical systems. In hydrogenic atom, they are termed as (a) free hydrogen atom (FHA) (b) confined hydrogen atom (CHA) (c) shell-confined hydrogen atom (SCHA) (d) left-confined hydrogen atom (LCHA). By placing $R_{a}, R_{b}$ at the location of radial nodes of respective \emph{free} $n,\ell$ states, a new kind of degeneracy may arise. At a given $n$ of FHA, there exists $\frac{n(n+1)(n+2)}{6}$ number of iso-energic states with energy $-\frac{Z^{2}}{2n^{2}}$. Furthermore, within a given $n$, the individual contribution of each of these four potentials has also been enumerated. This incidental degeneracy concept is further explored and analyzed in certain well-known \emph{plasma} (Debye and exponential cosine screened) systems. Multipole oscillator strength, $f^{(k)}$, and polarizability, $α^{(k)}$, are evaluated for (a)-(d) in some low-lying states $(k=1-4)$. In excited states, \emph{negative} polarizability is also observed. In this context, metallic behavior of H-like systems in SCHA is discussed and demonstrated. Additionally analytical closed-form expression of $f^{(k)}$ and $α^{(k)}$ are reported for $1s,2s,2p,3d,4f,5g$ states of FHA. Finally, Shannon entropy and Onicescu {\color{red}information} energies are investigated in ground state in SCHA and LCHA in both position and momentum spaces. Much of the results are reported here for first time.

preprint2021arXiv

Analysis of Compton profile through information theory in H-like atoms inside impenetrable sphere

Confinement of atoms inside various cavities has been studied for nearly eight decades. However, the Compton profile for such systems has not yet been investigated. Here we construct the Compton profile (CP) for a H atom radially confined inside a \emph{hard} spherical enclosure, as well as in \emph{free condition}. Some exact analytical relations for the CP's of circular or nodeless states of free atom is presented. By means of a scaling idea, this has been further extended to the study of an H-like atom trapped inside an impenetrable cavity. The accuracy of these constructed CP has been confirmed by computing various momentum moments. Apart from that, several information theoretical measures, like Shannon entropy ($S$) and Onicescu energy ($E$) have been exploited to characterize these profiles. Exact closed form expressions are derived for $S$ and $E$ using the ground state CP in free H-like atoms. A detailed study reveals that, increase in confinement inhibits the rate of dissipation of kinetic energy. At a fixed $\ell$, this rate diminishes with rise in $n$. However, at a certain $n$, this rate accelerates with progress in $\ell$. A similar analysis on the respective free counterpart displays an exactly opposite trend as that in confined system. However, in both free and confined environments, CP generally gets broadened with rise in $Z$. Representative calculations are done numerically for low-lying states of the confined systems, taking two forms of position-space wave functions: (a) exact (b) highly accurate eigenfunctions through a generalized pseudospectral method. In essence, CPs are reported for confined H atom (and isoelectronic series) and investigated adopting an information-theoretic framework.

preprint2021arXiv

Excitation energies through Becke's exciton model within a Cartesian-grid KS DFT

Photon-induced electronic excitations are ubiquitously observed in organic chromophore. In this context, we present a simple, alternative time-independent DFT procedure, for computation of single-particle excitation energies, in particular, the lower bound excited singlet states, which are of primary interest in photochemistry. This takes inspiration from recently developed Becke's exciton model, where a key step constitutes the accurate evaluation of correlated singlet-triplet splitting energy. It introduces a non-empirical model, both from "adiabatic connection theorem" and "virial theorem" to analyze the role of 2e$^-$ integral in such calculations. The latter quantity is efficiently mapped onto a real grid and computed accurately using a purely numerical strategy. Illustrative calculations are performed on 10 $π$-electron organic chromophores within a Cartesian-grid implementation of pseudopotential Kohn-Sham (KS) DFT, developed in our laboratory, taking SBKJC-type basis functions within B3LYP approximation. The triplet and singlet excitation energies corresponding to first singly excited configuration, are found to be in excellent agreement with TD-B3LYP calculations. Further, we perform the same for a set of larger molecular systems using the asymptotically corrected LC-BLYP, in addition to B3LYP. A systematic comparison with theoretical best estimates demonstrates the viability and suitability of current approach in determining optical gaps, combining predictive accuracy with moderate computational cost.

preprint2021arXiv

Shannon entropy in confined He-like ions within a density functional formalism

Shannon entropy in position ($S_{\rvec}$) and momentum ($S_{\pvec}$) spaces, along with their sum ($S_t$) are presented for unit-normalized densities of He, Li$^+$ and Be$^{2+}$ ions, spatially confined at the center of an impenetrable spherical enclosure defined by a radius $r_c$. Both ground as well as some selected low-lying singly excited states, \emph{viz.,} 1sns (n $=$ 2-4) $^3$S, 1snp (n $=$ 2-3) $^3$P, 1s3d $^3$D are considered within a density functional methodology that makes use of a work-function-based exchange potential along with two correlation potentials (local Wigner-type parametrized functional as well as the more involved non-linear gradient- and Laplacian-dependent Lee-Yang-Parr functional). The radial Kohn-Sham (KS) equation is solved using an optimal spatial discretization scheme via the generalized pseudospectral (GPS) method. A detailed systematic analysis of the confined system (relative to corresponding free system) has been performed for these quantities with respect to $r_c$ in tabular and graphical forms, \emph{with and without} electron correlation. Due to compression, the pattern of entropy in aforementioned states gets characterized by various crossovers at intermediate and lower $r_c$ regions. The impact of electron correlation is more pronounced in weaker confinement limit, and appears to decay with rise in confinement strength. The exchange-only results are quite good to provide a decent qualitative discussion. The lower-bounds provided by entropic uncertainty relation holds good in all cases. Several other new interesting features are observed.

preprint2014arXiv

Confinement in 3D polynomial oscillators through a generalized pseudospectral method

Spherical confinement in 3D harmonic, quartic and other higher oscillators of even order is studied. The generalized pseudospectral method is employed for accurate solution of relevant Schrödinger equation in an \emph{optimum, non-uniform} radial grid. Eigenvalues, eigenfunctions, position expectation values, radial densities in \emph{low and high-lying} states are presented in case of \emph{small, intermediate and large} confinement radius. The \emph{degeneracy breaking} in confined situation as well as correlation in its \emph{energy ordering} with respect to the respective unconfined counterpart is discussed. For all instances, current results agree excellently with best available literature results. Many new states are reported here for first time. In essence, a simple, efficient method is provided for accurate solution of 3D polynomial potentials enclosed within spherical impenetrable walls.

preprint2014arXiv

Ro-vibrational spectroscopy of molecules represented by a Tietz-Hua oscillator potential

Accurate low and high-lying bound states of Tietz-Hua oscillator potential are presented. The radial Schrödinger equation is solved efficiently by means of the generalized pseudospectral method that enables optimal spatial discretization. Both $\ell=0$ and rotational states are considered. Ro-vibrational levels of six diatomic molecules \emph{viz.}, H$_2$, HF, N$_2$, NO, O$_2$, O$_2^+$ are obtained with good accuracy. Most of the states are reported here for the first time. A detailed analysis of variation of eigenvalues with $n, \ell$ quantum numbers is made. Results are compared with literature data, wherever possible. These are also briefly contrasted with the Morse potential results.

preprint2014arXiv

Ro-vibrational studies of diatomic molecules in a shifted Deng-Fan oscillator potential

Bound-state spectra of shifted Deng-Fan oscillator potential are studied by means of a generalized pseudospectral method. Very accurate results are obtained for \emph{both low as well as high} states by a non-uniform optimal discretization of the radial Schrödinger equation. Excellent agreement with literature data is observed in \emph{both $s$-wave and rotational} states. Detailed variation of energies with respect to potential parameters is discussed. Application is made to the ro-vibrational levels of four representative diatomic molecules (H$_2$, LiH, HCl, CO). Nine states having $\{n,\ell\} =0,1,2$ are calculated with good accuracy along with 15 other higher states for each of these molecules. Variation of energies with respect to state indices $n$, $\ell$ show \emph{behavior} similar to that in the Morse potential. Many new states are reported here for the first time. In short, a simple, accurate and efficient method is presented for this and other similar potentials in molecular physics.

preprint2014arXiv

Studies on bound-state spectra of Manning-Rosen potential

Accurate ro-vibrational energies, eigenfunctions, radial densities, expectation values are presented for the exponential-type Manning-Rosen (MR) potential. Bound states accurate up to ten significant figure are obtained by employing a simple, reliable generalized pseudospectral method. \emph{All} 55 eigenstates with $n \leq 10$ are treated for arbitrary values of potential parameters, covering a wide range of interaction, through a \emph{non-uniform, optimal} spatial radial discretization. A detailed investigation has been made on energy changes with respect to \emph{screening and other} potential parameters. A systematic estimation of \emph{critical} screening parameters are given for these eigenstates. Special emphasis has been given to \emph{higher} states and in the vicinity of \emph{critical screening} region. A thorough comparison with literature results is made wherever possible. This \emph{surpasses} the accuracy of \emph{all} other existing methods currently available. Several \emph{new} states are reported for the first time. In short, a simple, efficient scheme for accurate calculation of this and other molecular potentials is offered.

preprint2014arXiv

Studies on the bound-state spectrum of hyperbolic potential

Bound states of hyperbolic potential is investigated by means of a generalized pseudospectral method. Significantly improved eigenvalues, eigenfunctions are obtained efficiently for arbitrary $n, \ell$ quantum states by solving the relevant non-relativistic Schrödinger equation allowing a non-uniform, optimal spatial discretization. Eigenvalues accurate up to tenth decimal place are reported for a large range of potential parameters; thus covering a wide range of interaction. Excellent agreement with available literature results is observed in all occasions. Special attention is paid for higher states. Some new states are given. Energy variations with respect to parameters in the potential are studied in considerable detail for the first time.

preprint2013arXiv

A new DFT method for atoms and molecules in Cartesian grid

Electronic structure calculation of atoms and molecules, in the past few decades has largely been dominated by density functional methods. This is primarily due to the fact that this can account for electron correlation effects in a rigorous, tractable manner keeping the computational cost at a manageable level. With recent advances in methodological development, algorithmic progress as well as computer technology, larger physical, chemical and biological systems are amenable to quantum mechanical calculations than ever before. Here we report the development of a new method for accurate reliable description of atoms, molecules within the Hohenberg-Kohn-Sham density functional theory (DFT). In a Cartesian grid, atom-centered localized basis set, electron density, molecular orbitals, two-body potentials are directly built on the grid. We employ a Fourier convolution method for classical Coulomb potentials by making an Ewald-type decomposition technique in terms of short- and long-range interactions. One-body matrix elements are obtained from standard recursion algorithms while two-body counterparts are done by direct numerical integration. A systematic analysis of our results obtained on various properties, such as component energy, total energy, ionization energy, potential energy curve, atomization energy, etc., clearly demonstrates that the method is capable of producing quite accurate and competitive (with those from other methods in the literature) results. In brief, a new variational DFT method is presented for atoms and molecules, \emph{completely} in Cartesian grid.

preprint2013arXiv

Accurate calculation of the bound states of Hellmann potential

Bound states of the Hellmann potential, which is a superposition of the attractive Coulomb ($-A/r$) and the Yukawa ($Be^{-Cr}/r$) potential, are calculated by using a generalized pseudospectral method. Energy eigenvalues accurate up to thirteen to fourteen significant figures, and densities are obtained through a nonuniform, optimal spatial discretization of the radial Schrödinger equation. Both ground and excited states are reported for arbitrary values of the potential parameters covering a wide range of interaction. Calculations have been made for higher states as well as for stronger couplings. Some new states are reported here for the first time, which could be useful for future works. The present results are significantly improved in accuracy over all other existing literature values and offers a simple, accurate and efficient scheme for these and other singular potentials in quantum mechanics.

preprint2013arXiv

Accurate ro-vibrational spectroscopy of diatomic molecules in a Morse oscillator potential

This work presents the bound-state spectra of Morse oscillator, which remains one of the oldest important model potentials for molecules. Accurate ro-vibrational energies are obtained by means of a generalized pseudospectral method that offers an optimal, non-uniform discretization of the radial grid. Both s-wave ($\ell=0$) and rotational ($\ell \neq 0$) states for low and high quantum numbers are calculated for four representative diatomic molecules, namely H$_2$, LiH, HCl and CO. First nine states belonging to a maximum of $n, \ell =2$ are computed with good accuracy, along with nine other high-lying states for each of these molecules. Present results \emph{surpass} the accuracy of \emph{all} hitherto published calculations found so far, except the tridiagonal J-matrix method, which produces similar accuracy as ours. Detailed variation of energies with respect to state indices $n,\ell$ show interesting behavior. A host of new states including the higher ones are reported as well. This offers a simple general efficient scheme for calculating these and other similar potentials in molecular physics.

preprint2013arXiv

Bound states of the generalized spiked harmonic oscillator

Bound states of the generalized spiked harmonic oscillator potential are calculated accurately by using the generalized pseudospectral method. Energy eigenvalues, various expectation values, radial densities are obtained through a nonuniform, optimal spatial discretization of the radial Schrödinger equation efficiently. Ground and excited states corresponding to arbitrary values of $n$ and $\ell$ are reported for potential parameters covering a wide range of interaction. Both weak and strong coupling is considered. The effect of potential parameters on eigenvalues and densities are discussed. Almost all of the results are reported here for the first time, which could be useful for future studies.

preprint2013arXiv

Calculation of the bound states of power-law and logarithmic potentials through a generalized pseudospectral method

Bound states of the power-law and logarithmic potentials are calculated using a generalized pseudospectral method. The solution of the single-particle Schrödinger equation in a nonuniform and optimal spatial discretization offers accurate eigenvalues, densities and expectation values. The calculations are carried out for states with arbitrary $n$ and $\ell$ quantum numbers. Comparisons are made with the available literature data and excellent agreement is observed. In all the cases, the present method yields considerably improved results over the other existing calculations. Some new states are reported.

preprint2013arXiv

Calculation of the spiked harmonic oscillators through a generalized pseudospectral method

The generalized pseudospectral method is employed for the accurate calculation of eigenvalues, densities and expectation values for the spiked harmonic oscillators. This allows \emph{nonuniform} and \emph{optimal} spatial discretization of the corresponding single-particle radial Schrödinger equation satisfying the Dirichlet boundary conditions leading to the standard diagonalization of the symmetric matrices. The present results for a large range of potential parameters are in excellent agreement with those from the other accurate methods available in the literature. The ground and excited states (both low as well as high angular momentum states) are obtained with equal ease and accuracy. Some new states including the higher excited states are reported here for the first time. This offers a simple, accurate and efficient method for the treatment of these and a wide variety of other singular potentials of physical and chemical interest in quantum mechanics.

preprint2013arXiv

Density functional studies on the hollow resonances in Li-isoelectronic sequence (Z=4--10

In this sequel to our work on triply excited hollow resonances in three-electron atomic systems, a density functional theory (DFT)-based formalism is employed to investigate similar resonances in Li-isoelectronic series (Z=4--10). A combination of the work-function-based local nonvariational exchange potential and the popular gradient plus Laplacian included Lee-Yang-Parr correlation energy functional is used. First, all the 8 n=2 intrashell states of B$^{2+}$, N$^{4+}$ and F$^{6+}$ are presented, which are relatively less studied in the literature compared to the remaining 4 members. Then calculations are performed for the 8 $2l2l'$n$l"$ (3$\leq$n$\leq$6) hollow resonance series; {\em viz.,} 2s$^2$ns $^2$S$^e$, 2s$^2$np $^2$P$^o$, 2s$^2$nd $^2$D$^e$, 2s2pns $^4$P$^o$, 2s2pnp $^4$D$^e$, 2p$^2$ns $^4$P$^e$, 2p$^2$np $^4$D$^o$ and 2p$^2$ns $^2$D$^e$, of all the 7 positive ions. Next, as an illustration, higher resonance positions of the 2s$^2$ns $^2$S$^e$ series are calculated for all the ions with a maximum of n=25. The excitation energies calculated from this single-determinantal approach are in excellent agreement with the available literature data (for the n=2 intrashell states the deviation is within 0.125% and excepting only one case, the same for the resonance series is well below 0.50%). With an increase in Z, the deviations tend to decrease. Radial densities are also presented for some of the selected states. The only result available in the literature for the lower resonances (corresponding to a maximum of n=17) have been reported very recently. The n$>$16 ($>17$ for F$^{6+}$) resonances are examined here for the first time. This gives a promising viable and general DFT scheme for the accurate calculation of these and other hollow resonances in many-electron atoms.

preprint2013arXiv

Studies on some exponential-screened Coulomb potentials

The generalized pseudospectral method is employed to study the bound-state spectra of some of the exponentially screened Coulomb potentials, \emph{viz.}, the exponential cosine screened Coulomb (ECSC) and general exponential screened Coulomb (GESC) potential, with special emphasis on \emph{higher} states and \emph{stronger} interaction. Eigenvalues accurate up to eleven significant figures are obtained through a non-uniform optimal spatial discretization of the radial Schrödinger equation. All the 55 eigenstates of ECSC potential with $n \le 10$ and 36 eigenstates of GESC potential with $n \le 8$ are considered for arbitrary values of the screening parameter, covering a wide range of interaction. Excited states as high as up to $n=18$ have been computed with high accuracy for the first time. Excellent agreement with the literature data has been observed in all cases. All the GESC eigenstates are calculated with much greater accuracy than the existing methods available in literature. Many $l \ne 0$ states of this potential are reported here. In both cases, a detailed variation of energies with respect to the parameters in potential is monitored.

preprint2013arXiv

Studies on some singular potentials in quantum mechanics

A simple methodology is suggested for the efficient calculation of certain central potentials having singularities. The generalized pseudospectral method used in this work facilitates {\em nonuniform} and optimal spatial discretization. Applications have been made to calculate the energies, densities and expectation values for two singular potentials of physical interest, {\em viz.,} (i) the harmonic potential plus inverse quartic and sextic perturbation and (ii) the Coulomb potential with a linear and quadratic term for a broad range of parameters. The first 10 states belonging to a maximum of $\ell=8$ and 5 for (i) and (ii) have been computed with good accuracy and compared with the most accurate available literature data. The calculated results are in excellent agreement, especially in the light of the difficulties encountered in these potentials. Some new states are reported here for the first time. This offers a general and efficient scheme for calculating these and other similar potentials of physical and mathematical interest in quantum mechanics accurately.

preprint2013arXiv

Studies on the hollow states of atomic lithium by a density functional approach

Density functional calculations are performed for twelve $2l2l'$n$l"$ (n$\geq$2) triply excited hollow resonance series of Li, {\em viz.,} 2s$^2$ns $^2$S$^e$, 2s$^2$np $^2$P$^o$, 2s$^2$nd $^2$D$^e$, 2p$^2$ns $^2$D$^e$,$^4$P$^e$, 2s2pns $^4$P$^o$, 2s2pnp $^4$D$^e$, 2p$^2$np $^2$F$^o$,$^4$D$^o$, 2p$^2$nd $^2$G$^e$, $^4$F$^e$ and 2s2pnd $^4$F$^o$, covering a total of about 270 low-, moderately high- and high-lying states, with n as high as up to 25. The work-function-based exchange potential and the nonlinear gradient plus Laplacian included Lee-Yang-Parr correlation energy functional is used. The relevant Kohn-Sham-type equation is solved numerically using the generalized pseudospectral method offering nonuniform, optimal spatial discretization to obtain the orbitals and densities. Except for the one state, the discrepancy in the calculated state energies remains well within 0.98%, whereas the excitation energies deviate by 0.02--0.58% compared to the available experimental and other theoretical results. Additionally companion calculations are also presented for the 37 $3l3l'$n$l"$ (n$\geq$3) doubly hollow states (seven are n=3 intrashell type) of Li with both K and L shells empty (up to n=6) in the photon energy range 175.63--180.51 eV, with varying symmetries and multiplicities. Our calculation shows good agreement with the recent literature data for the only two such doubly hollow states reported so far, {\em viz.,} 3s$^2$3p $^2$P$^o$ and 3s3p$^2$ $^4$P$^e$. The resonance series are found to be inextricably entangled to each other, leading to complicated behavior in their positions. Many new states are reported here for the first time. This provides a simple, efficient and general scheme for the accurate calculation of these and other multiply excited Rydberg series of many-electron atomic systems within density functional theory.

preprint2013arXiv

The generalized pseudospectral approach to the bound states of Hulthen and Yukawa potentials

The generalized pseudospectral method is employed to calculate the bound states of Hulthén and Yukawa potentials in quantum mechanics, with special emphases on higher excited states and stronger couplings. Accurate energy eigenvalues, expectation values and radial probability densities are obtained through a nonuniform and optimal spatial discretization of the radial Schrödinger equation. Results accurate up to thirteen to fourteen significant figures are reported for all the 55 eigenstates of both these potentials with $n\leq$10 for arbitrary values of the screening parameters covering a wide range of interaction. Furthermore, excited states as high as up to $n=17$ have been computed with good accuracy for both these potentials. Excellent agreement with the available literature data has been observed in all cases. The $n>6$ states of Yukawa potential has been considerably improved over all other existing results currently available, while the same for Hulthén potential are reported here for the first time. Excepting the $1s$ and $2s$ states of Yukawa potential, the present method surpasses in accuracy all other existing results in the stronger coupling region for all other states of both these systems. This offers a simple and efficient scheme for the accurate calculation of these and other screened Coulomb potentials.

preprint2010arXiv

Bound state spectra of the 3D rational potential

We present bound state spectra of the 3D rational potential, $V(r)=r^2 + λr^2/(1+gr^2)$, $g>0$, by means of the generalized pseudospectral method. All the thirty states corresponding to $n$=0--9 are considered for the first time for a broad range of coupling parameters. These results surpass the accuracy of \emph{all} other existing calculations published so far except the finite-difference method, which yields similar accuracy as ours. Variation of energies and radial distribution functions is followed with respect to the interaction parameters. Special emphasis has been laid on \emph{higher} excitations and \emph{negative} values of the interaction, where relatively less work has been reported. The energy sequence is found to be different for positive and negative interaction; numerically following a mirror-image relationship \emph{usually}, if not always. Additionally, twenty energy splittings arising from certain levels belonging to $n$=0--9 are systematically studied as functions of the potential parameters. Several new states (including the higher ones) are presented.

preprint2010arXiv

Density functional calculation of many-electron systems in cartesian coordinate grid

A recently developed density functional method, within Hohenberg-Kohn-Sham framework, is used for faithful description of atoms, molecules in Cartesian coordinate grid, by using an LCAO-MO ansatz. Classical Coulomb potential is obtained by means of a Fourier convolution technique. All two-body potentials (including exchange-correlation (XC)) are constructed directly on real grid, while their corresponding matrix elements are computed from numerical integration. Detailed systematic investigation is made for a representative set of atoms/molecules through a number of properties like total energies, component energies, ionization energies, orbital energies, etc. Two nonlocal XC functionals (FT97 and PBE) are considered for pseudopotential calculation of 35 species while preliminary all-electron results are reported for 6 atoms using the LDA XC density functional. Comparison with literature results, wherever possible, exhibits near-complete agreement. This offers a simple efficient route towards accurate reliable calculation of many-electron systems in the Cartesian grid. Future prospect of this method is also discussed.

preprint2010arXiv

Ground and excited states of Li$^-$, Be$^-$ through a density-based approach

Density functional calculations are performed for ground [He]2s$^2$ $^1$S$^e$, and three metastable bound excited states, 1s2s2p$^2$ $^5$P$^e$, 1s2p$^3$ $^5$S$^o$, 1s2s2p3p $^5$P$^e$ of Li$^-$ and [He]2s2p$^2$ $^4$P$^e$, [He]2p$^3$ $^4$S$^o$, 1s2s2p$^3$ $^6$S$^o$ of Be$^-$ each. The work-function-based exchange potential is used, while the correlation effects are included by employing the Lee-Yang-Parr potential. The relevant nonrelativistic KS equation is solved by means of a generalized pseudospectral discretization scheme offering nonuniform and optimal spatial grid. Computed total energies, radial densities, selected density moments, as well as two transition wavelengths (1s2s2p$^2$ $^5$P$^e \to$1s2p$^3$ $^5$S$^o$ of Li$^-$, [He]2s2p$^2$ $^4$P$^e \to$ [He]2p$^3$ $^4$S$^o$ of Be$^-$) show reasonably good agreement with the available theoretical and experimental data. The term energies show an absolute deviation of 0.007--0.171% with the largest deviation being observed for the even-parity $^5$P state of Li$^-$. The transition wavelengths of Li$^-$, Be$^-$ are calculated within 0.891 and 0.438% of the experimental values. This offers a simple practical route towards accurate reliable calculation of excited states of anions within density functional theory.

preprint2010arXiv

Pseudopotential density functional treament of atoms and molecules in cartesian coordinate grid

This is a follow-up of our recently proposed work on pseudopotential calculation (Ref. [21]) of atoms and molecules within DFT framework, using cartesian coordinate grid. Detailed results are presented to demonstrate the usefulness, applicability of the same for a larger set of species (5 atoms; 53 molecules) and exchange-correlation functionals (local, nonlocal). A thorough comparison on total, component, ionization, atomization energies, eigenvalues, potential energy curves with available literature data shows excellent agreement. Additionally, HOMO energies for a series of molecules show significant improvements by using the Leeuwen-Baerends exchange potential, compared to other functionals considered. Comparison with experiments has been made, wherever possible.

preprint2006arXiv

Spherically confined isotropic harmonic oscillator

The generalized pseudospectral Legendre method is used to carry out accurate calculations of eigenvalues of the spherically confined isotropic harmonic oscillator with impenetrable boundaries. The energy of the confined state is found to be equal to that of the unconfined state when the radius of confinement is suitably chosen as the location of the radial nodes in the unconfined state. This incidental degeneracy condition is numerically shown to be valid in general. Further, the full set of pairs of confined states defined by the quantum numbers [(n+1, \ell) ; (n, \ell+2)], n = 1,2,.., and with the radius of confinement {(2 \ell +3)/2}^{1/2} a.u., which represents the single node in the unconfined (1, \ell) state, is found to display a constant energy level separation exactly given by twice the oscillator frequency. The results of similar numerical studies on the confined Davidson oscillator with impenetrable boundary as well as the confined isotropic harmonic oscillator with finite potential barrier are also reported .The significance of the numerical results are discussed.