No description available.

Electronic structure methods aim to accurately describe the behaviour of the electrons in molecules and materials. To be applicable to arbitrary systems, these methods cannot depend on observations of specific chemical phenomena and must be derived solely from the fundamental physical constants and laws that govern all electrons. Such methods are called ab initio methods. Ab initio methods directly solve the electronic Schrödinger equation to obtain the electronic energy and wavefunction. For more than one electron, solving the electronic Schrödinger equation is impossible, so it is imperative to develop approximate methods that cater to the needs of their users, which can vary depending on the chemical systems under study, the available computational resources and time, and the desired level of accuracy. The most accessible ab initio approaches, including Hartree-Fock methods and Kohn-Sham density functional theory methods, assume that only one electronic configuration is needed to describe the system. While these single-reference methods are successful when describing systems where a single electron configuration dominates, like most closed-shell ground-state organic molecules in their equilibrium geometries, single-reference methods are unreliable for molecules in nonequilibrium geometries (e.g., transition states) and molecules containing unpaired electrons (e.g., transition metal complexes and radicals). For these types of multireference systems, accurate results can only be obtained if multiple electronic configurations are accounted for. Wavefunctions that incorporate many electronic configurations are called multideterminant wavefunctions. This thesis presents a systematic approach to developing multideterminant wavefunctions. First, we establish a framework that outlines the structural components of a multideterminant wavefunction and propose several novel wavefunction ansätze. Then, we present a software package that is designed to aid the development of new wavefunctions and algorithms. Using this approach, we develop an algorithm for evaluating the geminal wavefunctions, a class of multideterminant wavefunctions that are expressed with respect to electron pairs. Finally, we explore using machine learning to solve the Schrödinger equation by presenting a neural network wavefunction ansatz and optimizing its parameters using stochastic gradient descent. / Thesis / Doctor of Science (PhD)

<p dir="ltr">This thesis explores the intersection of quantum interference, entanglement, and quantum algorithms in the context of chemical reactions. The initial exploration delves into the constructive quantum interference in the photoassociation reaction of a 87Rb Bose Einstein condensate (BEC), where a coherent superposition of multiple bare spin states is achieved and it’s impact on photo-association (PA) was studied. Employing a quantum processor, the study illustrates that interferences can function as a resource for coherent control in photochemical reactions, presenting a universally applicable framework relevant to a spectrum of ultracold chemical reactions. The subsequent inquiry scrutinizes the entanglement dynamics between the spin and momentum degrees of freedom in an optically confined BEC of 87Rb atoms, induced by Raman and RF fields. Significantly, this study unveils substantial spin momentum entanglement under specific experimental conditions, indicating potential applications in the realm of quantum information processing. Finally, the third study advances a quantum algorithm for the computation of scattering matrix elements in chemical reactions, adeptly navigating the complexities of quantum interactions. This algorithm, rooted in the time-dependent method and Möller operator formulation, is applied to scenarios such as 1D semi-infinite square well potentials and co-linear hydrogen exchange reactions, showcasing its potential to enhance our comprehension of intricate quantum interactions within chemical systems.</p>

Conventional correlation methods suffer from the slow convergence of electron correlation energies with respect to the size of orbital expansions. This problem is due to the fact that orbital products alone cannot describe the behavior of the exact wave function at short inter-electronic distances. Explicitly correlated methods overcome this basis set problem by including the inter-electronic distances (rij) explicitly in wave function expansions. Here, the origin of the basis set problem of conventional wave function methods is reviewed, and a short history of explicitly correlated methods is presented. The F12 methods are the focus herein, as they are the most practical explicitly correlated methods to date. Moreover, some of the key developments in modern F12 technology, which have significantly improved the efficiency and accuracy of these methods, are also reviewed. In this work, the extension of the perturbative coupled-cluster F12 method, CCSD(T)F12, developed in our group for the treatment of high-spin open-shell molecules (J. Zhang and E. F. Valeev, J. Chem. Theory Comput., 2012, 8, 3175.), is also documented. Its performance is assessed for accurate prediction of chemical reactivity. The reference data include reaction barrier heights, electronic reaction energies, atomization energies, and enthalpies of formation from the following sources: (1) the DBH24/08 database of 22 reaction barriers (Truhlar et al., J. Chem. Theory Comput., 2007, 3, 569.), (2) the HJO12 set of isogyric reaction energies (Helgaker et al., Modern Electronic Structure Theory, Wiley, Chichester, first ed., 2000.), and (3) the HEAT set of atomization energies and heats of formation (Stanton et al., J. Chem. Phys., 2004, 121, 11599.). Two types of analyses were performed, which target the two distinct uses of explicitly correlated CCSD(T) models: as a replacement for the basis-set-extrapolated CCSD(T) in highly accurate composite methods like HEAT and as a distinct model chemistry for standalone applications. Hence, (1) the basis set error of each component of the CCSD(T)F12 contribution to the chemical energy difference in question and (2) the total error of the CCSD(T)F12 model chemistry relative to the benchmark values are analyzed in detail. Two basis set families were utilized in the calculations: the standard aug-cc-p(C)VXZ (X = D, T, Q) basis sets for the conventional correlation methods and the cc-p(C)VXZ-F12 (X = D, T, Q) basis sets of Peterson and co-workers that are specifically designed for explicitly correlated methods. The conclusion is that the performance of the two families for CCSD correlation contributions (which are the only components affected by the explicitly correlated terms in our formulation) are nearly identical with triple- and quadruple-ζ quality basis sets, with some differences at the double-ζ level. Chemical accuracy (~4.18 kJ/mol) for reaction barrier heights, electronic reaction energies, atomization energies, and enthalpies of formation is attained, on average, with the aug-cc-pVDZ, aug-cc-pVTZ, cc- pCVTZ-F12/aug-cc-pCVTZ, and cc-pCVDZ-F12 basis sets, respectively, at the CCSD(T)F12 level of theory. The corresponding mean unsigned errors are 1.72 kJ/ mol, 1.5 kJ/mol, ~ 2 kJ/mol, and 2.17 kJ/mol, and the corresponding maximum unsigned errors are 4.44 kJ/mol, 3.6 kJ/mol, ~ 5 kJ/mol, and 5.75 kJ/mol. In addition to accurate energy calculations, our studies were extended to the computation of molecular properties with the MP2-F12 method, and its performance was assessed for prediction of the electric dipole and quadrupole moments of the BH, CO, H2O, and HF molecules (J. Zhang and E. F. Valeev, in preparation for submission). First, various MP2- F12 contributions to the electric dipole and quadrupole moments were analyzed. It was found that the unrelaxed one-electron density contribution is much larger than the orbital response contribution in the CABS singles correction, while both contributions are important in the MP2 correlation contribution. In contrast, the majority of the F12 correction originates from orbital response effects. In the calculations, the two basis set families, the aug-cc-pVXZ (X = D, T, Q) and cc-pVXZ-F12 (X = D, T, Q) basis sets, were also employed. The two basis set series show noticeably different performances at the double-ζ level, though the difference is smaller at triple- and quadruple-ζ levels. In general, the F12 calculations with the aug-cc- pVXZ series give better results than those with the cc-pVXZ-F12 family. In addition, the contribution of the coupling from the MP2 and F12 corrections was investigated. Although the computational cost of the F12 calculations can be significantly reduced by neglecting the coupling terms, this does increase the errors in most cases. With the MP2-F12C/aug-cc-pVDZ calculations, dipole moments close to the basis set limits can be obtained; the errors are around 0.001 a.u. For quadrupole moments, the MP2-F12C/aug-cc-pVTZ calculations can accurately approximate the MP2 basis set limits (within 0.001 a.u.). / Ph. D.

Accurate correlated electronic structure methods involve a significant amount of computations and can be only employed to small molecular systems. For example, the coupled-cluster singles, doubles, and perturbative triples model (CCSD(T)), which is known as the ``gold standard" of quantum chemistry for its accuracy, usually can treat molecules with 20-30 atoms. To extend the reach of accurate correlated electronic structure methods to larger molecular systems, we work towards two directions: parallel computing and reduced-cost/scaling approaches. Parallel computing can utilize more computational resources to handle systems that demand more substantial computational efforts. Reduced-cost/scaling approaches, which introduce approximations to the existing electronic structure methods, can significantly reduce the amount of computation and storage requirements. In this work, we introduce a new distributed-memory massively parallel implementation of standard and explicitly correlated (F12) coupled-cluster singles and doubles (CCSD) with canonical bigO{N^6} computational complexity ( C. Peng, J. A. Calvin, F. Pavov{s}evi'c, J. Zhang, and E. F. Valeev, textit{J. Phys. Chem. A} 2016, textbf{120}, 10231.), based on the TiledArray tensor framework. Excellent strong scaling is demonstrated on a multi-core shared-memory computer, a commodity distributed-memory computer, and a national-scale supercomputer. We also present a distributed-memory implementation of the density-fitting (DF) based CCSD(T) method. (C. Peng, J. A. Calvin, and E. F. Valeev, textit{in preparation for submission}) An improved parallel DF-CCSD is presented utilizing lazy evaluation for tensors with more than two unoccupied indices, which makes the DF-CCSD storage requirements always smaller than those of the non-iterative triples correction (T). Excellent strong scaling is observed on both shared-memory and distributed-memory computers equipped with conventional Intel Xeon processors and the Intel Xeon Phi (Knights Landing) processors. With the new implementation, the CCSD(T) energies can be evaluated for systems containing 200 electrons and 1000 basis functions in a few days using a small size commodity cluster, with even more massive computations possible on leadership-class computing resources. The inclusion of F12 correction to the CCSD(T) method makes it converge to basis set limit much more rapidly. The large-scale parallel explicitly correlated coupled-cluster program makes the accurate estimation of the coupled-cluster basis set limit for molecules with 20 or more atoms a routine. Thus, it can be used rigorously to test the emerging reduced-scaling coupled-cluster approaches. Moreover, we extend the pair natural orbital (PNO) approach to excited states through the equation-of-motion coupled cluster singles and doubles (EOM-CCSD) method. (C. Peng, M. C. Clement, and E. F. Valeev, textit{submitted}) We simulate the PNO-EOM-CCSD method using an existing massively parallel canonical EOM-CCSD program. We propose the use of state-averaged PNOs, which are generated from the average of the pair density of excited states, to span the PNO space of all the excited states. The doubles amplitudes in the CIS(D) method are used to compute the state-averaged pair density of excited states. The issue of incorrect states in the state-averaged pair density, caused by an energy reordering of excited states between the CIS(D) and EOM-CCSD, is resolved by simply computing more states than desired. We find that with a truncation threshold of $10^{-7}$, the truncation error for the excitation energy is already below 0.02 eV for the systems tested, while the average number of PNOs is reduced to 50-70 per pair. The accuracy of the PNO-EOM-CCSD method on local, Rydberg and charge transfer states is also investigated. / Ph. D.

<p dir="ltr">Theoretical development and implementation of Partition Density Functional Theory, a quantum density embedding framework for electronic structure simulation.</p>

First, I introduce the concept of heavy fermion systems and discuss the ease of tuning their properties with external parameters including pressure, chemical doping, and magnetic fields to induce new quantum states such as unconventional superconductivity. I then delve into the limited use of dimensionality as a tuning knob for quantum criticality and highlight the new possibilities available if heavy fermion behavior is discovered in the single-layer limit. Chapter 1 establishes the van der Waals material, CeSiI, as a heavy fermion system and is the first material where heavy fermion behavior exists down to the few-layer limit. The chapter further explores the bulk magnetic properties and electronic structure of CeSiI at high magnetic fields. The quasi-two-dimensional electronic character of CeSiI leads to anisotropic hybridization between local moments and conduction electrons, a phenomenon previously only realized in theoretical calculations. With the heavy fermion properties of CeSiI established, Chapter 2 investigates the effects of pressure and La-doping on CeSiI, aiming to push it from the antiferromagnetic region of the Doniach phase diagram towards a quantum critical point. Preliminary evidence suggests that CeSiI is too distant from quantum criticality. Instead, La-doping is utilized to explore single-ion Kondo physics at the dilute Ce limit in CeSiI. Additionally, CeGaI, with a crystal structure similar to CeSiI, is examined. Although no Kondo physics is observed, the magnetic and electronic properties remain coupled to each other. Chapter 3 delves into a separate project focusing on the study of polymers composed of perylene diimide and various organic linkers. It explores how the structure of the polymer influences its pseudocapacitance properties. The chapter demonstrates the significance of contortion in device performance, aiming to provide insights for future endeavors in developing environmentally friendly energy storage systems.

Quantum chemistry methodologies can be used to address a wide variety of chemical problems. Key to the success of quantum chemistry methodologies, however, is the selection of suitable methodologies for specific problems of interest, which often requires significant assessment. To gauge a number of methodologies, the utility of density functionals (BLYP, B97D, TPSS, M06L, PBE0, B3LYP, M06, and TPSSh) in predicting reaction energetics was examined for model studies of C-O bond activation of methoxyethane and methanol. These species provide excellent representative examples of lignin degradation via C-O bond cleavage. PBE0, which performed better than other considered DFT functionals, was used to investigate late 3d (Fe, Co, and Ni), 4d (Ru, Rh, and Pd), and 5d (Re, Os, and Ir) transition metal atom mediated Cβ -O bond activation of the β–O–4 linkage of lignin. Additionally, the impact of the choice of DFT functionals, basis sets, implicit solvation models, and layered quantum chemical methods (i.e., ONIOM, Our Own N-layered Integrated molecular Orbital and molecular Mechanics) was investigated for the prediction of pKa for a set of Ni-group metal hydrides (M = Ni, Pd, and Pt) in acetonitrile. These investigations have provided insight about the utility of a number of theoretical methods in the computation of thermodynamic properties of transition metal hydrides in solution. As single reference wavefunction methods commonly perform poorly in describing molecular systems that involve bond-breaking and forming or electronic near-degeneracies and are typically best described with computationally costly multireference wavefunction-based methods, it is imperative to a priori analyze the multireference character for molecular systems so that the proper methodology choice is applied. In this work, diagnostic criteria for assessing the multireference character of 4d transition metal-containing molecules was investigated. Four diagnostics were considered in this work, including the weight of the leading configuration of the CASSCF wavefunction, C02; T1, the Frobenius norm of the coupled cluster amplitude vector related to single excitations and D1, the matrix norm of the coupled cluster amplitude vector arising from coupled cluster calculations; and the percent total atomization energy, %TAE. This work demonstrated the need to have different diagnostic criteria for 4d molecules than for main group molecules.

In the first two chapters, we provide a new way to think about the Green’s function-basedGW approximation and Bethe-Salpeter equation (BSE). The former is the most popular beyond-mean-field method for band structures of solids and an increasingly popular one for ionization potentials and electron affinities of molecules. The latter is widely used to compute neutral excitation energies and spectra for solids as well as, increasingly, molecules. Inspired by quantum chemistry approaches, we obtain a computational scaling reduction and avoid approximating certain dynamical quantities. The new formalism suggests further improvements to the GW and BSE methods. In chapters four and five, we derive and test a cheap, approximate version of the GW and BSE for large molecules and then extend the strategy to periodic systems. In chapter six, we assess another Green’s function-based method, the constrained random phase approximation with exact diagonalization, usually applied to solids. This method allows one to treat electron correlation within an active space of important orbitals while also including some of the external orbital space effects. In chapters seven and eight, we implement the BSE in the PySCF software package for periodic systems using Gaussian density fitting and then apply it to a challenging system, the superatomic solid Re₆Se₈Cl₂.

京都大学 / 新制・課程博士 / 博士(工学) / 甲第25306号 / 工博第5265号 / 新制||工||2001(附属図書館) / 京都大学大学院工学研究科分子工学専攻 / (主査)教授 佐藤 啓文, 教授 佐藤 徹, 教授 作花 哲夫 / 学位規則第4条第1項該当 / Doctor of Agricultural Science / Kyoto University / DGAM