An efficient and accurate decomposition of the Fermi operator.
We present a method to compute the Fermi function of the Hamiltonian for a system of independent fermions based on an exact decomposition of the grand-canonical potential. This scheme does not rely on the localization of the orbitals and is insensitive to ill-conditioned Hamiltonians. It lends itsel...
Main Authors: | , , |
---|---|
Formato: | Journal article |
Idioma: | English |
Publicado: |
2008
|
_version_ | 1826293056751534080 |
---|---|
author | Ceriotti, M Kühne, T Parrinello, M |
author_facet | Ceriotti, M Kühne, T Parrinello, M |
author_sort | Ceriotti, M |
collection | OXFORD |
description | We present a method to compute the Fermi function of the Hamiltonian for a system of independent fermions based on an exact decomposition of the grand-canonical potential. This scheme does not rely on the localization of the orbitals and is insensitive to ill-conditioned Hamiltonians. It lends itself naturally to linear scaling as soon as the sparsity of the system's density matrix is exploited. By using a combination of polynomial expansion and Newton-like iterative techniques, an arbitrarily large number of terms can be employed in the expansion, overcoming some of the difficulties encountered in previous papers. Moreover, this hybrid approach allows us to obtain a very favorable scaling of the computational cost with increasing inverse temperature, which makes the method competitive with other Fermi operator expansion techniques. After performing an in-depth theoretical analysis of computational cost and accuracy, we test our approach on the density functional theory Hamiltonian for the metallic phase of the LiAl alloy. |
first_indexed | 2024-03-07T03:24:12Z |
format | Journal article |
id | oxford-uuid:b8786ccc-e07e-42a9-99f3-e9b0dd7b6aa7 |
institution | University of Oxford |
language | English |
last_indexed | 2024-03-07T03:24:12Z |
publishDate | 2008 |
record_format | dspace |
spelling | oxford-uuid:b8786ccc-e07e-42a9-99f3-e9b0dd7b6aa72022-03-27T04:56:10ZAn efficient and accurate decomposition of the Fermi operator.Journal articlehttp://purl.org/coar/resource_type/c_dcae04bcuuid:b8786ccc-e07e-42a9-99f3-e9b0dd7b6aa7EnglishSymplectic Elements at Oxford2008Ceriotti, MKühne, TParrinello, MWe present a method to compute the Fermi function of the Hamiltonian for a system of independent fermions based on an exact decomposition of the grand-canonical potential. This scheme does not rely on the localization of the orbitals and is insensitive to ill-conditioned Hamiltonians. It lends itself naturally to linear scaling as soon as the sparsity of the system's density matrix is exploited. By using a combination of polynomial expansion and Newton-like iterative techniques, an arbitrarily large number of terms can be employed in the expansion, overcoming some of the difficulties encountered in previous papers. Moreover, this hybrid approach allows us to obtain a very favorable scaling of the computational cost with increasing inverse temperature, which makes the method competitive with other Fermi operator expansion techniques. After performing an in-depth theoretical analysis of computational cost and accuracy, we test our approach on the density functional theory Hamiltonian for the metallic phase of the LiAl alloy. |
spellingShingle | Ceriotti, M Kühne, T Parrinello, M An efficient and accurate decomposition of the Fermi operator. |
title | An efficient and accurate decomposition of the Fermi operator. |
title_full | An efficient and accurate decomposition of the Fermi operator. |
title_fullStr | An efficient and accurate decomposition of the Fermi operator. |
title_full_unstemmed | An efficient and accurate decomposition of the Fermi operator. |
title_short | An efficient and accurate decomposition of the Fermi operator. |
title_sort | efficient and accurate decomposition of the fermi operator |
work_keys_str_mv | AT ceriottim anefficientandaccuratedecompositionofthefermioperator AT kuhnet anefficientandaccuratedecompositionofthefermioperator AT parrinellom anefficientandaccuratedecompositionofthefermioperator AT ceriottim efficientandaccuratedecompositionofthefermioperator AT kuhnet efficientandaccuratedecompositionofthefermioperator AT parrinellom efficientandaccuratedecompositionofthefermioperator |