Source author record

Yingjin Ma

Yingjin Ma 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

8works
5topics
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

8 published item(s)

preprint2020arXiv

Multi-Reference Epstein-Nesbet Perturbation Theory with Density Matrix Renormalization Group Reference Wavefunction

The accurate electronic structure calculation for strongly correlated chemical systems requires an adequate description for both static and dynamic electron correlation, and is a persistent challenge for quantum chemistry. In order to account for static and dynamic electron correlations accurately and efficiently, in this work we propose a new method by integrating the density matrix renormalization group (DMRG) method and multi-reference second-order Epstein-Nesbet perturbation theory (ENPT2) with a selected configuration interaction (SCI) approximation. Compared with previous DMRG-based dynamic correlation methods, the DMRG-ENPT2 method extends the range of applicability, allowing us to efficiently calculate systems with very large active space beyond 30 orbitals. We demonstrate this by performing calculations on H$_2$S with an active space of (16e, 15o), hexacene with an active space of (26e, 26o) and 2D H$_{64}$ square lattice with an active space of (42e, 42o).

preprint2020arXiv

Portably parallel construction of a CI wave function from a matrix-product state using the Charm++ framework

The constructions of configuration interaction (CI) expansions from a matrix-product state (MPS) involves numerous matrix operations and the skillful sampling of important configurations when in a huge Hilbert space. In this work, we present an efficient procedure for constructing CI expansions from MPS using the Charm++ parallel programming framework, upon which automatic load balancing and object migration facilities can be employed. This procedure was employed in the MPS-to-CI utility (Moritz et al., J. Chem. Phys. 2007, 126, 224109), sampling-reconstructed complete active space algorithm (SR-CAS, Boguslawski et al., J. Chem. Phys. 2011, 134, 224101) and entanglement-driven genetic algorithm (EDGA, Luo et al., J. Chem. Theory Comput. 2017, 13, 4699-4710). It enhances productivity and allows the sampling programs evolve to their population-expansion versions (e.g., EDGA with population expansion [PE-EDGA]). Examples of 1,2-dioxetanone and firefly dioxetanone anion (FDO-) molecules demonstrated that 1) the procedure could be flexibly employed among various multi-core architectures; 2) the parallel efficiencies could be persistently improved simply by increasing the proportion of asynchronous executions; 3) PE-EDGA could construct a CAS-type CI wave function from a huge Hilbert space, with 0.9952 CI completeness and 96.7% correlation energy via ~1.66x10^6 configurations (only 0.0000028% of the total configurations) of a bi-radical state of FDO- molecule using the full valence active space within a few hours.

preprint2019arXiv

Approximate analytical gradients and nonadiabatic couplings for the state-average density matrix renormalization group self-consistent field method

We present an approximate scheme for analytical gradients and nonadiabatic couplings for calculating state-average density matrix renormalization group self-consistent-field wavefunction. Our formalism follows closely the state-average complete active space self-consistent-field (SA-CASSCF) \emph{ansatz}, which employs a Lagrangian, and the corresponding Lagrange multipliers are obtained from a solution of the coupled-perturbed CASSCF (CP-CASSCF) equations. We introduce a definition of the matrix product state (MPS) Lagrange multipliers based on a single-site tensor in a mixed-canonical form of the MPS, such that a sweep procedure is avoided in the solution of the CP-CASSCF equations. We apply our implementation to the optimization of a conical intersection in 1,2-dioxetanone, where we are able to fully reproduce the SA-CASSCF result up to arbitrary accuracy.

preprint2018arXiv

Externally-Contracted Multi-Reference Configuration Interaction Method Using a DMRG Reference Wave Function

The recent development of the density matrix renormalization group (DMRG) method in multireference quantum chemistry makes it practical to evaluate static correlation in a large active space, while dynamic correlation provides a critical correction to the DMRG reference for strong-correlated systems and is usually obtained using multi-reference perturbation (MRPT) or configuration interaction (MRCI) methods with internal contraction (ic) approximation. These methods can use active space scalable to relatively larger size references than has previously been possible. However, they are still hardly applicable to systems with active space larger than 30 orbitals because of high computation and storage costs of high-order reduced density matrices (RDMs) and the number of virtual orbitals are normally limited to few hundreds. In this work, we propose a new effective implementation of DMRG-MRCI, in which we use re-constructed CASCI-type configurations from DMRG wave function via the entropy-driving genetic algorithm (EDGA), and integrate with MRCI by an external contraction (ec) scheme. This bypasses the bottleneck of computing high-order RDMs in traditional DMRG dynamic correlation methods with ic approximation and the number of MRCI configurations is not dependent on the number of virtual orbitals. Therefore, DMRG-ec-MRCI method is promising for dealing with larger active space than 30 orbitals and large basis sets. We demonstrate the capability of our DMRG-ec-MRCI method in several benchmark applications, including the evaluation of potential energy curve of Cr$_{2}$, single-triplet gaps of higher n-acene molecules and the energy of Eu-BTBP(NO$_3$)$_3$ complex.

preprint2015arXiv

Density-Matrix Renormalization Group Algorithm with Multi-Level Active Space

The density-matrix renormalization group (DMRG) method, which can deal with a large active space composed of tens of orbitals, is nowadays widely used as an efficient addition to traditional complete active space (CAS)-based approaches. In this paper, we present the DMRG algorithm with a multi-level (ML) control of the active space based on chemical intuition-based hierarchical orbital ordering, which is called as ML-DMRG with its self-consistent field variant ML-DMRG-SCF. Ground and excited state calculations of H2O, N2, indole, and Cr2 with comparisons to DMRG references using fixed number of kept states (M) illustrate that MLtype DMRG calculations can obtain noticeable efficiency gains. It is also shown that the orbital re-ordering based on hierarchical multiple active subspaces may be beneficial for reducing computational time for not only ML-DMRG calculations but also DMRG ones with fixed M values.

preprint2014arXiv

Self-Consistent Embedding of Density-Matrix Renormalization Group Wavefunctions in a Density Functional Environment

We present the first implementation of a density matrix renormalization group algorithm embedded in an environment described by density functional theory. The frozen density embedding scheme is used with a freeze-and-thaw strategy for a self-consistent polarization of the orbital-optimized wavefunction and the environmental densities with respect to each other.

preprint2013arXiv

Assessment of various natural orbitals as the basis of large active space density matrix renormalization group calculations

It is well-known that not only the orbital ordering but also the choice of the orbitals themselves as the basis may significantly influence the computational efficiency of density-matrix renormalization group (DMRG) calculations. In this study, for assessing the efficiency of using various natural orbitals (NOs) as the DMRG basis, we performed benchmark DMRG calculations with different bases, which included the NOs obtained by various traditional electron correlation methods, as well as NOs acquired from preliminary moderate DMRG calculations (e.g., preserved states less than 500). The tested systems included N$_2$, transition metal Cr$_2$ systems, as well as 1-D hydrogen polyradical chain systems under equilibrium and dissociation conditions and 2-D hydrogen aggregates. The results indicate that a good compromise between the requirement for low computational costs of acquiring NOs and the demand for high efficiency of NOs as the basis of DMRG calculations may be very dependent on the studied systems' diverse electron correlation characteristics and the size of the active space. It is also shown that a DMRG-complete active space configuration interaction (DMRG-CASCI) calculation in a basis of carefully chosen NOs can provide a less expensive alternative to the standard DMRG-CASSCF calculation and avoid the convergence difficulties of orbital optimization for large active spaces. The effect of different NO ordering schemes on DMRG-CASCI calculations is also discussed.

preprint2013arXiv

Calculating excited states of molecular aggregates by the renormalized excitonic method

In this paper, we apply the recently developed ab initio renormalized excitonic method (REM) to the excitation energy calculations of various molecular aggregates, through the extension of REM to the time-dependent density functional theory (TDDFT). Tested molecular aggregate systems include one-dimensional hydrogen-bonded water chains, ring crystals with $π$-$π$ stacking or van-der Waals interactions and the general aqueous systems with polar and non-polar solutes. The basis set factor as well as the effect of the exchange-correlation functionals are also investigated. The results indicate that the REM-TDDFT method with suitable basis set and exchange-correlation functionals can give good descriptions of excitation energies and excitation area for lowest electronic excitations in the molecular aggregate systems with economic computational costs. It's shown that the deviations of REM-TDDFT excitation energies from those by standard TDDFT are much less than 0.1 eV and the computational time can be reduced by one order.