DEVELOPMENT AND APPLICATIONS OF TIME-DEPENDENT DENSITY MATRIX FUNCTIONAL THEORY

En esta tesis estudiamos las teorías sobre la Matriz Densidad Reducida (MDR) como un marco prometedor. Nos enfocamos sobre esta teorías desde dos aspectos: Primero, usamos algunos modelos sencillos hechos con dos partículas las cuales estan armónicamente confinadas como una base para ilustrar la utilidad de la matriz densidad. Para tales sistemas, usamos la MDR de un cuerpo para calcular algunas cantidades de interés tales como densidad de momentum. Posteriormente obtenemos los orbitales naturales y su número de ocupación para algunos de los modelos, y en uno de los casos expresamos la MDR de dos cuerpos de manera exacta en términos de la MDR de un cuerpo. También usamos el teorema diferencial del virial para establecer una descripción unificada de la familia entera de estos sistemas modelo en términos de la densidad. En la seguna parte cambiamos a casos fuera del equilibrio y analizamos la así llamada jerarquía BBGKY de ecuaciones para describir la evolución temporal de un sistema de muchos cuerpos en términos de sus MDRs (a todos los órdenes). Proveemos un exhaustivo estudio de los desafíos y problemas abiertos ligados a la truncación de tales jerarquías de ecuaciones para hacerlas aplicables. Restringimos nuestro análisis a la evolución acoplada de la MDR de uno y dos cuerpos, donde los efectos de correlación de alto orden estan embebidos dentro de la aproximación usada para cerrar las ecuaciones. Probamos que dentro de esta aproximación, el número de electrones y la energía total se conservan, sin importar la aproximación usada. Luego, demostramos que aplicando los esquemas de truncación de estado base para llevar los electrones a comportamientos indeseables y no físicos, tales como la violación e incluso la divergencia en la densidad electrónica local, tanto en regímenes correlacionados débiles y fuertes.

Coherence loss and recovery of an electron spin coupled inhomogeneously to a one-dimensional interacting spin bath: An adaptive time-dependent density-matrix renormalization group study

Coherence evolution and echo effect of an electron spin, which is coupled inhomogeneously to an interacting one-dimensional finite spin bath via hyperfine-type interaction, are studied using the adaptive time-dependent density-matrix renormalization group method. It is found that the interplay of the coupling inhomogeneity and the transverse intrabath interactions results in two qualitatively different coherence evolutions, namely, a coherence-preserving evolution characterized by periodic oscillation and a complete decoherence evolution. Correspondingly, the echo effects induced by an electron-spin flip at time tau exhibit stable recoherence pulse sequence for the periodic evolution and a single peak at root 2 tau for the decoherence evolution, respectively. With the diagonal intrabath interaction included, the specific feature of the periodic regime is kept, while the root 2 tau-type echo effect in the decoherence regime is significantly affected. To render the experimental verifications possible, the Hahn echo envelope as a function of tau is calculated, which eliminates the inhomogeneous broadening effect and serves for the identification of the different status of the dynamic coherence evolution, periodic versus decoherence.

Density matrix Loschmidt echo and quantum phase transitions

We introduce the concept of the Loschmidt echo (LE) to the space of the reduced density matrix of spin and fermionic systems to study the density matrix LEs (DMLEs) of the one-dimensional extended Hubbard model and the transverse field Ising model. Our results show that the DMLEs are remarkably influenced by the criticality of the system, and the method is a convenient way to study quantum phase transitions.

Correlation effects for semiconducting single-wall carbon nanotubes: A density matrix renormalization group study

In this paper we report the applicability of the density matrix renormalization group (DMRG) approach to the cylindrical single wall carbon nanotube (SWCN) for the purpose of its correlation effect. By applying the DMRG approach to the t+U+V model, with t and V being the hopping and Coulomb energies between the nearest neighboring sites, respectively, and U the on-site Coulomb energy, we calculate the phase diagram for the SWCN with chiral numbers (n(1)=3, n(2)=2), which reflects the competition between the correlation energy U and V. Within reasonable parameter ranges, we investigate possible correlated ground states, the lowest excitations, and the corresponding correlation functions in which the connection with the excitonic insulator is particularly addressed.

Resonant two-photon cross-section calculations for transitions to 4f5d state in Pr3+: YAG using density matrix theory

The density matrix resonant two-photon absorption (TPA) theory is applied to a rare-earth ion-doped laser crystal. TPA cross sections for transitions from the ground state to the first 4f5d state in Pr3+:YAG are calculated. The results indicate the density matrix TPA theory is attractive in studying TPA in laser crystals. (C) 2000 Elsevier Science B.V. All rights reserved.

Calculations of cross sections of resonant two-photon transitions to the second excited 4f5d state in Pr3+: Y3Al5O12 using density matrix theory

The density matrix resonant two-photon absorption (TPA) theory applicable to laser crystals doped with rare earth ions is described. Using this theory, resonant TPA cross sections for transitions from the ground state to the second excited state of the 4f5d configuration in cm(4)s Pr3+:Y3Al5O12 are calculated. The peak value of TPA cross section calculated is 2.75 x 10(-50) cm(4)s which is very close to the previous experimental value 4 x 10(-50) cm(4) s. The good agreement of calculated data with measured values demonstrates that the density matrix resonant TPA theory can predict resonant TPA intensity much better than the standard second-order perturbation TPA theory.

Tensor density matrix theory for determining population and alignment of symmetric top molecule using rotationally unresolved fluorescence

General expressions used for transforming raw laser-induced fluorescence (LIF) intensity into the population and alignment parameters of a symmetric top molecule are derived by employing the density matrix approach. The molecular population and alignment are described by molecular state multipoles. The results are presented for a general excitation-detection geometry and then applied to some special geometries. In general cases, the LIF intensity is a complex function of the initial molecular state multipoles, the dynamic factors and the excitation-detection geometrical factors. It contains a population and 14 alignment multipoles. How to extract all initial state multipoles from the rotationally unresolved emission LIF intensity is discussed in detail.

Iterative Linearized Density Matrix Propagation for Modeling Coherent Energy Transfer in Photosynthetic Light Harvesting

We present results of calculations [1] that employ a new mixed quantum classical iterative density matrix propagation approach (ILDM , or so called Is‐Landmap) [2] to explore the survival of coherence in different photo synthetic models. Our model studies confirm the long lived quantum coherence , while conventional theoretical tools (such as Redfield equation) fail to describe these phenomenon [3,4]. Our ILDM method is a numerical exactly propagation scheme and can be served as a bench mark calculation tools[2]. Result get from ILDM and from other recent methods have been compared and show agreement with each other[4,5]. Long lived coherence plateau has been attribute to the shift of harmonic potential due to the system bath interaction, and the harvesting efficiency is a balance between the coherence and dissipation[1]. We use this approach to investigate the excitation energy transfer dynamics in various light harvesting complex include Fenna‐Matthews‐Olsen light harvesting complex[1] and Cryptophyte Phycocyanin 645 [6]. [1] P.Huo and D.F.Coker ,J. Chem. Phys. 133, 184108 (2010) . [2] E.R. Dunkel, S. Bonella, and D.F. Coker, J. Chem. Phys. 129, 114106 (2008). [3] A. Ishizaki and G.R. Fleming, J. Chem. Phys. 130, 234111 (2009). [4] A. Ishizaki and G.R. Fleming, Proc. Natl. Acad. Sci. 106, 17255 (2009). [5] G. Tao and W.H. Miller, J. Phys. Chem. Lett. 1, 891 (2010). [6] P.Huo and D.F.Coker in preparation

Density Matrix Renormalization Group for Dummies

We describe the Density Matrix Renormalization Group algorithms for time dependent and time independent Hamiltonians. This paper is a brief but comprehensive introduction to the subject for anyone willing to enter in the field or write the program source code from scratch. An open source version of the code can be found at: http://www.dmrg.it.

Density matrix functional theory that includes pairing correlations

The extension of density functional theory (DFT) to include pairing correlations without formal violation of the particle-number conservation condition is described. This version of the theory can be considered as a foundation of the application of existing DFT plus pairing approaches to atoms, molecules, ultracooled and magnetically trapped atomic Fermi gases, and atomic nuclei where the number of particles is conserved exactly. The connection with Hartree-Fock-Bogoliubov (HFB) theory is discussed, and the method of quasilocal reduction of the nonlocal theory is also described. This quasilocal reduction allows equations of motion to be obtained which are much simpler for numerical solution than the equations corresponding to the nonlocal case. Our theory is applied to the study of some even Sn isotopes, and the results are compared with those obtained in the standard HFB theory and with the experimental ones.

Development and applications of density matrix functional theory of generalized Hubbard Models

In dieser Doktorarbeit wird eine akkurate Methode zur Bestimmung von Grundzustandseigenschaften stark korrelierter Elektronen im Rahmen von Gittermodellen entwickelt und angewandt. In der Dichtematrix-Funktional-Theorie (LDFT, vom englischen lattice density functional theory) ist die Ein-Teilchen-Dichtematrix γ die fundamentale Variable. Auf der Basis eines verallgemeinerten Hohenberg-Kohn-Theorems ergibt sich die Grundzustandsenergie Egs[γgs] = min° E[γ] durch die Minimierung des Energiefunktionals E[γ] bezüglich aller physikalischer bzw. repräsentativer γ. Das Energiefunktional kann in zwei Beiträge aufgeteilt werden: Das Funktional der kinetischen Energie T[γ], dessen lineare Abhängigkeit von γ genau bekannt ist, und das Funktional der Korrelationsenergie W[γ], dessen Abhängigkeit von γ nicht explizit bekannt ist. Das Auffinden präziser Näherungen für W[γ] stellt die tatsächliche Herausforderung dieser These dar. Einem Teil dieser Arbeit liegen vorausgegangene Studien zu Grunde, in denen eine Näherung des Funktionals W[γ] für das Hubbardmodell, basierend auf Skalierungshypothesen und exakten analytischen Ergebnissen für das Dimer, hergeleitet wird. Jedoch ist dieser Ansatz begrenzt auf spin-unabhängige und homogene Systeme. Um den Anwendungsbereich von LDFT zu erweitern, entwickeln wir drei verschiedene Ansätze zur Herleitung von W[γ], die das Studium von Systemen mit gebrochener Symmetrie ermöglichen. Zuerst wird das bisherige Skalierungsfunktional erweitert auf Systeme mit Ladungstransfer. Eine systematische Untersuchung der Abhängigkeit des Funktionals W[γ] von der Ladungsverteilung ergibt ähnliche Skalierungseigenschaften wie für den homogenen Fall. Daraufhin wird eine Erweiterung auf das Hubbardmodell auf bipartiten Gittern hergeleitet und an sowohl endlichen als auch unendlichen Systemen mit repulsiver und attraktiver Wechselwirkung angewandt. Die hohe Genauigkeit dieses Funktionals wird aufgezeigt. Es erweist sich jedoch als schwierig, diesen Ansatz auf komplexere Systeme zu übertragen, da bei der Berechnung von W[γ] das System als ganzes betrachtet wird. Um dieses Problem zu bewältigen, leiten wir eine weitere Näherung basierend auf lokalen Skalierungseigenschaften her. Dieses Funktional ist lokal bezüglich der Gitterplätze formuliert und ist daher anwendbar auf jede Art von geordneten oder ungeordneten Hamiltonoperatoren mit lokalen Wechselwirkungen. Als Anwendungen untersuchen wir den Metall-Isolator-Übergang sowohl im ionischen Hubbardmodell in einer und zwei Dimensionen als auch in eindimensionalen Hubbardketten mit nächsten und übernächsten Nachbarn. Schließlich entwickeln wir ein numerisches Verfahren zur Berechnung von W[γ], basierend auf exakten Diagonalisierungen eines effektiven Vielteilchen-Hamilton-Operators, welcher einen von einem effektiven Medium umgebenen Cluster beschreibt. Dieser effektive Hamiltonoperator hängt von der Dichtematrix γ ab und erlaubt die Herleitung von Näherungen an W[γ], dessen Qualität sich systematisch mit steigender Clustergröße verbessert. Die Formulierung ist spinabhängig und ermöglicht eine direkte Verallgemeinerung auf korrelierte Systeme mit mehreren Orbitalen, wie zum Beispiel auf den spd-Hamilton-Operator. Darüber hinaus berücksichtigt sie die Effekte kurzreichweitiger Ladungs- und Spinfluktuationen in dem Funktional. Für das Hubbardmodell wird die Genauigkeit der Methode durch Vergleich mit Bethe-Ansatz-Resultaten (1D) und Quanten-Monte-Carlo-Simulationen (2D) veranschaulicht. Zum Abschluss wird ein Ausblick auf relevante zukünftige Entwicklungen dieser Theorie gegeben.

Off-diagonal helicity density matrix elements for vector mesons produced in polarized e+e- processes

Final-state qq̄ interactions give origin to nonzero values of the off-diagonal element ρ1,-1 of the helicity density matrix of vector mesons produced in e+e- annihilations, as has been confirmed by recent OPAL data on φ, D*, and K*. New predictions are given for ρ1,-1 of several mesons produced at large XE and small pT - i.e., collinear with the parent jet - in the annihilation of polarized e+ and e-; the results depend strongly on the elementary dynamics and allow further nontrivial tests of the standard model.

Quantum master equation scheme of time-dependent density functional theory to time-dependent transport in nanoelectronic devices

In this work a practical scheme is developed for the first-principles study of time-dependent quantum transport. The basic idea is to combine the transport master equation with the well-known time-dependent density functional theory. The key ingredients of this paper include (i) the partitioning-free initial condition and the consideration of the time-dependent bias voltages which base our treatment on the Runge-Gross existence theorem; (ii) the non-Markovian master equation for the reduced (many-body) central system (i.e., the device); and (iii) the construction of Kohn-Sham master equations for the reduced single-particle density matrix, where a number of auxiliary functions are introduced and their equations of motion (EOMs) are established based on the technique of spectral decomposition. As a result, starting with a well-defined initial state, the time-dependent transport current can be calculated simultaneously along with the propagation of the Kohn-Sham master equation and the EOMs of the auxiliary functions.

Scattering matrix approach to electronic dephasing in long-range electron transfer

Based on the Buttiker dephasing model, we propose an analytical scattering matrix approach to the long-range electron transfer phenomena. The present efficient scheme smoothly interpolates between the superexchange and the sequential hopping mechanisms. Various properties such as the drastic dephasing-assisted enhancement and turnover behaviors are demonstrated in good agreement with those obtained via the dynamical reduced density-matrix methods. These properties are further elucidated as results of the interplay among the dephasing strength, the tunneling parameter, and the bridge length of the electron transfer system. (C) 2001 American Institute of Physics.