Double time window targeting technique: Real-time DMRG dynamics in the Pariser-Parr-Pople model

We present a generalized adaptive time-dependent density matrix renormalization-group (DMRG) scheme, called the double time window targeting (DTWT) technique, which gives accurate results with nominal computational resources, within reasonable computational time. This procedure originates from the amalgamation of the features of pace keeping DMRG algorithm, first proposed by Luo et al. [Phys. Rev. Lett. 91, 049701 (2003)] and the time-step targeting algorithm by Feiguin and White [Phys. Rev. B 72, 020404 (2005)]. Using the DTWT technique, we study the phenomena of spin-charge separation in conjugated polymers (materials for molecular electronics an spintronics), which have long-range electron-electron interactions and belong to the class of strongly correlated low-dimensional many-body systems. The issue of real-time dynamics within the Pariser-Parr-Pople (PPP) model which includes long-range electron correlations has not been addressed in the literature so far. The present study on PPP chains has revealed that, (i) long-range electron correlations enable both the charge and spin degree of freedom of the electron, to propagate faster in the PPP model compared to Hubbard model, (ii) for standard parameters of the PPP model as applied to conjugated polymers, the charge velocity is almost twice that of the spin velocity, and (iii) the simplistic interpretation of long-range correlations by merely renormalizing the U value of the Hubbard model fails to explain the dynamics of doped holes/electrons in the PPP model.

An electron paramagnetic resonance study of Pr0.6Ca0.4MnO3 across the charge-ordering transition

We report the first electron paramagnetic resonance studies of single crystals and powders of Pr0.6Ca0.4MnO3 in the 300-4.2 K range, covering the charge-ordering transition (Tco) at ~240 K and antiferromagnetic transition (TN) at ~170 K. The asymmetry parameter for the Dysonian single-crystal spectra shows an anomalous increase at Tco. Below Tco the g-value increases continuously, suggesting a gradual strengthening of the orbital ordering. The linewidth undergoes a sudden increase at Tco and continues to increase down to TN. The intensity increases as the temperature is decreased until Tco is reached, due to the renormalization of the magnetic susceptibility arising from the build-up of ferromagnetic correlations.

Isobaric electrical resistance along the critical line in nickel: An experimental test of universality

High-precision measurement of the electrical resistance of nickel along its critical line, a first attempt of this kind, as a function of pressure to 47.5 kbar is reported. Our analysis yields the values of the critical exponents α=α’=-0.115±0.005 and the amplitude ratios ‖A/A’‖=1.17±0.07 and ‖D/D’‖=1.2±0.1. These values are in close agreement with those predicted by renormalization-group (RG) theory. Moreover, this investigation provides an unambiguous experimental verification to one of the key consequences of RG theory that the critical exponents and amplitudes ratios are insensitive to pressure variation in nickel, a Heisenberg ferromagnet.

Ordering near the percolation threshold in models of two-dimensional interacting bosons with quenched dilution

Randomly diluted quantum boson and spin models in two dimensions combine the physics of classical percolation with the well-known dimensionality dependence of ordering in quantum lattice models. This combination is rather subtle for models that order in two dimensions but have no true order in one dimension, as the percolation cluster near threshold is a fractal of dimension between 1 and 2: two experimentally relevant examples are the O(2) quantum rotor and the Heisenberg antiferromagnet. We study two analytic descriptions of the O(2) quantum rotor near the percolation threshold. First a spin-wave expansion is shown to predict long-ranged order, but there are statistically rare points on the cluster that violate the standard assumptions of spin-wave theory. A real-space renormalization group (RSRG) approach is then used to understand how these rare points modify ordering of the O(2) rotor. A new class of fixed points of the RSRG equations for disordered one-dimensional bosons is identified and shown to support the existence of long-range order on the percolation backbone in two dimensions. These results are relevant to experiments on bosons in optical lattices and superconducting arrays, and also (qualitatively) for the diluted Heisenberg antiferromagnet La-2(Zn,Mg)(x)Cu1-xO4.

Mooij correlation in disordered metals

We point out that the Mooij correlation follows naturally from a dynamically disordered tight-binding Hamiltonian with random modulations of both the diagonal and the off-diagonal matrix elements which are known to act in opposition. The dynamic disorder is treated exactly while the static disorder is incorporated approximately as an effective additional time-dependent disorder affecting the diffusive electron. Such a time translation of static disorder is known to manifest itself in certain limits as a renormalization of the diffusion coefficient. The calculated conductivity exhibits the Mooij correlation at high temperatures, where quantum coherence associated with the static disorder can be ignored.

Spatially dependent correlation functions in the Anderson model

We report results from a first principles calculation of spatially dependent correlation functions around a magnetic impurity in metals described by the nondegenerate Anderson model. Our computations are based on a combination of perturbative scaling theory and numerical renormalization group methods. Results for the conduction election charge density around the impurity and correlation functions involving the conduction electron and impurity charge and spin densities will be presented. The behavior in various regimes including the mixed valent regime will be explored.

Metal-Insulator Transition in the Hubbard Model on a Triangular Lattice

We discuss the results of an extensive mean-field investigation of the half-filled Hubbard model on a triangular lattice at zero temperature. At intermediate U we find a first-order metal-insulator transition from an incommensurate spiral magnetic metal to a semiconducting state with a commensurate linear spin density wave ordering stabilized by the competition between the kinetic energy and the frustrated nature of the magnetic interaction. At large U the ground state is that of a classical triangular antiferromagnet within our approximation. In the incommensurate spiral metallic phase the Fermi surface has parts in which the wave function renormalization Z is extremely small. The evolution of the Fermi surface and the broadening of the quasi-particle band along with the variation of the plasma frequency and a charge stiffness constant with U/t are discussed.

Infrared divergences and confinement in massive (2+1)-dimensional QED

In (2+1)-dimensional quantum electrodynamics with massless photons and massive matter fields, it is shown that the mass renormalization of the latter is infrared divergent at one loop. This result remains unchanged at two loops. A simple argument based on a similar divergence of the Coulomb potential leads us to conjecture that charged states are not observable in this model. This argument holds in 1+1 dimensions also.

Probing Single and Bilayer Graphene Field Effect Transistors By Raman Spectroscopy

This article is a review of our work related to Raman studies of single layer and bilayer graphenes as a function Fermi level shift achieved by electrochemically top gating a field effect transistor. Combining the transport and in situ Raman studies of the field effect devices, a quantitative understanding is obtained of the phonon renormalization due to doping of graphene. Results are discussed in the light of time dependent perturbation theory, with electron phonon coupling parameter as an input from the density functional theory. It is seen that phonons near and Gamma and K points of the Brillouin zone are renormalized very differently by doping. Further, Gamma-phonon renormalization is different in bilayer graphene as compared to single layer, originating from their different electronic band structures near the zone boundary K-point. Thus Raman spectroscopy is not only a powerful probe to characterize the number of layers and their quality in a graphene sample, but also to quantitatively evaluate electron phonon coupling required to understand the performance of graphene devices.

Supersolid Phase in One-Dimensional Bose-Hubbard Model With Extended Range Interactions: DMRG And Field Theoretic Study at Different Densities

We use the Density Matrix Renormalization Group and the Abelian bosonization method to study the effect of density on quantum phases of one-dimensional extended Bose-Hubbard model. We predict the existence of supersolid phase and also other quantum phases for this system. We have analyzed the role of extended range interaction parameters on solitonic phase near half-filling. We discuss the effects of dimerization in nearest neighbor hopping and interaction as well as next nearest neighbor interaction on the plateau phase at half-filling.

Effect of dynamical asymmetry on the viscosity of a random copolymer melt

The variation of the viscosity as a function of the sequence distribution in an A-B random copolymer melt is determined. The parameters that characterize the random copolymer are the fraction of A monomers f, the parameter lambda which determines the correlation in the monomer identities along a chain and the Flory chi parameter chi(F) which determines the strength of the enthalpic repulsion between monomers of type A and B. For lambda>0, there is a greater probability of finding like monomers at adjacent positions along the chain, and for lambda<0 unlike monomers are more likely to be adjacent to each other. The traditional Markov model for the random copolymer melt is altered to remove ultraviolet divergences in the equations for the renormalized viscosity, and the phase diagram for the modified model has a binary fluid type transition for lambda>0 and does not exhibit a phase transition for lambda<0. A mode coupling analysis is used to determine the renormalization of the viscosity due to the dependence of the bare viscosity on the local concentration field. Due to the dissipative nature of the coupling. there are nonlinearities both in the transport equation and in the noise correlation. The concentration dependence of the transport coefficient presents additional difficulties in the formulation due to the Ito-Stratonovich dilemma, and there is some ambiguity about the choice of the concentration to be used while calculating the noise correlation. In the Appendix, it is shown using a diagrammatic perturbation analysis that the Ito prescription for the calculation of the transport coefficient, when coupled with a causal discretization scheme, provides a consistent formulation that satisfies stationarity and the fluctuation dissipation theorem. This functional integral formalism is used in the present analysis, and consistency is verified for the present problem as well. The upper critical dimension for this type of renormaliaation is 2, and so there is no divergence in the viscosity in the vicinity of a critical point. The results indicate that there is a systematic dependence of the viscosity on lambda and chi(F). The fluctuations tend to increase the viscosity for lambda<0, and decrease the viscosity for lambda>0, and an increase in chi(F) tends to decrease the viscosity. (C) 1996 American Institute of Physics.

Spatially dependent zero-frequency response functions and correlation functions in the Kondo model

We present the details of a formalism for calculating spatially varying zero-frequency response functions and equal-time correlation functions in models of magnetic and mixed-valence impurities of metals. The method is based on a combination of perturbative, thermodynamic scaling theory [H. R. Krishna-murthy and C. Jayaprakash, Phys. Rev. B 30, 2806 (1984)] and a nonperturbative technique such as the Wilson renormalization group. We illustrate the formalism for the spin-1/2 Kondo problem and present results for the conduction-spin-density�impurity-spin correlation function and conduction-electron charge density near the impurity. We also discuss qualitative features that emerge from our calculations and discuss how they can be carried over to the case of realistic models for transition-metal impurities.

The concentration dependent cooperative friction coefficient of dilute polymer solutions at the theta point

We combine multiple scattering and renormalization group methods to calculate the leading order dimensionless virial coefficient k(s) for the friction coefficient of dilute polymer solutions under conditions where the osmotic second virial coefficient vanishes (i.e., at the theta point T-theta). Our calculations are formulated in terms of coupled kinetic equations for the polymer and solvent, in which the polymers are modeled as continuous chains whose configurations evolve under the action of random forces in, the velocity field of the solvent. To lowest order in epsilon=4-d, we find that k(s) = 1.06. This result compares satisfactorily with existing experimental estimates of k(s), which are in the range 0.7-0.8. It is also in good agreement with other theoretical results on chains and suspensions at T-theta. Our calculated k(s) is also found to be identical to the leading order virial coefficient of the tracer friction coefficient at the theta point. We discuss possible reasons for the difficulties encountered when attempting to evaluate k(s) by extrapolating prior renormalization group calculations from semidilute concentrations to the infinitely dilute limit. (C) 1996 American Institute of Physics.

Quantum-confinement effects on the ordering of the lowest-lying excited states in conjugated chains

The symmetrized density-matrix renormalization-group approach is applied within the extended Hubbard-Peierls model (with parameters U/t, V/t, and bond alternation delta) to study the ordering of the lowest one-photon (1(1)B(u)(-)) and two-photon (2(1)A(g)(+)) states in one-dimensional conjugated systems with chain lengths N up to N = 80 sites. Three different types of crossovers are studied, as a function of U/t, delta, and N. The ''U crossover'' emphasizes the larger ionic character of the 2A(g) state compared to the lowest triplet excitation. The ''delta crossover'' shows strong dependence on both N and U/t. the ''N crossover'' illustrates the more localized nature of the 2A(g) excitation relative to the 1B(u) excitation at intermediate correlation strengths.

Spin distributions in antiferromagnetically coupled Mn2+Cu2+ systems: from the pair to the infinite chain

Probably the most informative description of the ground slate of a magnetic molecular species is provided by the spin density map. Such a map may be experimentally obtained from polarized neutron diffraction (PND) data or theoretically calculated using quantum chemical approaches. Density functional theory (DFT) methods have been proved to be well-adapted for this. Spin distributions in one-dimensional compounds may also be computed using the density matrix renormalization group (DMRG) formalism. These three approaches, PND, DFT, and DMRG, have been utilized to obtain new insights on the ground state of two antiferromagnetically coupled Mn2+Cu2+ compounds, namely [Mn(Me-6-[14]ane-N-4)Cu(oxpn)](CF3SO3)(2) and MnCu(pba)(H2O)(3) . 2H(2)O, with Me-6-[14]ane-N-4 = (+/-)-5,7,7,12,14,14-hexamethyl-1,4,8,11-tetraazacyclotetradecane, oxpn = N,N'-bis(3-aminopropyl)oxamido and pba = 1,3-propylenebis(oxamato). Three problems in particular have been investigated: the spin distribution in the mononuclear precursors [Cu(oxpn)] and [Cu(pba)](2-), the spin density maps in the two Mn2+Cu2+ compounds, and the evolution of the spin distributions on the Mn2+ and Cu2+ sites when passing from a pair to a one-dimensional ferrimagnet.