Quantum dynamics of the CH3 fragment: a curvilinear coordinate system and kinetic energy operators

A curvilinear coordinate system for AB(3) fragments is given. The corresponding exact kinetic energy operator is derived and a series of simpler, progressively more approximate kinetic energy operators are suggested. The operators are tailored for quantum dynamics simulations using the multiconfigurational time-dependent Hartree approach. It is outlined how these fragment coordinates can be utilized to set up coordinate systems for larger systems such as AB(3)C or AB(3)CD. Calculations of the vibrational levels of CH(3) and quantum dynamics studies investigate the accuracy of the different kinetic energy operators suggested.     

Information-theoretic approach to kinetic-energy functionals: the nearly uniform electron gas

We strengthen the connection between information theory and quantum mechanical systems using a recently developed dequantization procedure which results in a decomposition of the kinetic energy as the sum of a classical term and a purely quantum term. For the nearly uniform electron gas, we thereby approximate the noninteracting kinetic energy as the sum of the Thomas-Fermi term, which is exact for the uniform electron gas, and the Weizsäcker term, which is proportional to the Fisher information. Electron correlation is included via a nonlocal analytical expression which is a functional of the (N-1)-conditional probability density. This expression is evaluated via a statistically rigorous Monte-Carlo procedure to obtain the correlation energy as a functional of the electron density. We show that this functional is well aproximated by a term which is proportional to the Shannon entropy. Thus the kinetic energy is expressed as the standard Thomas-Fermi term plus terms which are proportional to two of the cornerstones of information theory: the Fisher information, which is a measure of localization, and the Shannon entropy, which is a measure of delocalization.     

Link between the kinetic‐ and exchange‐energy functionals in the generalized gradient approximation

An approximate kinetic‐energy functional of the generalized gradient approximation form was derived following the “conjointness conjecture” of Lee, Lee, and Parr. The functional shares the analytical form of its gradient dependency with the exchange‐energy functionals of Becke and Perdew, Burke, and Ernzerhof. The two free parameters of this functional were determined using the exact values of the kinetic energy of He and Xe atoms. A set of 12 closed‐shell atoms was used to test the accuracy of the proposed functional and more than 30 others taken from the literature. It is shown that the conjointness conjecture leads to a very good class of kinetic‐energy functionals. Moreover, the functional developed in this work is shown to be one of the most accurate despite its simple analytical form. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Quantum tunneling in an exactly solvable double-barrier potential: barrier transmission near the zero energy

The one-dimensional wave equation is solved in the presence of a symmetric double-barrier potential. An exact, analytical solution is obtained for the scattering states. The transmission and reflection amplitudes are calculated using the method of logarithmic derivative of the exact wave function. In the case of double-barrier potentials, perfect transmission (or zero reflectivity) at zero kinetic energy is non-intuitive. This phenomenon has been revealed and called the “threshold anomaly” in the previous investigations. Here we show that it is a critical phenomenon provided that the inter-barrier distance satisfies a resonance condition. When the resonance condition is fulfilled, the perfect transmission occurs at any energy, including the zero one.     

Kinetic Energy Partition Method for Competing Modes

We present a new basis set expansion method for quantum dynamics systems with two competing modes where the interaction potentials are equally dominant. The new idea introduced here is a kinetic energy partition scheme instead of the usual division of the potential energy. The partition results in two kinetic energy terms with their effective masses. By distributing each partial kinetic energy to the respective potential, the full Hamiltonian can be expressed as the sum of the two competing modes. The solution procedure is illustrated by using a system consisting of a particle under the action of two harmonic potentials with different equilibrium distances and force constants. Next we apply this method to obtain the potential energy curves for the prototype hydrogen molecule ion. This new expansion converges very fast to the exact solutions for both eigenvalues and eigenfunctions.     

Atoms-in-molecules in momentum space: a Hirshfeld partitioning of electron momentum densities

A direct application of the Hirshfeld atomic partitioning (HAP) scheme is implemented for molecular electron momentum densities (EMDs). The momentum density contributions of individual atoms in diverse molecular systems are analyzed along with their topographical features and the kinetic energies of the atomic partitions. The proposed p-space HAP-based charge scheme does seem to possess the desirable attributes expected of any atoms in molecules partitioning. In addition to this, the main strength of the p-space HAP is the exact knowledge of the kinetic energy functional and the inherent ease in computing the kinetic energy. The charges derived from HAP in momentum space are found to match chemical intuition and the generally known chemical characteristics such as electronegativity, etc.     

Kinetic analysis of solid-state reactions: Evaluation of approximations to temperature integral and their applications

The temperature integral, which has no exact analytical solution, is involved in the analysis of the experiment data obtained under nonisothermal conditions. Some approximations for the temperature integral have been proposed in the literature for the determination of the kinetic parameters, in particular the activation energy. Those approximations are classified into two categories, that is, exponential and rational approximations. The precision of them for estimating the temperature integral was evaluated within a certain continuous range rather than at several discrete points. Some applications of the approximations in the kinetic methods were presented. The relative errors of the activation energy and pre-exponential factor with four rational approximations by employing model-fitting method were calculated. The relative errors of the activation energy for a series of conversion rate with four rational and four exponential approximations by employing linear integral isoconversional methods were evaluated.     

The molecular potential energy surface and vibrational energy levels of methyl fluoride. Part II

New analytical bending and stretching, ground electronic state, potential energy surfaces for CH(3)F are reported. The surfaces are expressed in bond-length, bond-angle internal coordinates. The four-dimensional stretching surface is an accurate, least squares fit to over 2000 symmetrically unique ab initio points calculated at the CCSD(T) level. Similarly, the five-dimensional bending surface is a fit to over 1200 symmetrically unique ab initio points. This is an important first stage towards a full nine-dimensional potential energy surface for the prototype CH(3)F molecule. Using these surfaces, highly excited stretching and (separately) bending vibrational energy levels of CH(3)F are calculated variationally using a finite basis representation method. The method uses the exact vibrational kinetic energy operator derived for XY(3)Z systems by Manson and Law (preceding paper, Part I, Phys. Chem. Chem. Phys., 2006, 8, DOI: 10.1039/b603106d). We use the full C(3v) symmetry and the computer codes are designed to use an arbitrary potential energy function. Ultimately, these results will be used to design a compact basis for fully coupled stretch-bend calculations of the vibrational energy levels of the CH(3)F system.     

Numerical investigation of a lower bound to the Thomas—Fermi kinetic energy and the validity of the Lieb conjecture for H−, He,Li+, B3+, O6+, Ne8+, and Mg12+

The purpose of the present work was twofold: First, a lower bound of Gálvez and Dehesa to the Thomas-Fermi kinetic energy was investigated for the two-electron systems H⁻, He, Li⁺, B³⁺, O⁶⁺, Ne⁸⁺, and Mg¹²⁺. Second, a conjecture of Lieb, relating the kinetic energy to the Thomas-Fermi kinetic energy, was examined. For both investigations, the analytical approximations of Benesch to the radial electron densities of the two-electron systems were used. These approximations are based on the explicitly correlated (Hylleraas-type) 20 variational parameter wave functions of Hart and Herzberg. © 1996 John Wiley & Sons, Inc.     

Energy‐surfaces from the upper bound of the Pauli kinetic energy

Based on the Kohn–Sham Pauli potential and the Kohn–Sham electron density, the upper bound of the Pauli kinetic energy is tested as a suitable replacement for the exact Pauli kinetic energy for application in orbital‐free density functional calculations. It is found that bond lengths for strong and moderately bound systems can be qualitatively predicted, but with a systematic shift toward larger bond distances with a relative error of 6% up to 30%. Angular dependence of the energy‐surface cannot be modeled with the proposed functional. Therefore, the upper bound model is the first parameter‐free functional expression for the kinetic energy that is able to qualitatively reproduce binding curves with respect to bond distortions. © 2016 Wiley Periodicals, Inc.     

Kinetic analysis of solid‐state reactions: Precision of the activation energy calculated by integral methods

The integral methods are extensively used for the kinetic analysis of solid‐state reactions. As the Arrhenius integral function [p(x)] does not have an exact analytical solution, different approximated equations have been proposed in the literature for performing the kinetic analysis of experimental integral data. Since the first approximation of Van Krevelen, a large number of equations have been proposed with the objective of increasing the precision in the determination of the Arrhenius integral, as checked from the standard deviation of the approximated function with regard to the real exact value of the integral. However, the main application of these equations is the determination of the kinetic parameters, in particular activation energies, and not the computation of the Arrhenius integral. A systematic analysis of the errors involved in the determination of the activation energy from these integral methods is still missing. A comparative study of the precision of the activation energy as a function of x and T computed from the different integral methods has been carried out. © 2005 Wiley Periodicals, Inc. Int J Chem Kinet 37: 658–666, 2005     

Method for the design of broad energy range focusing reflectrons

A novel method for the design of reflections capable of focusing large kinetic energy ranges is presented. The design method itself is a numerical approach that provides a geometrically flexible alternative to traditional analytical design solutions. This design method has been used to produce a reflectron that provides unit mass resolution for product spectra in a tandem reflectron time-of-flight (TOF) mass spectrometer despite a kinetic energy range of 1950–2700 eV. In this application, the systematic progression of reflectron design results in a practical, nonlinear field reflectron with the use of only two grids. Design improvements are proposed for more flexible systems, although geometric constraints in the current instrument limit their experimental evaluation.     

Application of Exact Analytic Total Energy Functional for Hooke’s Atom to He, Li+ and Be++: An Examination of the Universality of the Energy Functional in DFT

We employ the recently generated energy density functional for Hooke’s atom [Ludeña EV et al (2004) Intern J Quantum Chem 99:297], to which we introduce a simplification for the kinetic energy term, to evaluate the total energy of the helium atom and of the two-electron ions Li⁺ and Be⁺⁺. Using accurate representations for the one-particle densities of these systems we show that the energy density functional for Hooke’s atom leads, in these cases, to energy values that are below the exact ones. We discuss the implication of this finding with respect to the existence of a universal energy functional in DFT.     

Numeric kinetic energy operators for molecules in polyspherical coordinates

Generalized curvilinear coordinates, as, e.g., polyspherical coordinates, are in general better adapted to the resolution of the nuclear Schro?dinger equation than rectilinear ones like the normal mode coordinates. However, analytical expressions of the kinetic energy operators (KEOs) for molecular systems in polyspherical coordinates may be prohibitively complicated for large systems. In this paper we propose a method to generate a KEO numerically and bring it to a form practicable for dynamical calculations. To examine the new method we calculated vibrational spectra and eigenenergies for nitrous acid (HONO) and compare it with results obtained with an exact analytical KEO derived previously [F. Richter, P. Rosmus, F. Gatti, and H.-D. Meyer, J. Chem. Phys. 120, 6072 (2004)]. In a second example we calculated π → π* photoabsorption spectrum and eigenenergies of ethene (C(2)H(4)) and compared it with previous work [M. R. Brill, F. Gatti, D. Lauvergnat, and H.-D. Meyer, Chem. Phys. 338, 186 (2007)]. In this ethene study the dimensionality was reduced from 12 to 6 by freezing six internal coordinates. Results for both molecules show that the proposed method for obtaining an approximate KEO is reliable for dynamical calculations. The error in eigenenergies was found to be below 1 cm(-1) for most states calculated.     

Density functional theory for molecular and periodic systems using density fitting and continuous fast multipole method: Stress tensor

A full implementation of the analytical stress tensor for periodic systems is reported in the TURBOMOLE program package within the framework of Kohn–Sham density functional theory using Gaussian-type orbitals as basis functions. It is the extension of the implementation of analytical energy gradients (Lazarski et al., Journal of Computational Chemistry 2016, 37, 2518–2526) to the stress tensor for the purpose of optimization of lattice vectors. Its key component is the efficient calculation of the Coulomb contribution by combining density fitting approximation and continuous fast multipole method. For the exchange-correlation (XC) part the hierarchical numerical integration scheme (Burow and Sierka, Journal of Chemical Theory and Computation 2011, 7, 3097–3104) is extended to XC weight derivatives and stress tensor. The computational efficiency and favorable scaling behavior of the stress tensor implementation are demonstrated for various model systems. The overall computational effort for energy gradient and stress tensor for the largest systems investigated is shown to be at most two and a half times the computational effort for the Kohn–Sham matrix formation. © 2019 Wiley Periodicals, Inc.     

Photokinetic methods: A mathematical analysis of the rate equations in photochromic systems

In this article the analytical integration of kinetic equations describing the dynamic behavior of reversible photoreactions has been undertaken. Photochemical systems giving rise to competing thermal and photochemical steps are kinetically analyzed; the rate law is set up and analytical solutions are obtained under precisely defined boundary conditions. More complicated kinetic systems, where several species interact both thermally and photochemically, are also investigated. The kinetic treatment leads to a system of coupled differential equations which are amenable to analytical solutions, under the appropriate experimental and boundary conditions. The usefulness of the equations developed is illustrated by their application to some spirooxazine and chromene photochromic systems: two examples are described in detail. © 1999 John Wiley & Sons, Inc., Int J Chem Kinet 31: 303–313, 1999     

Positive equilibria of weakly reversible power law kinetic systems with linear independent interactions

In this paper, we extend our study of power law kinetic systems whose kinetic order vectors (which we call “interactions”) are reactant-determined (i.e. reactions with the same reactant complex have identical vectors) and are linear independent per linkage class. In particular, we consider PL-TLK systems, i.e. such whose T-matrix (the matrix with the interactions as columns indexed by the reactant complexes), when augmented with the rows of characteristic vectors of the linkage classes, has maximal column rank. Our main result states that any weakly reversible PL-TLK system has a complex balanced equilibrium. On the one hand, we consider this result as a “Higher Deficiency Theorem” for such systems since in our previous work, we derived analogues of the Deficiency Zero and the Deficiency One Theorems for mass action kinetics (MAK) systems for them, thus covering the “Low Deficiency” case. On the other hand, our result can also be viewed as a “Weak Reversibility Theorem” (WRT) in the sense that the statement “any weakly reversible system with a kinetics from the given set has a positive equilibrium” holds. According to the work of Deng et al. and more recently of Boros, such a WRT holds for MAK systems. However, we show that a WRT does not hold for two proper MAK supersets: the set PL-NIK of non-inhibitory power law kinetics (i.e. all kinetic orders are non-negative) and the set PL-FSK of factor span surjective power law kinetics (i.e. different reactants imply different interactions).