Researcher profile

Himadri Pathak

Himadri Pathak contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
3topics
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

10 published item(s)

preprint2026arXiv

Extending the Handover-Iterative VQE to Challenging Strongly Correlated Systems: $N_2$ and Fe-S Cluster

Accurately describing strongly correlated electronic systems remains a central challenge in quantum chemistry, as electron-electron interactions give rise to complex many-body wavefunctions that are difficult to capture with conventional approximations. Classical wavefunction-based approaches, such as the Semistochastic Heat-bath Configuration Interaction (SHCI) and the Density Matrix Renormalization Group (DMRG), currently define the state of the art, systematically converging toward the Full Configuration Interaction (FCI) limit, but at a rapidly increasing computational cost. Quantum computing algorithms promise to alleviate this scaling bottleneck by leveraging entanglement and superposition to represent correlated states more compactly. We introduced the Handover-Iterative Variational Quantum Eigensolver (HI-VQE) as a practical quantum computing algorithm with an iterative "handover" mechanism that dynamically exchanges information between quantum and classical computers, even using Noisy Intermediate-Scale Quantum (NISQ) computers. In this work, we extend the HI-VQE to benchmark two prototypical strongly correlated systems, the nitrogen molecule $N_2$ and iron-sulfur (Fe-S) cluster, which serve as stringent tests for both classical and quantum electronic-structure methods. By comparing HI-VQE results against Heat-bath Configuration Interaction (HCI) benchmarks, we assess its accuracy, scalability, and ability to capture multireference correlation effects. Achieving quantitative agreement on these canonical systems demonstrates a viable pathway toward quantum-enhanced simulations of complex bioinorganic molecules, catalytic mechanisms, and correlated materials.

preprint2026arXiv

Time-dependent Hole States in Multiconfigurational Time-Dependent Hartree-Fock Approaches: A Time-Domain Generalization of Extended Koopmans' Theorem

We introduce a framework for resolving electron-hole dynamics within wavefunction-based multiconfigurational time-dependent Hartree-Fock (MCTDHF) theory. Central to this framework is a time-domain generalization of the extended Koopmans' theorem, which rigorously defines time-dependent hole states through single-electron removal. From this foundation, we prove the existence of exact equations of motion for time-dependent Dyson orbitals, enabling instantaneous construction of photofragments' reduced density matrices. The formalism further yields a systematic procedure to extract hole-resolved observables, such as channel-resolved photoelectron momentum distributions, directly from time-dependent \textit{ab initio} wavefunctions. As a demonstration, we employ an attosecond $ω-2ω$ laser strategy to control hole dynamics, thereby resolving a long-standing challenge in MCTDHF simulations. This advance opens a pathway for exploring correlated multielectron dynamics in atoms and molecules under ultrafast laser fields.

preprint2025arXiv

Time-dependent Hole States in Multiconfigurational Time-Dependent Hartree-Fock Approaches: Applications in Photoionization of Water Molecule

By simulating the real-time multielectron wavefunction with the multi-configurational time-dependent Hartree-Fock (MCTDHF) approach, we conduct an \textit{ab initio} study of the single-photon ionization process of a body-fixed water molecule ($\mathrm{H_2O}$) driven by attosecond pulses. To this end, we present a full-dimensional implementation of the MCTDHF method based on one-center expansions, allowing for the simulation of arbitrarily polarized lasers and multi-center polyatomic potentials. With a rigorous definition of the time-dependent hole state (TDHS) using the time-domain generalization of extended Koopmans' theorem (TD-EKT), we derive the reduced ion density matrix within the MCTDHF framework, which inherently encodes the total and channel-resolved photoionization cross sections of $\mathrm{H_2O}$. The cross sections obtained are benchmarked against existing experimental and theoretical results, validating the TDHS formalism. Furthermore, by adjusting the phase delay and intensity ratio of a pair of orthogonally polarized attosecond pulses, we explore the ultrafast control of attosecond coherence between electronic states of $\mathrm{H_2O^+}$.

preprint2023arXiv

Real-time Equation-of-Motion Coupled-Cluster Cumulant Green's Function Method: Heterogeneous Parallel Implementation Based on the Tensor Algebra for Many-body Methods Infrastructure

We report the implementation of the real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green's function method [J. Chem. Phys. 152, 174113 (2020)] within the Tensor Algebra for Many-body Methods (TAMM) infrastructure. TAMM is a massively parallel heterogeneous tensor library designed for utilizing forthcoming exascale computing resources. The two-body electron repulsion matrix elements are Cholesky-decomposed, and we imposed spin-explicit forms of the various operators when evaluating the tensor contractions. Unlike our previous real algebra Tensor Contraction Engine (TCE) implementation, the TAMM implementation supports fully complex algebra. The RT-EOM-CC singles (S) and doubles (D) time-dependent amplitudes are propagated using a first-order Adams--Moulton method. This new implementation shows excellent scalability tested up to 500 GPUs using the Zn-porphyrin molecule with 655 basis functions, with parallel efficiencies above 90\% up to 400 GPUs. The TAMM RT-EOM-CCSD was used to study core photo-emission spectra in the formaldehyde and ethyl trifluoroacetate (ESCA) molecules. Simulations of the latter involve as many as 71 occupied and 649 virtual orbitals. The relative quasiparticle ionization energies and overall spectral functions agree well with available experimental results.

preprint2022arXiv

Real-time equation-of-motion CC cumulant and CC Green's function simulations of photoemission spectra of water and water dimer

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses results obtained with the real-time equation-of-motion CC cumulant approach (RT-EOM-CC), and CC Green's function (CCGF) approaches in applications to the water and water dimer molecules. We compare the ionization potentials obtained with these methods for the valence region with the results obtained with the CCSD(T) formulation as a difference of energies for N and N-1 electron systems. All methods show good agreement with each other. They also agree well with experiment, with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.

preprint2022arXiv

Time-dependent optimized coupled-cluster method with doubles and perturbative triples [TD-OCCD(T)] for first principles simulation of multielectron dynamics

We report the formulation of a new, cost-effective approximation method in the time-dependent optimized coupled-cluster (TD-OCC) framework [T. Sato et al., J. Chem. Phys. 148, 051101 (2018)] for first-principles simulations of multielectron dynamics in an intense laser field. The method, designated as TD-OCCD(T), is a time-dependent, orbital-optimized extension of the "gold-standard" CCSD(T) method in the ground-state electronic structure theory. The equations of motion for the orbital functions and the coupled-cluster amplitudes are derived based on the real-valued time-dependent variational principle using the fourth-order Lagrangian. The TD-OCCD(T) is size extensive and gauge invariant, and scales as O(N^7) with respect to the number of active orbitals N. The pilot application of the TD-OCCD(T) method to the strong-field ionization and high-order harmonic generation from a Kr atom is reported in comparison with the results of the previously developed methods, such as the time-dependent complete-active-space self-consistent field (TD-CASSCF), TD-OCC with double and triple excitations (TD-OCCDT), TD-OCC with double excitations (TD-OCCD), and the time-dependent Hartree-Fock (TDHF) methods.

preprint2020arXiv

Relativistic double-ionization equation-of-motion coupled-cluster method: Application to low-lying doubly ionized states

The article deals with the extension of the relativistic double-ionization equation-of-motion coupled-cluster (DI-EOMCC) method [H. Pathak et al. Phys. Rev. A 90, 010501(R) (2014)] for the molecular systems. The Dirac-Coulomb (DC) Hamiltonian with four-component spinors is considered to take care of the relativistic effects. The implemented method is employed to compute a few low-lying doubly ionized states of noble gas atoms (Ar, Kr, Xe, and Rn) and Cl 2 , Br 2 , HBr, and HI. Additionally, we presented results with two intermediate schemes in the four-component relativistic DI-EOMCC framework to understand the role of electron correlation. The computed double ionization spectra for the atomic systems are compared with the values from the non-relativistic DI-EOMCC method with spin-orbit coupling (SOC) [Z. Wang et al. J. Chem. Phys. 142, 144109 (2015)] and the values from the National Institute of Science and Technology (NIST) database. Our atomic results are found to be in good agreement with the NIST values. Further, the obtained results for the molecular systems agree well with the available experimental values.

preprint2020arXiv

Study of laser-driven multielectron dynamics of Ne atom using time-dependent optimized second-order many-body perturbation theory

We calculate the high-harmonic generation (HHG) spectra, strong-field ionization, and time-dependent dipole-moment of Ne using explicitly time-dependent optimized second-order many-body perturbation method (TD-OMP2) where both orbitals and amplitudes are time-dependent. We consider near-infrared (800 nm) and mid-infrared (1200 nm) laser pulses with very high intensities ($5\times10^{14}$, $8\times10^{14}$ , and $1\times10^{15}$ W/cm$^2$), required for strong-field experiments with the high-ionization potential (21.6 eV) atom. We compare the result of the TD-OMP2 method with the time-dependent complete-active-space self-consistent field method and the time-dependent Hartree-Fock method. Further, we report the implementation of the TD-CC2 method within the chosen active space, which is also a second-order approximation to the TD-CCSD method, and present results of time-dependent dipole-moment and HHG spectra with an intensity of $5\times10^{13}$ W/cm$^2$ at a wavelength of 800 nm. It is found that the TD-CC2 method is not stable in the case with a higher laser intensity, and it does not provide a gauge-invariant description of the physical properties, which makes TD-OMP2 a superior choice to reach out to larger chemical systems, especially for the study of strong-field dynamics. The obtained results indicate that the TD-OMP2 method shows moderate performance, overestimating the response of Ne, while TDHF underestimates it. Nevertheless, it is remarkable that stable computation of such highly nonlinear nonperturbative phenomena is possible within the framework of time-dependent perturbation method, by virtue of the nonperturbative inclusion of the laser-electron interaction and time-dependent optimization of orbitals.

preprint2020arXiv

Time-dependent optimized coupled-cluster method for multielectron dynamics II. A coupled electron-pair approximation

We report the implementation of a cost-effective approximation method within the framework of time-dependent optimized coupled-cluster (TD-OCC) method [J. Chem. Phys. 148, 051101 (2018)] for real-time simulations of intense laser-driven multielectron dynamics. The method, designated as TD-OCEPA0, is a time-dependent extension of the simplest version of the coupled-electron pair approximation with optimized orbitals [J. Chem. Phys. 139, 054104 (2013)]. It is size extensive, gauge invariant, and computationally much more efficient than the TD-OCC with double excitations (TD-OCCD). We employed this method to simulate the electron dynamics in Ne and Ar atoms exposed to intense near infrared laser pulses with various intensities. The computed results, including high-harmonic generation spectra and ionization yields, are compared with those of various other methods ranging from uncorrelated time-dependent Hartree-Fock (TDHF) to fully-correlated (within the active orbital space) time-dependent complete-active-space self-consistent-field (TD-CASSCF). The TD-OCEPA0 results show a good agreement with TD-CASSCF ones for moderate laser intensities. For higher intensities, however, TD-OCEPA0 tends to overestimate the correlation effect, as occasionally observed for CEPA0 in the ground-state correlation energy calculations.

preprint2020arXiv

Time-dependent optimized coupled-cluster method for multielectron dynamics III: A second-order many-body perturbation approximation

We report successful implementation of the time-dependent second-order many-body perturbation theory using optimized orthonormal orbital functions called time-dependent optimized second-order many-body perturbation theory [TD-OMP2] to reach out to relatively larger chemical systems for the study of intense-laser-driven multielectron dynamics. We apply this method to strong-field ionization and high-order harmonic generation (HHG) of Ar. The calculation results are benchmarked against ab initio time-dependent complete-active-space self-consistent field (TD-CASSCF), time-dependent optimized coupled-cluster double (TD-OCCD), and time-dependent Hartree-Fock (TDHF) methods, as well as a single active electron (SAE) model to explore the role of electron correlation.