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

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.

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.