Search

In order to simulate the exact universal density functional, approximations are nowadays constructed by permitting more flexibility in its ansatz. In view of the difficulty of defining a systematically improvable form for it, this paper argues that an alternative way could be considered. It falls within the class of hybrid functionals with multi-determinant wave functions. The parameter controlling the hybridization is considered as variable. The invariance of the exact result with respect to changes in this variable is used to introduce information about the system under consideration, and to correct the density functional result. The construction considered in this paper accelerates convergence from the model system to the physical one, in the vicinity of the latter. The method, at the present level of implementation, should be seen as a starting point for further development, and not necessarily as a computationally advantageous tool. [ABSTRACT FROM AUTHOR]

We present a new procedure for computing a rovibrational spectrum of a polyatomic molecule and apply it to methane. The Schrödinger equation is solved, numerically exactly, by using a nested contracted basis. Rovibrational wavefunctions are computed in a |v>|JKM> basis, where |v> is a vibrational wavefunction and |JKM> is a symmetric top wavefunction. In turn, the |v> are obtained by solving a vibrational Schrödinger equation with basis functions that are products of contracted bend and stretch functions. At all stages of the calculation we exploit parity symmetry. The calculations are done in internal coordinates that facilitate the treatment of large amplitude motion. An Eckart molecule-fixed frame is used by numerically computing coefficients of the kinetic energy operator. The efficacy of the method is demonstrated by calculating a large number of converged J = 10 methane rovibrational levels in the Tetradecad polyad. No previous calculation of rovibrational levels of methane includes as many levels as we report in this paper. [ABSTRACT FROM AUTHOR]

We show that by using a rectangular collocation method with a small basis of parameterized functions, it is possible to compute a vibrational spectrum by solving the Schrödinger equation in 7D from a small number of ab initio calculations without a potential surface. The method is ideal for spectra of molecules adsorbed on a surface. In this paper, it is applied to calculate experimentally relevant energy levels of acetic acid adsorbed on the (101) surface of anatase TiO2. In this case, to obtain levels of experimental accuracy, increasing the number of dimensions from 4 to 7 increases the number of required potential points from about 1000 to about 10 000 and the number of basis functions from 126 to 792: the scaling is very attractive. [ABSTRACT FROM AUTHOR]

We study, analytically as well as numerically, the dynamics that arises from the interaction of a polar polarizable rigid rotor with single unipolar electromagnetic pulses of varying length, Δτ, with respect to the rotational period of the rotor, τr. In the sudden, non-adiabatic limit, Δτ ≪ τr, we derive analytic expressions for the rotor's wavefunctions, kinetic energies, and field-free evolution of orientation and alignment. We verify the analytic results by solving the corresponding time-dependent Schrödinger equation numerically and extend the temporal range of the interactions considered all the way to the adiabatic limit, Δτ > τr, where general analytic solutions beyond the field-free case are no longer available. The effects of the orienting and aligning interactions as well as of their combination on the post-pulse populations of the rotational states are visualized as functions of the orienting and aligning kick strengths in terms of population quilts. Quantum carpets that encapsulate the evolution of the rotational wavepackets provide the space-time portraits of the resulting dynamics. The population quilts and quantum carpets reveal that purely orienting, purely aligning, or even-break combined interactions each exhibit sui generis dynamics. In the intermediate temporal regime, we find that the wavepackets as functions of the orienting and aligning kick strengths show resonances that correspond to diminished kinetic energies at particular values of the pulse duration. [ABSTRACT FROM AUTHOR]

A vibrational SCF model is presented in which the functions forming the single-mode functions in the product wavefunction are expressed in terms of internal coordinates and the coordinates used for each mode are optimized variationally. This model involves no approximations to the kinetic energy operator and does not require a Taylor-series expansion of the potential. The non-linear optimization of coordinates is found to give much better product wavefunctions than the limited variations considered in most previous applications of SCF methods to vibrational problems. The approach is tested using published potential energy surfaces for water, ammonia, and formaldehyde. Variational flexibility allowed in the current ansätze results in excellent zero-point energies expressed through single-product states and accurate fundamental transition frequencies realized by short configuration-interaction expansions. Fully variational optimization of single-product states for excited vibrational levels also is discussed. The highlighted methodology constitutes an excellent starting point for more sophisticated treatments, as the bulk characteristics of many-mode coupling are accounted for efficiently in terms of compact wavefunctions (as evident from the accurate prediction of transition frequencies). [ABSTRACT FROM AUTHOR]

Hylleraas-CI is a superposition-of-configurations method in which each configuration is constructed from a Slater-type orbital (STO) product to which is appended (linearly) at most one interelectron distance rij. Computations of the kinetic energy for atoms by this method have been difficult due to the lack of formulas expressing these matrix elements for general angular momentum in terms of overlap and potential-energy integrals. It is shown here that a strategic application of angular-momentum theory, including the use of vector spherical harmonics, enables the reduction of all atomic kinetic-energy integrals to overlap and potential-energy matrix elements. The new formulas are validated by showing that they yield correct results for a large number of integrals published by other investigators. [ABSTRACT FROM AUTHOR]

Accurate densities were determined from configuration interaction wave functions for atoms and ions of Li, Be, and B with up to four electrons. Exchange-correlation potentials, Vxc(r), and functional derivatives of the noninteracting kinetic energy, δK[p]/δp(r), obtained from these densities were used to examine their discontinuities as the number of electrons N increases across integer boundaries for N = 1, N = 2, and N = 3. These numerical results are consistent with conclusions that the discontinuities are characterized by a jump in the chemical potential while the shape of Vxc(r) varies continuously as an integer boundary is crossed. The discontinuity of the Vxc(r) is positive, depends on the ionization potential, electron affinity, and orbital energy differences, and the discontinuity in δK[p]/δp(r) depends on the difference between the energies of the highest occupied and lowest unoccupied orbitals. The noninteracting kinetic energy and the exchange correlation energy have been computed for integer and noninteger values of N between 1 and 4. [ABSTRACT FROM AUTHOR]

Similarity of equations of motion for the classical and quantum trajectories is used to introduce a friction term dependent on the wavefunction phase into the time-dependent Schrödinger equation. The term describes irreversible energy loss by the quantum system. The force of friction is proportional to the velocity of a quantum trajectory. The resulting Schrödinger equation is nonlinear, conserves wavefunction normalization, and evolves an arbitrary wavefunction into the ground state of the system (of appropriate symmetry if applicable). Decrease in energy is proportional to the average kinetic energy of the quantum trajectory ensemble. Dynamics in the high friction regime is suitable for simple models of reactions proceeding with energy transfer from the system to the environment. Examples of dynamics are given for single and symmetric and asymmetric double well potentials. [ABSTRACT FROM AUTHOR]