Search

We propose a method for calculating the non-radiative decay rates for polyatomic molecules including anharmonic effects of the potential energy surface (PES) in the Franck-Condon region. The method combines the n-mode representation method to construct the ab initio PES and the nearly exact time-dependent density matrix renormalization group method (TD-DMRG) to simulate quantum dynamics. In addition, in the framework of TD-DMRG, we further develop an algorithm to calculate the final-state-resolved rate coefficient which is very useful to analyze the contribution from each vibrational mode to the transition process. We use this method to study the internal conversion (IC) process of azulene after taking into account the anharmonicity of the ground state PES. The results show that even for this semi-rigid molecule, the intramode anharmonicity enhances the IC rate significantly, and after considering the two-mode coupling effect, the rate increases even further. The reason is that the anharmonicity enables the C–H vibrations to receive electronic energy while C–H vibrations do not contribute on the harmonic PES as the Huang-Rhys factor is close to 0.

A class of frustrated one-dimensional periodic Heisenberg spin systems formed either by triangular unit cells with spin 1/2 or by composite unit cells formed by two different structural units, triangles and small linear segments formed by an odd number of spin-1/2 is investigated. Based on perturbative processing and numerical calculations of the density matrix renormalization group method, the gapless character of the exact energy spectrum of excitation for these systems was found. Their instability with respect to regular (Peierls) oscillations of interactions between structural units is demonstrated. The corresponding critical exponents for the energies of the ground state are estimated numerically. For some frustrated systems, a quantum phase transition associated with the spin symmetry of the ground state, caused by frustration, has been discovered.

In this work, we propose a new method to calculate the molecular nonradiative electronic relaxation rates based on the numerically exact time-dependent density matrix renormalization group theory (TD-DMRG). This method could go beyond the existing frameworks under the harmonic approximation (HA) of the potential energy surface (PES) so that the anharmonic effect could be considered, which is of vital importance when the electronic energy gap is much larger than the vibrational frequency. We calculate the internal conversion (IC) rates in a two-mode model with Morse potential to investigate the validity of HA. We find that HA is unsatisfactory unless only the lowest several vibrational states of the lower electronic state are involved in the transition process when the adiabatic excitation energy is relatively low. As the excitation energy increases, HA first underestimates and then overestimates the IC rates when the excited state PES shifts towards the dissociative side of the ground state PES. On the contrary, HA slightly overestimates the IC rates when the excited state PES shifts towards the repulsive side. In both cases, higher temperature enlarges the error of HA. As a real example to demonstrate the effectiveness and scalability of the method, we calculate the IC rates of azulene from $S_1$ to $S_0$ on the ab initio anharmonic PES approximated by 1-mode representation. The calculated IC rates of azulene under HA are consistent with the analytically exact results. The rates on anharmonic PES are 30%-40% higher than the rates under HA.

The magnetization process of the S = 1 antiferromagnetic chain with the biquadratic interaction is investigated using the numerical diagonalization. As a result, it is found that the quantum spin nematic liquid phase appears below the saturation magnetization for sufficiently large negative biquadratic interaction. The ground state phase diagram is also presented. This result is consistent with the previous density matrix renormalization group study. In addition the present study confirms that the two magnon excitation is gapless, while the single magnon one is gapped in the spin nematic liquid phase.

We introduce the transcorrelated Density Matrix Renormalization Group (tcDMRG) theory for the efficient approximation of the energy for strongly correlated systems. tcDMRG encodes the wave function as a product of a fixed Jastrow or Gutzwiller correlator and a matrix product state. The latter is optimized by applying the imaginary-time variant of time-dependent (TD) DMRG to the non-Hermitian transcorrelated Hamiltonian. We demonstrate the efficiency of tcDMRG at the example of the two-dimensional Fermi-Hubbard Hamiltonian, a notoriously difficult target for the DMRG algorithm, for different sizes, occupation numbers, and interaction strengths. We demonstrate fast energy convergence of tcDMRG, which indicates that tcDMRG could increase the efficiency of standard DMRG beyond quasi-monodimensional systems and provides a generally powerful approach toward the dynamic correlation problem of DMRG., 31 pages, 4 figures, 4 tables

Marcus theory has been successfully applied to molecular design for organic semiconductors with the aid of quantum chemistry calculations for the molecular parameters: the intermolecular electronic coupling V and the intramolecular charge reorganization energy λ. The assumption behind this is the localized nature of the electronic state for representing the charge carriers, being holes or electrons. As far as the quantitative description of carrier mobility is concerned, the direct application of Marcus semiclassical theory usually led to underestimation of the experimental data. A number of effects going beyond such a semiclassical description will be introduced here, including the quantum nuclear effect, dynamic disorder, and delocalization effects. The recently developed quantum dynamics simulation at the time-dependent density matrix renormalization group theory is briefly discussed. The latter was shown to be a quickly emerging efficient quantum dynamics method for the complex system.

We propose a multireference linearized coupled cluster theory using matrix product states (MPS-LCC) which provides remarkably accurate ground-state energies, at a computational cost that has the same scaling as multireference configuration interaction singles and doubles (MRCISD), for a wide variety of electronic Hamiltonians. These range from first-row dimers at equilibrium and stretched geometries, to highly multireference systems such as the chromium dimer and lattice models such as periodic two-dimensional 1-band and 3-band Hubbard models. The MPS-LCC theory shows a speed up of several orders of magnitude over the usual DMRG algorithm while delivering energies in excellent agreement with converged DMRG calculations. Also, in all the benchmark calculations presented here MPS-LCC outperformed the commonly used multi-reference quantum chemistry methods in some cases giving energies in excess of an order of magnitude more accurate. As a size-extensive method that can treat large active spaces, MPS-LCC opens up the use of multireference quantum chemical techniques in strongly-correlated \emph{ab-initio} Hamiltonians, including two and three-dimensional solids., Comment: Journal of Chemical Physics (2015)

We present how to compute vibrational eigenstates with tree tensor network states (TTNSs), the underlying ansatz behind the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method. The eigenstates are computed with an algorithm that is based on the density matrix renormalization group (DMRG). We apply this to compute the vibrational spectrum of acetonitrile (CH$_3$CN) to high accuracy and compare TTNSs with matrix product states (MPSs), the ansatz behind the DMRG. The presented optimization scheme converges much faster than ML-MCTDH-based optimization. For this particular system, we found no major advantage of the more general TTNS over MPS. We highlight that for both TTNS and MPS, the usage of an adaptive bond dimension significantly reduces the amount of required parameters. We furthermore propose a procedure to find good trees.

Configuration-interaction-type calculations on electronic and vibrational structure are often the method of choice for the reliable approximation of many-particle wave functions and energies. The exponential scaling, however, limits their application range. An efficient approximation to the full configuration interaction solution can be obtained with the density matrix renormalization group (DMRG) algorithm without a restriction to a predefined excitation level. In a standard DMRG implementation, however, excited states are calculated with a ground-state optimization in the space orthogonal to all lower lying wave function solutions. A trivial parallelization is therefore not possible and the calculation of highly excited states becomes prohibitively expensive, especially in regions with a high density of states. Here, we introduce two variants of the density matrix renormalization group algorithm that allow us to target directly specific energy regions and therefore highly excited states. The first one, based on shift-and-invert techniques, is particularly efficient for low-lying states, but is not stable in regions with a high density of states. The second one, based on the folded auxiliary operator, is less efficient, but more accurate in targeting high-energy states. We apply the algorithm to the solution of the nuclear Schroedinger equation, but emphasize that it can be applied to the diagonalization of general Hamiltonians as well, such as the electronic Coulomb Hamiltonian to address X-ray spectra. In combination with several root-homing algorithms and a stochastic sampling of the determinant space, excited states of interest can be adequately tracked and analyzed during the optimization. We demonstrate that we can accurately calculate prominent spectral features of large molecules such as the sarcosyn-glycine dipeptide., 26 pages, 7 figures, 3 tables

The Auxiliary-Field Quantum Monte Carlo (AFQMC) algorithm is a powerful quantum many-body method that can be used successfully as an alternative to standard quantum chemistry approaches to compute the ground state of many body systems, such as molecules and solids, with high accuracy. In this article we use AFQMC with trial wave-functions built from non-orthogonal multi Slater determinant expansions to study the energetics of molecular systems, including the 55 molecules of the G1 test set and the isomerization path of the $[Cu_{2}O_{2}]^{2+}$ molecule. The main goal of this study is to show the ability of non-orthogonal multi Slater determinant expansions to produce high-quality, compact trial wave-functions for quantum Monte Carlo methods. We obtain systematically improvable results as the number of determinants is increased, with high accuracy typically obtained with tens of determinants. Great reduction in the average error and traditional statistical indicators are observed in the total and absorption energies of the molecules in the G1 test set with as few as 10-20 determinants. In the case of the relative energies along the isomerization path of the $[Cu_{2}O_{2}]^{2+}$, our results compare favorably with other advanced quantum many-body methods, including DMRG and complete-renormalized CCSD(T). Discrepancies in previous studies for this molecular problem are identified and attributed to the differences in the number of electrons and active spaces considered in such calculations.

We have carried out density matrix renormalization group calculations on the T1 state of linear polyenes applying the Pariser-Parr-Pople (PPP) Model. The geometry optimization for the polyene oligomers C(2n)H(2n+2) (n = 4,5,6,...,15) shows that the S0 to T1 excitation region is composed of a soliton-antisoliton pair located symmetrically away from the center of the chain and leads to single- and double-bond interconversions in between. The distance between the soliton and antisoliton centers in T1 state changes with the length of the chain, contradictory to earlier conclusions obtained with PPP-SDCI or ab initio SCI methods. The inconsistency most possibly comes from the insufficient consideration of the electron correlations in small-scale CI methods.

We present a method for the approximate solution of the coupled-cluster (CC) equation, based upon the singular value decomposition (SVD). The key idea of this method is that we use SVD for the cluster amplitudes to exploit the physically important states in the CC equation. This method enables us to significantly reduce the computational requirements for CC calculations without losing the essence of the method. Relationships to the density matrix renormalization group theory and the local correlation methods are mentioned. We perform pilot calculations on some atoms and molecules to investigate the applicability of the method.

Using a newly developed density matrix renormalization group algorithm, we have computed exact solutions of the Schrödinger equation for water at two geometries in a basis of 41 orbitals. Calculations of this size cannot be carried out using any other method. We compare our results with high-order coupled cluster and configuration interaction calculations.

In earlier work [A. Y. Sokolov and G. K.-L. Chan, J. Chem. Phys. 144, 064102 (2016)], we introduced a time-dependent formulation of the second-order N-electron valence perturbation theory (t-NEVPT2) which (i) had a lower computational scaling than the usual internally contracted perturbation formulation and (ii) yielded the fully uncontracted NEVPT2 energy. Here, we present a combination of t-NEVPT2 with a matrix product state (MPS) reference wavefunction (t-MPS-NEVPT2) that allows us to compute uncontracted dynamic correlation energies for large active spaces and basis sets, using the time-dependent density matrix renormalization group algorithm. In addition, we report a low-scaling MPS-based implementation of strongly contracted NEVPT2 (sc-MPS-NEVPT2) that avoids computation of the four-particle reduced density matrix. We use these new methods to compute the dissociation energy of the chromium dimer and to study the low-lying excited states in all-trans polyenes (C_4H_6 to C_(24)H_(26)), incorporating dynamic correlation for reference wavefunctionswith up to 24 active electrons and orbitals.

This paper describes an interface between the density matrix renormalization group (DMRG) method and the complete active-space self-consistent field (CASSCF) method and its analytical gradient, as well as an extension to the second-order perturbation theory (CASPT2) method. This interfacing allows large active-space multi-reference computations to be easily performed. The interface and its extension are both implemented in terms of reduced density matrices (RDMs) which can be efficiently computed via the DMRG sweep algorithm. We also present benchmark results showing that, in practice, the DMRG-CASSCF calculations scale with active-space size in a polynomial manner in the case of quasi-1D systems. Geometry optimization of a binuclear iron-sulfur cluster using the DMRG-CASSCF analytical gradient is demonstrated, indicating that the inclusion of the valence p-orbitals of sulfur and double-shell d-orbitals of iron lead to non-negligible changes in the geometry compared to the results of small active-space calculations. With the exception of the selection of M values, many computational settings in these practical DMRG calculations have been tuned and black-boxed in our interface, and so the resulting DMRG-CASSCF and DMRG-CASPT2 calculations are now available to novice users as a common tool to compute strongly correlated electronic wavefunctions.

The present work studies the quasi one-dimensional $Ni(dmit)_2$-based compounds within a correlated model. More specifically, we focus our attention on the composed influence of the electronic dimerization-factor and the repulsion, on the transport properties and the localization of the electronic density in the ground-state. Those properties are studied through the computation of the charge gaps (difference between the ionization potential and the electro-affinity: IP-EA) and the long- and short-bond orders of an infinite quarter-filled chain within a $t-J(t,U)$ model. The comparison between the computed gaps and the experimental activation energy of the semiconductor $NH_2Me_2 [Ni(dmit)_2]_2$ allows us to estimate the on-site electronic repulsion of the $Ni(dmit)_2$ molecule to $1.16eV$.

We present selected results for the Holstein molecular crystal model in one space dimension as determined by the Global–Local variational method, including complete polaron energy bands, ground state energies, and effective masses. We juxtapose our results with specific comparable results of numerous other methodologies of current interest, including quantum Monte Carlo, cluster diagonalization, dynamical mean field theory, density matrix renormalization group, semiclassical analysis, weak-coupling perturbation theory, and strong-coupling perturbation theory. Taken as a whole, these methodologies are mutually confirming and provide a comprehensive and quantitatively accurate description of polaron properties in essentially any regime. In particular, this comparison confirms the Global–Local variational method as being highly accurate over a wide range of the polaron parameter space, from the nonadiabatic limit to the extremes of high adiabaticity, from weak coupling through intermediate coupling to strong coupling.

Spin-density maps, deduced from polarized neutron diffraction experiments, for both the pair and chain compounds of the system Mn2+Cu2+ have been reported recently. These results have motivated us to investigate theoretically the spin populations in such alternant mixed-spin systems. In this paper, we report our studies on the one-dimensional ferrimagnetic systems (SA,SB)N where N is the number of AB pairs. We have considered all cases in which the spin SA takes on allowed values in the range 1 to 72 while the spin SB is held fixed at 12. The theoretical studies have been carried out on the isotropic Heisenberg model, using the density matrix renormalization group method. The effect of the magnitude of the larger spin SA on the quantum fluctuations in both A and B sublattices has been studied as a function of the system size N. We have investigated systems with both periodic and open boundary conditions, the latter with a view to understanding end-of-chain effects. The spin populations have been followed ...

Using a symmetry adapted density matrix renormalization group (DMRG) formulation, we have calculated various optical spectra, including linear absorption, electroabsorption, third-harmonic generation, and two-photon absorption, for a series of linear polyene molecules (up to 40 carbon sites) modeled by an extended Hubbard–Peierls Hamiltonian within the intermediate correlation regime. The theoretical two-photon absorption spectrum for trans-octatetraene is in good agreement with experiment. By comparing the experimental linear absorption spectra and theoretical higher Ag state energies, it is suggested that the so-called cis-band actually corresponds to the threshold of the mAg band.

By analogy with inorganic semiconductors such as GaAs or Si, electron-hole models may be expected to provide a useful description of the excited states of conjugated polymers. Here, these models are tested against density matrix renormalization group (DMRG) calculations. The DMRG method is used to generate nearly-exact descriptions of the ground state, 11Bu optical gap state, and the band gap of the Pariser-Parr-Pople (PPP) Hamiltonian of polyenes with between 2 and 40 carbon atoms. These are compared with both bare electron-hole (singles configuration interaction theory and the random phase approximation) and dressed electron-hole (second and third order Green’s function) methods. For the optical gap, only second-order Green’s function results were obtained. When an unscreened (Ohno) electron-electron interaction potential is used, the dressed electron-hole methods work well for the band gap. The difference between the band gap predicted by bare and dressed electron-hole methods increases with chain leng...

We report the symmetrized density matrix renormalization group (DMRG) study of neutral and doped oligomers of poly-para-phenylene (PPP) system within an extended Hubbard model. Model parameters are determined by comparing the existing results for an interacting small system. We compute a number of properties in the ground state as well as in the one-photon, two-photon and triplet states to completely characterize these states. Bond-order studies show that the lowest two-photon state corresponds to a localized excitation while one-photon and triplet excitations are extended in nature. The bipolaronic state shows clear evidence for charge separation and disproportionation into two polarons. We find that the extended nature of one-photon and triplet states of the neutral system are very similar to those of the bipolaronic ground states.

Faceting diagrams for the step-faceting zone, the step droplet zone, and the Gruber-Mullins-Pokrovsky-Talapov (GMPT) zone for a crystal surface are obtained by using the density matrix renormalization group method to calculate the surface tension. The model based on these calculations is the restricted solid-on-solid (RSOS) model with a point-contact-type step-step attraction (p-RSOS model) on a square lattice. The point-contact-type step-step attraction represents the energy gain obtained by forming a bonding state with orbital overlap at the meeting point of the neighboring steps. In the step-faceting zone, disconnectedness in the surface tension leads to the formation of a faceted macrostep on a vicinal surface at equilibrium. The disconnectedness in the surface tension also causes the first-order shape transition for the equilibrium shape of a crystal droplet. The lower zone boundary line (ZBL), which separates the step-faceting zone and the step droplet zone, is obtained by the condition γ 1 = lim n → ∞ γ n / n , where γn is the step tension of the n-th merged step. The upper ZBL, which separates the GMPT zone and the step droplet zone, is obtained by the condition Aq,eff = 0 and Bq,eff = 0, where Aq,eff and Bq,eff represent the coefficients for the | q → | 2 term and the | q → | 3 term, respectively, in the | q → | -expanded form of the surface free energy f eff ( q → ) . Here, q → is the surface gradient relative to the (111) surface. The reason why the vicinal surface inclined in the 〈101〉 direction does not exhibit step-faceting is explained in terms of the one-dimensional spinless quasi-impenetrable attractive bosons at absolute zero.

Renormalization is a powerful concept in the many-body problem. Inspired by the highly successful density matrix renormalization group (DMRG) algorithm, and the quantum chemical graphical representation of configuration space, we introduce a new theoretical tool: Hilbert space renormalization, to describe many-electron correlations. While in DMRG, the many-body states in nested Fock subspaces are successively renormalized, in Hilbert space renormalization, many-body states in nested Hilbert subspaces undergo renormalization. This provides a new way to classify and combine configurations. The underlying wavefunction ansatz, namely the Hilbert space matrix product state (HS-MPS), has a very rich and flexible mathematical structure. It provides low-rank tensor approximations to any configuration interaction (CI) space through restricting either the 'physical indices' or the coupling rules in the HS-MPS. Alternatively, simply truncating the 'virtual dimension' of the HS-MPS leads to a family of size-extensive wave function ansaetze that can be used efficiently in variational calculations. We make formal and numerical comparisons between the HS-MPS, the traditional Fock-space MPS used in DMRG, and traditional CI approximations. The analysis and results shed light on fundamental aspects of the efficient representation of many-electron wavefunctions through the renormalization of many-body states., 23 pages, 14 figures, The following article has been submitted to The Journal of Chemical Physics

Density matrix theory is implemented in a variational quantum Monte Carlo computation of electronic properties of atoms and molecules. Differences between electronic densities from conventional and density matrix methods are detected. However, calculated properties present similar behavior and partial antisymmetry can be ignored in the cases studied.

The invariance of the renormalization group equation under the action of the prolonged Lie field is examined and it is shown that the system exhibits infinite dimensional symmetry. From the characteristics, detailed solutions of the renormalization group equation is derived. As an illustration, the method is applied to QED and QCD.

For scalar quantum electrodynamics on a three-dimensional toroidal lattice with a fine lattice spacing, we consider the renormalization problem of choosing counter terms depending on the lattice spacing, so that the theory stays finite as the spacing goes to zero. We employ a renormalization group method which analyzes the flow of the mass and the vacuum energy as a problem in discrete dynamical systems. The main result is that counter terms can be chosen so that at the end of the iteration these quantities take preassigned values. No use is made of perturbation theory. The renormalization group transformations are defined with bounded fields, an approximation which can be justified in Balaban’s approach to the renormalization group.

Polyacenes in their armchair geometry (phenacenes) have recently been found to possess appealing electronic and optical properties with higher chemical stability and comparatively larger band gap as compared to linear polyacenes. They also behave as high-temperature superconductors upon alkali metal doping. Moreover, the optical properties of crystalline picene can be finely tuned by applying external pressure. We investigated the variation of optical gap as a function of altering the interplanar distances between parallel cofacial phenacene dimers. We employed both time-dependent density functional theory and density matrix renormalization group (DMRG) technique to investigate the lowest singlet excitations in phenacene dimer. Our study showed that the lowest singlet excitation in these systems evolved as a function of interplanar separation. The optical excitation energy gap decreases as a function of inverse interplanar separation of the phenacene dimer. The distant dependent variation of optical absorption at the dimer level may be comparable with experimental observation in picene crystal under pressure. DMRG study also demonstrates that besides picene, electronic properties of higher phenacenes can also be tunable by altering interplanar separation.

Spin Orbit Coupling (SOC) is introduced to molecular ab initio density matrix renormalization group (DMRG) calculations. In the presented scheme, one first approximates the electronic ground state and a number of excited states of the Born-Oppenheimer (BO) Hamiltonian with the aid of the DMRG algorithm. Owing to the spin-adaptation of the algorithm, the total spin S is a good quantum number for these states. After the non-relativistic DMRG calculation is finished, all magnetic sublevels of the calculated states are constructed explicitly, and the SOC operator is expanded in the resulting basis. To this end, spin orbit coupled energies and wavefunctions are obtained as eigenvalues and eigenfunctions of the full Hamiltonian matrix which is composed of the SOC operator matrix and the BO Hamiltonian matrix. This treatment corresponds to a quasi-degenerate perturbation theory approach and can be regarded as the molecular equivalent to atomic Russell-Saunders coupling. For the evaluation of SOC matrix elements, the full Breit-Pauli SOC Hamiltonian is approximated by the widely used spin-orbit mean field operator. This operator allows for an efficient use of the second quantized triplet replacement operators that are readily generated during the non-relativistic DMRG algorithm, together with the Wigner-Eckart theorem. With a set of spin-orbit coupled wavefunctions at hand, the molecular g-tensors are calculated following the scheme proposed by Gerloch and McMeeking. It interprets the effective molecular g-values as the slope of the energy difference between the lowest Kramers pair with respect to the strength of the applied magnetic field. Test calculations on a chemically relevant Mo complex demonstrate the capabilities of the presented method.

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., 4 pages, 3 figures

The ab-initio density matrix renormalization group (DMRG) is a tool that can be applied to a wide variety of interesting problems in quantum chemistry. Here, we examine the density matrix renormalization group from the vantage point of the quantum chemistry user. What kinds of problems is the DMRG well-suited to? What are the largest systems that can be treated at practical cost? What sort of accuracies can be obtained, and how do we reason about the computational difficulty in different molecules? By examining a diverse benchmark set of molecules: π-electron systems, benchmark main-group and transition metal dimers, and the Mn-oxo-salen and Fe-porphine organometallic compounds, we provide some answers to these questions, and show how the density matrix renormalization group is used in practice.

Properly combined with the Lanczos iteration algorithm, the density matrix renormalization group enables us to directly calculate the dipole-allowed singlet states and the transition moment elements between these states and the ground state for longer chains. Our new scheme does not resort to reflection symmetry and electron-hole symmetry, nor to the construction of different superblocks for the ground state and these singlets separately. As an example, we calculate the ground state and first-excited state in the extended Hubbard model. Our calculations show that for a large U, the transition matrix elements decrease with the chain length; for V around U/2, the transition matrix elements can increase with the chain length, which is consistent with recent experimental results.

We extend the active space decomposition method, recently developed by us, to more than two active sites using the density matrix renormalization group algorithm. The fragment wave functions are described by complete or restricted active-space wave functions. Numerical results are shown on a benzene pentamer and a perylene diimide trimer. It is found that the truncation errors in our method decrease almost exponentially with respect to the number of renormalization states M, allowing for numerically exact calculations (to a few microhartrees or less) with M = 128 in both cases. This rapid convergence is because the renormalization steps are used only for the interfragment electron correlation.

We report an extension of our previous development that incorporated quantum-chemical density matrix renormalization group (DMRG) into the complete active space second-order perturbation theory (CASPT2) [Y. Kurashige and T. Yanai, J. Chem. Phys. 135, 094104 (2011)]. In the previous study, the combined theory, referred to as DMRG-CASPT2, was built upon the use of pseudo-canonical molecular orbitals (PCMOs) for one-electron basis. Within the PCMO basis, the construction of the four-particle reduced density matrix (4-RDM) using DMRG can be greatly facilitated because of simplicity in the multiplication of 4-RDM and diagonal Fock matrix in the CASPT2 equation. In this work, we develop an approach to use more suited orbital basis in DMRG-CASPT2 calculations, e.g., localized molecular orbitals, in order to extend the domain of applicability. Because the multiplication of 4-RDM and generalized Fock matrix is no longer simple in general orbitals, an approximation is made to it using the cumulant reconstruction neglecting higher-particle cumulants. Also, we present the details of the algorithm to compute 3-RDM of the DMRG wavefunction as an extension of the 2-RDM algorithm of Zgid et al. [J. Chem. Phys. 128, 144115 (2008)] and Chan et al. [J. Chem. Phys. 128, 144117 (2008)]. The performance of the extended DMRG-CASPT2 approach was examined for large-scale multireference systems, such as low-lying excited states of long-chain polyenes and isomerization potential of {[Cu(NH3)3]2O2}(2+).

A double quantum dot device, connected to two channels that only interact through interdot Coulomb repulsion, is analyzed using the numerical renormalization group technique. Using a two-impurity Anderson model, and realistic parameter values [S. Amasha, A. J. Keller, I. G. Rau, A. Carmi, J. A. Katine, H. Shtrikman, Y. Oreg, and D. Goldhaber-Gordon, Phys. Rev. Lett. 110, 046604 (2013)], it is shown that, by applying a moderate magnetic field and independently adjusting the gate potential of each quantum dot at half-filling, a spin-orbital SU(2) Kondo state can be achieved where the Kondo resonance originates from spatially separated parts of the device. Our results clearly link this spatial separation effect to currents with opposing spin polarizations in each channel, i.e., the device acts as a spin filter. In addition, an experimental probe of this polarization effect is suggested, pointing to the exciting possibility of experimentally probing the internal structure of an SU(2) Kondo state.

The applicability of ab initio multireference wavefunction-based methods to the study of magnetic complexes has been restricted by the quickly rising active-space requirements of oligonuclear systems and dinuclear complexes with S1 spin centers. Ab initio density matrix renormalization group (DMRG) methods built upon an efficient parameterization of the correlation network enable the use of much larger active spaces, and therefore may offer a way forward. Here, we apply DMRG-CASSCF to the dinuclear complexes [Fe2OCl6](2-) and [Cr2O(NH3)10](4+). After developing the methodology through systematic basis set and DMRG M testing, we explore the effects of extended active spaces that are beyond the limit of conventional methods. We find that DMRG-CASSCF with active spaces including the metal d orbitals, occupied bridging-ligand orbitals, and their virtual double shells already capture a major portion of the dynamic correlation effects, accurately reproducing the experimental magnetic coupling constant (J) of [Fe2OCl6](2-) with (16e,26o), and considerably improving the smaller active space results for [Cr2O(NH3)10](4+) with (12e,32o). For comparison, we perform conventional MRCI+Q calculations and find the J values to be consistent with those from DMRG-CASSCF. In contrast to previous studies, the higher spin states of the two systems show similar deviations from the Heisenberg spectrum, regardless of the computational method.

We present the first implementation of the relativistic quantum chemical two- and four-component density matrix renormalization group algorithm that includes a variational description of scalar-relativistic effects and spin-orbit coupling. Numerical results based on the four-component Dirac-Coulomb Hamiltonian are presented for the standard reference molecule for correlated relativistic benchmarks: thallium hydride.

We present a second-order perturbation theory based on a density matrix renormalization group self-consistent field (DMRG-SCF) reference function. The method reproduces the solution of the complete active space with second-order perturbation theory (CASPT2) when the DMRG reference function is represented by a sufficiently large number of renormalized many-body basis, thereby being named DMRG-CASPT2 method. The DMRG-SCF is able to describe non-dynamical correlation with large active space that is insurmountable to the conventional CASSCF method, while the second-order perturbation theory provides an efficient description of dynamical correlation effects. The capability of our implementation is demonstrated for an application to the potential energy curve of the chromium dimer, which is one of the most demanding multireference systems that require best electronic structure treatment for non-dynamical and dynamical correlation as well as large basis sets. The DMRG-CASPT2/cc-pwCV5Z calculations were performed with a large (3d double-shell) active space consisting of 28 orbitals. Our approach using large-size DMRG reference addressed the problems of why the dissociation energy is largely overestimated by CASPT2 with the small active space consisting of 12 orbitals (3d4s), and also is oversensitive to the choice of the zeroth-order Hamiltonian.

By using the properties of orthogonal polynomials, we present an exact unitary transformation that maps the Hamiltonian of a quantum system coupled linearly to a continuum of bosonic or fermionic modes to a Hamiltonian that describes a one-dimensional chain with only nearest-neighbour interactions. This analytical transformation predicts a simple set of relations between the parameters of the chain and the recurrence coefficients of the orthogonal polynomials used in the transformation, and allows the chain parameters to be computed using numerically stable algorithms that have been developed to compute recurrence coefficients. We then prove some general properties of this chain system for a wide range of spectral functions, and give examples drawn from physical systems where exact analytic expressions for the chain properties can be obtained. Crucially, the short range interactions of the effective chain system permits these open quantum systems to be efficiently simulated by the density matrix renormalization group methods., 24 pages, 1 figure, final version

In previous work we have shown that the Density Matrix Renormalization Group (DMRG) enables near-exact calculations in active spaces much larger than are possible with traditional Complete Active Space algorithms. Here, we implement orbital optimisation with the Density Matrix Renormalization Group to further allow the self-consistent improvement of the active orbitals, as is done in the Complete Active Space Self-Consistent Field (CASSCF) method. We use our resulting DMRGCASSCF method to study the low-lying excited states of the all-trans polyenes up to C24H26 as well as \beta-carotene, correlating with near-exact accuracy the optimised complete \pi-valence space with up to 24 active electrons and orbitals, and analyse our results in the light of the recent discovery from Resonance Raman experiments of new optically dark states in the spectrum., Comment: 16 pages, 8 figures

We present the density matrix renormalization group self-consistent field (DMRG-SCF) approach that is analogous to the complete active space self-consisted field (CASSCF) method but instead of using for the description of the active space the full configuration interaction (FCI) method, the DMRG-SCF uses the density matrix renormalization group (DMRG) method. The DMRG-SCF approach, similarly to CASSCF, properly describes the multiconfigurational character of the wave function but avoids the exponential scaling of the FCI method and replaces it with a polynomial scaling. Hence, calculations for a larger number of orbitals and electrons in the active space are possible since the DMRG method provides an efficient tool to automatically select from the full Hilbert space the many-body contracted basis states that are the most important for the description of the wave function.

We present a spin-adapted density matrix renormalization group (DMRG) algorithm designed to target spin and spatial symmetry states that can be difficult to obtain while using a non-spin-adapted algorithm. The algorithmic modifications that have to be introduced into the usual density matrix renormalization group scheme in order to spin adapt it are discussed, and it is demonstrated that the introduced modifications do not change the overall scaling of the method. The new approach is tested on HNCO, a model system, that has a singlet-triplet curve crossing between states of the same symmetry. The advantages of the spin-adapted DMRG scheme are discussed, and it is concluded that the spin-adapted DMRG method converges better in almost all cases and gives more parallel curves to the full configuration interaction result than the non-spin-adapted method. It is shown that the spin-adapted DMRG energies can be lower than the ones obtained from the non-spin-adapted scheme. Such a counterintuitive result is explained by noting that the spin-adapted method is not a special case of the non-spin-adapted one; consequently, the spin-adapted result is not an upper bound for the non-spin-adapted energy.

Over the past few years, it has been shown in various studies on small molecules with only a few electrons that the density-matrix renormalization group (DMRG) method converges to results close to the full configuration-interaction limit for the total electronic energy. In order to test the capabilities of the method for molecules with complex electronic structures, we performed a study on the potential-energy curves of the ground state and the first excited state of 1sigma+ symmetry of the cesium hydride molecule. For cesium relativistic effects cannot be neglected, therefore we have used the generalized arbitrary-order Douglas-Kroll-Hess protocol up to tenth order, which allows for a complete decoupling of the Dirac Hamiltonian. Scalar-relativistic effects are thus fully incorporated in the calculations. The potential curves of the cesium hydride molecule feature an avoided crossing between the ground state and the first excited state, which is shown to be very well described by the DMRG method. Compared to multireference configuration-interaction results, the potential curves hardly differ in shape, for both the ground state and the excited state, but the total energies from the DMRG calculations are in general consistently lower. However, the DMRG energies are as accurate as corresponding coupled cluster energies at the equilibrium distance, but convergence to the full configuration-interaction limit is not achieved.

The density-matrix renormalization group algorithm has emerged as a promising new method in ab initio quantum chemistry. However, many problems still need to be solved before this method can be applied routinely. At the start of such a calculation, the orbitals originating from a preceding quantum chemical calculation must be placed in a specific order on a one-dimensional lattice. This ordering affects the convergence of the density-matrix renormalization group iterations significantly. In this paper, we present two approaches to obtain optimized orderings of the orbitals. First, we use a genetic algorithm to optimize the ordering with respect to a low total electronic energy obtained at a predefined stage of the density-matrix renormalization group algorithm with a given number of total states kept. In addition to that, we derive orderings from the one- and two-electron integrals of our test system. This test molecule is the chromium dimer, which is known to possess a complicated electronic structure. For this molecule, we have carried out calculations for the various orbital orderings obtained. The convergence behavior of the density-matrix renormalization group iterations is discussed in detail.

The renormalization group for a field theory on a lattice, or equivalently, in the language of statistical mechanics, a system of continuous spins on a lattice, described by Wilson, is considered in detail. The conditions under which a class of transformations are renormalization group transformations are studied.

A new renormalization procedure is introduced via a method that assigns finite values to divergent expressions in quantum field theory. It is shown that this procedure works at the two‐loop level in Feynman diagrams with overlapping divergences. Applications to λφ4 theory and quantum electrodynamics are made. Also, the Casimir effect is evaluated. Comparison of these results with those from standard renormalization methods proves them to be coincident.

We introduce a correlation function description of the renormalization group approach to critical phenomena. Our work is based on treating the renormalization group operator as a linear mapping on the set of Ursell functions, rather than as a nonlinear mapping on the space of Hamiltonians. We mainly consider the ’’mean‐spin’’ renormalization group, but the closely related ’’decimation’’ transformation is also considered. Using this approach, we demonstrate for a suitable class of system that the spectrum of the renormalization group operator is bounded and countably infinitely degenerate. We give counterexamples to the notion that there must be convergence to a renormalization group fixed point. Our formulation of the renormalization group is sufficiently general so that convergence to a fixed point does not necessarily imply hyperscaling, i.e., the vanishing of the anomalous dimension of the vacuum, ω*. In the case of convergence to a fixed point we find δ= (d+σ−ω*)/(d−σ−ω*), with ω*?0 necessarily. Above...

Complex‐dimensional renormalization is defined for an arbitrary Feynman amplitude and shown to be equivalent to BPH renormalization. Using quantum electrodynamics as an example, Ward identities are proved; here Carlson's theorem extends the identities from integer to complex dimension. Both complex dimensional and analytic regularization are necessary at intermediate stages.