Test for basis-set errors in relativistic Dirac-Fock-Slater calculations with a numerical AO-DFS-basis

A fully relativistic four-component Dirac-Fock-Slater program for diatomics, with numerically given AO's as basis functions is presented. We discuss the problem of the errors due to the finite basis-set, and due to the influence of the negative energy solutions of the Dirac Hamiltonian. The negative continuum contributions are found to be very small.

Electron-Nuclear Correlation in Laser-Driven Diatomic Molecules

The interaction of short intense laser pulses with atoms/molecules produces a multitude of highly nonlinear processes requiring a non-perturbative treatment. Detailed study of these highly nonlinear processes by numerically solving the time-dependent Schrodinger equation becomes a daunting task when the number of degrees of freedom is large. Also the coupling between the electronic and nuclear degrees of freedom further aggravates the computational problems. In the present work we show that the time-dependent Hartree (TDH) approximation, which neglects the correlation effects, gives unreliable description of the system dynamics both in the absence and presence of an external field. A theoretical framework is required that treats the electrons and nuclei on equal footing and fully quantum mechanically. To address this issue we discuss two approaches, namely the multicomponent density functional theory (MCDFT) and the multiconfiguration time-dependent Hartree (MCTDH) method, that go beyond the TDH approximation and describe the correlated electron-nuclear dynamics accurately. In the MCDFT framework, where the time-dependent electronic and nuclear densities are the basic variables, we discuss an algorithm to calculate the exact Kohn-Sham (KS) potentials for small model systems. By simulating the photodissociation process in a model hydrogen molecular ion, we show that the exact KS potentials contain all the many-body effects and give an insight into the system dynamics. In the MCTDH approach, the wave function is expanded as a sum of products of single-particle functions (SPFs). The MCTDH method is able to describe the electron-nuclear correlation effects as the SPFs and the expansion coefficients evolve in time and give an accurate description of the system dynamics. We show that the MCTDH method is suitable to study a variety of processes such as the fragmentation of molecules, high-order harmonic generation, the two-center interference effect, and the lochfrass effect. We discuss these phenomena in a model hydrogen molecular ion and a model hydrogen molecule. Inclusion of absorbing boundaries in the mean-field approximation and its consequences are discussed using the model hydrogen molecular ion. To this end, two types of calculations are considered: (i) a variational approach with a complex absorbing potential included in the full many-particle Hamiltonian and (ii) an approach in the spirit of time-dependent density functional theory (TDDFT), including complex absorbing potentials in the single-particle equations. It is elucidated that for small grids the TDDFT approach is superior to the variational approach.

A potential energy surface for the ground state of acetylene, C2H2

A potential energy function has been derived for the ground state surface of C2H2 as a many-body expansion. The 2- and 3-body terms have been obtained by preliminary investigation of the ground state surfaces of CH2( 3B1) and C2H( 2Σ+). A 4-body term has been derived which reproduces the energy, geometry and harmonic force field of C2H2. The potential has a secondary minimum corresponding to the vinylidene structure and the geometry and energy of this are in close agreement with predictions from ab initio calculations. The saddle point for the HCCH-H2CC rearrangement is predicted to lie 2•530 eV above the acetylene minimum.

Potential energy surface and MULTIMODE vibrational analysis of C2H3+

A full dimensional, ab initio-based semiglobal potential energy surface for C2H3+ is reported. The ab initio electronic energies for this molecule are calculated using the spin-restricted, coupled cluster method restricted to single and double excitations with triples corrections [RCCSD(T)]. The RCCSD(T) method is used with the correlation-consistent polarized valence triple-zeta basis augmented with diffuse functions (aug-cc-pVTZ). The ab initio potential energy surface is represented by a many-body (cluster) expansion, each term of which uses functions that are fully invariant under permutations of like nuclei. The fitted potential energy surface is validated by comparing normal mode frequencies at the global minimum and secondary minimum with previous and new direct ab initio frequencies. The potential surface is used in vibrational analysis using the "single-reference" and "reaction-path" versions of the code MULTIMODE. (c) 2006 American Institute of Physics.

SiO(2) in density functional theory and beyond

We present the first-principle electronic structure calculation on an amorphous material including many-body corrections within the GW approximation. We show that the inclusion of the local field effects in the exchange-correlation potential is crucial to quantitatively describe amorphous systems and defect states. We show that the mobility gap of amorphous silica coincides with the band gap of quartz, contrary to the traditional picture and the densityfunctional theory results. (C) 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

Reentrant quantum Hall effect in bilayer system at high filling factors

We report on the measurements of the quantum Hall effect states in double quantum well structures at the filling factors v = 4N + 1 and 4N + 3, where N is the Landau index number, in the presence of the in-plane magnetic field. The quantum Hall states at these filling factors vanish and reappear several times. Repeated reentrance of the transport gap occurs due to the periodic vanishing of the tunneling amplitude in the presence of the in-plane field. When the gap vanishes, the transport becomes anisotropic. The anisotropy persist at half-odd filling factors, when bilayer quantum Hall states are recovered with increase of the tilt angle. (C) 2007 Elsevier B.V. All rights reserved.

EXACT AND APPROXIMATE RELATIONS FOR THE SPIN-DEPENDENCE OF THE EXCHANGE ENERGY IN HIGH MAGNETIC FIELDS

The exchange energy of an arbitrary collinear-spin many-body system in an external magnetic field is a functional of the spin-resolved charge and current densities, E(x)[n(up arrow), n(down arrow), j(up arrow), j(down arrow)]. Within the framework of density-functional theory (DFT), we show that the dependence of this functional on the four densities can be fully reconstructed from either of two extreme limits: a fully polarized system or a completely unpolarized system. Reconstruction from the limit of an unpolarized system yields a generalization of the Oliver-Perdew spin scaling relations from spin-DFT to current-DFT. Reconstruction from the limit of a fully polarized system is used to derive the high-field form of the local-spin-density approximation to current-DFT and to magnetic-field DFT.

Friedel oscillations in one-dimensional metals: From Luttinger`s theorem to the Luttinger liquid

Charge density and magnetization density profiles of one-dimensional metals are investigated by two complementary many-body methods: numerically exact (Lanczos) diagonalization, and the Bethe-Ansatz local-density approximation with and without a simple self-interaction correction. Depending on the magnetization of the system, local approximations reproduce different Fourier components of the exact Friedel oscillations. (C) 2008 Elsevier B.V. All rights reserved.

Effective theory for trapped few-fermion systems

We apply the general principles of effective field theories to the construction of effective interactions suitable for few- and many-body calculations in a no-core shell model framework. We calculate the spectrum of systems with three and four two-component fermions in a harmonic trap. In the unitary limit, we find that three-particle results are within 10% of known semianalytical values even in small model spaces. The method is very general, and can be readily extended to other regimes, more particles, different species (e.g., protons and neutrons in nuclear physics), or more-component fermions (as well as bosons). As an illustration, we present calculations of the lowest-energy three-fermion states away from the unitary limit and find a possible inversion of parity in the ground state in the limit of trap size large compared to the scattering length. Furthermore, we investigate the lowest positive-parity states for four fermions, although we are limited by the dimensions we can currently handle in this case.

Transport properties and equation of state of quark-nuclear matter

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Intersubband excitations at finite temperatures and their roles in different quasi-one-dimensional systems

We theoretically study many-body excitations in three different quasi-one-dimensional (Q1D) electron systems: (i) those formed on the surface of liquid Helium; (ii) in two coupled semiconductor quantum wires; and (iii) Q1D electrons embedded in polar semiconductor-based quantum wires. Our results show intersubband coupling between higher subbands and the two lowest subbands affecting even the lower energy intersubband plasmons on the liquid Helium surface. Concerning the second system, we show a pronounced extra peak appearing in the intersubband impurity spectral function for temperatures as high as 20 K. We finally show coupled intersubband plasmon-phonon modes surviving for temperatures up to 300 K.

Influence of nearest-neighbor Coulomb interactions on the phase diagram of the ferromagnetic Kondo model

The influence of a nearest-neighbor Coulomb repulsion of strength V on the properties of the ferromagnetic Kondo model is analyzed using computational techniques. The Hamiltonian studied here is defined on a chain using localized S = 1/2 spins, and one orbital per site. Special emphasis is given to the influence of the Coulomb repulsion on the regions of phase separation recently discovered in this family of models, as well as on the double-exchange-induced ferromagnetic ground state. When phase separation dominates at V= 0, the Coulomb interaction breaks the large domains of the two competing phases into small islands of one phase embedded into the other. This is in agreement with several experimental results, as discussed in the text. Vestiges of the original phase separation regime are found in the spin structure factor as incommensurate peaks, even at large values of V. In the ferromagnetic regime close to density n = 0.5, the Coulomb interaction induces tendencies to charge ordering without altering the fully polarized character of the state. This regime of charge-ordered ferromagnetism may be related with experimental observations of a similar phase by Chen and Cheong [Phys. Rev. Lett. 76, 4042 (1996)]. Our results reinforce the recently introduced notion [see, e.g., S. Yunoki et al., Phys. Rev. Lett. 80, 845 (1998)] that in realistic models for manganites analyzed with unbiased many-body techniques, the ground state properties arise from a competition between ferromagnetism and phase-separation - charge-ordering tendencies. ©1999 The American Physical Society.

Dimensional effects on the momentum distribution of bosonic trimer states

The momentum distribution is a powerful probe of strongly interacting systems that are expected to display universal behavior. This is contained in the contact parameters which relate few- and many-body properties. Here we consider a Bose gas in two dimensions and explicitly show that the two-body contact parameter is universal and then demonstrate that the momentum distribution at next-to-leading order has a logarithmic dependence on momentum which is vastly different from the three-dimensional case. Based on this, we propose a scheme for measuring the effective dimensionality of a quantum many-body system by exploiting the functional form of the momentum distribution. © 2013 American Physical Society.

Excitonic instabilities and spontaneous time-reversal symmetry breaking on the honeycomb lattice

We elucidate the close relationship between spontaneous time-reversal symmetry breaking and the physics of excitonic instabilities in strongly correlated multiband systems. The underlying mechanism responsible for the spontaneous breaking of time-reversal symmetry in a many-body system is closely related to the Cooper-like pairing instability of interband particle-hole pairs involving higher-order symmetries. Studies of such pairing instabilities have, however, mainly focused on the mean-field aspects of the virtual exciton condensate, which ignores the presence of the underlying collective Fermi-liquid excitations. We show that this relationship can be exploited to systematically derive the coupling of the condensate order parameter to the intraband Fermi-liquid particle-hole excitations. Surprisingly, we find that the static susceptibility is negative in the ordered phase when the coupling to the Fermi-liquid collective excitations are included, suggesting that a uniform condensate of virtual excitons, with or without time-reversal breaking, is an unstable phase at T = 0.

Accurate ab initio potential energy surfaces for the (3)A '' and (3)A ' electronic states of the O(P-3) plus HBr system

In this work, we report the construction of potential energy surfaces for the (3)A '' and (3)A' states of the system O(P-3) + HBr. These surfaces are based on extensive ab initio calculations employing the MRCI+Q/CBS+SO level of theory. The complete basis set energies were estimated from extrapolation of MRCI+Q/aug-cc-VnZ(-PP) (n = Q, 5) results and corrections due to spin-orbit effects obtained at the CASSCF/aug-cc-pVTZ(-PP) level of theory. These energies, calculated over a region of the configuration space relevant to the study of the reaction O(P-3) + HBr -> OH + Br, were used to generate functions based on the many-body expansion. The three-body potentials were interpolated using the reproducing kernel Hilbert space method. The resulting surface for the (3)A '' electronic state contains van der Waals minima on the entrance and exit channels and a transition state 6.55 kcal/mol higher than the reactants. This barrier height was then scaled to reproduce the value of 5.01 kcal/mol, which was estimated from coupled cluster benchmark calculations performed to include high-order and core-valence correlation, as well as scalar relativistic effects. The (3)A' surface was also scaled, based on the fact that in the collinear saddle point geometry these two electronic states are degenerate. The vibrationally adiabatic barrier heights are 3.44 kcal/mol for the (3)A '' and 4.16 kcal/mol for the (3)A' state. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4705428]