Search

We present our successful implementation of the quantum electrodynamics coupled-cluster method with single and double excitations (QED-CCSD) for electronic and bosonic amplitudes, covering both individual and mixed excitation processes within the ExaChem program package, which relies on the Tensor Algebra for Many-body Methods (TAMM) infrastructure. TAMM is a parallel heterogeneous tensor library designed for utilizing modern computing platforms, from laptops to leadership-class computing resources. This developed computational framework extends the traditional CCSD method to incorporate the intricate interplay between electronic and bosonic degrees of freedom, providing a comprehensive description of quantum phenomena. We discuss theoretical foundations, algorithmic details, and numerical benchmarks to demonstrate how the integration of bosonic degrees of freedom alters the electronic ground state. The interactions between electrons and photons within an optical cavity are modeled using the Pauli-Fierz Hamiltonian within the dipole approximation in the length gauge. The integration of QED effects within the CCSD framework contributes to a more accurate and versatile model for simulating complex quantum systems, thereby opening avenues for a better understanding, prediction, and manipulation of various physical phenomena.

In this paper, we analyze the properties of the recently proposed real-time equation-of-motion coupled-cluster (RT-EOM-CC) cumulant Green's function approach [J. Chem. Phys. 2020, 152, 174113]. We specifically focus on identifying the limitations of the original time-dependent coupled cluster (TDCC) ansatz and propose an enhanced double TDCC ansatz ensuring the exactness in the expansion limit. Additionally, we introduce a practical cluster-analysis-based approach for characterizing the peaks in the computed spectral function from the RT-EOM-CC cumulant Green's function approach, which is particularly useful for the assignments of satellite peaks when many-body effects dominate the spectra. Our preliminary numerical tests focus on reproducing, approximating, and characterizing the exact impurity Green's function of the three-site and four-site single impurity Anderson models using the RT-EOM-CC cumulant Green's function approach. The numerical tests allow us to have a direct comparison between the RT-EOM-CC cumulant Green's function approach and other Green's function approaches in the numerical exact limit.

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.

We review time-dependent multiconfiguration self-consistent-field (TD-MCSCF) method and time-dependent optimized coupled-cluster (TD-OCC) method for first-principles simulations of high-field phenomena such as tunneling ionization and high-order harmonic generation in atoms and molecules irradiated by a strong laser field. These methods provide a flexible and systematically improvable description of the multielectron dynamics by expressing the all-electron wavefunction by configuration interaction expansion or coupled-cluster expansion, using time-dependent one-electron orbital functions. The time-dependent variational principle plays a key role to derive these methods satisfying gauge invariance and Ehrenfest theorem. The real-time/real-space implementation with an absorbing boundary condition enables the simulation of high-field processes involving multiple excitation and ionization. We present a detailed, comprehensive discussion of such features of TD-MCSCF and TD-OCC methods.

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.

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.

We present a cost-effective treatment of the triple excitation amplitudes in the time-dependent optimized coupled-cluster (TD-OCC) framework called TD-OCCDT(4) for studying intense laser-driven multielectron dynamics. It considers triple excitation amplitudes correct up to fourth-order in many-body perturbation theory and achieves a computational scaling of O(N7), with N being the number of active orbital functions. This method is applied to the electron dynamics in Ne and Ar atoms exposed to an intense near-infrared laser pulse with various intensities. We benchmark our results against the time-dependent complete-active-space self-consistent field (TD-CASSCF), time-dependent optimized coupled-cluster with double and triple excitations (TD-OCCDT), time-dependent optimized coupled-cluster with double excitations (TD-OCCD), and the time-dependent Hartree-Fock (TDHF) methods to understand how this approximate scheme performs in describing nonperturbatively nonlinear phenomena, such as field-induced ionization and high-harmonic generation. We find that the TD-OCCDT(4) method performs equally well as the TD-OCCDT method, almost perfectly reproducing the results of fully-correlated TD-CASSCF with a more favorable computational scaling., Comment: arXiv admin note: text overlap with arXiv:2001.02206

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., Comment: 19 pages, 13 figures

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.

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., Comment: 8 page

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.

Nowadays, companies such as Amazon, Alibaba, and even pizza chains are pushing forward to use drones, also called UAVs (Unmanned Aerial Vehicles), for service provision, such as package and food delivery. As governments intend to use these immense economic benefits that UAVs have to offer, urban planners are moving forward to incorporate so-called UAV flight zones and UAV highways in their smart city designs. However, the high-speed mobility and behavior dynamics of UAVs need to be monitored to detect and, subsequently, to deal with intruders, rogue drones, and UAVs with a malicious intent. This paper proposes a UAV defense system for the purpose of intercepting and escorting a malicious UAV outside the flight zone. The proposed UAV defense system consists of a defense UAV swarm, which is capable to self-organize its defense formation in the event of intruder detection, and chase the malicious UAV as a networked swarm. Modular design principles have been used for our fully localized approach. We developed an innovative auto-balanced clustering process to realize the intercept- and capture-formation. As it turned out, the resulting networked defense UAV swarm is resilient against communication losses. Finally, a prototype UAV simulator has been implemented. Through extensive simulations, we show the feasibility and performance of our approach., Comment: IEEE Conference on Local Computer Networks (LCN), 2017

Time-dependent coupled-cluster method with time-varying orbital functions, called time-dependent optimized coupled- cluster (TD-OCC) method, is formulated for multielectron dynamics in an intense laser field. We have successfully derived equations of motion for CC amplitudes and orthonormal orbital functions based on the real action functional, and implemented the method including double excitations (TD-OCCD) and double and triple excitations (TD-OCCDT) within the optimized active orbitals. The present method is size extensive and gauge invariant, a polynomial cost-scaling alternative to the time-dependent multiconfiguration self-consistent-field method. The first application of the TD-OCC method to intense-laser driven correlated electron dynamics in Ar atom is reported.

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., Comment: 8 pages, 3 figures

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., Comment: 10 pages

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., Comment: main text : 5 pages

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., Comment: 6 pages, 1 figure

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 %)., Comment: 26 Pages, 3 figures, 6 Tables. Comments are welcome

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., Comment: 5 pages, 1 figure, Comments are welcome

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., Comment: 14 pages, 6 figures

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., Comment: Accepted in Physical Review A

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., Comment: 8 pages, 3 figures

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., Comment: 5 pages, 4 figure

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 .