-
Machine learning configuration interaction for ab initio potential energy curves
Authors:
J. P. Coe
Abstract:
The concept of machine learning configuration interaction (MLCI) [J. Chem. Theory Comput. 2018, 14, 5739], where an artificial neural network (ANN) learns on the fly to select important configurations, is further developed so that accurate ab initio potential energy curves can be efficiently calculated. This development includes employing the artificial neural network also as a hash function for t…
▽ More
The concept of machine learning configuration interaction (MLCI) [J. Chem. Theory Comput. 2018, 14, 5739], where an artificial neural network (ANN) learns on the fly to select important configurations, is further developed so that accurate ab initio potential energy curves can be efficiently calculated. This development includes employing the artificial neural network also as a hash function for the efficient deletion of duplicates on the fly so that the singles and doubles space does not need to be stored and this barrier to scalability is removed. In addition configuration state functions are introduced into the approach so that pure spin states are guaranteed, and the transferability of data between geometries is exploited. This improved approach is demonstrated on potential energy curves for the nitrogen molecule, water, and carbon monoxide. The results are compared with full configuration interaction values, when available, and different transfer protocols are investigated. It is shown that, for all of the considered systems, accurate potential energy curves can now be efficiently computed with MLCI. For the potential curves of N$_{2}$ and CO, MLCI can achieve lower errors than stochastically selecting configurations while also using substantially less processor hours.
△ Less
Submitted 30 October, 2019; v1 submitted 20 August, 2019;
originally announced August 2019.
-
Lattice Density-Functional Theory for Quantum Chemistry
Authors:
J. P. Coe
Abstract:
We propose a lattice density-functional theory for {\it ab initio} quantum chemistry or physics as a route to an efficient approach that approximates the full configuration interaction energy and orbital occupations for molecules with strongly-correlated electrons. We build on lattice density-functional theory for the Hubbard model by deriving Kohn-Sham equations for a reduced then full quantum ch…
▽ More
We propose a lattice density-functional theory for {\it ab initio} quantum chemistry or physics as a route to an efficient approach that approximates the full configuration interaction energy and orbital occupations for molecules with strongly-correlated electrons. We build on lattice density-functional theory for the Hubbard model by deriving Kohn-Sham equations for a reduced then full quantum chemistry Hamiltonian, and demonstrate the method on the potential energy curves for the challenging problem of modelling elongating bonds in a linear chain of six hydrogen atoms. Here the accuracy of the Bethe-ansatz local-density approximation is tested for this quantum chemistry system and we find that, despite this approximate functional being designed for the Hubbard model, the shapes of the potential curves generally agree with the full configuration interaction results. Although there is a discrepancy for very stretched bonds, this is lower than when using standard density-functional theory with the local-density approximation.
△ Less
Submitted 18 April, 2019; v1 submitted 22 February, 2019;
originally announced February 2019.
-
Machine Learning Configuration Interaction
Authors:
J. P. Coe
Abstract:
We propose the concept of machine learning configuration interaction (MLCI) whereby an artificial neural network is trained on-the-fly to predict important new configurations in an iterative selected configuration interaction procedure. We demonstrate that the neural network can discriminate between important and unimportant configurations, that it has not been trained on, much better than by chan…
▽ More
We propose the concept of machine learning configuration interaction (MLCI) whereby an artificial neural network is trained on-the-fly to predict important new configurations in an iterative selected configuration interaction procedure. We demonstrate that the neural network can discriminate between important and unimportant configurations, that it has not been trained on, much better than by chance. MLCI is then used to find compact wavefunctions for carbon monoxide at both stretched and equilibrium geometries. We also consider the multireference problem of the water molecule with elongated bonds. Results are contrasted with those from other ways of selecting configurations: first-order perturbation, random selection and Monte Carlo configuration interaction. Compared with these other serial calculations, this prototype MLCI is competitive in its accuracy, converges in significantly fewer iterations than the stochastic approaches, and requires less time for the higher-accuracy computations.
△ Less
Submitted 17 August, 2018;
originally announced August 2018.
-
Testing density-functional approximations on a lattice and the limits of the related Hohenberg-Kohn-type theorem
Authors:
V. V. França,
J. P. Coe,
I. D'Amico
Abstract:
We present a metric-space approach to quantify the performance of density-functional approximations for interacting many-body systems and to explore the validity of the Hohenberg-Kohn-type theorem on fermionic lattices. This theorem demonstrates the existence of one-to-one mappings between particle densities, wave functions and external potentials. We then focus on these quantities, and quantify h…
▽ More
We present a metric-space approach to quantify the performance of density-functional approximations for interacting many-body systems and to explore the validity of the Hohenberg-Kohn-type theorem on fermionic lattices. This theorem demonstrates the existence of one-to-one mappings between particle densities, wave functions and external potentials. We then focus on these quantities, and quantify how far apart in metric space the approximated and exact ones are. We apply our method to the one-dimensional Hubbard model for different types of external potentials, and assess its validity on one of the most used approximations in density-functional theory, the local density approximation (LDA). We find that the potential distance may have a very different behaviour from the density and wave function distances, in some cases even providing the wrong assessments of the LDA performance trends. We attribute this to the systems reaching behaviours which are borderline for the applicability of the one-to-one correspondence between density and external potential. On the contrary the wave function and density distances behave similarly and are always sensitive to system variations. Our metric-based method correctly predicts the regimes where the LDA performs fairly well and the regimes where it fails. This suggests that our method could be a practical tool for testing the efficiency of density-functional approximations.
△ Less
Submitted 29 August, 2017;
originally announced August 2017.
-
Uniqueness of density-to-potential mapping for fermionic lattice systems
Authors:
J. P. Coe,
I. D'Amico,
V. V. França
Abstract:
We demonstrate that, for a fermionic lattice system, the ground-state particle density uniquely determines the external potential except for the sites corresponding to nodes of the wave function, and the limiting case where the Pauli exclusion principle completely determines the occupation of all sites. Our fundamental finding completes, for this general class of systems, the one-to-one correspond…
▽ More
We demonstrate that, for a fermionic lattice system, the ground-state particle density uniquely determines the external potential except for the sites corresponding to nodes of the wave function, and the limiting case where the Pauli exclusion principle completely determines the occupation of all sites. Our fundamental finding completes, for this general class of systems, the one-to-one correspondence between ground states, their densities, and the external potential at the base of the Hohenberg-Kohn theorem. Moreover we demonstrate that the mapping from wave function to potential is unique not just for the ground state, but also for excited states. To illustrate our findings, we develop a practical inversion scheme to determine the external potential from a given density. Our results hold for a general class of lattice models, which includes the Hubbard model.
△ Less
Submitted 8 July, 2016;
originally announced July 2016.
-
Multireference X-Ray Emission and Absorption Spectroscopy calculations from Monte Carlo Configuration Interaction
Authors:
J. P. Coe,
M. J. Paterson
Abstract:
We adapt the method of Monte Carlo configuration interaction to calculate core-hole states and use this for the computation of X-ray emission and absorption values. We consider CO, CH$_{4}$, NH$_{3}$, H$_{2}$O, HF, HCN, CH$_{3}$OH, CH$_{3}$F, HCl and NO using a 6-311G** basis. We also look at carbon monoxide with a stretched geometry and discuss the dependence of its results on the cutoff used. Th…
▽ More
We adapt the method of Monte Carlo configuration interaction to calculate core-hole states and use this for the computation of X-ray emission and absorption values. We consider CO, CH$_{4}$, NH$_{3}$, H$_{2}$O, HF, HCN, CH$_{3}$OH, CH$_{3}$F, HCl and NO using a 6-311G** basis. We also look at carbon monoxide with a stretched geometry and discuss the dependence of its results on the cutoff used. The Monte Carlo configuration interaction results are compared with EOM-CCSD values for X-ray emission and with experiment for X-ray absorption. Oscillator strengths are also computed and we quantify the multireference nature of the wavefunctions to suggest when approaches based on a single reference would be expected to be successful.
△ Less
Submitted 6 May, 2015;
originally announced May 2015.
-
Approaching exact hyperpolarizabilities via sum-over-states Monte Carlo configuration interaction
Authors:
J. P. Coe,
M. J. Paterson
Abstract:
We propose using sum-over-states calculations with the compact wavefunctions of Monte Carlo configuration interaction to approach accurate values for higher-order dipole properties up to second hyperpolarizabilities in a controlled way. We apply the approach to small systems that can generally be compared with full configuration interaction (FCI) results. We consider hydrogen fluoride with a 6-31g…
▽ More
We propose using sum-over-states calculations with the compact wavefunctions of Monte Carlo configuration interaction to approach accurate values for higher-order dipole properties up to second hyperpolarizabilities in a controlled way. We apply the approach to small systems that can generally be compared with full configuration interaction (FCI) results. We consider hydrogen fluoride with a 6-31g basis and then look at results, including frequency dependent properties, in an aug-cc-pVDZ basis. We extend one calculation beyond FCI by using an aug-cc-pVTZ basis. The properties of an H$_{4}$ molecule with multireference character are calculated in an aug-cc-pVDZ basis. We then investigate this method on a strongly multireference system with a larger FCI space by modelling the properties of carbon monoxide with a stretched geometry. The behavior of the approach with increasing basis size is considered by calculating results for the neon atom using aug-cc-pVDZ to aug-cc-pVQZ. We finally test if the unusual change in polarizability between the first two states of molecular oxygen can be reproduced by this method in a 6-31g basis.
△ Less
Submitted 25 September, 2014;
originally announced September 2014.
-
Applying Monte Carlo configuration interaction to transition metal dimers: exploring the balance between static and dynamic correlation
Authors:
J. P. Coe,
P. Murphy,
M. J. Paterson
Abstract:
We calculate potential curves for transition metal dimers using Monte Carlo configuration interaction (MCCI). These results, and their associated spectroscopic values, are compared with experimental and computational studies. The multireference nature of the MCCI wavefunction is quantified and we estimate the important orbitals. We initially consider the ground state of the chromium dimer. Next we…
▽ More
We calculate potential curves for transition metal dimers using Monte Carlo configuration interaction (MCCI). These results, and their associated spectroscopic values, are compared with experimental and computational studies. The multireference nature of the MCCI wavefunction is quantified and we estimate the important orbitals. We initially consider the ground state of the chromium dimer. Next we calculate potential curves for Sc$_{2}$ where we contrast the lowest triplet and quintet states. We look at the molybdenum dimer where we compare non-relativistic results with the partial inclusion of relativistic effects via effective core potentials, and report results for scandium nickel.
△ Less
Submitted 22 May, 2014;
originally announced May 2014.
-
Characterizing a configuration interaction excited state using natural transition geminals
Authors:
J. P. Coe,
M. J. Paterson
Abstract:
We introduce natural transition geminals as a means to qualitatively understand a transition where double excitations are important. The first two $A_{1}$ singlet states of the CH cation are used as an initial example. We calculate these states with configuration interaction singles (CIS) and state-averaged Monte Carlo configuration interaction (SA-MCCI). For each method we compare the important n…
▽ More
We introduce natural transition geminals as a means to qualitatively understand a transition where double excitations are important. The first two $A_{1}$ singlet states of the CH cation are used as an initial example. We calculate these states with configuration interaction singles (CIS) and state-averaged Monte Carlo configuration interaction (SA-MCCI). For each method we compare the important natural transition geminals with the dominant natural transition orbitals. We then compare SA-MCCI and full configuration interaction (FCI) with regards to the natural transition geminals using the beryllium atom. We compare using the natural transition geminals with analyzing the important configurations in the CI expansion to give the dominant transition for the beryllium atom and the carbon dimer. Finally we calculate the natural transition geminals for two electronic excitations of formamide.
△ Less
Submitted 28 February, 2014;
originally announced February 2014.
-
State-averaged Monte Carlo configuration interaction applied to electronically excited states
Authors:
J. P. Coe,
M. J. Paterson
Abstract:
We introduce state-averaging into the method of Monte Carlo configuration interaction (SA-MCCI) to allow the stable and efficient calculation of excited states. We show that excited potential curves for H$_{3}$, including a crossing with the ground state, can be accurately reproduced using a small fraction of the FCI space. A recently introduced error measure for potential curves [J. P. Coe and M.…
▽ More
We introduce state-averaging into the method of Monte Carlo configuration interaction (SA-MCCI) to allow the stable and efficient calculation of excited states. We show that excited potential curves for H$_{3}$, including a crossing with the ground state, can be accurately reproduced using a small fraction of the FCI space. A recently introduced error measure for potential curves [J. P. Coe and M. J. Paterson, J. Chem. Phys., 137, 204108 (2012)] is shown to also be a fair approach when considering potential curves for multiple states. We demonstrate that potential curves for LiF using SA-MCCI agree well with the FCI results and the avoided crossing occurs correctly. The seam of conical intersections for CH$_{2}$ found by Yarkony [J. Chem. Phys., 104, 2932 (1996)] is used as a test for SA-MCCI and we compare potential curves from SA-MCCI with FCI results for this system for the first three triplet states. We then demonstrate the improvement from using SA-MCCI on the dipole of the $2$ $^{1}A_{1}$ state of carbon monoxide. We then look at vertical excitations for small organic molecules up to the size of butadiene where the SA-MCCI energies and oscillator strengths are compared with CASPT2 values [M. Schreiber, M. R. Silva-Junior, S. P. A. Sauer, and W. Thiel, J. Chem. Phys., 128, 134110 (2008)]. We finally see if the SA-MCCI results for these excitation energies can be improved by using MCCIPT2 with approximate natural orbitals when the PT2 space is not onerously large.
△ Less
Submitted 28 February, 2014;
originally announced February 2014.
-
Development of Monte Carlo configuration interaction: Natural orbitals and second-order perturbation theory
Authors:
J. P. Coe,
M. J. Paterson
Abstract:
Approximate natural orbitals are investigated as a way to improve a Monte Carlo configuration interaction (MCCI) calculation. We introduce a way to approximate the natural orbitals in MCCI and test these and approximate natural orbitals from MP2 and QCISD in MCCI calculations of single-point energies. The efficiency and accuracy of approximate natural orbitals in MCCI potential curve calculations…
▽ More
Approximate natural orbitals are investigated as a way to improve a Monte Carlo configuration interaction (MCCI) calculation. We introduce a way to approximate the natural orbitals in MCCI and test these and approximate natural orbitals from MP2 and QCISD in MCCI calculations of single-point energies. The efficiency and accuracy of approximate natural orbitals in MCCI potential curve calculations for the double hydrogen dissociation of water, the dissociation of carbon monoxide and the dissociation of the nitrogen molecule are then considered in comparison with standard MCCI when using full configuration interaction as a benchmark. We also use the method to produce a potential curve for water in an aug-cc-pVTZ basis. A new way to quantify the accuracy of a potential curve is put forward that takes into account all of the points and that the curve can be shifted by a constant. We adapt a second-order perturbation scheme to work with MCCI (MCCIPT2) and improve the efficiency of the removal of duplicate states in the method. MCCIPT2 is tested in the calculation of a potential curve for the dissociation of nitrogen using both Slater determinants and configuration state functions.
△ Less
Submitted 11 April, 2013;
originally announced April 2013.
-
Monte Carlo configuration interaction applied to multipole moments, ionisation energies and electron affinities
Authors:
J. P. Coe,
D. J. Taylor,
M. J. Paterson
Abstract:
The method of Monte Carlo configuration interaction (MCCI) [1,2] is applied to the calculation of multipole moments. We look at the ground and excited state dipole moments in carbon monoxide. We then consider the dipole of NO, the quadrupole of the nitrogen molecule and of BH. An octupole of methane is also calculated. We consider experimental geometries and also stretched bonds. We show that thes…
▽ More
The method of Monte Carlo configuration interaction (MCCI) [1,2] is applied to the calculation of multipole moments. We look at the ground and excited state dipole moments in carbon monoxide. We then consider the dipole of NO, the quadrupole of the nitrogen molecule and of BH. An octupole of methane is also calculated. We consider experimental geometries and also stretched bonds. We show that these non-variational quantities may be found to relatively good accuracy when compared with FCI results, yet using only a small fraction of the full configuration interaction space. MCCI results in the aug-cc-pVDZ basis are seen to generally have reasonably good agreement with experiment. We also investigate the performance of MCCI when applied to ionisation energies and electron affinities of atoms in an aug-cc-pVQZ basis. We compare the MCCI results with full configuration-interaction quantum Monte Carlo [3,4] and `exact' non-relativistic results [3,4]. We show that MCCI could be a useful alternative for the calculation of atomic ionisation energies however electron affinities appear much more challenging for MCCI. Due to the small magnitude of the electron affinities their percentage errors can be high, but with regards to absolute errors MCCI performs similarly for ionisation energies and electron affinities.
△ Less
Submitted 21 January, 2013;
originally announced January 2013.
-
Calculations of Potential Energy Surfaces Using Monte Carlo Configuration Interaction
Authors:
J. P. Coe,
D. J. Taylor,
M. J. Paterson
Abstract:
We apply the method of Monte Carlo configuration interaction (MCCI) to calculate ground-state potential energy curves for a range of small molecules and compare the results with full configuration interaction. We show that the MCCI potential energy curve can be calculated to relatively good accuracy, as quantified using the non-parallelity error, using only a very small fraction of the FCI space.…
▽ More
We apply the method of Monte Carlo configuration interaction (MCCI) to calculate ground-state potential energy curves for a range of small molecules and compare the results with full configuration interaction. We show that the MCCI potential energy curve can be calculated to relatively good accuracy, as quantified using the non-parallelity error, using only a very small fraction of the FCI space. In most cases the potential curve is of better accuracy than its constituent single-point energies. We finally test the MCCI program on systems with basis sets beyond full configuration interaction: a lattice of fifty hydrogen atoms and ethylene. The results for ethylene agree fairly well with other computational work while for the lattice of fifty hydrogens we find that the fraction of the full configuration interaction space we were able to consider appears to be too small as, although some qualitative features are reproduced, the potential curve is less accurate.
△ Less
Submitted 11 December, 2012;
originally announced December 2012.
-
Present status and first results of the final focus beam line at the KEK Accelerator Test Facility
Authors:
P. Bambade,
M. Alabau Pons,
J. Amann,
D. Angal-Kalinin,
R. Apsimon,
S. Araki,
A. Aryshev,
S. Bai,
P. Bellomo,
D. Bett,
G. Blair,
B. Bolzon,
S. Boogert,
G. Boorman,
P. N. Burrows,
G. Christian,
P. Coe,
B. Constance,
Jean-Pierre Delahaye,
L. Deacon,
E. Elsen,
A. Faus-Golfe,
M. Fukuda,
J. Gao,
N. Geffroy
, et al. (69 additional authors not shown)
Abstract:
ATF2 is a final-focus test beam line which aims to focus the low emittance beam from the ATF damping ring to a vertical size of about 37 nm and to demonstrate nanometer level beam stability. Several advanced beam diagnostics and feedback tools are used. In December 2008, construction and installation were completed and beam commissioning started, supported by an international team of Asian, Europe…
▽ More
ATF2 is a final-focus test beam line which aims to focus the low emittance beam from the ATF damping ring to a vertical size of about 37 nm and to demonstrate nanometer level beam stability. Several advanced beam diagnostics and feedback tools are used. In December 2008, construction and installation were completed and beam commissioning started, supported by an international team of Asian, European, and U.S. scientists. The present status and first results are described.
△ Less
Submitted 5 July, 2012;
originally announced July 2012.
-
Pinning effect and QPT-like behavior for two particles confined by a core-shell potential
Authors:
Pier Paolo Marchisio,
Jeremy P. Coe,
Irene D'Amico
Abstract:
We study the ground state entanglement, energy and fidelities of a two-electron system bounded by a core-shell potential, where the core width is varied continuously until it eventually vanishes. This simple system displays a rich and complex behavior: as the core width is varied, this system is characterized by two peculiar transitions where, for different reasons, it displays characteristics sim…
▽ More
We study the ground state entanglement, energy and fidelities of a two-electron system bounded by a core-shell potential, where the core width is varied continuously until it eventually vanishes. This simple system displays a rich and complex behavior: as the core width is varied, this system is characterized by two peculiar transitions where, for different reasons, it displays characteristics similar to a few-particle quantum phase transition. The first occurrence corresponds to something akin to a second order quantum phase transition, while the second transition is marked by a discontinuity, with respect to the driving parameter, in the first derivatives of quantities like energy and entanglement. The study of this system allows to shed light on the sudden variation of entanglement and energy observed in S. Abdullah et al. [Phys. Rev. B 80, 235302 (2009)]. We also compare the core-shell system with a system where a core well is absent: this shows that, even when extremely narrow, the core well has a relevant `pinning' effect. Interestingly, depending on the potential symmetry, the pinning of the wavefunction may either halve or double the system entanglement (with respect to the no-core-well system) when the ground state is already bounded to the outer (shell) well. In the process we discuss the system fidelity and show the usefulness of considering the particle density fidelity as opposed to the more commonly used -- but much more difficult to access -- wavefunction fidelity. In particular we demonstrate that -- for ground-states with nodeless spatial wavefunctions -- the particle density fidelity is zero if and only if the wavefunction fidelity is zero.
△ Less
Submitted 17 January, 2012;
originally announced January 2012.
-
Quantum mechanics in metric space: wave functions and their densities
Authors:
I. D'Amico,
J. P. Coe,
V. V. Franca,
K. Capelle
Abstract:
Hilbert space combines the properties of two fundamentally different types of mathematical spaces: vector space and metric space. While the vector-space aspects of Hilbert space, such as formation of linear combinations of state vectors, are routinely used in quantum mechanics, the metric-space aspects of Hilbert space are much less exploited. Here we show that a suitable metric stratifies Fock sp…
▽ More
Hilbert space combines the properties of two fundamentally different types of mathematical spaces: vector space and metric space. While the vector-space aspects of Hilbert space, such as formation of linear combinations of state vectors, are routinely used in quantum mechanics, the metric-space aspects of Hilbert space are much less exploited. Here we show that a suitable metric stratifies Fock space into concentric spheres. Maximum and minimum distances between wave functions are derived and geometrically interpreted in terms of this metric. Unlike the usual Hilbert-space analysis, our results apply also to the reduced space of only ground-state wave functions and to that of particle densities, each of which forms a metric space but not a Hilbert space. The Hohenberg-Kohn mapping between densities and ground-state wave functions, which is highly complex and nonlocal in coordinate description, is found, for three different model systems, to be very simple in metric space, where it is represented by a monotonic mapping of vicinities onto vicinities. Surprisingly, it is also found to be nearly linear over a wide range of parameters.
△ Less
Submitted 11 February, 2011;
originally announced February 2011.
-
Feasibility of approximating spatial and local entanglement in long-range interacting systems using the extended Hubbard model
Authors:
J. P. Coe,
V. V. França,
I. D'Amico
Abstract:
We investigate the extended Hubbard model as an approximation to the local and spatial entanglement of a one-dimensional chain of nanostructures where the particles interact via a long range interaction represented by a `soft' Coulomb potential. In the process we design a protocol to calculate the particle-particle spatial entanglement for the Hubbard model and show that, in striking contrast with…
▽ More
We investigate the extended Hubbard model as an approximation to the local and spatial entanglement of a one-dimensional chain of nanostructures where the particles interact via a long range interaction represented by a `soft' Coulomb potential. In the process we design a protocol to calculate the particle-particle spatial entanglement for the Hubbard model and show that, in striking contrast with the loss of spatial degrees of freedom, the predictions are reasonably accurate. We also compare results for the local entanglement with previous results found using a contact interaction (PRA, 81 (2010) 052321) and show that while the extended Hubbard model recovers a better agreement with the entanglement of a long-range interacting system, there remain realistic parameter regions where it fails to predict the quantitative and qualitative behaviour of the entanglement in the nanostructure system.
△ Less
Submitted 5 January, 2011; v1 submitted 30 July, 2010;
originally announced July 2010.
-
Geometry induced entanglement transitions in nanostructures
Authors:
J. P. Coe,
S. Abdullah,
I. D'Amico
Abstract:
We model quantum dot nanostructures using a one-dimensional system of two interacting electrons. We show that strong and rapid variations may be induced in the spatial entanglement by varying the nanostructure geometry. We investigate the position-space information entropy as an indicator of the entanglement in this system. We also consider the expectation value of the Coulomb interaction and t…
▽ More
We model quantum dot nanostructures using a one-dimensional system of two interacting electrons. We show that strong and rapid variations may be induced in the spatial entanglement by varying the nanostructure geometry. We investigate the position-space information entropy as an indicator of the entanglement in this system. We also consider the expectation value of the Coulomb interaction and the ratio of this expectation to the expectation of the confining potential and their link to the entanglement. We look at the first derivative of the entanglement and the position-space information entropy to infer information about a possible quantum phase transition.
△ Less
Submitted 19 March, 2010;
originally announced March 2010.
-
The Hubbard model as an approximation to the entanglement in nanostructures
Authors:
J. P. Coe,
V. V. França,
I. D'Amico
Abstract:
We investigate how well the one-dimensional Hubbard model describes the entanglement of particles trapped in a string of quantum wells. We calculate the average single-site entanglement for two particles interacting via a contact interaction and consider the effect of varying the interaction strength and the interwell distance. We compare the results with the ones obtained within the one-dimens…
▽ More
We investigate how well the one-dimensional Hubbard model describes the entanglement of particles trapped in a string of quantum wells. We calculate the average single-site entanglement for two particles interacting via a contact interaction and consider the effect of varying the interaction strength and the interwell distance. We compare the results with the ones obtained within the one-dimensional Hubbard model with on-site interaction. We suggest an upper bound for the average single-site entanglement for two electrons in M wells and discuss analytical limits for very large repulsive and attractive interactions. We investigate how the interplay between interaction and potential shape in the quantum well system dictates the position and size of the entanglement maxima and the agreement with the theoretical limits. Finally we calculate the spatial entanglement for the quantum well system and compare it to its average single-site entanglement.
△ Less
Submitted 18 March, 2010;
originally announced March 2010.
-
The entanglement of few-particle systems when using the local-density approximation
Authors:
J. P. Coe,
I. D'Amico
Abstract:
In this chapter we discuss methods to calculate the entanglement of a system using density-functional theory. We firstly introduce density-functional theory and the local-density approximation (LDA). We then discuss the concept of the `interacting LDA system'. This is characterised by an interacting many-body Hamiltonian which reproduces, uniquely and exactly, the ground state density obtained f…
▽ More
In this chapter we discuss methods to calculate the entanglement of a system using density-functional theory. We firstly introduce density-functional theory and the local-density approximation (LDA). We then discuss the concept of the `interacting LDA system'. This is characterised by an interacting many-body Hamiltonian which reproduces, uniquely and exactly, the ground state density obtained from the single-particle Kohn-Sham equations of density-functional theory when the local-density approximation is used. We motivate why this idea can be useful for appraising the local-density approximation in many-body physics particularly with regards to entanglement and related quantum information applications. Using an iterative scheme, we find the Hamiltonian characterising the interacting LDA system in relation to the test systems of Hooke's atom and helium-like atoms. The interacting LDA system ground state wavefunction is then used to calculate the spatial entanglement and the results are compared and contrasted with the exact entanglement for the two test systems. For Hooke's atom we also compare the entanglement to our previous estimates of an LDA entanglement. These were obtained using a combination of evolutionary algorithm and gradient descent, and using an LDA-based perturbative approach. We finally discuss if the position-space information entropy of the density---which can be obtained directly from the system density and hence easily from density-functional theory methods---can be considered as a proxy measure for the spatial entanglement for the test systems.
△ Less
Submitted 10 March, 2010;
originally announced March 2010.
-
Effect of confinement potential geometry on entanglement in quantum dot-based nanostructures
Authors:
S. Abdullah,
J. P. Coe,
I. D'Amico
Abstract:
We calculate the spatial entanglement between two electrons trapped in a nanostructure for a broad class of confinement potentials, including single and double quantum dots, and core-shell quantum dot structures.
By using a parametrized confinement potential, we are able to switch from one structure to the others with continuity and to analyze how the entanglement is influenced by the changes…
▽ More
We calculate the spatial entanglement between two electrons trapped in a nanostructure for a broad class of confinement potentials, including single and double quantum dots, and core-shell quantum dot structures.
By using a parametrized confinement potential, we are able to switch from one structure to the others with continuity and to analyze how the entanglement is influenced by the changes in the confinement geometry. We calculate the many-body wave function by `exact' diagonalization of the time independent Schrödinger equation. We discuss the relationship between the entanglement and specific cuts of the wave function, and show that the wave function at a single highly symmetric point could be a good indicator for the entanglement content of the system. We analyze the counterintuitive relationship between spatial entanglement and Coulomb interaction, which connects maxima (minima) of the first to minima (maxima) of the latter. We introduce a potential quantum phase transition which relates quantum states characterized by different spatial topology. Finally we show that by varying shape, range and strength of the confinement potential, it is possible to induce strong and rapid variations of the entanglement between the two electrons. This property may be used to tailor nanostructures according to the level of entanglement required by a specific application.
△ Less
Submitted 2 December, 2009; v1 submitted 19 August, 2009;
originally announced August 2009.
-
Exact and LDA entanglement of tailored densities in an interacting one-dimensional electron system
Authors:
J. P. Coe,
I. D'Amico
Abstract:
We calculate the `exact' potential corresponding to a one-dimensional interacting system of two electrons with a specific, tailored density. We use one-dimensional density-functional theory with a local-density approximation (LDA) on the same system and calculate densities and energies, which are compared with the `exact' ones. The `interacting-LDA system' corresponding to the LDA density is the…
▽ More
We calculate the `exact' potential corresponding to a one-dimensional interacting system of two electrons with a specific, tailored density. We use one-dimensional density-functional theory with a local-density approximation (LDA) on the same system and calculate densities and energies, which are compared with the `exact' ones. The `interacting-LDA system' corresponding to the LDA density is then found and its potential compared with the original one. Finally we calculate and compare the spatial entanglement of the electronic systems corresponding to the interacting-LDA and original interacting system.
△ Less
Submitted 28 May, 2009;
originally announced May 2009.
-
Reverse engineering in many-body quantum physics: What many-body system corresponds to an effective single-particle equation?
Authors:
J. P. Coe,
K. Capelle,
I. D'Amico
Abstract:
The mapping, exact or approximate, of a many-body problem onto an effective single-body problem is one of the most widely used conceptual and computational tools of physics. Here, we propose and investigate the inverse map of effective approximate single-particle equations onto the corresponding many-particle system. This approach allows us to understand which interacting system a given single-p…
▽ More
The mapping, exact or approximate, of a many-body problem onto an effective single-body problem is one of the most widely used conceptual and computational tools of physics. Here, we propose and investigate the inverse map of effective approximate single-particle equations onto the corresponding many-particle system. This approach allows us to understand which interacting system a given single-particle approximation is actually describing, and how far this is from the original physical many-body system. We illustrate the resulting reverse engineering process by means of the Kohn-Sham equations of density-functional theory. In this application, our procedure sheds light on the non-locality of the density-potential mapping of DFT, and on the self-interaction error inherent in approximate density functionals.
△ Less
Submitted 19 January, 2009; v1 submitted 3 September, 2008;
originally announced September 2008.
-
Entanglement and density-functional theory: testing approximations on Hooke's atom
Authors:
J. P. Coe,
A. Sudbery,
I. D'Amico
Abstract:
We present two methods of calculating the spatial entanglement of an interacting electron system within the framework of density-functional theory. These methods are tested on the model system of Hooke's atom for which the spatial entanglement can be calculated exactly. We analyse how the strength of the confining potential affects the spatial entanglement and how accurately the methods that we in…
▽ More
We present two methods of calculating the spatial entanglement of an interacting electron system within the framework of density-functional theory. These methods are tested on the model system of Hooke's atom for which the spatial entanglement can be calculated exactly. We analyse how the strength of the confining potential affects the spatial entanglement and how accurately the methods that we introduced reproduce the exact trends. We also compare the results with the outcomes of standard first-order perturbation methods. The accuracies of energies and densities when using these methods are also considered.
△ Less
Submitted 3 August, 2010; v1 submitted 21 December, 2007;
originally announced December 2007.