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.

Scaling of the entanglement spectrum near quantum phase transitions

The entanglement spectrum describing quantum correlations in many-body systems has been recently recognized as a key tool to characterize different quantum phases, including topological ones. Here we derive its analytically scaling properties in the vicinity of some integrable quantum phase transitions and extend our studies also to nonintegrable quantum phase transitions in one-dimensional spin models numerically. Our analysis shows that, in all studied cases, the scaling of the difference between the two largest nondegenerate Schmidt eigenvalues yields with good accuracy critical points and mass scaling exponents.

Entanglement entropy dynamics of Heisenberg chains

By means of the time dependent density matrix renormalization group algorithm we study the zero-temperature dynamics of the Von Neumann entropy of a block of spins in a Heisenberg chain after a sudden quench in the anisotropy parameter. In the absence of any disorder the block entropy increases linearly with time and then saturates. We analyse the velocity of propagation of the entanglement as a function of the initial and final anisotropies and compare our results, wherever possible, with those obtained by means of conformal field theory. In the disordered case we find a slower ( logarithmic) evolution which may signal the onset of entanglement localization.

Phase diagram of spin-1 bosons on one-dimensional lattices

Spinor Bose condensates loaded in optical lattices have a rich phase diagram characterized by different magnetic order. Here we apply the density matrix renormalization group to accurately determine the phase diagram for spin-1 bosons loaded on a one-dimensional lattice. The Mott lobes present an even or odd asymmetry associated to the boson filling. We show that for odd fillings the insulating phase is always in a dimerized state. The results obtained in this work are also relevant for the determination of the ground state phase diagram of the S=1 Heisenberg model with biquadratic interaction.

Bilinear-biquadratic spin-1 chain undergoing quadratic Zeeman effect

The Heisenberg model for spin-1 bosons in one dimension presents many different quantum phases, including the famous topological Haldane phase. Here we study the robustness of such phases in front of a SU(2) symmetry-breaking field as well as the emergence of unique phases. Previous studies have analyzed the effect of such uniaxial anisotropy in some restricted relevant points of the phase diagram. Here we extend those studies and present the complete phase diagram of the spin-1 chain with uniaxial anisotropy. To this aim, we employ the density-matrix renormalization group together with analytical approaches. The complete phase diagram can be realized using ultracold spinor gases in the Mott insulator regime under a quadratic Zeeman effect.

Strong electronic correlation in the hydrogen chain: A variational Monte Carlo study

In this paper, we report a fully ab initio variational Monte Carlo study of the linear and periodic chain of hydrogen atoms, a prototype system providing the simplest example of strong electronic correlation in low dimensions. In particular, we prove that numerical accuracy comparable to that of benchmark density-matrix renormalization-group calculations can be achieved by using a highly correlated Jastrow-antisymmetrized geminal power variational wave function. Furthermore, by using the so-called "modern theory of polarization" and by studying the spin-spin and dimer-dimer correlations functions, we have characterized in detail the crossover between the weakly and strongly correlated regimes of this atomic chain. Our results show that variational Monte Carlo provides an accurate and flexible alternative to highly correlated methods of quantum chemistry which, at variance with these methods, can be also applied to a strongly correlated solid in low dimensions close to a crossover or a phase transition.

Full characterization of the quantum linear-zigzag transition in atomic chains

A string of repulsively interacting particles exhibits a phase transition to a zigzag structure, by reducing the transverse trap potential or the interparticle distance. Based on the emergent symmetry Z2 it has been argued that this instability is a quantum phase transition, which can be mapped to an Ising model in transverse field. An extensive Density Matrix Renormalization Group analysis is performed, resulting in an high-precision evaluation of the critical exponents and of the central charge of the system, confirming that the quantum linear-zigzag transition belongs to the critical Ising model universality class. Quantum corrections to the classical phase diagram are computed, and the range of experimental parameters where quantum effects play a role is provided. These results show that structural instabilities of one-dimensional interacting atomic arrays can simulate quantum critical phenomena typical of ferromagnetic systems.

Nonequilibrium critical scaling from quantum thermodynamics

The emerging field of quantum thermodynamics is contributing important results and insights into archetypal many-body problems, including quantum phase transitions. Still, the question whether out-of-equilibrium quantities, such as fluctuations of work, exhibit critical scaling after a sudden quench in a closed system has remained elusive. Here, we take a novel approach to the problem by studying a quench across an impurity quantum critical point. By performing density matrix renormalization group computations on the two-impurity Kondo model, we are able to establish that the irreversible work produced in a quench exhibits finite-size scaling at quantum criticality. This scaling faithfully predicts the equilibrium critical exponents for the crossover length and the order parameter of the model, and, moreover, implies a new exponent for the rescaled irreversible work. By connecting the irreversible work to the two-impurity spin correlation function, our findings can be tested experimentally.

Multi-zone TAP-reactors theory and application - IV. Ideal and non-ideal boundary conditions

We study the influence of non-ideal boundary and initial conditions (BIC) of a temporal analysis of products (TAP) reactor model on the data (observed exit flux) analysis. The general theory of multi-response state-defining experiments for a multi-zone TAP reactor is extended and applied to model several alternative boundary and initial conditions proposed in the literature. The method used is based on the Laplace transform and the transfer matrix formalism for multi-response experiments. Two non-idealities are studied: (1) the inlet pulse not being narrow enough (gas pulse not entering the reactor in Dirac delta function shape) and (2) the outlet non-ideality due to imperfect vacuum. The effect of these non-idealities is analyzed to the first and second order of approximation. The corresponding corrections were obtained and discussed in detail. It was found that they are negligible. Therefore, the model with ideal boundary conditions is proven to be completely adequate to the description and interpretation of transport-reaction data obtained with TAP-2 reactors.

Cavity cooling of atoms: Within and without a cavity

We compare the efficiencies of two optical cooling schemes, where a single particle is either inside or outside an optical cavity, under experimentally-realisable conditions. We evaluate the cooling forces using the general solution of a transfer matrix method for a moving scatterer inside a general one-dimensional system composed of immobile optical elements. Assuming the same atomic saturation parameter, we find that the two cooling schemes provide cooling forces and equilibrium temperatures of comparable magnitude.

Scattering theory of multilevel atoms interacting with arbitrary radiation fields

We present a generic transfer matrix approach for the description of the interaction of atoms possessing multiple ground state and excited state sublevels with light fields. This model allows us to treat multi-level atoms as classical scatterers in light fields modified by, in principle, arbitrarily complex optical components such as mirrors, resonators, dispersive or dichroic elements, or filters. We verify our formalism for two prototypical sub-Doppler cooling mechanisms and show that it agrees with the standard literature.

Optical Cooling Using the Dipole Force

The term `laser cooling' is applied to the use of optical means to cool the motional energies of either atoms and molecules, or micromirrors. In the literature, these two strands are kept largely separate; both, however suffer from severe limitations. Laser cooling of atoms and molecules largely relies on the internal level structure of the species being cooled. As a result, only a small number of elements and a tiny number of molecules can be cooled this way. In the case of micromirrors, the problem lies in the engineering of micromirrors that need to satisfy a large number of constraints---these include a high mechanical Q-factor, high reflectivity and very good optical quality, weak coupling to the substrate, etc.---in order to enable efficient cooling. During the course of this thesis, I will draw these two sides of laser cooling closer together by means of a single, generically applicable scattering theory that can be used to explain the interaction between light and matter at a very general level. I use this `transfer matrix' formalism to explore the use of the retarded dipole--dipole interaction as a means of both enhancing the efficiency of micromirror cooling systems and rendering the laser cooling of atoms and molecules less species selective. In particular, I identify the `external cavity cooling' mechanism, whereby the use of an optical memory in the form of a resonant element (such as a cavity), outside which the object to be cooled sits, can potentially lead to the construction of fully integrated optomechanical systems and even two-dimensional arrays of translationally cold atoms, molecules or even micromirrors.

Ultrafast nonlinearities of minibands in metallodielectric Bragg resonators

The nonlinear properties of metallodielectric DBRs are investigated via optical pump-probe techniques using a widely tunable, dual-colour, high-repetition rate, ultrafast setup. As a consequence of the Bragg-arranged multilayers, the electric field penetrates to different depths of the nanostructure at different excitation resonances, strongly enhancing the intrinsic nonlinear response of the metal in comparison with bulk films. The analyzed spectral response of these structures reveals how their nonlinear behavior is dominated by the pump-induced modification of the metal dielectric function. Fitting the simulated changes of the optical resonances using transfer-matrix methods matches experiment well, and shows the key effects of the spectral dependence of the spatial mode profiles.

Dynamics of the entanglement spectrum in spin chains

We study the dynamics of the entanglement spectrum, that is the time evolution of the eigenvalues of the reduced density matrices after a bipartition of a one-dimensional spin chain. Starting from the ground state of an initial Hamiltonian, the state of the system is evolved in time with a new Hamiltonian. We consider both instantaneous and quasi adiabatic quenches of the system Hamiltonian across a quantum phase transition. We analyse the Ising model that can be exactly solved and the XXZ for which we employ the time-dependent density matrix renormalisation group algorithm. Our results show once more a connection between the Schmidt gap, i.e. the difference of the two largest eigenvalues of the reduced density matrix and order parameters, in this case the spontaneous magnetisation.

Nonclassicality and criticality in symmetry protected magnetic phases

Quantum and global discord in a spin-1 Heisenberg chain subject to single-ion anisotropy (uniaxial field) are studied using exact diagonalisation and the density matrix renormalisation group (DMRG). We find that these measures of quantum non-classicality are able to detect the quantum phase transitions confining the symmetry protected Haldane phase and show critical scaling with universal exponents. Moreover, in the case of thermal states, we find that quantum discord can increase with increasing temperature.