Entanglement as Measure of Electron-Electron Correlation in Quantum Chemistry Calculations

In quantum chemistry calculations, the correlation energy is defined as the difference between the Hartree-Fock limit energy and the exact solution of the nonrelativistic Schrodinger equation. With this definition, the electron correlation effects are not directly observable. In this report, we show that the entanglement can be used as an alternative measure of the electron correlation in quantum chemistry calculations. Entanglement is directly observable and it is one of the most striking properties of quantum mechanics. As an example we calculate the entanglement for He atom and H2 molecule with different basis sets.Comment: 12 pages, 2 figure

Correlation energy, pair-distribution functions and static structure factors of jellium

We discuss and clarify a simple and accurate interpolation scheme for the spin-resolved electron static structure factor (and corresponding pair correlation function) of the 3D unpolarized homogeneous electron gas which, along with some analytic properties of the spin-resolved pair-correlation functions, we have just published. We compare our results with the very recent spin-resolved scheme by Schmidt et al., and focus our attention on the spin-resolved correlation energies and the high-density limit of the correlation functions.Comment: 8 pages, 3 figures. Proceedings of the conference on Statistical Mechanics and Strongly Correlated Systems (Bachelet, Parisi & Vulpiani Eds.) to appear as a special issue of Physica A (Elsevier, Amsterdam 2000

Computational Complexity in Electronic Structure

In quantum chemistry, the price paid by all known efficient model chemistries is either the truncation of the Hilbert space or uncontrolled approximations. Theoretical computer science suggests that these restrictions are not mere shortcomings of the algorithm designers and programmers but could stem from the inherent difficulty of simulating quantum systems. Extensions of computer science and information processing exploiting quantum mechanics has led to new ways of understanding the ultimate limitations of computational power. Interestingly, this perspective helps us understand widely used model chemistries in a new light. In this article, the fundamentals of computational complexity will be reviewed and motivated from the vantage point of chemistry. Then recent results from the computational complexity literature regarding common model chemistries including Hartree-Fock and density functional theory are discussed.Comment: 14 pages, 2 figures, 1 table. Comments welcom

Role of Micellar Entanglement Density on Kinetics of Shear Banding Flow Formation

We investigate the effects of micellar entanglement number on the kinetics of shear banding flow formation in a Taylor-Couette flow. Three sets of wormlike micellar solutions, each set with a similar fluid elasticity and zero-shear-rate viscosity, but with varying entanglement densities, are studied under start-up of steady shear. Our experiments indicate that in the set with the low fluid elasticity, the transient shear banding flow is characterized by the formation of a transient flow reversal in a range of entanglement densities. Outside of this range, the transient flow reversal is not observed. For the sets of medium and high elasticities, the transient flow reversals exist for relatively small entanglement densities, and disappear for large entanglement densities. Our analysis shows that wall slip and elastic instabilities do not affect the transient flow feature. We identify a correlation between micellar entanglement number, the width of the stress plateau, and the extent of the transient flow reversal. As the micellar entanglement number increases, the width of the stress plateau first increases, then, at a higher micellar entanglement number, plateau width decreases. Therefore, we hypothesize that the transient flow reversal is connected to the micellar entanglement number through the width of the stress plateau

Semiclassical Dynamics with Quantum Trajectories: Formulation and Comparison with the Semiclassical Initial Value Representation Propagator

We present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths

Introduction to Quantum Algorithms for Physics and Chemistry

In this introductory review, we focus on applications of quantum computation to problems of interest in physics and chemistry. We describe quantum simulation algorithms that have been developed for electronic-structure problems, thermal-state preparation, simulation of time dynamics, adiabatic quantum simulation, and density functional theory.Comment: 44 pages, 5 figures; comments or suggestions for improvement are welcom

Energy Conserving Approximations to the Quantum Potential: Dynamics with Linearized Quantum Force

Solution of the Schrödinger equation within the de Broglie–Bohm formulation is based on propagation of trajectories in the presence of a nonlocal quantum potential. We present a new strategy for defining approximate quantum potentials within a restricted trial function by performing the optimal fit to the log-derivatives of the wave function density. This procedure results in the energy-conserving dynamics for a closed system. For one particular form of the trial function leading to the linear quantum force, the optimization problem is solved analytically in terms of the first and second moments of the weighted trajectory distribution. This approach gives exact time-evolution of a correlated Gaussian wave function in a locally quadratic potential. The method is computationally cheap in many dimensions, conserves total energy and satisfies the criterion on the average quantum force. Expectation values are readily found by summing over trajectory weights. Efficient extraction of the phase-dependent quantities is discussed. We illustrate the efficiency and accuracy of the linear quantum force approximation by examining a one-dimensional scattering problem and by computing the wavepacket reaction probability for the hydrogen exchange reaction and the photodissociation spectrum of ICN in two dimensions

Bohmian Dynamics on Subspaces Using Linearized Quantum Force

In the de Broglie–Bohm formulation of quantum mechanics the time-dependent Schrödinger equation is solved in terms of quantum trajectories evolving under the influence of quantum and classical potentials. For a practical implementation that scales favorably with system size and is accurate for semiclassical systems, we use approximate quantum potentials. Recently, we have shown that optimization of the nonclassical component of the momentum operator in terms of fitting functions leads to the energy-conserving approximate quantum potential. In particular, linear fitting functions give the exact time evolution of a Gaussian wave packet in a locally quadratic potential and can describe the dominant quantum-mechanical effects in the semiclassical scattering problems of nuclear dynamics. In this paper we formulate the Bohmian dynamics on subspaces and define the energy-conserving approximate quantum potential in terms of optimized nonclassical momentum, extended to include the domain boundary functions. This generalization allows a better description of the non-Gaussian wave packets and general potentials in terms of simple fitting functions. The optimization is performed independently for each domain and each dimension. For linear fitting functions optimal parameters are expressed in terms of the first and second moments of the trajectory distribution. Examples are given for one-dimensional anharmonic systems and for the collinear hydrogen exchange reaction

Stable Long-Time Semiclassical Description of Zero-Point Energy in High-Dimensional Molecular Systems

Semiclassical implementation of the quantum trajectory formalism [J. Chem. Phys. 120, 1181 (2004)] is further developed to give a stable long-time description of zero-point energy in anharmonic systems of high dimensionality. The method is based on a numerically cheap linearized quantum force approach; stabilizing terms compensating for the linearization errors are added into the time-evolution equations for the classical and nonclassical components of the momentum operator. The wave function normalization and energy are rigorously conserved. Numerical tests are performed for model systems of up to 40 degrees of freedom