@Article{bednarz_jcp_122_164302_2005,
Title = {Non-Born--Oppenheimer variational calculations of HT$^+$ bound states with zero angular momentum},
Author = {Eugeniusz Bednarz and Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2005},
Number = {16},
Pages = {164302},
Volume = {122},
Abstract = {We report fully nonadiabatic calculations of all rotationless bound states of HT$^+$ molecular ion ($t^+ p^+ e^-$) carried out in the framework of the variational method. We show that, in all the states, except the two highest ones, the bond in the system can be described as covalent. In the highest two states the bond becomes essentially ionic and HT$^+$ can be described as a T + H$^+$ complex. The wave function of the system was expanded in terms of spherically symmetric, explicitly correlated Gaussian functions with preexponential multipliers consisting of powers of the internuclear distance. Apart from the total energies of the states, we have calculated the expectation values of the $t-p$, $t-e$, and $p-e$ interparticle distances, their squares, and the nucleus-nucleus correlation functions.},
Doi = {10.1063/1.1884602},
Keywords = {variational techniques; bound states; Gaussian processes; wave functions; hydrogen ions; functional analysis; positive ions; bonds (chemical)},
Numpages = {10},
Publisher = {AIP}
}
@Article{bednarz_mp_103_1169_2005,
Title = {Integrals for non-Born--Oppenheimer calculations of molecules with three nuclei},
Author = {Bednarz, Eugeniusz and Bubin, Sergiy and Adamowicz, Ludwik},
Journal = {Mol. Phys.},
Year = {2005},
Pages = {1169-1182},
Volume = {103},
Abstract = {Variational molecular calculations without assuming the Born--Oppenheimer approximation of states corresponding to zero total angular momentum require the use of spherically symmetric basis functions explicitly depending on the electron--electron, electron--nucleus, and nucleus--nucleus distances. In our calculations, such functions have been the explicitly correlated Gaussians. For molecules with three nuclei, the Gaussians have to be multiplied by powers of all three internuclear distances to describe the highly correlated motion of the nuclei. In this work we have derived formulae for the overlap and the Hamiltonian matrix elements for such basis functions. Implementation of the formulae presents some unconventional numerical difficulties related to maintaining the required precision of the calculation. The implementation problems are discussed.},
Doi = {10.1080/00268970512331339314},
Publisher = {Taylor \& Francis}
}
@Article{bubin_jcp_118_3079_2003,
Title = {Variational calculations of excited states with zero total angular momentum (vibrational spectrum) of H$_2$ without use of the Born--Oppenheimer approximation},
Author = {Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2003},
Number = {7},
Pages = {3079-3082},
Volume = {118},
Abstract = {Very accurate, rigorous and fully variational, all-particle, non-Born--Oppenheimer calculations of the vibrational spectrum of the H$_2$ molecule have been performed. Very high accuracy has been achieved by expanding the wave functions in terms of explicitly correlated Gaussian functions with preexponential powers of the internuclear distance. An indicator of the high accuracy of the calculations is the new upper bound for the H$_2$ nonrelativistic nonadiabatic ground state energy equal to -1.1640250300 hartree.},
Doi = {10.1063/1.1537719},
Keywords = {hydrogen neutral molecules; wave functions; excited states; vibrational states; angular momentum; variational techniques},
Publisher = {AIP}
}
@Article{bubin_jcp_120_6051_2004,
Title = {Non-Born--Oppenheimer study of positronic molecular systems: $e^+$LiH},
Author = {Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2004},
Number = {13},
Pages = {6051-6055},
Volume = {120},
Abstract = {Very accurate non-Born--Oppenheimer variational calculations of the ground state of $e^+$LiH have been performed using explicitly correlated Gaussian functions with preexponential factors dependent on powers of the internuclear distance. In order to determine the positron detachment energy of $e^+$LiH and the dissociation energy corresponding to the $e^+$LiH fragmentation into HPs and Li$^+$ we also calculated non-BO energies of HPs, LiH, and Li$^+$. For all the systems the calculations provided the lowest ever-reported variational upper-bounds to the ground state energies. Annihilation rates of HPs and $e^+$LiH were also computed. The dissociation energy of $e^+$LiH into HPs and Li$^+$ was determined to be 0.036548 hartree.},
Doi = {10.1063/1.1651056},
Keywords = {lithium compounds; dissociation energies; ground states; variational techniques; positron annihilation; molecule-positron collisions},
Publisher = {AIP}
}
@Article{bubin_jcp_121_6249_2004,
Title = {Nonrelativistic molecular quantum mechanics without approximations: Electron affinities of LiH and LiD},
Author = {Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2004},
Number = {13},
Pages = {6249-6253},
Volume = {121},
Abstract = {We took the complete nonrelativistic Hamiltonians for the LiH and LiH$^-$ systems, as well as their deuterated isotopomers, we separated the kinetic energy of the center of mass motion from the Hamiltonians, and with the use of the variational method we optimized the ground-state nonadiabatic wave functions for the systems expanding them in terms of $n$-particle explicitly correlated Gaussian functions. With 3600 functions in the expansions we obtained the lowest ever ground-state energies of LiH, LiD, LiH$^-$, and LiD$^-$ and these values were used to determine LiH and LiD electrons affinities (EAs) yielding 0.33030 and 0.32713 eV, respectively. The present are the first high-accuracy ab initio quantum mechanical calculations of the LiH and LiD EAs that do not assume the Born-Oppenheimer approximation. The obtained EAs fall within the uncertainty brackets of the experimental results.},
Doi = {10.1063/1.1786580},
Keywords = {lithium compounds; electron affinity; negative ions; ab initio calculations; wave functions; ground states; isomerism; isotope effects; variational techniques; functional analysis},
Publisher = {AIP}
}
@Article{bubin_cpl_403_185_2005,
Title = {Nucleus-nucleus correlation function in non-Born-Oppenheimer molecular calculations: vibrationally excited states of HD$^+$},
Author = {Sergiy Bubin and Ludwik Adamowicz},
Journal = {Chem. Phys. Lett.},
Year = {2005},
Number = {1-3},
Pages = {185 - 191},
Volume = {403},
Abstract = {In the present work, we studied HD$^+$ ($dpe$) molecular ion in the framework of variational method without assuming the Born-Oppenheimer approximation. The non-adiabatic wave function was expanded in terms of explicitly correlated Gaussian basis functions. An algorithm for calculating the nucleus-nucleus correlation function (i.e., the probability density of one nucleus in the reference frame where the other one is at the origin) was derived, implemented, and used to depict all bound states of HD$^+$ with zero rotational energy.},
Doi = {10.1016/j.cplett.2005.01.009},
}
@Article{bubin_jcp_124_224317_2006,
Title = {Matrix elements of $N$-particle explicitly correlated Gaussian basis functions with complex exponential parameters},
Author = {Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2006},
Number = {22},
Pages = {224317},
Volume = {124},
Abstract = {In this work we present analytical expressions for Hamiltonian matrix elements with spherically symmetric, explicitly correlated Gaussian basis functions with complex exponential parameters for an arbitrary number of particles. The expressions are derived using the formalism of matrix differential calculus. In addition, we present expressions for the energy gradient that includes derivatives of the Hamiltonian integrals with respect to the exponential parameters. The gradient is used in the variational optimization of the parameters. All the expressions are presented in the matrix form suitable for both numerical implementation and theoretical analysis. The energy and gradient formulas have been programed and used to calculate ground and excited states of the He atom using an approach that does not involve the Born-Oppenheimer approximation.},
Doi = {10.1063/1.2204605},
Keywords = {helium neutral atoms; variational techniques; ground states; excited states; matrix algebra; Gaussian processes},
Numpages = {14},
Publisher = {AIP}
}
@Article{bubin_jcp_125_064309_2006,
Title = {Non-Born-Oppenheimer variational calculation of the ground-state vibrational spectrum of LiH$^+$},
Author = {Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2006},
Number = {6},
Pages = {064309},
Volume = {125},
Abstract = {Very accurate, rigorous, variational, non-Born-Oppenheimer (non-BO) calculations have been performed for the fully symmetric, bound states of the LiH$^+$ ion. These states correspond to the ground and excited vibrational states of LiH$^+$ in the ground electronic state. The non-BO wave functions of the states have been expanded in terms of spherical $N$-particle explicitly correlated Gaussian functions multiplied by even powers of the internuclear distance and 5600 Gaussians were used for each state. The calculations that, to our knowledge, are the most accurate ever performed for a diatomic system with three electrons have yielded six bound states. Average interparticle distances and nucleus-nucleus correlation function plots are presented.},
Doi = {10.1063/1.2244563},
Eid = {064309},
Keywords = {positive ions; vibrational states; variational techniques; ground states; excited states; bound states; lithium compounds; wave functions},
Numpages = {4},
Publisher = {AIP}
}
@Article{bubin_pra_74_052502_2006,
Title = {Nonrelativistic variational calculations of the positronium molecule and the positronium hydride},
Author = {Bubin, Sergiy and Adamowicz, Ludwik },
Journal = {Phys. Rev. A},
Year = {2006},
Number = {5},
Pages = {052502},
Volume = {74},
Abstract = {Highly accurate ground-state nonrelativistic variational calculations of Ps$_2$, $^1$HPs, and $^\infty$HPs are reported. The calculations have been performed using 5000 explicitly correlated Gaussian basis functions and yielded the lowest variational energy upper bounds for these systems to date. The relative accuracies of the energies obtained are of the order of 4$\times$10$^{-10}$ a.u. for Ps$_2$ and 2$\times$10$^{-9}$ a.u. for HP. Several expectation values have also been computed for each system, as well as electron-positron annihilation rates.},
Doi = {10.1103/PhysRevA.74.052502},
Numpages = {5},
Publisher = {American Physical Society}
}
@Article{bubin_jcp_126_214305_2007,
Title = {Calculations of the ground states of BeH and BeH$^+$ without the Born-Oppenheimer approximation},
Author = {Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2007},
Number = {21},
Pages = {214305},
Volume = {126},
Abstract = {Non-Born-Oppenheimer variational calculations employing explicitly correlated Gaussian basis functions have been performed for the ground states of the beryllium monohydride molecule (BeH) and its ion (BeH$^+$), as well as for the beryllium atom (Be) and its ion (Be$^+$). An approach based on the analytical energy gradient calculated with respect to the Gaussian exponential parameters was employed. The calculated energies were used to determine the ionization potential of BeH and the dissociation energies of BeH and BeH$^+$. Also, the generated wave functions were used to compute various expectation values, such as the average interparticle distances and the nucleus-nucleus correlation functions.},
Doi = {10.1063/1.2736699},
Keywords = {beryllium compounds; positive ions; variational techniques; Gaussian distribution; ionisation potential; dissociation energies},
Numpages = {5},
Publisher = {AIP}
}
@Article{bubin_jcp_128_114107_2008,
Title = {Energy and energy gradient matrix elements with $N$-particle explicitly correlated complex Gaussian basis functions with $L$ = 1},
Author = {Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2008},
Number = {11},
Pages = {114107},
Volume = {128},
Abstract = {In this work we consider explicitly correlated complex Gaussian basis functions for expanding the wave function of an $N$-particle system with the $L$=1 total orbital angular momentum. We derive analytical expressions for various matrix elements with these basis functions including the overlap, kinetic energy, and potential energy (Coulomb interaction) matrix elements, as well as matrix elements of other quantities. The derivatives of the overlap, kinetic, and potential energy integrals with respect to the Gaussian exponential parameters are also derived and used to calculate the energy gradient. All the derivations are performed using the formalism of the matrix differential calculus that facilitates a way of expressing the integrals in an elegant matrix form, which is convenient for the theoretical analysis and the computer implementation. The new method is tested in calculations of two systems: the lowest $P$ state of the beryllium atom and the bound $P$ state of the positronium molecule (with the negative parity). Both calculations yielded new, lowest-to-date, variational upper bounds, while the number of basis functions used was significantly smaller than in previous studies. It was possible to accomplish this due to the use of the analytic energy gradient in the minimization of the variational energy.},
Doi = {10.1063/1.2894866},
Keywords = {angular momentum; beryllium; Gaussian distribution; positronium; potential energy functions; variational techniques; wave functions},
Numpages = {15},
Publisher = {AIP}
}
@Article{bubin_pra_79_022501_2009,
Title = {Calculations of low-lying $^1 P$ states of the beryllium atom},
Author = {Bubin, Sergiy and Adamowicz, Ludwik},
Journal = {Phys. Rev. A},
Year = {2009},
Number = {2},
Pages = {022501},
Volume = {79},
Abstract = {High-accuracy nonrelativistic variational calculations employing explicitly correlated Gaussian basis functions have been performed to determine the energies and the expectation values of some operators for the lowest four $^1 P_1$ states of the beryllium atom. The states correspond to the electron configurations $1 s^2 2 s^1 n p^1$, where $n$=2, 3, 4, and 5. The calculations were performed for both finite and infinite mass of the Be nucleus. The basis set for each state was grown to the level of 5000 Gaussians. With that many functions we achieved a tight energy convergence. The reported values, to the best of our knowledge, are the most accurate ever obtained for the four states.},
Doi = {10.1103/PhysRevA.79.022501},
Numpages = {6},
Publisher = {American Physical Society}
}
@Article{bubin_jcp_135_214104_2011,
Title = {Accurate variational calculations of the ground $^2 P^{o} (1s^2 2s^2 2p)$ and excited $^2 S (1s^2 2s 2p^2)$ and $^2 P^{o} (1s^2 2s^2 3p)$ states of singly ionized carbon atom},
Author = {Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2011},
Pages = {214104},
Volume = {135},
Abstract = {In this article we report accurate nonrelativistic variational calculations of the ground and two excited states of C$^+$ ion. We employ extended and well optimized basis sets of all-electron explicitly correlated Gaussians to represent the wave functions of the states. The optimization of the basis functions is performed with a procedure employing the analytic gradient of the energy with respect to the nonlinear parameters of the Gaussians. The calculations explicitly include the effects due to the finite nuclear mass. The calculated transition energies between the three states are compared to the experimentally derived values. Finally, we present expectation values of some small positive and negative powers of the interparticle distances and contact densities.},
Doi = {10.1063/1.3664900},
}
@Article{bubin_pra_83_022505_2011,
Title = {Correlated-Gaussian calculations of the ground and low-lying excited states of the boron atom},
Author = {Bubin, Sergiy and Adamowicz, Ludwik },
Journal = {Phys. Rev. A},
Year = {2011},
Number = {2},
Pages = {022505},
Volume = {83},
Abstract = {Benchmark variational calculations of the four lowest $^2 P$ and $^2 S$ states of the boron atom (including the ground state) have been performed. The wave functions of the states have been expanded in terms of all-particle explicitly correlated Gaussian basis functions and the finite mass of the nucleus has been explicitly accounted for. Variational upper bounds for the nonrelativistic finite- and infinite-nuclear-mass energies of all considered states have been obtained with the relative convergence of the order of 10$^{-7}$--10$^{-8}$. Expectation values of the powers of the inter-particle distances and Dirac $\delta$ functions depending on those distances have also been computed. These calculations provide reference values that can be used to test other high-level quantum chemistry methods.},
Doi = {10.1103/PhysRevA.83.022505},
Numpages = {5},
Publisher = {American Physical Society}
}
@Article{bubin_jcp_136_134305_2012,
Title = {Explicitly correlated Gaussian calculations of the $^2 P^o$ Rydberg spectrum of the lithium atom},
Author = {Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2012},
Number = {13},
Pages = {134305},
Volume = {136},
Abstract = {Accurate quantum-mechanical nonrelativistic variational calculations are performed for the nine lowest members of the $^2 P^o$ Rydberg series ($1s^2 np^1$, $n$=2,...,10) of the lithium atom. The effect of the finite nuclear mass is included in the calculations allowing for determining the isotopic shifts of the energy levels. The wave functions of the states are expanded in terms of all-electron explicitly correlated Gaussian functions. The exponential parameters of the Gaussians are variationally optimized with the aid of the analytical energy gradient determined with respect to those parameters. The calculated state energies are compared with the available experimental data.},
Doi = {10.1063/1.3698584},
Keywords = {Gaussian processes; isotope shifts; lithium; Rydberg states; variational techniques; wave functions},
Numpages = {6},
Publisher = {AIP}
}
@Article{bubin_jcp_137_104315_2012,
Title = {Assessment of the accuracy the experimental energies of the $^1 P^o$ $1s^2 2s 6p$ and $1s^2 2s7p$ states of $^9$Be based on variational calculations with explicitly correlated Gaussians},
Author = {Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2012},
Number = {10},
Pages = {104315},
Volume = {137},
Abstract = {Benchmark variational calculations are performed for the six lowest states of the $^1 P^o$ $1s^2 2s np$ state series of the $^9$Be atom. The wave functions of the states are expanded in terms of all-particle, explicitly correlated Gaussian basis functions and the effect of the finite nuclear mass is directly included in the calculations. The exponential parameters of the Gaussians are variationally optimized using the analytical energy gradient determined with respect to those parameters. Besides providing reference non-relativistic energies for the considered states, the calculations also allow to assess the accuracy of the experimental energies of the $^1 P^o$ $1s^2 2s 6p$ and $1s^2 2s 7p$ states and suggest their refinement.},
Doi = {10.1063/1.4742764},
Eid = {104315},
Keywords = {beryllium; ground states; nuclear mass; variational techniques; wave functions},
Numpages = {6},
Publisher = {AIP}
}
@Article{bubin_pra_87_042510_2013,
Title = {Prediction of $^2 S$ Rydberg energy levels of $^6$Li and $^7$Li based on quantum-mechanical calculations performed with explicitly correlated Gaussian functions},
Author = {Bubin, Sergiy and Adamowicz, Ludwik},
Journal = {Phys. Rev. A},
Year = {2013},
Pages = {042510},
Volume = {87},
Abstract = {Accurate variational quantum-mechanical calculations are performed for the nine lowest $^2 S$ ($1s^2 ns$), $n=2,..,10$ states of the lithium atom. The effect of the finite nuclear mass is explicitly included in the calculations allowing for the determination of the isotopic shifts of the energy levels. The wave functions of the states are expanded in terms of all-electron explicitly correlated Gaussian functions and their exponential parameters are variationally optimized with the aid of the analytical energy gradient determined with respect to those parameters. The experimental results for the lower states ($n=3,...,6$) and the calculated results for the higher states ($n=7,...,10$) fitted with quantum-defect-like formulas are used to predict the energies of $^2 S$ $1s^2 ns$ states for $^7$Li and $^6$Li with $n$ up to 30.},
Doi = {10.1103/PhysRevA.87.042510},
Issue = {4},
Numpages = {7},
Publisher = {American Physical Society}
}
@Article{bubin_jcp_140_024301_2014,
Title = {Prediction of $^1 P$ Rydberg energy levels of beryllium based on calculations with explicitly correlated Gaussians},
Author = {Bubin, Sergiy and Adamowicz, Ludwik},
Journal = {J. Chem. Phys.},
Year = {2014},
Number = {2},
Pages = {024301},
Volume = {140},
Abstract = {Benchmark variational calculations are performed for the seven lowest $1s^2 2s np (^1 P)$, $n = 2...8$, states of the beryllium atom. The calculations explicitly include the effect of finite mass of $^9$Be nucleus and account perturbatively for the mass-velocity, Darwin, and spin-spin relativistic corrections. The wave functions of the states are expanded in terms of all-electron explicitly correlated Gaussian functions. Basis sets of up to 12500 optimized Gaussians are used. The maximum discrepancy between the calculated nonrelativistic and experimental energies of $1s^2 2s np (^1 P) \rightarrow 1s^2 2s^2 (^1 S)$ transition is about 12 cm$^{-1}$. The inclusion of the relativistic corrections reduces the discrepancy to bellow 0.8 cm$^{-1}$.},
Doi = {10.1063/1.4858275},
}
@Article{bubin_prl_118_043001_2017,
Title = {Lowest $^{2}S$ Electronic Excitations of the Boron Atom},
Author = {Bubin, Sergiy and Adamowicz, Ludwik},
Journal = {Phys. Rev. Lett.},
Year = {2017},
Pages = {043001},
Volume = {118},
Abstract = {A theoretical ab initio approach for calculating bound states of small atoms is developed and implemented. The approach is based on finite-nuclear-mass [non-Born-Oppenheimer (non-BO)] nonrelativistic variational calculations performed with all-particle explicitly correlated Gaussian functions and includes the leading relativistic and quantum electrodynamics energy corrections determined using the non-BO wave functions. The approach is applied to determine the total and transition energies for the lowest four $^2 S$ electronic excitations of the boron atom. The transition energies agree with the available experimental values within 0.2-0.3 cm$^{-1}$. Previously, such accuracy was achieved for three- and four-electron systems.},
Doi = {10.1103/PhysRevLett.118.043001},
Issue = {4},
Numpages = {5},
Publisher = {American Physical Society}
}
@Article{bubin_jcp_123_134310_2005,
Title = {An accurate non-Born--Oppenheimer calculation of the first purely vibrational transition in LiH molecule},
Author = {Sergiy Bubin and Ludwik Adamowicz and Marcin Molski},
Journal = {J. Chem. Phys.},
Year = {2005},
Number = {13},
Pages = {134310},
Volume = {123},
Abstract = {In this work we study the ground and the first vibrationally excited states of LiH molecule. We performed an extensive nonrelativistic variational calculations of the two states without using the Born--Oppenheimer approximation. The results are analyzed and compared with the data extracted from recent experiments. The 0$\leftarrow$1 transition energy obtained in the calculations converged to a value which is less than a wave number above the transition energy estimated from the available experimental data concerning the LiH rovibrational transitions. We discuss the remaining discrepancy and the procedure used to determine the ``experimental'' transition frequencies.},
Doi = {10.1063/1.2047487},
Keywords = {lithium compounds; vibrational states; rotational-vibrational states; variational techniques},
Numpages = {5},
Publisher = {AIP}
}
@Article{bubin_jcp_122_041102_2005,
Title = {Charge asymmetry in HD$^+$},
Author = {Sergiy Bubin and Eugeniusz Bednarz and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2005},
Number = {4},
Pages = {041102},
Volume = {122},
Abstract = {Expanding the wave functions of the ground and excited states of HD$^+$ (or $pde$) in terms of spherically symmetric explicitly correlated Gaussian functions with preexponential multipliers consisting of powers of the internuclear distance, and using the variational method, we performed very accurate nonadiabatic calculations of all bound states of this system corresponding to the zero total angular momentum quantum number (vibrational states; $v$=0-22). The total and the transition energies obtained agree with the best available calculations. For each state we computed the expectation values of the $d-p$, $d-e$, and $p-e$ interparticle distances. This is the first time these quantities were computed for HD$^+$ using rigorous nonadiabatic wave functions. While up to the $v$=20 state some asymmetry is showing in the $d-e$ and $p-e$ distances, for $v$=21 and $v$=22 we observe a complete breakdown of the Born--Oppenheimer approximation and localization of the electron almost entirely at the deuteron.},
Doi = {10.1063/1.1850905},
Keywords = {positive ions; hydrogen ions; deuterium; wave functions; ground states; excited states; variational techniques; bound states; molecular configurations},
Numpages = {4},
Publisher = {AIP}
}
@InCollection{bubin_fwqc_2004,
Title = {Quantum mechanical calculations on molecules containing positrons},
Author = {S. Bubin and M. Cafiero and L. Adamowicz},
Booktitle = {Fundamental World of Quantum Chemistry},
Publisher = {Kluwer},
Year = {2004},
Editor = {E. J. Br\"{a}ndas and E. S. Kryachko},
Pages = {521},
Volume = {3}
}
@Article{bubin_acp_131_377_2005,
Title = {Non-Born-Oppenheimer variational calculations of atoms and molecules with explicitly correlated Gaussian basis functions},
Author = {Bubin, Sergiy and Cafiero, Mauricio and Adamowicz, Ludwik},
Journal = {Adv. Chem. Phys.},
Year = {2005},
Pages = {377-475},
Volume = {131},
Doi = {10.1002/0471739464.ch6},
Keywords = {Nonrelativistic Hamiltonian;dipole approximation;energy gradients},
Publisher = {Wiley}
}
@Article{bubin_cpl_647_122_2016,
Title = {Universal all-particle explicitly-correlated Gaussians for non-Born--Oppenheimer calculations of molecular rotationless states},
Author = {Sergiy Bubin and Martin Formanek and Ludwik Adamowicz},
Journal = {Chem. Phys. Lett. },
Year = {2016},
Pages = {122-126},
Volume = {647},
Abstract = {Complex explicitly correlated one-center all-particle Gaussian functions (CECGs) are tested as basis functions for molecular non-Born--Oppenheimer (non-BO) calculations. The tests concern the complete spectrum of the 23 bound pure vibrational states of the HD$^+$ ion. The tests show that CECGs are very effective in representing both the nuclear-nuclear correlation and the vibrational oscillation of the wave functions in terms of the internuclear distance. This finding is important because it paves the way for high-accuracy non-BO calculations of bound vibrational states of molecules with more than two nuclei that have never been performed before.},
Doi = {10.1016/j.cplett.2016.01.056},
}
@Article{bubin_jcp_131_234112_2009,
Title = {Isotope shift in the electron affinity of lithium},
Author = {Sergiy Bubin and Jacek Komasa and Monika Stanke and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2009},
Number = {23},
Pages = {234112},
Volume = {131},
Abstract = {Very accurate electron affinity (EA) calculations of $^6$Li and $^7$Li (and $^\infty$Li) have been performed using explicitly correlated Gaussian functions and a variational approach that explicitly includes the nuclear motion in the calculations (i.e., the approach that does not assume the Born--Oppenheimer approximation). The leading relativistic and quantum electrodynamics corrections to the electron affinities were also calculated. The results are the most accurate theoretical values obtained for the studied systems to date. Our best estimates of the $^7$Li and $^6$Li EAs are 4984.9842(30) and 4984.9015(30) cm$^{-1}$, respectively, and of the $^7$Li/$^6$Li EA isotope shift is 0.0827 cm$^{-1}$.},
Doi = {10.1063/1.3275804},
Keywords = {electron affinity; isotope shifts; lithium; quantum electrodynamics; relativistic corrections; variational techniques},
Numpages = {5},
Publisher = {AIP}
}
@Article{bubin_jcp_132_114109_2010,
Title = {Isotope shifts of the three lowest $^1 S$ states of the B$^+$ ion calculated with a finite-nuclear-mass approach and with relativistic and quantum electrodynamics corrections},
Author = {Sergiy Bubin and Jacek Komasa and Monika Stanke and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2010},
Number = {11},
Pages = {114109},
Volume = {132},
Abstract = {We present very accurate quantum mechanical calculations of the three lowest $S$-states [$1s^2 2s^2 (^1 S_0)$, $1s^2 2p^2 (^1 S_0)$, and $1s^2 2s 3s (^1 S_0)$] of the two stable isotopes of the boron ion, $^{10}$B$^+$ and $^{11}$B$^+$. At the nonrelativistic level the calculations have been performed with the Hamiltonian that explicitly includes the finite mass of the nucleus as it was obtained by a rigorous separation of the center-of-mass motion from the laboratory frame Hamiltonian. The spatial part of the nonrelativistic wave function for each state was expanded in terms of 10000 all-electron explicitly correlated Gaussian functions. The nonlinear parameters of the Gaussians were variationally optimized using a procedure involving the analytical energy gradient determined with respect to the nonlinear parameters. The nonrelativistic wave functions of the three states were subsequently used to calculate the leading $\alpha^2$ relativistic corrections ($\alpha$ is the fine structure constant; $\alpha = 1/c$, where $c$ is the speed of light) and the $\alpha^3$ quantum electrodynamics (QED) correction. We also estimated the $\alpha^4$ QED correction by calculating its dominant component. A comparison of the experimental transition frequencies with the frequencies obtained based on the energies calculated in this work shows an excellent agreement. The discrepancy is smaller than 0.4 cm$^{-1}$.},
Doi = {10.1063/1.3358999},
Keywords = {boron; fine structure; Gaussian distribution; isotope shifts; quantum electrodynamics; relativistic corrections; wave functions},
Numpages = {6},
Publisher = {AIP}
}
@Article{bubin_pra_81_052504_2010,
Title = {Isotope shifts of the $1s^2 2s^2 (^1 S_0) \rightarrow 1 s^2 2p^2 (^1 S_0)$ transition in the doubly ionized carbon ion C$^{2+}$},
Author = {Bubin, Sergiy and Komasa, Jacek and Stanke, Monika and Adamowicz, Ludwik},
Journal = {Phys. Rev. A},
Year = {2010},
Number = {5},
Pages = {052504},
Volume = {81},
Abstract = {Highly accurate quantum mechanical calculations are performed for the $1s^2 2s^2 (^1 S_0) \rightarrow 1 s^2 2p^2 (^1 S_0)$ transition energy in the isotopomers of C$^{2+}$ ion to determine the isotope shifts. Explicitly correlated Gaussian functions and a variational approach that explicitly includes the nuclear motion are employed in the calculations. The leading relativistic and quantum electrodynamics corrections to the transition energy are also calculated using the perturbation theory with the nonrelativistic wave function as the zero-order approximation. It is determined that the $^{12}$C$^{2+}$ transitions energy, which is obtained from the calculations to be 182519.031 cm$^{-1}$ (vs the experimental value of 182519.88 cm$^{-1}$ , an excellent sub-wave-number agreement) up-shifts by 1.755 cm$^{-1}$ for $^{13}$C$^{2+}$ and by additional 1.498 cm$^{-1}$ for $^{14}$C$^{2+}$. Those shifts are sufficiently large to be measured experimentally.},
Doi = {10.1103/PhysRevA.81.052504},
Numpages = {7},
Publisher = {American Physical Society}
}
@Article{bubin_cpl_477_12_2009,
Title = {Non-adiabatic corrections to the energies of the pure vibrational states of H$_2$},
Author = {Sergiy Bubin and Filip Leonarski and Monika Stanke and Ludwik Adamowicz},
Journal = {Chem. Phys. Lett.},
Year = {2009},
Number = {1-3},
Pages = {12 - 16},
Volume = {477},
Abstract = {Nonrelativistic energies of all fifteen pure vibrational states of the H$_2$ molecule have been recalculated with much higher accuracy than before. In the calculations we employed explicitly correlated Gaussian functions and an approach where the Born-Oppenheimer (BO) approximation is not assumed. The wave function of each state was expanded in terms of 10000 Gaussians whose nonlinear parameters were optimized using a procedure involving the analytical energy gradient. The obtained non-BO energies combined with the recent BO adiabatic energies of Pachucki and Komasa [K. Pachucki, J. Komasa, J. Chem. Phys. 129 (2008) 034102] allowed us to determine new improved values of the non-adiabatic corrections for the considered states.},
Doi = {10.1016/j.cplett.2009.06.060},
}
@Article{bubin_jcp_130_124120_2009,
Title = {Charge asymmetry in pure vibrational states of the HD molecule},
Author = {Sergiy Bubin and Filip Leonarski and Monika Stanke and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2009},
Number = {12},
Pages = {124120},
Volume = {130},
Abstract = {Very accurate variational calculations of all rotationless states (also called pure vibrational states) of the HD molecule have been performed within the framework that does not assume the Born--Oppenheimer (BO) approximation. The non-BO wave functions of the states describing the internal motion of the proton, the deuteron, and the two electrons were expanded in terms of one-center explicitly correlated Gaussian functions multiplied by even powers of the internuclear distance. Up to 6000 functions were used for each state. Both linear and nonlinear parameters of the wave functions of all 18 states were optimized with a procedure that employs the analytical gradient of the energy with respect to the nonlinear parameters of the Gaussians. These wave functions were used to calculate expectation values of the interparticle distances and some other related quantities. The results allow elucidation of the charge asymmetry in HD as a function of the vibrational excitation.},
Doi = {10.1063/1.3094047},
Eid = {124120},
Keywords = {deuterium; hydrogen neutral molecules; variational techniques; vibrational states; wave functions},
Numpages = {6},
Publisher = {AIP}
}
@Article{bubin_cr_113_36_2013,
Title = {Born--Oppenheimer and Non-Born--Oppenheimer, Atomic and Molecular Calculations with Explicitly Correlated Gaussians},
Author = {Bubin, Sergiy and Pavanello, Michele and Tung, Wei-Cheng and Sharkey, Keeper L. and Adamowicz, Ludwik},
Journal = {Chem. Rev.},
Year = {2013},
Number = {0},
Pages = {36--79},
Volume = {113},
Doi = {10.1021/cr200419d},
}
@Article{bubin_prl_111_193401_2013,
Title = {Excited States of Positronic Lithium and Beryllium},
Author = {Bubin, Sergiy and Prezhdo, Oleg V.},
Journal = {Phys. Rev. Lett.},
Year = {2013},
Pages = {193401},
Volume = {111},
Abstract = {Using a variational method with an explicitly correlated Gaussian basis, we study the $e^+$-Li and $e^+$-Be complexes in the ground and lowest excited states with higher spin multiplicity. Our calculations provide rigorous theoretical confirmation that a positron can be attached to the excited states: $1s 2s 2p \, ^4 P^o$ and $1s^2 2s 2p \, ^3 P^o$ for $e^+$-Li and $e^+$-Be, respectively. The result is particularly notable for the $e^+$-Be complex, as the excited $^3 P^o$ state lies below the autoionization threshold. We report accurate binding energies, annihilation rates and structural properties of these positron-atom systems. The existence of the ground and metastable excited states with bound positron opens up a new route to the presently lacking experimental verification of stability of a positron binding to any neutral atom.},
Doi = {10.1103/PhysRevLett.111.193401},
Issue = {19},
Numpages = {5},
Publisher = {American Physical Society}
}
@Article{bubin_pra_87_054501_2013,
Title = {Instability of tripositronium},
Author = {Bubin, Sergiy and Prezhdo, Oleg V. and Varga, K\'alm\'an},
Journal = {Phys. Rev. A},
Year = {2013},
Pages = {054501},
Volume = {87},
Abstract = {The stability of tripositronium, a system consisting of three electrons and three positrons, has been investigated systematically by varying the repulsion strength between like-charged particles. The possibility of the existence of a Ps$_3$ bound state that is stable against dissociation appears utterly unlikely based on the results of variational calculations employing all-particle explicitly correlated Gaussian basis functions.},
Doi = {10.1103/PhysRevA.87.054501},
Issue = {5},
Numpages = {4},
Publisher = {American Physical Society}
}
@Article{bubin_jcp_138_164308_2013,
Title = {Prediction of $^2 D$ Rydberg energy levels of $^6$Li and $^7$Li based on very accurate quantum mechanical calculations performed with explicitly correlated Gaussian functions},
Author = {Sergiy Bubin and Keeper L. Sharkey and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2013},
Number = {16},
Pages = {164308},
Volume = {138},
Abstract = {Very accurate variational nonrelativistic finite-nuclear-mass calculations employing all-electron explicitly correlated Gaussian basis functions are carried out for six Rydberg $^2 D$ states ($1s^2 nd$, $n=6, ..., 11$) of the $^7$Li and $^6$Li isotopes. The exponential parameters of the Gaussian functions are optimized using the variational method with the aid of the analytical energy gradient determined with respect to these parameters. The experimental results for the lower states ($n=3, ..., 6$) and the calculated results for the higher states ($n=7, ..., 11$) fitted with quantum-defect-like formulas are used to predict the energies of $^2 D$ $1s^2 nd$ states for $^7$Li and $^6$Li with $n$ up to 30.},
Doi = {10.1063/1.4801855},
Keywords = {isotopes; lithium; Rydberg states; variational techniques},
Numpages = {6},
Publisher = {AIP}
}
@Article{bubin_jcp_131_044128_2009,
Title = {Non-Born--Oppenheimer calculations of the BH molecule},
Author = {Sergiy Bubin and Monika Stanke and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2009},
Number = {4},
Pages = {044128},
Volume = {131},
Abstract = {Variational calculations employing explicitly correlated Gaussian basis functions have been performed for the ground state of the boron monohydride molecule (BH) and for the boron atom (B). Up to 2000 Gaussians were used for each system. The calculations did not assume the Born--Oppenheimer (BO) approximation. In the optimization of the wave function, we employed the analytical energy gradient with respect to the Gaussian exponential parameters. In addition to the total nonrelativistic energies, we computed scalar relativistic corrections (mass-velocity and Darwin). With those added to the total energies, we estimated the dissociation energy of BH. The non-BO wave functions were also used to compute some expectation values involving operators dependent on the interparticle distances.},
Doi = {10.1063/1.3195061},
Keywords = {boron compounds; dissociation energies; ground states; relativistic corrections; variational techniques; wave functions},
Numpages = {5},
Publisher = {AIP}
}
@Article{bubin_cpl_500_229_2010,
Title = {Lower vibrational transitions of the $^3$He$^4$He$^+$ ion calculated without the Born-Oppenheimer approximation and with leading relativistic corrections},
Author = {Sergiy Bubin and Monika Stanke and Ludwik Adamowicz},
Journal = {Chem. Phys. Lett.},
Year = {2010},
Number = {4-6},
Pages = {229 - 231},
Volume = {500},
Abstract = {Very accurate variational calculations of the five lowest vibrational states of the $^3$He$^4$He$^+$ ion are carried out within a framework that does not assume the Born-Oppenheimer (BO) approximation, i.e., treating the two nuclei and three electrons forming the system on an equal footing. The non-BO wave functions are expanded in terms of one-center explicitly correlated Gaussian functions multiplied by even powers of the internuclear distance. The wave functions are used to calculate the leading relativistic corrections. The approach reproduces the experimental $^3$He$^4$He$^+$ fundamental transition within 0.055 cm$^{-1}$ and similar accuracy is expected for the higher yet unmeasured vibrational transitions determined in the present calculations.},
Doi = {10.1016/j.cplett.2010.10.021},
}
@Article{bubin_jcp_134_024103_2011,
Title = {Vibrational transitions of the $^7$LiH$^+$ ion calculated without the Born--Oppenheimer approximation and with leading relativistic corrections},
Author = {Sergiy Bubin and Monika Stanke and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2011},
Number = {2},
Pages = {024103},
Volume = {134},
Abstract = {We recently presented very accurate calculations of the fundamental vibrational frequency of the $^7$LiH$^+$ and $^3$He$^4$He$^+$ ions [Stanke et al. Phys. Rev. A 79, 060501(R) (2009)] performed without the Born--Oppenheimer approximation and included leading relativistic corrections. The accuracy of those calculations was estimated to be of the order of 0.06 cm$^{-1}$. In the present work we extend the calculations to the remaining pure vibrational states of $^7$LiH$^+$ and similarly accurate results are generated. They may lead to the experimental search for still unidentified lines corresponding to those transitions.},
Doi = {10.1063/1.3525679},
Keywords = {lithium compounds; positive ions; relativistic corrections; vibrational states},
Numpages = {4},
Publisher = {AIP}
}
@Article{bubin_jcp_135_74110_2011,
Title = {Accurate non-Born-Oppenheimer calculations of the complete pure vibrational spectrum of D$_2$ with including relativistic corrections},
Author = {Sergiy Bubin and Monika Stanke and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2011},
Number = {7},
Pages = {074110},
Volume = {135},
Abstract = {In this work we report very accurate variational calculations of the complete pure vibrational spectrum of the D$_2$ molecule performed within the framework where the Born-Oppenheimer (BO) approximation is not assumed. After the elimination of the center-of-mass motion, D$_2$ becomes a three-particle problem in this framework. As the considered states correspond to the zero total angular momentum, their wave functions are expanded in terms of all-particle, one-center, spherically symmetric explicitly correlated Gaussian functions multiplied by even non-negative powers of the internuclear distance. The nonrelativistic energies of the states obtained in the non-BO calculations are corrected for the relativistic effects of the order of $\alpha^2$ (where $\alpha = 1/c$ is the fine structure constant) calculated as expectation values of the operators representing these effects.},
Doi = {10.1063/1.3625955},
Keywords = {angular momentum; deuterium; fine structure; relativistic corrections; variational techniques; vibrational states; wave functions; },
Publisher = {AIP}
}
@Article{bubin_pra_83_042520_2011,
Title = {Complete pure vibrational spectrum of HD calculated without the Born-Oppenheimer approximation and including relativistic corrections},
Author = {Bubin, Sergiy and Stanke, Monika and Adamowicz, Ludwik},
Journal = {Phys. Rev. A},
Year = {2011},
Number = {4},
Pages = {042520},
Volume = {83},
Abstract = {All 18 bound pure vibrational levels of the HD molecule have been calculated within the framework that does not assume the Born-Oppenheimer (BO) approximation. The nonrelativistic energies of the states have been corrected for the relativistic effects of the order of $\alpha^2$ (where $\alpha$ is the fine structure constant), calculated using the perturbation theory with the nonrelativistic non-BO wave functions being the zero-order approximation. The calculations were performed by expanding the non-BO wave functions in terms of one-center explicitly correlated Gaussian functions multiplied by even powers of the internuclear distance and by performing extensive optimization of the Gaussian nonlinear parameters. Up to 10000 basis functions were used for each state.},
Doi = {10.1103/PhysRevA.83.042520},
Numpages = {5},
Publisher = {American Physical Society}
}
@Article{bubin_jcp_140_154303_2014,
Title = {Accurate non-Born-Oppenheimer calculations of the complete pure vibrational spectrum of ditritium using all-particle explicitly correlated Gaussian functions},
Author = {Bubin, Sergiy and Stanke, Monika and Adamowicz, Ludwik},
Journal = {J. Chem. Phys.},
Year = {2014},
Number = {15},
Pages = {154303},
Volume = {140},
Abstract = {Very accurate variational calculations of the complete pure vibrational spectrum of the ditritium (T$_2$) molecule are performed within the framework where the Born-Oppenheimer approximation is not assumed. After separating out the center-of-mass motion from the total laboratory-frame Hamiltonian, T$_2$ becomes a three-particle problem. States corresponding to the zero total angular momentum, which are pure vibrational states, are spherically symmetric in this framework. The wave functions of these states are expanded in terms of all-particle, one-center, spherically symmetric explicitly correlated Gaussian functions multiplied by even non-negative powers of the internuclear distance. In the calculations the total energies, the dissociation energies, and expectation values of some operators dependent on interparticle distances are determined.},
Doi = {10.1063/1.4870935},
}
@Article{bubin_pra_75_062504_2007,
Title = {Relativistic corrections to the ground-state energy of the positronium molecule},
Author = {Bubin, Sergiy and Stanke, Monika and K\c{e}dziera, Dariusz and Adamowicz, Ludwik },
Journal = {Phys. Rev. A},
Year = {2007},
Number = {6},
Pages = {062504},
Volume = {75},
Abstract = {The leading-order relativistic corrections to the ground-state energy of the positronium molecule (Ps$_2$) have been computed within the framework of perturbation theory. As the zero-order wave function we used a highly accurate nonrelativistic variational expansion in terms of 6000 explicitly correlated Gaussians that yielded the lowest variational upper bound for this system to date. We also report some expectation values representing the properties of Ps$_2$.},
Doi = {10.1103/PhysRevA.75.062504},
Numpages = {9},
Publisher = {American Physical Society}
}
@Article{bubin_pra_76_022512_2007,
Title = {Improved calculations of the lowest vibrational transitions in HeH$^{+}$ },
Author = {Bubin, Sergiy and Stanke, Monika and K\c{e}dziera, Dariusz and Adamowicz, Ludwik },
Journal = {Phys. Rev. A},
Year = {2007},
Number = {2},
Pages = {022512},
Volume = {76},
Abstract = {More accurate variational calculations of the lowest three pure vibrational states ($v$=0,1,2) of the $^4$HeH$^+$ molecular ion have been carried out without assuming the Born-Oppenheimer approximation. In the calculations we included the complete set of $\alpha^2$ relativistic corrections, i.e., mass-velocity, Darwin, spin-spin, and orbit-orbit. This allowed us to improve the agreement between the theory and the experiment for the vibrational frequencies of the 1$\leftarrow$0 and 2$\leftarrow$1 transitions as compared to our previous calculations [Stanke et al., Phys. Rev. Lett. 96, 233002 (2006)].},
Doi = {10.1103/PhysRevA.76.022512},
Numpages = {7},
Publisher = {American Physical Society}
}
@Article{bubin_cpl_494_21_2010,
Title = {Accurate non-Born-Oppenheimer calculations of the lowest vibrational energies of D$_2$ and T$_2$ with including relativistic corrections},
Author = {Sergiy Bubin and Monika Stanke and Marcin Molski and Ludwik Adamowicz},
Journal = {Chem. Phys. Lett.},
Year = {2010},
Number = {1-3},
Pages = {21 - 25},
Volume = {494},
Abstract = {In this work we report very accurate variational calculations of the two lowest vibrational states of the D$_2$ and T$_2$ molecules within the framework that does not assume the Born-Oppenheimer approximation. The non-relativistic energies of the states obtained in the non-BO calculations are corrected for the relativistic effects of the order of $\alpha^2$ calculated as expectation values of the operators representing these effects. The $v=0 \rightarrow 1$ transition energy of D$_2$ obtained in the calculations is compared with the transition frequency obtained from the experimental spectra. The comparison shows the need to include corrections higher than second-order in $\alpha$ to further improve the agreement between the theory and the experiment.},
Doi = {10.1016/j.cplett.2010.05.081},
}
@Article{bubin_pra_84_012509_2011,
Title = {Ground-state energy and relativistic corrections for positronium hydride},
Author = {Bubin, Sergiy and Varga, K\'{a}lm\'{a}n},
Journal = {Phys. Rev. A},
Year = {2011},
Number = {1},
Pages = {012509},
Volume = {84},
Abstract = {Variational calculations of the ground state of positronium hydride (HPs) are reported, including various expectation values, electron-positron annihilation rates, and leading relativistic corrections to the total and dissociation energies. The calculations have been performed using a basis set of 4000 thoroughly optimized explicitly correlated Gaussian basis functions. The relative accuracy of the variational energy upper bound is estimated to be of the order of $2 \times 10^{-10}$, which is a significant improvement over previous nonrelativistic results.},
Doi = {10.1103/PhysRevA.84.012509},
Numpages = {8},
Publisher = {American Physical Society}
}
@Article{cafiero_pccp_5_1491_2003,
Title = {Non-Born-Oppenheimer calculations of atoms and molecules},
Author = {Cafiero, M and Bubin, S and Adamowicz, L},
Journal = {Phys. Chem. Chem. Phys.},
Year = {2003},
Number = {8},
Pages = {1491-1501},
Volume = {5},
Abstract = {We review a recent development in high-accuracy non-Born--Oppenheimer calculations of atomic and molecular systems in a basis of explicitly correlated Gaussian functions. Much of the recent progress in this area is due to the derivation and implementation of analytical gradients of the energy functional with respect to variational linear and non-linear parameters of the basis functions. This method has been used to obtain atomic and molecular ground and excited state energies and the corresponding wave functions with accuracy that exceeds previous calculations. Further, we have performed the first calculations of non-linear electrical properties of molecules without the Born--Oppenheimer approximation for systems with more than one electron. The results for the dipole moments of such systems as HD and LiH agree very well with experiment. After reviewing our non-Born--Oppenheimer results we will discuss ways this method can be extended to deal with larger molecular systems with and without an external perturbation.},
Doi = {10.1039/b211193d},
}
@Article{kedziera_jcp_125_014318_2006,
Title = {Darwin and mass-velocity relativistic corrections in the non-Born-Oppenheimer calculations of pure vibrational states of H$_2$},
Author = {Dariusz K\c{e}dziera and Monika Stanke and Sergiy Bubin and Maria Barysz and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2006},
Number = {1},
Pages = {014318},
Volume = {125},
Abstract = {The Darwin and mass-velocity relativistic corrections have been calculated for all pure vibrational states of the H$_2$ using the perturbation theory and very accurate variational wave functions obtained without assuming the Born-Oppenheimer (BO) approximation. Expansions in terms of explicitly correlated Gaussians with premultipliers in the form of even powers of the internuclear distance were used for the wave functions. With the inclusion of the two relativistic corrections to the non-BO energies the transition energies for the highest states agree more with the experimental results.},
Doi = {10.1063/1.2209691},
Keywords = {hydrogen neutral molecules; relativistic corrections; vibrational states; perturbation theory; variational techniques; wave functions},
Numpages = {7},
Publisher = {AIP}
}
@Article{kedziera_jcp_125_084303_2006,
Title = {Darwin and mass-velocity relativistic corrections in non-Born-Oppenheimer variational calculations},
Author = {Dariusz K\c{e}dziera and Monika Stanke and Sergiy Bubin and Maria Barysz and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2006},
Number = {8},
Pages = {084303},
Volume = {125},
Abstract = {The Pauli approach to account for the mass-velocity and Darwin relativistic corrections has been applied to the formalism for quantum mechanical molecular calculations that does not assume the Born-Oppenheimer (BO) approximation regarding separability of the electronic and nuclear motions in molecular systems. The corrections are determined using the first order perturbation theory and are derived for the non-BO wave function of a diatomic system expressed in terms of explicitly correlated Gaussian functions with premultipliers in the form of even powers of the internuclear distance. As a numerical example we used calculations of the transition energies for pure vibrational states of the HD$^+$ ion.},
Doi = {10.1063/1.2236113},
Keywords = {relativistic corrections; vibrational states; bond lengths; quantum chemistry; molecular electronic states; perturbation theory; deuterium; hydrogen ions},
Numpages = {12},
Publisher = {AIP}
}
@Article{mitroy_rmp_85_693_2013,
Title = {Theory and application of explicitly correlated Gaussians},
Author = {Mitroy, Jim and Bubin, Sergiy and Horiuchi, Wataru and Suzuki, Yasuyuki and Adamowicz, Ludwik and Cencek, Wojciech and Szalewicz, Krzysztof and Komasa, Jacek and Blume, D. and Varga, K\'alm\'an},
Journal = {Rev. Mod. Phys.},
Year = {2013},
Pages = {693--749},
Volume = {85},
Abstract = {The variational method complemented with the use of explicitly correlated Gaussian basis functions is one of the most powerful approaches currently used for calculating the properties of few-body systems. Despite its conceptual simplicity, the method offers great flexibility, high accuracy, and can be used to study diverse quantum systems, ranging from small atoms and molecules to light nuclei, hadrons, quantum dots, and Efimov systems. The basic theoretical foundations are discussed, recent advances in the applications of explicitly correlated Gaussians in physics and chemistry are reviewed, and the strengths and weaknesses of the explicitly correlated Gaussians approach are compared with other few-body techniques.},
Doi = {10.1103/RevModPhys.85.693},
Issue = {2},
Publisher = {American Physical Society}
}
@Article{pavanello_jcp_123_104306_2005,
Title = {Non-Born--Oppenheimer calculations of the pure vibrational spectrum of HeH$^+$},
Author = {Michele Pavanello and Sergiy Bubin and Marcin Molski and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2005},
Number = {10},
Pages = {104306},
Volume = {123},
Abstract = {Very accurate calculations of the pure vibrational spectrum of the HeH$^+$ ion are reported. The method used does not assume the Born--Oppenheimer approximation, and the motion of both the electrons and the nuclei are treated on equal footing. In such an approach the vibrational motion cannot be decoupled from the motion of electrons, and thus the pure vibrational states are calculated as the states of the system with zero total angular momentum. The wave functions of the states are expanded in terms of explicitly correlated Gaussian basis functions multipled by even powers of the internuclear distance. The calculations yielded twelve bound states and corresponding eleven transition energies. Those are compared with the pure vibrational transition energies extracted from the experimental rovibrational spectrum.},
Doi = {10.1063/1.2012332},
Keywords = {vibrational states; helium compounds; wave functions; positive ions},
Numpages = {7},
Publisher = {AIP}
}
@Article{pavanello_ijqc_108_2291_2008,
Title = {Accurate Born--Oppenheimer calculations of the low-lying $c ^3 \Sigma_g^+$ and $a ^3 \Sigma_u^+$ excited states of helium dimer},
Author = {Michele Pavanello and Mauricio Cafiero and Sergiy Bubin and Ludwik Adamowicz},
Journal = {Int. J. Quantum Chem.},
Year = {2008},
Pages = {2291-2298},
Volume = {108},
Abstract = {Accurate variational Born--Oppenheimer calculations of the two lowest excited $\Sigma$ triplet states of the helium dimer at their respective equilibrium geometries are reported. The wave functions of the states are expanded in terms of explicitly correlated Gaussian functions with shifted centers. The obtained energies are the best variational estimates ever calculated for these states. One-electron densities are also presented and discussed. The results are compared with the experimental values and previous calculations.},
Doi = {10.1002/qua.21757},
Keywords = {helium dimmer; Born--Oppenheimer; nonrelativistic; excited states; correlated Gaussians},
Publisher = {Wiley}
}
@Article{sharkey_jcp_132_184106_2010,
Title = {Analytical energy gradient in variational calculations of the two lowest $^3 P$ states of the carbon atom with explicitly correlated Gaussian basis functions},
Author = {Keeper L. Sharkey and Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2010},
Number = {18},
Pages = {184106},
Volume = {132},
Abstract = {Variational calculations of ground and excited bound states on atomic and molecular systems performed with basis functions that explicitly depend on the interparticle distances can generate very accurate results provided that the basis function parameters are thoroughly optimized by the minimization of the energy. In this work we have derived the algorithm for the gradient of the energy determined with respect to the nonlinear exponential parameters of explicitly correlated Gaussian functions used in calculating $n$-electron atomic systems with two $p$-electrons and ($n-2$) $s$-electrons. The atomic Hamiltonian we used was obtained by rigorously separating out the kinetic energy of the center of mass motion from the laboratory-frame Hamiltonian and explicitly depends on the finite mass of the nucleus. The advantage of having the gradient available in the variational minimization of the energy is demonstrated in the calculations of the ground and the first excited $^3 P$ state of the carbon atom. For the former the lowest energy upper bound ever obtained is reported.},
Doi = {10.1063/1.3419931},
Keywords = {bound states; carbon; excited states; Gaussian distribution; ground states; variational techniques},
Numpages = {11},
Publisher = {AIP}
}
@Article{sharkey_jcp_134_044120_2011,
Title = {An algorithm for calculating atomic $D$ states with explicitly correlated Gaussian functions},
Author = {Keeper L. Sharkey and Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2011},
Number = {4},
Pages = {044120},
Volume = {134},
Abstract = {An algorithm for the variational calculation of atomic $D$ states employing $n$-electron explicitly correlated Gaussians is developed and implemented. The algorithm includes formulas for the first derivatives of the Hamiltonian and overlap matrix elements determined with respect to the Gaussian nonlinear exponential parameters. The derivatives are used to form the energy gradient which is employed in the variational energy minimization. The algorithm is tested in the calculations of the two lowest $D$ states of the lithium and beryllium atoms. For the lowest $D$ state of Li the present result is lower than the best previously reported result.},
Doi = {10.1063/1.3523348},
Keywords = {beryllium; Gaussian processes; lithium; variational techniques},
Numpages = {9},
Publisher = {AIP}
}
@Article{sharkey_jcp_134_194114_2011,
Title = {Refinement of the experimental energy levels of higher $^2 D$ Rydberg states of the lithium atom with very accurate quantum mechanical calculations},
Author = {Keeper L. Sharkey and Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2011},
Number = {19},
Pages = {194114},
Volume = {134},
Abstract = {Very accurate variational non-relativistic calculations are performed for four higher Rydberg $^2 D$ states ($1s^2 nd^1$, $n$ = 8, ..., 11) of the lithium atom ($^7$Li). The wave functions of the states are expanded in terms of all-electron explicitly correlated Gaussian functions and finite nuclear mass is used. The exponential parameters of the Gaussians are optimized using the variational method with the aid of the analytical energy gradient determined with respect to those parameters. The results of the calculations allow for refining the experimental energy levels determined with respect to the $^2 S$ $1s^2 2s^1$ ground state.},
Doi = {10.1063/1.3591836},
Keywords = {atomic mass; atomic structure; ground states; lithium; Rydberg states; variational techniques},
Numpages = {4},
Publisher = {AIP}
}
@Article{sharkey_pra_83_012506_2011,
Title = {Lower Rydberg $^{2}D$ states of the lithium atom: Finite-nuclear-mass calculations with explicitly correlated Gaussian functions},
Author = {Sharkey, Keeper L. and Bubin, Sergiy and Adamowicz, Ludwik},
Journal = {Phys. Rev. A},
Year = {2011},
Number = {1},
Pages = {012506},
Volume = {83},
Abstract = {Very accurate variational nonrelativistic calculations are performed for the five lowest Rydberg $^2 D$ states ($1s^2 nd^1$, $n$=3,...,7) of the lithium atom ($^7$Li). The finite-nuclear-mass approach is employed and the wave functions of the states are expanded in terms of all-electron explicitly correlated Gaussian function. Four thousand Gaussians are used for each state. The calculated relative energies of the states determined with respect to the $^2 S$ $1s^2 2s^1$ ground state are systematically lower than the experimental values by about 2.5 cm$^{-1}$. As this value is about the same as the difference between the experimental relative energy between $^7$Li$^+$ and $^7$Li in their ground-state energy and the corresponding calculated nonrelativistic relative energy, we attribute it to the relativistic effects not included in the present calculations.},
Doi = {10.1103/PhysRevA.83.012506},
Numpages = {5},
Publisher = {American Physical Society}
}
@Article{sharkey_pra_84_044503_2011,
Title = {$^{1}D$ states of the beryllium atom: Quantum mechanical nonrelativistic calculations employing explicitly correlated Gaussian functions},
Author = {Sharkey, Keeper L. and Bubin, Sergiy and Adamowicz, Ludwik},
Journal = {Phys. Rev. A},
Year = {2011},
Pages = {044503},
Volume = {84},
Abstract = {Very accurate finite-nuclear-mass variational nonrelativistic calculations are performed for the lowest five $^1 D$ states ($1s^2 2p^2$, $1s^2 2s^1 3d^1$, $1s^2 2s^1 4d^1$, $1s^2 2s^1 5d^1$, and $1s^2 2s^1 6d^1$) of the beryllium atom ($^9$Be). The wave functions of the states are expanded in terms of all-electron explicitly correlated Gaussian functions. The exponential parameters of the Gaussians are optimized using the variational method with the aid of the analytical energy gradient determined with respect to those parameters. The calculations exemplify the level of accuracy that is now possible with Gaussians in describing bound states of a four-electron system where some of the electrons are excited into higher angular states.},
Doi = {10.1103/PhysRevA.84.044503},
Issue = {4},
Numpages = {3},
Publisher = {American Physical Society}
}
@Article{sharkey_cpl_616_254_2014,
Title = {Singlet--triplet energy splitting between $^1 D$ and $^3 D$ ($1s^2 2s nd$) $n$=3,4,5, and 6, Rydberg states of the beryllium atom ($^9 Be$) calculated with all-electron explicitly correlated Gaussian functions},
Author = {Keeper L. Sharkey and Sergiy Bubin and Ludwik Adamowicz},
Journal = {Chem. Phys. Lett.},
Year = {2014},
Pages = {254-258},
Volume = {616-617},
Abstract = {Abstract Accurate variational nonrelativistic quantum-mechanical calculations are performed for the five lowest $^1 D$ and four lowest $^3 D$ states of the $^9 Be$ isotope of the beryllium atom. All-electron explicitly correlated Gaussian (ECG) functions are used in the calculations and their nonlinear parameters are optimized with the aid of the analytical energy gradient determined with respect to these parameters. The effect of the finite nuclear mass is directly included in the Hamiltonian used in the calculations. The singlet--triplet energy gaps between the corresponding $^1 D$ and $^3 D$ states, are reported.},
Doi = {10.1016/j.cplett.2014.10.012},
}
@Article{sharkey_pra_80_062510_2009,
Title = {Algorithm for quantum-mechanical finite-nuclear-mass variational calculations of atoms with two $p$ electrons using all-electron explicitly correlated Gaussian basis functions},
Author = {Sharkey, Keeper L. and Pavanello, Michele and Bubin, Sergiy and Adamowicz, Ludwik },
Journal = {Phys. Rev. A},
Year = {2009},
Number = {6},
Pages = {062510},
Volume = {80},
Abstract = {A new algorithm for calculating the Hamiltonian matrix elements with all-electron explicitly correlated Gaussian functions for quantum-mechanical calculations of atoms with two $p$ electrons or a single $d$ electron have been derived and implemented. The Hamiltonian used in the approach was obtained by rigorously separating the center-of-mass motion and it explicitly depends on the finite mass of the nucleus. The approach was employed to perform test calculations on the isotopes of the carbon atom in their ground electronic states and to determine the finite-nuclear-mass corrections for these states.},
Doi = {10.1103/PhysRevA.80.062510},
Numpages = {10},
Publisher = {American Physical Society}
}
@Article{stanke_pra_79_060501_2009,
Title = {Fundamental vibrational transitions of the $^3$He$^4$He$^+$ and $^7$LiH$^+$ ions calculated without assuming the Born-Oppenheimer approximation and with including leading relativistic corrections},
Author = {Stanke, Monica and Bubin, Sergiy and Adamowicz, Ludwik},
Journal = {Phys. Rev. A},
Year = {2009},
Number = {6},
Pages = {060501},
Volume = {79},
Abstract = {Very accurate variational calculations of the fundamental pure vibrational transitions of the $^3He$^4He$^+$ and $^7$LiH$^+$ ions are performed within the framework that does not assume the Born-Oppenheimer (BO) approximation. The non-BO wave functions expanded in terms of one-center explicitly correlated Gaussian functions multiplied by even powers of the internuclear distance are used to calculate the leading relativistic corrections. Up to 10000 Gaussian functions are used for each state. It is shown that the experimental $^3He$^4He$^+$ fundamental transitions is reproduced within 0.06 cm$^{-1}$ by the calculations. A similar precision is expected for the calculated, but still unmeasured, fundamental transition of $^7$LiH$^+$. Thus, three-electron diatomic systems are calculated with a similar accuracy as two-electron systems.},
Doi = {10.1103/PhysRevA.79.060501},
Numpages = {4},
Publisher = {American Physical Society}
}
@Article{stanke_pra_79_032507_2009,
Title = {Non-Born-Oppenheimer calculations of the lowest vibrational energy of HD including relativistic corrections},
Author = {Stanke, Monika and Bubin, Sergiy and Molski, Marcin and Adamowicz, Ludwik},
Journal = {Phys. Rev. A},
Year = {2009},
Number = {3},
Pages = {032507},
Volume = {79},
Abstract = {In this work we report variational calculations of the two lowest vibrational states of the HD molecule within the framework that does not assume the Born-Oppenheimer (BO) approximation. The nonrelativistic energies of the states were corrected for the relativistic effects of the order of $\alpha^2$ (where $\alpha=1/c$), calculated as expectation values of the operators representing these effects with the nonrelativistic non-BO wave functions. The non-BO wave functions were expanded in terms of the one-center explicitly correlated Gaussian functions multiplied by even powers of the internuclear distance. The $v=0 \rightarrow 1$ transition energy obtained in the calculations is compared with the previous calculations, as well as with the transition frequency obtained from the experimental spectra. The comparison shows the need to include corrections higher than second order in $\alpha$ to further improve the agreement between the theory and the experiment.},
Doi = {10.1103/PhysRevA.79.032507},
Numpages = {8},
Publisher = {American Physical Society}
}
@Article{stanke_jcp_126_194312_2007,
Title = {Relativistic corrections to the non-Born-Oppenheimer energies of the lowest singlet Rydberg states of $^3$He and $^4$He},
Author = {Monika Stanke and Dariusz K\c{e}dziera and Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2007},
Number = {19},
Pages = {194312},
Volume = {126},
Abstract = {In this work the authors present an approach to calculate the leading-order relativistic corrections for ground and excited states of helium isotopomers. In the calculations they used variational wave functions expanded in terms of explicitly correlated Gaussians obtained without assuming the Born-Oppenheimer approximation.},
Doi = {10.1063/1.2735305},
Keywords = {helium neutral atoms; relativistic corrections; Rydberg states; ground states; isomerism; isotope effects; wave functions},
Numpages = {6},
Publisher = {AIP}
}
@Article{stanke_jcp_127_134107_2007,
Title = {Electron affinity of $^7$Li calculated with the inclusion of nuclear motion and relativistic corrections},
Author = {Monika Stanke and Dariusz K\c{e}dziera and Sergiy Bubin and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2007},
Number = {13},
Pages = {134107},
Volume = {127},
Abstract = {Explicitly correlated Gaussian functions have been used to perform very accurate variational calculations for the ground states of $^7$Li and $^7$Li$^-$. The nuclear motion has been explicitly included in the calculations (i.e., they have been done without assuming the Born-Oppenheimer (BO) approximation). An approach based on the analytical energy gradient calculated with respect to the Gaussian exponential parameters was employed. This led to a noticeable improvement of the previously determined variational upper bound to the nonrelativistic energy of Li$^-$. The Li energy obtained in the calculations matches those of the most accurate results obtained with Hylleraas functions. The finite-mass (non-BO) wave functions were used to calculate the $\alpha^2$ relativistic corrections ($\alpha = 1/c$). With those corrections and the $\alpha^3$ and $\alpha^4$ corrections taken from Pachucki and Komasa [J. Chem. Phys. 125, 204304 (2006)] , the electron affinity (EA) of $^7$Li was determined. It agrees very well with the most recent experimental EA.},
Doi = {10.1063/1.2755767},
Eid = {134107},
Keywords = {electron affinity; ground states; lithium; negative ions; relativistic corrections; variational techniques; wave functions},
Numpages = {5},
Publisher = {AIP}
}
@Article{stanke_pra_75_052510_2007,
Title = {Ionization potential of $^{9}$Be calculated including nuclear motion and relativistic corrections},
Author = {Stanke, Monika and K\c{e}dziera, Dariusz and Bubin, Sergiy and Adamowicz, Ludwik},
Journal = {Phys. Rev. A},
Year = {2007},
Number = {5},
Pages = {052510},
Volume = {75},
Abstract = {Variational calculations employing explicitly correlated Gaussian functions have been performed for the ground states of $^9$Be and $^9$Be$^+$ including the nuclear motion [i.e., without assuming the Born-Oppenheimer (BO) approximation]. An approach based on the analytical energy gradient calculated with respect to the Gaussian exponential parameters was employed, leading to energies of the two systems noticeably improved over those found in the recent paper of Pachucki and Komasa [Physical Review A 73, 052502 (2006)]. The non-BO wave functions were used to calculate the $\alpha^2$ relativistic corrections ($\alpha=e^2/ \hbar c$). With those corrections and the $\alpha^3$ and $\alpha^4$ corrections taken from Pachucki and Komasa, a new value of the ionization potential (IP) of 9Be was determined. It agrees very well with the most recent experimental IP.},
Doi = {10.1103/PhysRevA.75.052510},
Numpages = {5},
Publisher = {American Physical Society}
}
@Article{stanke_prl_99_043001_2007,
Title = {Lowest Excitation Energy of $^{9}$Be},
Author = {Stanke, Monika and K\c{e}dziera, Dariusz and Bubin, Sergiy and Adamowicz, Ludwik},
Journal = {Phys. Rev. Lett.},
Year = {2007},
Number = {4},
Pages = {043001},
Volume = {99},
Abstract = {Variational calculations employing explicitly correlated Gaussian functions and explicitly including the nuclear motion [i.e., without assuming the Born-Oppenheimer (BO) approximation] have been performed to determine the lowest singlet transition energy in the 9Be atom. The non-BO wave functions were used to calculate the $\alpha^2$ relativistic corrections ($\alpha$=1/137.035999679). With those corrections and with the $\alpha^3$ and $\alpha^4$ QED corrections determined previously by others, we obtained 54677.35 cm$^{-1}$ for the $3 ^1 S \rightarrow 2 ^1 S$ transition energy. This result falls within the error bracket for the experimental transition of 54677.26(10) cm${-}1$. This is the first time an electronic transition of Be has been calculated from first principles with the experimental accuracy.},
Doi = {10.1103/PhysRevLett.99.043001},
Numpages = {4},
Publisher = {American Physical Society}
}
@Article{stanke_pra_77_022506_2008,
Title = {Complete $\alpha^2$ relativistic corrections to the pure vibrational non-Born-Oppenheimer energies of HeH$^+$ },
Author = {Stanke, Monika and K\c{e}dziera, Dariusz and Bubin, Sergiy and Adamowicz, Ludwik },
Journal = {Phys. Rev. A},
Year = {2008},
Number = {2},
Pages = {022506},
Volume = {77},
Abstract = {We report the implementation of the complete set of the lowest-order relativistic corrections of the order of $\alpha^2$ (where $\alpha$ is the fine structure constant) for calculating vibrational states of diatomic molecular systems within the framework that does not assume the Born-Oppenheimer approximation. To test the accuracy of the approach we have performed calculations for all rotationless vibrational states (also called pure vibrational states or $S$ states) of the HeH$^+$ ion in the ground electronic state. For the lowest transitions, where very precise experimental results are available, an excellent agreement with the experimental values has been achieved.},
Doi = {10.1103/PhysRevA.77.022506},
Numpages = {5},
Publisher = {American Physical Society}
}
@Article{stanke_pra_76_052506_2007,
Title = {Lowest vibrational states of $^{4}$He$^{3}$He$^{+}$: Non-Born-Oppenheimer calculations},
Author = {Stanke, Monika and K\c{e}dziera, Dariusz and Bubin, Sergiy and Molski, Marcin and Adamowicz, Ludwik },
Journal = {Phys. Rev. A},
Year = {2007},
Number = {5},
Pages = {052506},
Volume = {76},
Abstract = {Very accurate quantum mechanical calculations of the first five vibrational states of the $^4$He$^3$He$^+$ molecular ion are reported. The calculations have been performed explicitly including the coupling of the electronic and nuclear motions [i.e., without assuming the Born-Oppenheimer (BO) approximation]. The nonrelativistic non-BO wave functions were used to calculate the $\alpha^2$ relativistic mass velocity, Darwin, and spin-spin interaction corrections. For the lowest vibrational transition, whose experimental energy is established with high precision, the calculated and the experimental results differ by only 0.16 cm$^{-1}$.},
Doi = {10.1103/PhysRevA.76.052506},
Numpages = {7},
Publisher = {American Physical Society}
}
@Article{stanke_jcp_128_114313_2008,
Title = {Orbit-orbit relativistic corrections to the pure vibrational non-Born-Oppenheimer energies of H$_2$},
Author = {Monika Stanke and Dariusz K\c{e}dziera and Sergiy Bubin and Marcin Molski and Ludwik Adamowicz},
Journal = {J. Chem. Phys.},
Year = {2008},
Number = {11},
Pages = {114313},
Volume = {128},
Abstract = {We report the derivation of the orbit-orbit relativistic correction for calculating pure vibrational states of diatomic molecular systems with $\sigma$ electrons within the framework that does not assume the Born-Oppenheimer (BO) approximation. The correction is calculated as the expectation value of the orbit-orbit interaction operator with the non-BO wave function expressed in terms of explicitly correlated Gaussian functions multiplied by even powers of the internuclear distance. With that we can now calculate the complete relativistic correction of the order of $\alpha^2$ (where $\alpha = 1/c$). The new algorithm is applied to determine the full set of the rotationless vibrational levels and the corresponding transition frequencies of the H$_2$ molecule. The results are compared with the previous calculations, as well as with the frequencies obtained from the experimental spectra. The comparison shows the need to include corrections higher than second order in $\alpha$ to further improve the agreement between the theory and the experiment.},
Doi = {10.1063/1.2834926},
Keywords = {hydrogen neutral molecules; relativistic corrections; vibrational states},
Numpages = {15},
Publisher = {AIP}
}
@Article{stanke_prl_96_233002_2006,
Title = {Convergence of Experiment and Theory on the Pure Vibrational Spectrum of HeH$^{+}$},
Author = {Stanke, Monika and K\c{e}dziera, Dariusz and Molski, Marcin and Bubin, Sergiy and Barysz, Maria and Adamowicz, Ludwik},
Journal = {Phys. Rev. Lett.},
Year = {2006},
Number = {23},
Pages = {233002},
Volume = {96},
Abstract = {Very accurate quantum mechanical calculations of the pure vibrational spectrum of the HeH$^+$ molecular ion are reported and compared with newly obtained pure vibrational transitions extracted from the available experimental data. The calculations are performed without assuming the Born-Oppenheimer approximation regarding separability of the nuclear and electronic motions and include the first order relativistic mass-velocity and Darwin corrections. For the two lowest transitions, whose experimental energies are established with the highest precision, the calculated and the experimental results show very good agreement.},
Doi = {10.1103/PhysRevLett.96.233002},
Numpages = {4},
Publisher = {American Physical Society}
}
@Article{stanke_pra_80_022514_2009,
Title = {Five lowest $^1S$ states of the Be atom calculated with a finite-nuclear-mass approach and with relativistic and QED corrections},
Author = {Stanke, Monika and Komasa, Jacek and Bubin, Sergiy and Adamowicz, Ludwik},
Journal = {Phys. Rev. A},
Year = {2009},
Number = {2},
Pages = {022514},
Volume = {80},
Abstract = {We have performed very accurate quantum mechanical calculations of the five lowest $S$ states of the beryllium atom. In the nonrelativistic part of the calculations we used the variational method and we explicitly included the nuclear motion in the Schr\"{o}dinger equation. The nonrelativistic wave functions of the five states were expanded in terms of explicitly correlated Gaussian functions. These wave functions were used to calculate the leading $\alpha^2$ relativistic correction ($\alpha$ is the fine structure constant) and the $\alpha^3$ quantum electrodynamics (QED) correction. We also estimated the $\alpha^4$ QED correction by calculating its dominant component. A comparison of the experimental transition frequencies with the frequencies obtained based on the energies calculated in this work shows an excellent agreement.},
Doi = {10.1103/PhysRevA.80.022514},
Numpages = {6},
Publisher = {American Physical Society}
}
@Article{stanke_pra_77_062509_2008,
Title = {Three lowest $S$ states of $^{9}$Be$^{+}$ calculated with including nuclear motion and relativistic and QED corrections},
Author = {Stanke, Monika and Komasa, Jacek and K\c{e}dziera, Dariusz and Bubin, Sergiy and Adamowicz, Ludwik},
Journal = {Phys. Rev. A},
Year = {2008},
Number = {6},
Pages = {062509},
Volume = {77},
Abstract = {We have performed high-accuracy quantum mechanical calculations for the three lowest $S$ states of the beryllium ion ($^9$Be$^+$). The nonrelativistic part of the calculations was done with the variational approach and explicitly included the nuclear motion (i.e., the finite-nuclear-mass approach). The nonrelativistic wave functions were expanded in terms of explicitly correlated Gaussian functions. These nonrelativistic functions were subsequently used to calculate the leading $\alpha^2$ relativistic corrections ($\alpha=1/c$) and the $\alpha^3$ and $\alpha^4$ QED (quantum electrodynamics) corrections. In the $\alpha^4$ QED correction we only accounted for its dominant component typically contributing about 80\% of the correction. With those the present results are the most accurate ever obtained for $^9$Be$^+$. They also agree with the experimentally measured transitions within less than 0.1 cm$^{-1}$.},
Doi = {10.1103/PhysRevA.77.062509},
Numpages = {7},
Publisher = {American Physical Society}
}
@Article{stanke_pra_78_052507_2008,
Title = {Accuracy limits on the description of the lowest $S$ excitation in the Li atom using explicitly correlated Gaussian basis functions},
Author = {Stanke, Monika and Komasa, Jacek and K\c{e}dziera, Dariusz and Bubin, Sergiy and Adamowicz, Ludwik},
Journal = {Phys. Rev. A},
Year = {2008},
Number = {5},
Pages = {052507},
Volume = {78},
Abstract = {We have performed very accurate quantum-mechanical calculations for the two lowest $S$ states of the lithium atom in order to determine the transition energy. In the nonrelativistic part of the calculations performed using the variational method, we explicitly included the nuclear motion. The nonrelativistic wave function was expanded in terms of explicitly correlated Gaussian functions. Next, this wave function was used to calculate the leading $\alpha^2$ relativistic correction ($\alpha$ is the fine-structure constant) and the $\alpha^3$ QED correction. We also estimated the $\alpha^4$ QED correction by calculating its dominating component. The results obtained with Gaussians are compared with the most accurate results obtained recently with the Hylleraas-type basis functions.},
Doi = {10.1103/PhysRevA.78.052507},
Numpages = {8},
Publisher = {American Physical Society}
}
@Article{bubin_pra_95_062509_2017,
author = {Bubin, Sergiy and Stanke, Monika and Adamowicz, Ludwik},
title = {Relativistic corrections for non-Born-Oppenheimer molecular wave functions expanded in terms of complex explicitly correlated Gaussian functions},
journal = {Phys. Rev. A},
year = {2017},
volume = {95},
pages = {062509},
abstract = {In our previous work S. Bubin et al., Chem. Phys. Lett. 647, 122 (2016), it was established that complex explicitly correlated one-center all-particle Gaussian functions (CECGs) provide effective basis functions for very accurate nonrelativistic molecular non-Born-Oppenheimer calculations. In this work, we advance the molecular CECGs approach further by deriving and implementing algorithms for calculating the leading relativistic corrections within this approach. The algorithms are tested in the calculations of the corrections for all 23 bound pure vibrational states of the HD$^+$ ion.},
doi = {10.1103/PhysRevA.95.062509},
issue = {6},
numpages = {8},
publisher = {American Physical Society},
}
@Article{bubin_pra_86_043407_2012,
Title = {Strong laser-pulse-driven ionization and Coulomb explosion of hydrocarbon molecules},
Author = {Bubin, Sergiy and Atkinson, Mackenzie and Varga, K\'alm\'an and Xie, Xinhua and Roither, Stefan and Kartashov, Daniil and Baltu\v{s}ka, Andrius and Kitzler, Markus},
Journal = {Phys. Rev. A},
Year = {2012},
Pages = {043407},
Volume = {86},
Abstract = {Field ionization and Coulomb explosion of small hydrocarbon molecules driven by intense laser pulses are studied in a combined theoretical and experimental framework. The spectra of ejected protons calculated by the time-dependent density functional approach are in good agreement with the experimental data. The results of the simulations give detailed insight into the correlated electron and nuclear dynamics and complement the experiment with a time-dependent physical picture. It is demonstrated that the Coulomb explosion in the studied molecular systems is a sudden, all-at-once fragmentation where the ionization step is followed by a simultaneous ejection of the charged fragments.},
Doi = {10.1103/PhysRevA.86.043407},
Issue = {4},
Numpages = {5},
Publisher = {American Physical Society}
}
@Article{bubin_jpcs_436_12084_2013,
Title = {Interaction of electromagnetic fields and atomic clusters},
Author = {S Bubin and A G Russakoff and K Varga},
Journal = {J. Phys. Conf. Ser.},
Year = {2013},
Number = {1},
Pages = {012084},
Volume = {436},
Abstract = {In the framework of real-time real-space time-dependent density functional theory complemented with Ehrenfest molecular dynamics we have studied the response of nanostructures to intense femtosecond laser pulses. Examples of applications include laser desorption of hydrogen from graphene and Coulomb explosion of hydrocarbon molecules.},
Doi = {10.1088/1742-6596/436/1/012084},
}
@Article{bubin_jpcm_22_465306_2010,
Title = {Calculation of transmission probability by solving an eigenvalue problem},
Author = {Sergiy Bubin and K\'{a}lm\'{a}n Varga},
Journal = {J. Phys. Condens. Matter},
Year = {2010},
Number = {46},
Pages = {465306},
Volume = {22},
Abstract = {The electron transmission probability in nanodevices is calculated by solving an eigenvalue problem. The eigenvalues are the transmission probabilities and the number of nonzero eigenvalues is equal to the number of open quantum transmission eigenchannels. The number of open eigenchannels is typically a few dozen at most, thus the computational cost amounts to the calculation of a few outer eigenvalues of a complex Hermitian matrix (the transmission matrix). The method is implemented on a real space grid basis providing an alternative to localized atomic orbital based quantum transport calculations. Numerical examples are presented to illustrate the efficiency of the method.},
Doi = {10.1088/0953-8984/22/46/465306},
}
@Article{bubin_apl_98_154101_2011,
Title = {First-principles time-dependent simulation of laser assisted desorption of hydrogen atoms from H--Si(111) surface},
Author = {Sergiy Bubin and K\'{a}lm\'{a}n Varga},
Journal = {Appl. Phys. Lett.},
Year = {2011},
Number = {15},
Pages = {154101},
Volume = {98},
Abstract = {The dynamics of hydrogen desorption from H-terminated silicon surface clusters has been simulated in the framework of real space time-dependent density functional theory complemented with molecular dynamics for ions. It has been demonstrated that by choosing an appropriate frequency and intensity of the laser it is possible to remove the hydrogen layer from the surface without destroying the structure of underlying silicon. At the laser field intensities used in the current study (0.5-2.0 V/\AA) the desorption process is notably nonlinear.},
Doi = {10.1063/1.3580563},
Keywords = {ab initio calculations; density functional theory; desorption; hydrogen; hydrogenation; laser beam effects; molecular dynamics method},
Numpages = {3},
Publisher = {AIP}
}
@Article{bubin_jap_110_064905_2011,
Title = {Electron-ion dynamics in laser-assisted desorption of hydrogen atoms from H-Si(111) surface},
Author = {Sergiy Bubin and K\'{a}lm\'{a}n Varga},
Journal = {J. Appl. Phys.},
Year = {2011},
Number = {6},
Pages = {064905},
Volume = {110},
Abstract = {In the framework of real time real space time-dependent density functional theory we have studied the electron-ion dynamics of a hydrogen-terminated silicon surface H-Si(111) subjected to intense laser irradiation. Two surface fragments of different sizes have been used in the simulations. When the intensity and duration of the laser exceed certain levels (which depend on the wavelength) we observe the desorption of the hydrogen atoms, while the underlying silicon layer remains essentially undamaged. Upon further increase of the laser intensity, the chemical bonds between silicon atoms break as well. The results of the simulations suggest that with an appropriate choice of laser parameters it should be possible to remove the hydrogen layer from the H-Si(111) surface in a matter of a few tens of femtoseconds. We have also observed that at high laser field intensities (2-4 V/A in this work) the desorption occurs even when the laser frequency is smaller than the optical gap of the silicon surface fragments. Therefore, nonlinear phenomena must play an essential role in such desorption processes.},
Coden = {JAPIAU},
Doi = {10.1063/1.3638064},
Keywords = {bonds (chemical); density functional theory; desorption; elemental semiconductors; energy gap; hydrogen; laser beam effects; optical constants; silicon; },
Publisher = {AIP}
}
@Article{bubin_prb_85_205441_2012,
Title = {Electron and ion dynamics in graphene and graphane fragments subjected to high-intensity laser pulses},
Author = {Bubin, Sergiy and Varga, K\'alm\'an},
Journal = {Phys. Rev. B},
Year = {2012},
Pages = {205441},
Volume = {85},
Abstract = {In the framework of real-time real-space time-dependent density functional theory complemented with Ehrenfest molecular dynamics we have studied the response of small graphene and graphane fragments to intense femtosecond laser pulses. In particular, we have investigated how the response changes with laser pulses of different frequency (near IR, visible, and UV). The results of our simulations show that graphene has a very high immediate (i.e., within laser-pulse duration) damage threshold. They also suggest that, similar to the case of other hydrogenated surfaces, it should be possible to selectively desorb hydrogens from graphane without destroying the underlying carbon structure provided that the laser pulse parameters are properly chosen.},
Doi = {10.1103/PhysRevB.85.205441},
Issue = {20},
Numpages = {7},
Publisher = {American Physical Society}
}
@Article{bubin_prb_85_235435_2012,
Title = {Simulation of high-energy ion collisions with graphene fragments},
Author = {Bubin, Sergiy and Wang, Bin and Pantelides, Sokrates and Varga, K\'alm\'an},
Journal = {Phys. Rev. B},
Year = {2012},
Pages = {235435},
Volume = {85},
Abstract = {The collision of energetic ions and graphene fragments is studied in the framework of real-space finite-difference time-dependent density functional theory (TDDFT) coupled with classical molecular dynamics for nuclei. The amount of energy transferred from the projectile to the target is calculated to explore the defect formation mechanisms as a function of the projectile's energy. It is found that creation of defects in graphene due to the interaction of a fast proton with valence electrons is unlikely. In the case of projectiles with higher charges, the transferred energy increases significantly, leading to higher probability of bond breaking.},
Doi = {10.1103/PhysRevB.85.235435},
Issue = {23},
Numpages = {7},
Publisher = {American Physical Society}
}
@Article{driscoll_n_22_285702_2011,
Title = {Time-dependent density functional study of field emission from nanotubes composed of C, BN, SiC, Si, and GaN},
Author = {Joseph A Driscoll and Sergiy Bubin and William R French and Kalman Varga},
Journal = {Nanotechnology},
Year = {2011},
Number = {28},
Pages = {285702},
Volume = {22},
Abstract = {Field emission from various types of nanotubes is studied by propagating the electronic density in real space and time using time-dependent density functional theory. Capped (5, 5) C, BN, SiC, Si, and GaN nanotubes are considered. The GaN, SiC, and Si nanotubes were found to be significantly better field emitters than C and BN nanotubes, both in terms of current magnitude and sharpness of peaks in the energy spectra. By analyzing the electronic structure of the various systems it is seen that the nanotubes with the highest currents have electron densities that extend significantly from the nanotube in the emission direction.},
Doi = {10.1088/0957-4484/22/28/285702},
}
@Article{driscoll_prb_83_233405_2011,
Title = {Laser-induced electron emission from nanostructures: A first-principles study},
Author = {Driscoll, Joseph A. and Bubin, Sergiy and Varga, K\'alm\'an },
Journal = {Phys. Rev. B},
Year = {2011},
Number = {23},
Pages = {233405},
Volume = {83},
Abstract = {Time-dependent density functional theory simulation of laser-induced ionization is presented. Various test systems including a small wire-like molecule, C$_{12}$H$_{14}$, as well as carbon nanotubes with varying diameter are studied. It has been demonstrated that significant ionization electron current is produced when a laser pulse is applied. Moreover, pulse-like patterns of the current have been observed, which suggests that short laser pulses can be used to create spatially and temporally localized electron sources.},
Doi = {10.1103/PhysRevB.83.233405},
Numpages = {4},
Publisher = {American Physical Society}
}
@Article{driscoll_jap_110_024304_2011,
Title = {First-principles study of field emission from carbon nanotubes and graphene nanoribbons},
Author = {Joseph A. Driscoll and Brandon Cook and Sergiy Bubin and K\'{a}lm\'{a}n Varga},
Journal = {J. Appl. Phys.},
Year = {2011},
Pages = {024304},
Volume = {110},
Abstract = {A real-space, real-time implementation of time-dependent density functional theory is used to study electron field emission from nanostructures. Carbon nanotubes and graphene nanoribbons are used as model systems. The calculations show that carbon nanotubes with iron adsorbates have spin-polarized emission currents. Graphene nanoribbons are shown to be good field emitters with spatial variation of the emission current influenced by the presence of passivating hydrogen.},
Doi = {10.1063/1.3610511},
}
@Article{erattupuzha_jpb_50_125601_2017,
author = {Sonia Erattupuzha and Cody L Covington and Arthur Russakoff and Erik L\"otstedt and Seyedreza Larimian and V\'aclav Hanus and Sergiy Bubin and Markus Koch and Stefanie Gr\"afe and Andrius Baltu\v{s}ka and Xinhua Xie and Kaoru Yamanouchi and K\'alm\'an Varga and Markus Kitzler},
title = {Enhanced ionisation of polyatomic molecules in intense laser pulses is due to energy upshift and field coupling of multiple orbitals},
journal = {J. Phys. B},
year = {2017},
volume = {50},
number = {12},
pages = {125601},
abstract = {We present the results of a combined experimental and numerical study on strong-field ionisation of acetylene performed with the aim of identifying the mechanism behind the previously reported surprisingly large multi-electron ionisation probabilities of polyatomic molecules. Using coincidence momentum imaging techniques and time-dependent density functional simulations, we show that the reported efficient ionisation is due to the combined action of a significant geometrically induced energy upshift of the most relevant valence orbitals as the C--H distance stretches beyond about two times the equilibrium distance, and a strong increase in the coupling between multiple molecular orbitals concomitant with this stretch motion. The identified enhanced ionisation mechanism, which we refer to as EIC-MOUSE, is only effective for molecules aligned close to parallel to the laser polarisation direction, and is inhibited for perpendicularly aligned molecules because of a suppression of the C--H stretch motion during the onset of ionisation.},
doi = {10.1088/1361-6455/aa7098},
}
@Article{russakoff_pra_91_023422_2015,
Title = {Time-dependent density-functional study of the alignment-dependent ionization of acetylene and ethylene by strong laser pulses},
Author = {Russakoff, Arthur and Bubin, Sergiy and Xie, Xinhua and Erattupuzha, Sonia and Kitzler, Markus and Varga, K\'{a}lm\'{a}n},
Journal = {Phys. Rev. A},
Year = {2015},
Pages = {023422},
Volume = {91},
Doi = {10.1103/PhysRevA.91.023422},
Numpages = {11},
Publisher = {American Physical Society}
}
@Article{xie_pra_89_023429_2014,
Title = {Role of proton dynamics in efficient photoionization of hydrocarbon molecules},
Author = {Xie, Xinhua and Roither, Stefan and Sch\"offler, Markus and Xu, Huailiang and Bubin, Sergiy and L\"otstedt, Erik and Erattuphuza, Sonia and Iwasaki, Atsushi and Kartashov, Daniil and Varga, K\'alm\'an and G.~Paulus, Gerhard and Baltu\v{s}ka, Andrius and Yamanouchi, Kaoru and Kitzler, Markus},
Journal = {Phys. Rev. A},
Year = {2014},
Month = {Feb},
Pages = {023429},
Volume = {89},
Abstract = {We experimentally investigate the ionization mechanism behind the formation of remarkably high charge states observed in the laser-pulse-induced fragmentation of different hydrocarbon molecules by Roither et al. [Phys. Rev. Lett. 106, 163001 (2011)], who suggested enhanced ionization occurring at multiple C-H bonds as the underlying ionization mechanism. Using multiparticle coincidence momentum imaging we measure the yield of multiply charged fragmenting ethylene and acetylene molecules at several intensities and pulse durations ranging from the few-cycle regime to 25 fs. We observe, at constant intensity, a strong increase of the proton energy with increasing laser pulse duration. It is shown that this is caused by a strong increase in the yield of highly charged parent molecular ions with pulse duration. Based on experimental evidence we explain this increase by the necessary population of precursor states in the parent ion that feature fast C-H stretch dynamics to the critical internuclear distance, where efficient ionization via enhanced ionization takes place. For increasing pulse duration these precursor ionic states are more efficiently populated, which leads in turn to a higher enhanced-ionization probability for longer pulses. Our work provides experimental evidence for the existence of a multiple-bond version of enhanced ionization in polyatomic molecules.},
Doi = {10.1103/PhysRevA.89.023429},
Issue = {2},
Numpages = {9},
Publisher = {American Physical Society}
}
@Article{stanke_pra_97_012513_2018,
author = {Stanke, Monika and Bralin, Amir and Bubin, Sergiy and Adamowicz, Ludwik},
title = {Leading relativistic corrections for atomic $P$ states calculated with a finite-nuclear-mass approach and all-electron explicitly correlated Gaussian functions},
journal = {Phys. Rev. A},
year = {2018},
volume = {97},
pages = {012513},
month = {Jan},
abstract = {In this work we report progress in the development and implementation of quantum-mechanical methods for calculating bound ground and excited states of small atomic systems. The work concerns singlet states with the $L=1$ total orbital angular momentum ($P$ states). The method is based on the finite-nuclear-mass (non-Born-Oppenheimer; non-BO) approach and the use of all-particle explicitly correlated Gaussian functions for expanding the nonrelativistic wave function of the system. The development presented here includes derivation and implementation of algorithms for calculating the leading relativistic corrections for singlet states. The corrections are determined in the framework of the perturbation theory as expectation values of the corresponding effective operators using the non-BO wave functions. The method is tested in the calculations of the ten lowest $^1 P$ states of the helium atom and the four lowest $^1 P$ states of the beryllium atom.},
doi = {10.1103/PhysRevA.97.012513},
issue = {1},
numpages = {15},
publisher = {American Physical Society},
}
@Article{chavez_cpl_717_147_2019,
author = {Erik M. Chavez and Sergiy Bubin and Ludwik Adamowicz},
title = {Implementation of explicitly correlated complex Gaussian functions in calculations of molecular rovibrational $J=1$ states without Born-Oppenheimer approximation},
journal = {Chem. Phys. Lett.},
year = {2019},
volume = {717},
pages = {147-151},
abstract = {In our previous work (Bubin et al., 2016) it was established that complex explicitly-correlated one-center all-particle Gaussian functions (CECGs) provide an effective basis set for very accurate non-relativistic molecular non-Born-Oppenheimer (non-BO) calculations for vibrational ground and excited states corresponding to the rotational ground state. In this work we advance the molecular CECGs approach further by implementing and testing the algorithms for calculating the vibrational states corresponding to the first rotational excited state (the $J=1$ state). The test concerns all bound $J=1$ rovibrational states of the HD$^+$ ion.},
doi = {10.1016/j.cplett.2019.01.013},
}