Source author record

Himadri Pathak

Himadri Pathak appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

21 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.

preprint2016arXiv

Relativistic coupled-cluster study of RaF as a candidate for parity and time reversal violating interaction

We have employed both Z-vector method and the expectation value approach in the relativistic coupled-cluster framework to calculate the scalar-pseudoscalar (S-PS) P, T -odd interaction constant (W_s) and the effective electric field (Eeff) experienced by the unpaired electron in the ground electronic state of RaF. Further, the magnetic hyperfine structure constants of ^{223}Ra in RaF and ^{223}Ra+ are also calculated and compared with the experimental values wherever available to judge the extent of accuracy obtained in the employed methods. The outcome of our study reveals that the Z-vector method is superior than the expectation value approach in terms of accuracy obtained for the calculation of ground state property. The Z-vector calculation shows that RaF has a high E_eff (52.5 GV/cm) and W_s (141.2 kHz) which makes it a potential candidate for the eEDM experiment.

preprint2016arXiv

Relativistic equation-of-motion coupled-cluster method for the electron attachment problem

The article considers the successful implementation of relativistic equation-of-motion coupled cluster method for the electron attachment problem (EA-EOMCC) at the level of single- and double- excitation approximation. The implemented relativistic EA-EOMCC method is employed to calculate ionization potential values of alkali metal atoms (Li, Na, K, Rb, Cs, Fr) and the vertical electron affinity values of LiX (X = H, F, Cl, Br), NaY (Y = H, F, Cl) starting from their closed-shell configuration. Both four-component and exact two-component calculations are done for all the opted systems. Further, we have shown the effect of spin-orbit interaction considering the atomic systems. The results of our atomic calculations are compared with the values from the NIST database and the results are found to be very accurate (< 1 %).

preprint2016arXiv

Relativistic equation-of-motion coupled-cluster method using open-shell reference wavefunction

The open-shell reference relativistic equation-of-motion coupled-cluster method within its four-component description is successfully implemented with the consideration of single- and double- excitation approximation. The one-body and two-body matrix elements required for the correlation calculation are generated using Dirac-Coulomb Hamiltonian. As a first attempt, the implemented method is employed to calculate a few of the low-lying ionized states of heavy atomic (Ag, Cs, Au, Fr, Lr) and valence ionization potential of molecular (HgH, PbF) systems, where the effect of relativity does really matter to obtain highly accurate results. Not only the relativistic effect, but also the effect of electron correlation is crucial in these heavy atomic and molecular systems. To justify the fact, we have taken two further approximations in the four-component relativistic equation-of-motion framework to quantify how the effect of electron correlation plays a role in the calculated values at different level of the approximation. All these calculated results are compared with the available experimental data as well as with other theoretically calculated values to judge the extent of accuracy obtained in our calculations. A nice agreement is achieved for all the considered systems which validates the performance of the open-shell reference relativistic equation-of-motion coupled-cluster method.

preprint2016arXiv

Search for parity and time reversal violating effects in HgH: Relativistic coupled-cluster study

The high effective electric field ($E_\mathrm{eff}$) experienced by the unpaired electron in an atom or a molecule is one of the key ingredients in the success of electron electric dipole moment (eEDM) experiment and its precise calculation require a very accurate theory. We, therefore, employed the Z-vector method in the relativistic coupled-cluster framework and found that HgH has a very large $E_\mathrm{eff}$ value (123.2 GV/cm) which makes it a potential candidate for the next generation eEDM experiment. Our study also reveals that it has a large scalar-pseudoscalar ${\mathcal{P,T}}$-violating interaction constant, $W_\mathrm{s}$ = 284.2 kHz. To judge the accuracy of the obtained results we have calculated parallel and perpendicular magnetic HFS constants and compared with the available experimental values. The results of our calculation are found to be in nice agreement with the experimental values. Therefore, by looking at the HFS results we can say that both $E_\mathrm{eff}$ and $W_\mathrm{s}$ values are also very accurate. Further, We have derived the relationship between these quantities and the ratio which will help to get model independent value of eEDM and S-PS interaction constant.

preprint2015arXiv

A Relativistic Equation-of-motion Coupled-Cluster Investigation of the Trends of Single and Double Ionization Potentials in the He and Be Isoelectronic Systems

We employ four-component spinor relativistic equation-of-motion coupled-cluster (EOMCC) method within the single- and double- excitation approximation to calculate the single ionization potentials (IPs) and double ionization potentials (DIPs) of the He and Be isoelectronic sequences up to Ne. The obtained results are compared with the available results from the National Institute of Standards and Technology (NIST) database to test the performance of the EOMCC method. We also present intermediate results at different level of approximations in the EOMCC framework to gain insight of the effect of electron correlation. Furthermore, we investigate the dependence of the IPs and DIPs of these ions on the ionic charge and observe that these follow parabolic trends. Comparison between the trends of IPs and DIPs in both the classes of considered systems are categorically demonstrated.

preprint2015arXiv

Calculation of P,T-odd interaction constant of PbF using Z-vector method in the relativistic coupled-cluster framework

The effective electric field experienced by the unpaired electron in the ground state of PbF, which is a potential candidate in the search of electron electric dipole moment due to some special characteristics, is calculated using Z-vector method in the coupled cluster single- and double- excitation approximation with four component Dirac spinor. This is an important quantity to set the upper bound limit of the electron electric dipole moment. Further, we have calculated molecular dipole moment and parallel magnetic hyperfine structure constant (A$_\|$) of $^{207}$Pb in PbF to test the accuracy of the wave function obtained in the Z-vector method. The outcome of our calculations clearly suggests that the core electrons have significant contribution to the "atom in compound (AIC)" properties.

preprint2015arXiv

Implementation of Z-vector method in the relativistic coupled cluster framework to calculate first order energy derivatives: Application to SrF molecule

The molecular dipole moment and magnetic hyperfine structure constant demand an accurate wavefunction far from the nucleus and in near nuclear region, respectively. We, therefore, employ the so-called Z-vector method in the domain of relativistic coupled cluster theory to calculate the first order property of molecular systems in their open-shell ground state configuration. The implemented method is applied to calculate molecular dipole moment and parallel component of the magnetic hyperfine structure constant of SrF molecule. The results of our calculation are compared with the experimental and other available theoretically calculated values. We are successful in achieving good accordance with the experimental results. The result of our calculation of molecular dipole moment is in the accuracy of ~? 0.5 %, which is clearly an improvement over the previous calculation based on the expectation value method in the four component coupled cluster framework [V. S. Prasannaa et al, Phys. Rev. A 90, 052507 (2014)] and it is the best calculated value till date. Thus, it can be inferred that the Z vector method can provide an accurate wavefunction in both near and far nuclear region, which is evident from our calculated results.

preprint2015arXiv

Relativistic extended coupled cluster method for magnetic hyperfine structure constant

This article deals with the general implementation of 4-component spinor relativistic extended coupled cluster (ECC) method to calculate first order property of atoms and molecules in their open-shell ground state configuration. The implemented relativistic ECC is employed to calculate hyperfine structure (HFS) constant of alkali metals (Li, Na, K, Rb and Cs), singly charged alkaline earth metal atoms (Be+, Mg+, Ca+ and Sr+) and molecules (BeH, MgF and CaH). We have compared our ECC results with the calculations based on restricted active space configuration interaction (RAS-CI) method. Our results are in better agreement with the available experimental values than those of the RAS-CI values.

preprint2014arXiv

Relativistic Equation of Motion Coupled-Cluster Method: Application to the closed-shell atomic systems

We report our successful implementation of the full fledged relativistic equation of motion coupled cluster (EOMCC) method. This method is employed to compute the principal ionization potentials (IPs) of closed-shell rare gas atoms, He-like ions, Be-like ions along with Na+, Al+, K+, Be, and Mg. Four component Dirac spinors are used in the calculations and the one and two electron integrals are evaluated using the Dirac Coulomb Hamiltonian. Our results are in excellent agreement with those available measurements, which are taken from the National Institute of Science and Technology database (NIST). We also present results using the second order many-body perturbation theory (MBPT(2)) and random phase approximation (RPA) in the EOMCC framework. These results are compared with those of EOMCC at the level of single and double excitations in order to assess the role of the electron correlation effects in the intermediate schemes considered in our calculations .

preprint2014arXiv

Relativistic equation-of-motion coupled-cluster method for the double ionization potentials of the closed-shell atoms

We report the general implementation of the relativistic equation-of-motion coupled-cluster method to calculate double ionization spectra (DI-EOMCC) of atomic and molecular systems. As a first application, this method is employed to calculate the principal valence double ionization potential values of He and alkaline earth metal (Be, Mg, Ca, Sr and Ba) atoms. Our results are compared with the results available from the national institute of standards and technology (NIST) database and other ab initio calculations. We have achieved an accuracy of ~ 0.1%, which is an improvement over the first principles T-matrix calculations [J. Chem. Phys. 123, 144112 (2005)]. We also present results using the second-order many-body perturbation theory and the random- phase approximation in the equation-of-motion framework and these results are compared with the DI-EOMCC results.

preprint2014arXiv

Relativistic equation-of-motion coupled-cluster method for the ionization problem: Application to molecules

We report the implementation of 4-component spinor relativistic equation-of-motion coupledcluster method within the single- and double- excitation approximation to calculate ionization potential (EOM-CCSD) of molecules. We have applied this method to calculate vertical ionization potentials of the molecules, XH(X=F, Cl, Br, I) along with Cl2 and Br2 in their closed-shell configuration. We have also presented intermediate results using 2-nd order many-body perturbation theory level in the EOM framework (EOM-MBPT(2)) to understand the role of electron correlation. All the calculated values are compared with the available experimental results. Our results are found to be in well agreement with the sophisticated experiments and relative deviation of less than 1% achieved for all the considered systems.