Engineering Papers⌕ Search

Engineering topics

Chan, Garnet Kin-Lic

Publications and source records attributed to Chan, Garnet Kin-Lic.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Hamiltonian symmetries in auxiliary-field quantum Monte Carlo calculations for electronic structure

Here, we describe how to incorporate symmetries of the Hamiltonian into auxiliary-field quantum Monte Carlo (AFQMC) calculations. Focusing on the case of Abelian symmetries, we show that the computational cost of most steps of an AFQMC calculation is reduced by $N^{–1}_k$ , where $N_k$ is the number of irreducible representations of the symmetry group. We apply the formalism to a molecular system as well as to several crystalline solids. In the latter case, the lattice translational group provides increasing savings as the number of k points is increased, which is important in enabling calculations that approach the thermodynamic limit. The extension to non-Abelian symmetries is briefly discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Efficient Ab Initio Auxiliary-Field Quantum Monte Carlo Calculations in Gaussian Bases via Low-Rank Tensor Decomposition

We describe an algorithm to reduce the cost of auxiliary-field quantum Monte Carlo (AFQMC) calculations for the electronic structure problem. The technique uses a nested low-rank factorization of the electron repulsion integral (ERI). While the cost of conventional AFQMC calculations in Gaussian bases scales as O(N 4 ), where N is the size of the basis, we show that ground-state energies can be computed through tensor decomposition with reduced memory requirements and subquartic scaling. Here, the algorithm is applied to hydrogen chains and square grids, water clusters, and hexagonal BN. In all cases, we observe significant memory savings and, for larger systems, reduced, subquartic simulation time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Exact Fluctuations of Nonequilibrium Steady States from Approximate Auxiliary Dynamics

We describe a framework to reduce the computational effort to evaluate large deviation functions of time integrated observables within nonequilibrium steady states. We do this by incorporating an auxiliary dynamics into trajectory based Monte Carlo calculations, through a transformation of the system's propagator using an approximate guiding function. This procedure importance samples the trajectories that most contribute to the large deviation function, mitigating the exponential complexity of such calculations. Here, we illustrate the method by studying driven diffusion and interacting lattice models in one and two spatial dimensions. Finally, our work offers an avenue to calculate large deviation functions for high dimensional systems driven far from equilibrium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Importance sampling large deviations in nonequilibrium steady states. I

Large deviation functions hold information on the stability and response of systems driven into nonequilibrium steady states and in such a way are similar to free energies for systems at equilibrium. As with equilibrium free energies, evaluating large deviation functions numerically for all but the simplest systems is difficult because by construction they depend on exponentially rare events. In this first paper of a series, we assess different trajectory-based sampling methods capable of computing large deviation functions of time integrated observables within nonequilibrium steady states. We highlight some convergence criteria and best practices using a number of different models, including a biased Brownian walker, a driven lattice gas, and a model of self-assembly. We show how two popular methods for sampling trajectory ensembles, transition path sampling and diffusion Monte Carlo, suffer from exponentially diverging correlations in trajectory space as a function of the bias parameter when estimating large deviation functions. Improving the efficiencies of these algorithms requires introducing guiding functions for the trajectories.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Stripe order in the underdoped region of the two-dimensional Hubbard model

Competing inhomogeneous orders are a central feature of correlated electron materials, including the high-temperature superconductors. The two-dimensional Hubbard model serves as the canonical microscopic physical model for such systems. Multiple orders have been proposed in the underdoped part of the phase diagram, which corresponds to a regime of maximum numerical difficulty. By combining the latest numerical methods in exhaustive simulations, we uncover the ordering in the underdoped ground state. We find a stripe order that has a highly compressible wavelength on an energy scale of a few kelvin, with wavelength fluctuations coupled to pairing order. The favored filled stripe order is different from that seen in real materials. Our results demonstrate the power of modern numerical methods to solve microscopic models, even in challenging settings.

Zheng, Bo-Xiao↗
DOE OSTI2017

Gaussian and plane-wave mixed density fitting for periodic systems

We present a mixed density fitting scheme that utilizes both a Gaussian and a plane-wave fitting basis to accurately evaluate electron repulsion integrals in crystalline systems. We use this scheme to enable efficient all-electron Gaussian based periodic density functional and Hartree-Fock calculations

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Time-dependent N -electron valence perturbation theory with matrix product state reference wavefunctions for large active spaces and basis sets: Applications to the chromium dimer and all-trans polyenes

In earlier work [A. Y. Sokolov and G. K.-L. Chan, J. Chem. Phys. 144, 064102 (2016)], we introduced a time-dependent formulation of the second-order N -electron valence perturbation theory (t-NEVPT2) which (i) had a lower computational scaling than the usual internally contracted perturbation formulation and (ii) yielded the fully uncontracted NEVPT2 energy. We present a combination of t-NEVPT2 with a matrix product state (MPS) reference wavefunction (t-MPS-NEVPT2) that allows us to compute uncontracted dynamic correlation energies for large active spaces and basis sets, using the time-dependent density matrix renormalization group algorithm. In addition, we report a low-scaling MPS-based implementation of strongly contracted NEVPT2 (sc-MPS-NEVPT2) that avoids computation of the four-particle reduced density matrix. We use these new methods to compute the dissociation energy of the chromium dimer and to study the low-lying excited states in all-trans polyenes (C 4 H 6 to C 24 H 26 ), incorporating dynamic correlation for reference wavefunctions with up to 24 active electrons and orbitals

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Final Technical Report for Quantum Embedding for Correlated Electronic Structure in Large Systems and the Condensed Phase

This is the final technical report. We briefly describe some selected results below. Developments in density matrix embedding. DMET is a quantum embedding theory that we introduced at the beginning of the last funding period, around 2012-2013. Since the first DMET papers, which demonstrated proof-of- principle calculations on the Hubbard model and hydrogen rings, we have carried out a number of different developments, including: Extending the DMET technology to compute broken symmetry phases, including magnetic phases and super- conductivity (Pub. 13); Calibrating the accuracy of DMET and its cluster size convergence against other methods, and formulation of a dynamical cluster analog (Pubs. 4, 10) (see Fig. 1); Implementing DMET for ab-initio molecular calculations, and exploring different self-consistency criteria (Pubs. 9, 14); Using embedding to defi ne quantum classical interfaces Pub. 2; Formulating DMET for spectral functions (Pub. 7) (see Fig. 1); Extending DMET to coupled fermion-boson problems (Pub. 12). Together with these embedding developments, we have also implemented a wide variety of impurity solvers within our DMET framework, including DMRG (Pub. 3), AFQMC (Pub. 10), and coupled cluster theory (CC) (Pub. 9).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Gaussian-Based Coupled-Cluster Theory for the Ground-State and Band Structure of Solids

We present the results of Gaussian-based ground-state and excited-state equation-of-motion coupled-cluster theory with single and double excitations for three-dimensional solids. We focus on diamond and silicon, which are paradigmatic covalent semiconductors. In addition to ground-state properties (the lattice constant, bulk modulus, and cohesive energy), we compute the quasiparticle band structure and band gap. We sample the Brillouin zone with up to 64 k-points using norm-conserving pseudopotentials and polarized double- and triple-ζ basis sets, leading to canonical coupled-cluster calculations with as many as 256 electrons in 2176 orbitals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Cluster size convergence of the density matrix embedding theory and its dynamical cluster formulation: A study with an auxiliary-field quantum Monte Carlo solver

In this work, we investigate the cluster size convergence of the energy and observables using two forms of density matrix embedding theory (DMET): the original cluster form (CDMET) and a new formulation motivated by the dynamical cluster approximation (DCA-DMET). Both methods are applied to the half-filled one- and two-dimensional Hubbard models using a sign-problem free auxiliary-field quantum Monte Carlo impurity solver, which allows for the treatment of large impurity clusters of up to 100 sites. While CDMET is more accurate at smaller impurity cluster sizes, DCA-DMET exhibits faster asymptotic convergence towards the thermodynamic limit. We use our two formulations to produce new accurate estimates for the energy and local moment of the two-dimensional Hubbard model for U/t=2,4,6. These results compare favorably with the best data available in the literature, and help resolve earlier uncertainties in the moment for U/t=2.

36 MATERIALS SCIENCE↗
DOE OSTI2016

From plane waves to local Gaussians for the simulation of correlated periodic systems

Here, we present a simple, robust, and black-box approach to the implementation and use of local, periodic, atom-centered Gaussian basis functions within a plane wave code, in a computationally efficient manner. The procedure outlined is based on the representation of the Gaussians within a finite bandwidth by their underlying plane wave coefficients. The core region is handled within the projected augment wave framework, by pseudizing the Gaussian functions within a cutoff radius around each nucleus, smoothing the functions so that they are faithfully represented by a plane wave basis with only moderate kinetic energy cutoff. To mitigate the effects of the basis set superposition error and incompleteness at the mean-field level introduced by the Gaussian basis, we also propose a hybrid approach, whereby the complete occupied space is first converged within a large plane wave basis, and the Gaussian basis used to construct a complementary virtual space for the application of correlated methods. We demonstrate that these pseudized Gaussians yield compact and systematically improvable spaces with an accuracy comparable to their non-pseudized Gaussian counterparts. A key advantage of the described method is its ability to efficiently capture and describe electronic correlation effects of weakly bound and low-dimensional systems, where plane waves are not sufficiently compact or able to be truncated without unphysical artifacts. We investigate the accuracy of the pseudized Gaussians for the water dimer interaction, neon solid, and water adsorption on a LiH surface, at the level of second-order Møller–Plesset perturbation theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Matrix product operators, matrix product states, and ab initio density matrix renormalization group algorithms

Current descriptions of the ab initio density matrix renormalization group (DMRG) algorithm use two superficially different languages: an older language of the renormalization group and renormalized operators, and a more recent language of matrix product states and matrix product operators. The same algorithm can appear dramatically different when written in the two different vocabularies. Here, we carefully describe the translation between the two languages in several contexts. First, we describe how to efficiently implement the ab initio DMRG sweep using a matrix product operator based code, and the equivalence to the original renormalized operator implementation. Next we describe how to implement the general matrix product operator/matrix product state algebra within a pure renormalized operator-based DMRG code. Finally, we discuss two improvements of the ab initio DMRG sweep algorithm motivated by matrix product operator language: Hamiltonian compression, and a sum over operators representation that allows for perfect computational parallelism. The connections and correspondences described here serve to link the future developments with the past and are important in the efficient implementation of continuing advances in ab initio DMRG and related algorithms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

A state interaction spin-orbit coupling density matrix renormalization group method

We describe a state interaction spin-orbit (SISO) coupling method using density matrix renormalization group (DMRG) wavefunctions and the spin-orbit mean-field (SOMF) operator. We implement our DMRG-SISO scheme using a spin-adapted algorithm that computes transition density matrices between arbitrary matrix product states. To demonstrate the potential of the DMRG-SISO scheme we present accurate benchmark calculations for the zero-field splitting of the copper and gold atoms, comparing to earlier complete active space self-consistent-field and second-order complete active space perturbation theory results in the same basis. We also compute the effects of spin-orbit coupling on the spin-ladder of the iron-sulfur dimer complex [Fe{sub 2}S{sub 2}(SCH{sub 3}){sub 4}]{sup 3−}, determining the splitting of the lowest quartet and sextet states. We find that the magnitude of the zero-field splitting for the higher quartet and sextet states approaches a significant fraction of the Heisenberg exchange parameter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗