Source author record

Haibo Ma

Haibo Ma appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.chem-ph cond-mat.str-el cond-mat.mtrl-sci cond-mat.other

Catalog footprint

What is connected

9works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2021arXiv

Stochastic Adaptive Single-Site Time-Dependent Variational Principle

In recent years, the time-dependent variational principle (TDVP) method based on the matrix product state (MPS) wave function formulation has shown its great power in performing large-scale quantum dynamics simulations for realistic chemical systems with strong electron-vibration interactions. In this work, we propose a new stochastic adaptive single-site TDVP (SA-1TDVP) scheme to evolve the bond-dimension adaptively, which can integrate the tra-ditional advantages of both the high efficiency of single-site TDVP (1TDVP) variant and the high accuracy of the two-site TDVP (2TDVP) variant. Based on the assumption that the level statistics of entanglement Hamiltonians, which originate from the reduced density matrices of the MPS method, follows a Poisson or Wigner distribution, as generically predicted by random matrix theory, addi-tional random singular values are generated to expand the bond-dimension automatically. Tests on simulating the vibrationally-resolved quantum dynamics and absorption spectra in the pyrazine molecule and perylene bisimide (PBI) J-aggregate trimer as well as a spin-1/2 Heisenberg chain show that it can be automatic and as accurate as 2TDVP but reduce the computational time remarkably.

preprint2020arXiv

Charge transfer via deep hole in the J51/N2200 blend

In recently developed non-fullerene acceptor (NFA) based organic solar cells (OSCs), both the donor and acceptor parts can be excited by absorbing light photons. Therefore, both electron transfer and hole transfer channels could occur at the donor/acceptor interface for generating free charge carriers in NFA based OSCs. However, in many molecular and DNA systems, recent studies revealed the high charge transfer (CT) efficiency cannot be reasonably explained by a CT model with only highest occupied molecular orbitals (HOMOs) and lowest unoccupied molecular orbitals (LUMOs) of donor and acceptor molecules. In this work, taking an example of a full-polymer blend consisting of benzodithiophenealt-benzotriazole copolymers (J51) as donor and naphthalene diimide-bithiophene (N2200) as acceptor, in which the ultrafast hole transfer has been recently reported, we investigate its CT process and examine the different roles of various frontier molecular orbitals. Through a joint study of quantum mechanics electronic structure calculation and nonadiabatic dynamics simulation, we find the hole transfer between HOMOs of J51 and N2200 can hardly happen but the hole transfer from HOMO of N2200 to HOMO-1 of J51 is much more efficient. This points out the underlying importance of deep hole channel in CT process and indicates that including frontier molecular orbitals (FMOs) other than HOMOs and LUMOs is highly necessary to build a robust physical model for studying CT process in molecular optoelectronic materials.

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).

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.

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.

preprint2010arXiv

Effect of Electron-Electron Interactions on the Charge Carrier Transitions in trans-Polyacetylene

By employing a newly developed dynamical simulation method, which is a combination of classical molecular dynamics (MD) and the adaptive time-dependent density matrix renormalization group (TDDMRG), we investigate the dynamics of charge carrier transitions in trans-polyacetylene (PA) with the inclusion of both electron-phonon and electron -electron interactions. The calculations are performed within a modified Su-Schrieffer-Heeger (SSH) model in which electron-electron interactions are taken into account via the combination with extended Hubbard model (EHM). We find that removing an electron from a trans-PA chain bearing a positively charged polaron leads to the formation of a pair of charged solitons. Furthermore, we study the effect of electron-electron interactions on such charge carrier transitions in trans-PA. Our results show that increasing the on-site Coulomb interaction U and the nearest-neighbor Coulomb repulsion V will not change the qualitative behavior of the transition from a polaron to a soliton pair in the evolution process but will quantitatively reduce the moving velocities of the both formed solitons significantly and change the conditions for the soliton collisions.

preprint2008arXiv

Dynamical simulations of polaron transport in conjugated polymers with the inclusion of electron-electron interactions

Dynamical simulations of polaron transport in conjugated polymers in the presence of an external time-dependent electric field have been performed within a combined extended Hubbard model (EHM) and Su-Schrieffer-Heeger (SSH) model. Nearly all relevant electron-phonon and electron-electron interactions are fully taken into account by solving the time-dependent Schrödinger equation for the $π$-electrons and the Newton's equation of motion for the backbone monomer displacements by virtue of the combination of the adaptive time-dependent density matrix renormalization group (TDDMRG) and classical molecular dynamics (MD). We find that after a smooth turn-on of the external electric field the polaron is accelerated at first and then moves with a nearly constant velocity as one entity consisting of both the charge and the lattice deformation. An ohmic region (3 mV/$\textÅ$ $\leq E_0\leq$ 9 mV/$\textÅ$) where the stationary velocity increases linearly with the electric field strength is observed for the case of $U$=2.0 eV and $V$=1.0 eV. The maximal velocity is well above the speed of sound. Below 3 mV/$\textÅ$ the polaron velocity increases nonlinearly and in high electric fields with strength $E_0\geq$ 10.0 mV/$\textÅ$ the polaron will become unstable and dissociate. The relationship between electron-electron interaction strengths and polaron transport is also studied in detail. We find that the the on-site Coulomb interactions $U$ will suppress the polaron transport and small nearest-neighbor interactions $V$ values are also not beneficial to the polaronic motion while large $V$ values favor the polaron transport.