Researcher profile

Ankit Mahajan

Ankit Mahajan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2022arXiv

Selected configuration interaction wave functions in phaseless auxiliary field quantum Monte Carlo

We present efficient algorithms for using selected configuration interaction (sCI) trial wave functions in phaseless auxiliary field quantum Monte Carlo (ph-AFQMC). These advancements, geared towards optimizing computational performance for longer CI expansions, allow us to use up to a million configurations feasibly in the trial state for ph-AFQMC. In one example, we found the cost of ph-AFQMC per sample to increase only by a factor of about $3$ for a calculation with $10^4$ configurations compared to that with a single one, demonstrating the tiny computational overhead due to a longer expansion. This favorable scaling allows us to study the systematic convergence of the phaseless bias in AFQMC calculations with an increasing number of configurations and provides a means to gauge the accuracy of ph-AFQMC with other trial states. We also show how the scalability issues of sCI trial states for large system sizes could be mitigated by restricting them to a moderately sized orbital active space and leveraging the near-cancellation of out of active space phaseless errors.

preprint2020arXiv

Efficient local energy evaluation for multi-Slater wave functions in orbital space quantum Monte Carlo

Recent developments in selected configuration interaction methods have led to increased interest in using multi-Slater trial wave functions in various quantum Monte Carlo (QMC) methods. Here we present an algorithm for calculating the local energy of a multi-Slater wave function in orbital space QMC. For an ab initio Hamiltonian, our algorithm has a cost scaling of $O(n^5 + n_c)$, as opposed to the $O(n^4n_c)$ scaling of existing orbital space algorithms, where $n$ is the system size, and $n_c$ is the number of configurations in the wave function. We present our method using variational Monte Carlo calculations with the Jastrow multi-Slater wave function, although the formalism should be applicable for auxiliary field quantum Monte Carlo. We apply it to polyacetylene and demonstrate the possibility of using a much larger number of configurations than possible using existing methods.

preprint2020arXiv

Efficient multireference perturbation theory without high-order reduced density matrices

We present a stochastic approach to perform strongly contracted n-electron valence state perturbation theory (SC-NEVPT), which only requires one- and two-body reduced density matrices, without introducing approximations. We use this method to perform SC-NEVPT2 for CASSCF wave functions obtained from selected configuration interaction, although the approach is applicable to a larger class of wave functions, including those from variational Monte Carlo (VMC). The accuracy of this approach is demonstrated for small test systems, and the scaling is investigated with the number of virtual orbitals and the molecule size. We also find the SC-NEVPT2 energy to be relatively insensitive to the quality of the reference wave function. Finally, the method is applied to the Fe(II)-Porphyrin system with a (32e, 29o) active space, and to the isomerization of Cu2O2 in a (28e, 32o) active space.

preprint2019arXiv

An accelerated linear method for optimizing non-linear wavefunctions in variational Monte Carlo

Although the linear method is one of the most robust algorithms for optimizing non-linearly parametrized wavefunctions in variational Monte Carlo, it suffers from a memory bottleneck due to the fact at each optimization step a generalized eigenvalue problem is solved in which the Hamiltonian and overlap matrices are stored in memory. Here we demonstrate that by applying the Jacobi-Davidson algorithm, one can solve the generalized eigenvalue problem iteratively without having to build and store the matrices in question. The resulting direct linear method greatly lowers the cost and improves the scaling of the algorithm with respect to the number of parameters. To further improve the efficiency of optimization for wavefunctions with a large number of parameters, we use the first order method AMSGrad far from the minimum as it is very inexpensive, and only switch to the direct linear method near the end of the optimization where methods such as AMSGrad have long convergence tails. We apply this improved optimizer to various wavefunctions with both real and orbital space Jastrow factors for atomic systems such as Beryllium and Neon, molecular systems such as the Carbon dimer and Iron(II) Porphyrin, and model systems such as the Hubbard model and Hydrogen chains.

preprint2019arXiv

Multireference configuration interaction and perturbation theory without reduced density matrices

The computationally expensive evaluation and storage of high-rank reduced density matrices (RDMs) has been the bottleneck in the calculation of dynamic correlation for multireference wave functions in large active spaces. We present a stochastic formulation of multireference configuration interaction (MRCI) and perturbation theory (MRPT) that avoids the need for these expensive RDMs. The algorithm presented here is flexible enough to incorporate a wide variety of active space reference wave functions, including selected configuration interaction, matrix product states, and symmetry-projected Jastrow mean field wave functions. It enjoys the usual attractive features of Monte Carlo methods, such as embarrassing parallelizability and low memory costs. We find that the stochastic algorithm is already competitive with the deterministic algorithm for small active spaces, containing as few as 14 orbitals. We illustrate the utility of our stochastic formulation using benchmark applications.