Excitation configuration analysis for divide-and-conquer excited-state calculation method using dynamical polarizability

The authors previously developed a divide-and-conquer (DC)-based non-local excited-state calculation method for large systems using dynamical polarizability [Nakai and Yoshikawa, J. Chem. Phys. 146, 124123 (2017)]. This method evaluates the excitation energies and oscillator strengths using information on the dynamical polarizability poles. This article proposes a novel analysis of the previously developed method to obtain further configuration information on excited states, including excitation and de-excitation coefficients of each excitation configuration. Numerical applications to simple molecules, such as ethylene, hydrogen molecule, ammonia, and pyridazine, confirmed that the proposed analysis could accurately reproduce the excitation and de-excitation coefficients. The combination with the DC scheme enables both the local and non-local excited states of large systems with an excited nature to be treated.

Toward Accurate Prediction of Ion Mobility in Organic Semiconductors by Atomistic Simulation

A multiscale scheme (MLMS: Multi-Level Multi-Scale) to predict the ion mobility (μ) of amorphous organic semiconductors is proposed, which was successfully applied to the hole mobility predictions of 14 organic systems. An inverse relationship between μ and reorganization energy is observed due to local polaronic distortions. Another moderate inverse correlation between μ and distribution of site energy change exists, representing the effects of geometric flexibility. The former and the latter represent the intramolecular and intermolecular geometric effects, respectively. In addition, a linear correlation between transfer coupling and μ is observed, showing the importance of orbital overlaps between monomers. Especially, the highest hole mobility of C6-2TTN is due to its large transfer coupling. On the other hand, another high hole mobility of CBP turned out to come from the high first neighbor density (ρ_FND) of its first self-solvation, emphasizing the proper description of amorphous structural configurations with a sufficiently large number of monomers. In general, systems with either unusually high transfer coupling or high first neighbor density can potentially have high μ regardless of geometric effects. Especially, the newly suggested design parameter, ρ_FND, is pointing to a new direction as opposed to the traditional π-conjugation strategy.

Divide-and-Conquer Linear-Scaling Quantum Chemical Computations

Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.

Local electronic structure analysis by ab initio elongation method: A benchmark using DNA block polymers

The ab initio elongation (ELG) method based on a polymerization concept is a feasible way to perform linear-scaling electronic structure calculations for huge aperiodic molecules while maintaining computational accuracy. In the method, the electronic structures are sequentially elongated by repeating (1) the conversion of canonical molecular orbitals (CMOs) to region-localized MOs (RLMOs), that is, active RLMOs localized onto a region close to an attacking monomer or frozen RLMOs localized onto the remaining region, and the subsequent (2) partial self-consistent-field calculations for an interaction space composed of the active RLMOs and the attacking monomer. For each ELG process, one can obtain local CMOs for the interaction space and the corresponding local orbital energies. Local site information, such as the local highest-occupied/lowest-unoccupied MOs, can be acquired with linear-scaling efficiency by correctly including electronic effects from the frozen region. In this study, we performed a local electronic structure analysis using the ELG method for various DNA block polymers with different sequential patterns. This benchmark aimed to confirm the effectiveness of the method toward the efficient detection of a singular local electronic structure in unknown systems as a future practical application. We discussed the high-throughput efficiency of our method and proposed a strategy to detect singular electronic structures by combining with a machine learning technique.

Ab initio multi-level layered elongation method and its application to local interaction analysis between DNA bulge and ligand molecules

A multi-level layered elongation method was developed for efficiently analyzing the electronic states of local structures in large bio/nano-systems at the full ab initio level of theory. The original elongation method developed during the last three decades in our group has focused on the system in one direction from one terminal to the other terminal to sequentially construct the electronic states of a polymer, called a theoretical synthesis of polymers. In this study, an important region termed the central (C) part is targeted in a large polymer and the remainder are terminal (T) parts. The electronic structures along with polymer elongation are calculated repeatedly from both end T parts to the C central part at the same time. The important C part is treated with large basis sets (high level) and the other regions are treated with small basis sets (low level) in the ab initio theoretical framework. The electronic structures besides the C part can be reused for other systems with different structures at the C part, which renders the method computationally efficient. This multi-level layered elongation method was applied to the investigation on DNA single bulge recognition of small molecules (ligands). The reliability and validity of our approach were examined in comparison with the results obtained by direct calculations using a conventional quantum chemical method for the entire system. Furthermore, stabilization energies by the formation of the complex of bulge DNA and a ligand were estimated with basis set superposition error corrections incorporated into the elongation method.

iOI: An Iterative Orbital Interaction Approach for Solving the Self-Consistent Field Problem

An iterative orbital interaction (iOI) approach is proposed to solve, in a bottom-up fashion, the self-consistent field problem in quantum chemistry. While it belongs grossly to the family of fragment-based quantum chemical methods, iOI is distinctive in that (1) it divides and conquers not only the energy but also the wave function and that (2) the subsystem sizes are automatically determined by successively merging neighboring small subsystems until they are just enough for converging the wave function to a given accuracy. Orthonormal occupied and virtual localized molecular orbitals are obtained in a natural manner, which can be used for all post-SCF purposes.

Quantum chemistry-machine learning approach for predicting and elucidating molecular hyperpolarizability: Application to [2.2]paracyclophane-containing push-pull polymers

Nonlinear optical properties of organic chromophores are of great interest in diverse photonic and optoelectronic applications. To elucidate general trends in the behaviors of molecules, large amounts of data are required. Therefore, both an accurate and a rapid computational approach can significantly promote the theoretical design of molecules. In this work, we combined quantum chemistry and machine learning (ML) to study the first hyperpolarizability (β) in [2.2]paracyclophane-containing push-pull compounds with various terminal donor/acceptor pairs and molecular lengths. To generate reference β values for ML, the ab initio elongation finite-field method was used, allowing us to treat long polymer chains with linear scale efficiency and high computational accuracy. A neural network (NN) model was built for β prediction, and the relevant molecular descriptors were selected using a genetic algorithm. The established NN model accurately reproduced the β values (R² > 0.99) of long molecules based on the input quantum chemical properties (dipole moment, frontier molecular orbitals, etc.) of only the shortest systems and additional information about the actual system length. To obtain general trends in molecular descriptor-target property relationships learned by the NN, three approaches for explaining the ML decisions (i.e., partial dependence, accumulated local effects, and permutation feature importance) were used. The effect of donor/acceptor alternation on β in the studied systems was examined. The asymmetric extension of molecular regions end-capped with donors and acceptors produced unequal β responses. The results revealed how the electronic properties originating from the nature of substituents on the microscale controlled the magnitude of β according to the NN approximation. The applied approach facilitates the conceptual discoveries in chemistry by using ML to both (i) efficiently generate data and (ii) provide a source of information about causal correlations among system properties.

Computational Methods for Biochemical Simulations Implemented in GAMESS

Computational methods for modeling biochemical processes implemented in GAMESS package are reviewed; in particular, quantum mechanics combined with molecular mechanics (QM/MM), semi-empirical, and fragmentation approaches. A detailed summary of capabilities is provided for the QM/MM implementation in QuanPol program and the fragment molecular orbital (FMO) method. Molecular modeling and visualization packages useful for biochemical simulations with GAMESS are described. GAMESS capabilities with corresponding references are tabulated for reader's convenience.

Elongation method with intermediate mechanical and electrostatic embedding for geometry optimizations of polymers

The elongation method with intermediate mechanical and electrostatic embedding (ELG-IMEE) is proposed. The electrostatic embedding uses atomic charges generated by a charge sensitivity analysis (CSA) method and parameterized for three different population analyses, namely, the Merz-Singh-Kollman scheme, the charge model 5, and the atomic polar tensor. The obtained CSA models were tested on two model systems. Test calculations show that the electrostatic embedding provides several times of decrease in the difference of energies of testing and reference calculations in comparison with the conventional elongation approach (ELG). The mechanical embedding is implemented in a combination of the conventional elongation method and the ONIOM approach. Moreover, it was demonstrated that the geometry optimization with the ELG-IMEE reduces the errors in the optimized structures by about one order in root-mean-square deviation, when compared to ELG.

Solving the Schrödinger equation of atoms and molecules: Chemical-formula theory, free-complement chemical-formula theory, and intermediate variational theory

Chemistry is governed by the principle of quantum mechanics as expressed by the Schrödinger equation (SE) and Dirac equation (DE). The exact general theory for solving these fundamental equations is therefore a key for formulating accurately predictive theory in chemical science. The free-complement (FC) theory for solving the SE of atoms and molecules proposed by one of the authors is such a general theory. On the other hand, the working theory most widely used in chemistry is the chemical formula that refers to the molecular structural formula and chemical reaction formula, collectively. There, the central concepts are the local atomic concept, transferability, and from-atoms-to-molecule concept. Since the chemical formula is the most successful working theory in chemistry ever existed, we formulate our FC theory to have the structure reflecting the chemical formula. Our basic postulate is that as far as the SE is the principle of chemistry, its solutions for chemistry should have the structure that can be related to the chemical formulas. So, in this paper, we first formulate a theory that designs the wave function to reflect the structure of the chemical formula. We call this theory chemical formula theory (CFT). In the CFT, we place the valence ground and excited states of each atom at each position of the chemical formula of the molecule and let them interact using their free valences to form the ground and excited states of the molecule. The principle there is the variational principle so that the ground and excited states obtained satisfy the orthogonality and Hamiltonian-orthogonality relations. Then, we formulate the exact FC theory starting from the initial functions produced by the CFT. This FC theory is referred to as free-complement chemical-formula theory (FC-CFT), which is expected to describe efficiently the solution of the SE by the above reason. The FC-CFT wave function is modified from that of CFT. Since this modification is done by the exact SE, its analysis may give some insights to chemists that assist their chemistry. Thus, this theory would be not only exact but also conceptually useful. Furthermore, the intermediate theory between CFT and FC-CFT would also be useful. There, we use only integratable functions and apply the variational principle so that we refer to this theory as FC-CFT-variational (FC-CFT-V). It is an advanced theory of CFT. Since the variational method is straightforward and powerful, we can do extensive chemical studies in a reasonable accuracy. After finishing such studies, if we still need an exact level of solutions, we add the remaining functions of the FC-CFT and perform the exact calculations. Furthermore, when we deal with large and even giant molecules, the inter-exchange (iExg) theory for the antisymmetry rule introduced previously leads to a large simplification. There, the inter-exchanges between distant electron pairs fade away so that only Coulombic interactions survive. Further in giant systems, even an electrostatic description becomes possible. Then, the FC-CFT for exactly solving the SE would behave essentially to order N for large and giant molecular systems, though the pre-factor should be very large and must be minimized.

Time-Dependent Coupled Perturbed Hartree-Fock and Density-Functional-Theory Approach for Calculating Frequency-Dependent (Hyper)Polarizabilities with Nonorthogonal Localized Molecular Orbitals

The time-dependent coupled perturbed Hartree-Fock/density-functional-theory (TDHF/TDDFT) approach has been reformulated based on nonorthogonal localized molecular orbitals (NOLMOs). Based on the NOLMO Fock equation, we have derived the corresponding NOLMO-TDHF/TDDFT equations up to the third order, and the formula for the frequency-dependent (hyper)polarizabilities has been given. Our approach has been applied to calculate both static and dynamic (hyper)polarizabilities of molecules varying from small molecules to large molecules. The NOLMO-TDHF/TDDFT approach can reproduce the reference canonical molecular orbital (CMO) results for all of our testing calculations. With the help of ongoing development of optimized local virtual molecular orbitals, the NOLMO-TDHF/TDDFT approach would be a very efficient method for large system calculations and tp achieve linear scaling.

Automated property optimization via ab initio O(N) elongation method: Application to (hyper-)polarizability in DNA

An automated property optimization method was developed based on the ab initio O(N) elongation (ELG) method and applied to the optimization of nonlinear optical (NLO) properties in DNA as a first test. The ELG method mimics a polymerization reaction on a computer, and the reaction terminal of a starting cluster is attacked by monomers sequentially to elongate the electronic structure of the system by solving in each step a limited space including the terminal (localized molecular orbitals at the terminal) and monomer. The ELG-finite field (ELG-FF) method for calculating (hyper-)polarizabilities was used as the engine program of the optimization method, and it was found to show linear scaling efficiency while maintaining high computational accuracy for a random sequenced DNA model. Furthermore, the self-consistent field convergence was significantly improved by using the ELG-FF method compared with a conventional method, and it can lead to more feasible NLO property values in the FF treatment. The automated optimization method successfully chose an appropriate base pair from four base pairs (A, T, G, and C) for each elongation step according to an evaluation function. From test optimizations for the first order hyper-polarizability (β) in DNA, a substantial difference was observed depending on optimization conditions between "choose-maximum" (choose a base pair giving the maximum β for each step) and "choose-minimum" (choose a base pair giving the minimum β). In contrast, there was an ambiguous difference between these conditions for optimizing the second order hyper-polarizability (γ) because of the small absolute value of γ and the limitation of numerical differential calculations in the FF method. It can be concluded that the ab initio level property optimization method introduced here can be an effective step towards an advanced computer aided material design method as long as the numerical limitation of the FF method is taken into account.

The fragment molecular orbital method combined with density-functional tight-binding and the polarizable continuum model

The electronic gap in proteins is analyzed in detail, and it is shown that FMO-DFTB/PCM is efficient and accurate in describing the molecular structure of proteins in solution.

Intermediate electrostatic field for the generalized elongation method

An intermediate electrostatic field is introduced to improve the accuracy of fragment-based quantum-chemical computational methods by including long-range polarizations of biomolecules. The point charge distribution of the intermediate field is generated by a charge sensitivity analysis that is parameterized for five different population analyses, namely, atoms-in-molecules, Hirshfeld, Mulliken, natural orbital, and Voronoi population analysis. Two model systems are chosen to demonstrate the performance of the generalized elongation method (ELG) combined with the intermediate electrostatic field. The calculations are performed for the STO-3G, 6-31G, and 6-31G(d) basis sets and compared with reference Hartree-Fock calculations. It is shown that the error in the total energy is reduced by one order of magnitude, independently of the population analyses used. This demonstrates the importance of long-range polarization in electronic-structure calculations by fragmentation techniques.

Electron-correlated fragment-molecular-orbital calculations for biomolecular and nano systems

Recent developments in the fragment molecular orbital (FMO) method for theoretical formulation, implementation, and application to nano and biomolecular systems are reviewed. The FMO method has enabled ab initio quantum-mechanical calculations for large molecular systems such as protein-ligand complexes at a reasonable computational cost in a parallelized way. There have been a wealth of application outcomes from the FMO method in the fields of biochemistry, medicinal chemistry and nanotechnology, in which the electron correlation effects play vital roles. With the aid of the advances in high-performance computing, the FMO method promises larger, faster, and more accurate simulations of biomolecular and related systems, including the descriptions of dynamical behaviors in solvent environments. The current status and future prospects of the FMO scheme are addressed in these contexts.

Coupled-Perturbed SCF Approach for Calculating Static Polarizabilities and Hyperpolarizabilities with Nonorthogonal Localized Molecular Orbitals

Coupled-perturbed self-consistent-field (CPSCF) approach has been broadly used for polarizabilities and hyperpolarizabilities computation. To extend this application to large systems, we have reformulated the CPSCF equations with nonorthogonal localized molecular orbitals (NOLMOs). NOLMOs are the most localized representation of electronic degrees of freedom. Methods based on NOLMOs are potentially ideal for investigating large systems. In atomic orbital representation, with a static external electric field added, the wave function and SCF operator of unperturbed NOLMO-SCF wave function/orbitals are expanded to different orders of perturbations. We have derived the corresponding equations up to the third order, which are significantly different from those of a conventional CPSCF method because of the release of the orthogonal restrictions on MOs. The solution to these equations has been implemented. Several chemical systems are used to verify our method. This work represents the first step toward efficient calculations of molecular response and excitation properties with NOLMOs.

Localization of molecular orbitals: from fragments to molecule

Conspectus Localized molecular orbitals (LMO) not only serve as an important bridge between chemical intuition and molecular wave functions but also can be employed to reduce the computational cost of many-body methods for electron correlation and excitation. Therefore, how to localize the usually completely delocalized canonical molecular orbitals (CMO) into confined physical spaces has long been an important topic: It has a long history but still remains active to date. While the known LMOs can be classified into (exact) orthonormal and nonorthogonal, as well as (approximate) absolutely localized MOs, the ways for achieving these can be classified into two categories, a posteriori top-down and a priori bottom-up, depending on whether they invoke the global CMOs (or equivalently the molecular density matrix). While the top-down approaches have to face heavy tasks of minimizing or maximizing a given localization functional typically of many adjacent local extrema, the bottom-up ones have to invoke some tedious procedures for first generating a local basis composed of well-defined occupied and unoccupied subsets and then maintaining or resuming the locality when solving the Hartree-Fock/Kohn-Sham (HF/KS) optimization condition. It is shown here that the good of these kinds of approaches can be combined together to form a very efficient hybrid approach that can generate the desired LMOs for any kind of gapped molecules. Specifically, a top-down localization functional, applied to individual small subsystems only, is minimized to generate an orthonormal local basis composed of functions centered on the preset chemical fragments. The familiar notion for atomic cores, lone pairs, and chemical bonds emerges here automatically. Such a local basis is then employed in the global HF/KS calculation, after which a least action is taken toward the final orthonormal localized molecular orbitals (LMO), both occupied and virtual. This last step is very cheap, implying that, after the CMOs, the LMOs can be obtained essentially for free. Because molecular fragments are taken as the basic elements, the approach is in the spirit of "from fragments to molecule". Two representatives of highly conjugated molecules, that is, C12H2 and C60, are taken as showcases for demonstrating the success of the proposed approach. The use of the so-obtained LMOs will lead naturally to low-order scaling post-HF/KS methods for electron correlation or excitation. In addition, the underlying fragment picture allows for easy and pictorial interpretations of the correlation/excitation dynamics.

Intermediate electrostatic field for the elongation method

A simple way to improve the accuracy of the fragmentation methods is proposed. The formalism was applied to the elongation (ELG) method at restricted open-shell Hartree-Fock (ROHF) level of theory. The α-helix conformer of polyglycine was taken as a model system. The modified ELG method includes a simplified electrostatic field resulting from point-charge distribution of the system's environment. In this way the long-distance polarization is approximately taken into account. The field attenuates during the ELG process to eventually disappear when the final structure is reached. The point-charge distributions for each ELG step are obtained from charge sensitivity analysis (CSA) in force-field atoms resolution. The presence of the intermediate field improves the accuracy of ELG calculations. The errors in total energy and its kinetic and potential contributions are reduced by at least one-order of magnitude. In addition the SCF convergence of ROHF scheme is improved.

Effective preconditioning for ab initio ground state energy minimization with non-orthogonal localized molecular orbitals

The non-orthogonal localized molecular orbital (NOLMO) is the most localized representation of electronic degrees of freedom. As such, NOLMOs are thus potentially the most efficient for linear-scaling calculations of electronic structures for large systems. However, direct ab initio calculations with NOLMO have not been fully implemented and widely used, partly because of the slow convergence issue in the optimization of NOLMO. Towards realizing the potential of NOLMO for large systems, we applied an energy minimum variational principle for carrying out ab initio self-consistent-field (SCF) calculations with NOLMOs. We developed an effective preconditioning approach using the diagonal part of the second order derivatives and show that the convergence of the energy optimization is significantly improved. The speed of convergence of the energy and density are comparable with that of the conventional SCF approach, thus paving the way for the optimization of NOLMO in linear scaling calculations for large systems.

An Efficient Local Molecular Dynamics Polymerization Simulation Combined with an Ab Initio MO Method

A new local ab initio molecular dynamics method, namely elongation molecular dynamics (ELG-MD) is proposed for highly efficient simulations of aperiodic polymer systems. ELG-MD combines the elongation method (ELG) with the Gear predictor corrector (GPC) algorithm of molecular dynamics simulation. In this method, the local gradients acting on the atom’s nucleus in the active region are calculated by the ELG method while the equations of the nucleus’s motion are solved by the GPC algorithm. In this work, the first application of this ELG-MD method is described to investigate the stable conformation of polyglycine with surrounding water molecules. The water effects on the structure of polyglycine are examined. The ELG-MD simulations show that the formation of the polyglycine helix is strongly induced by the hydrogen bonds observed in two types of H-bond rings.