Showing 1–50 of 50 results for author: Turkowski, V

Thickness dependence of superconductivity in FeSe films

Authors: Jia Shi, Duy Le, Volodymyr Turkowski, Naseem Ud Din, Tao Jiang, Qiang Gu, Talat S. Rahman

Abstract: The films of FeSe on substrates have attracted attention because of their unusually high-temperature (Tc) superconducting properties whose origins continue to be debated. To disentangle the competing effects of the substrate and interlayer and intralayer processes, we present here results of density functional theory (DFT)-based analysis of the electronic structure of unsupported FeSe films consis… ▽ More The films of FeSe on substrates have attracted attention because of their unusually high-temperature (Tc) superconducting properties whose origins continue to be debated. To disentangle the competing effects of the substrate and interlayer and intralayer processes, we present here results of density functional theory (DFT)-based analysis of the electronic structure of unsupported FeSe films consisting of 1 to 5 layers (1L-5L). Furthermore, by solving the Bardeen-Schrieffer-Cooper (BCS) equation with spin-wave exchange attraction derived from the Hubbard model, we find the superconducting critical temperature Tc for 1L-5L and bulk FeSe systems in reasonable agreement with experimental data. Our results point to the importance of correlation effects in superconducting properties of single- and multi-layer FeSe films, independently of the role of substrate. △ Less

Submitted 19 February, 2023; originally announced February 2023.

Dark exciton energy splitting in monolayer WSe2: insights from time-dependent density-functional theory

Authors: Jia Shi, Volodymyr Turkowski, Duy Le, Talat S. Rahman

Abstract: We present here a formalism based on time-dependent density-functional theory (TDDFT) to describe characteristics of both intra- and inter-valley excitons in semiconductors, the latter of which had remained a challenge. Through the usage of an appropriate exchange-correlation kernel (nanoquanta), we trace the energy difference between the intra- and inter-valley dark excitons in monolayer (1L) WSe… ▽ More We present here a formalism based on time-dependent density-functional theory (TDDFT) to describe characteristics of both intra- and inter-valley excitons in semiconductors, the latter of which had remained a challenge. Through the usage of an appropriate exchange-correlation kernel (nanoquanta), we trace the energy difference between the intra- and inter-valley dark excitons in monolayer (1L) WSe2 to the domination of the exchange part in the exchange-correlation energies of these states. Furthermore, our calculated transition contribution maps establish the momentum resolved weights of the electron-hole excitations in both bright and dark excitons thereby providing a comprehensive understanding of excitonic properties of 1L WSe2. We find that the states consist of hybridized excitations around the corresponding valleys which leads to brightening of the dark excitons, i.e., significantly decreasing their lifetime which is reflected in the PL spectrum. Using many-body perturbation theory, we calculate the phonon contribution to the energy bandgap and the linewidths of the excited electrons, holes and (bright) exciton to find that as the temperature increases the bandgap significantly decreases, while the linewidths increase. Our work paves for describing the ultrafast charge dynamics of transition metal dichalcogenide within an ab initio framework. △ Less

Submitted 16 February, 2023; originally announced February 2023.

Mapping Spin Interactions from Conductance Peak Splitting in Coulomb Blockade

Authors: Eric D. Switzer, Xiao-Guang Zhang, Volodymyr Turkowski, Talat S. Rahman

Abstract: We investigate the transport properties of a quantum dot coupled to leads interacting with a multi-spin system using the generalized master equation within the Coulomb blockade regime. We find that if two states for each scattering region electron manifold are included, several signatures of the interacting spin system appear in steady-state transport properties. We provide a theoretical mapping o… ▽ More We investigate the transport properties of a quantum dot coupled to leads interacting with a multi-spin system using the generalized master equation within the Coulomb blockade regime. We find that if two states for each scattering region electron manifold are included, several signatures of the interacting spin system appear in steady-state transport properties. We provide a theoretical mapping of differential conductance peak signatures and all spin Hamiltonian parameters related to the inclusion of excited state transitions between uncharged and charged electron manifolds. Our predictions describe a scheme of only using a quantum dot and differential conductance to measure magnetic anisotropy, inter-spin exchange coupling, exchange coupling between the spin system and itinerant electron, and applied magnetic field response. △ Less

Submitted 30 June, 2023; v1 submitted 12 January, 2023; originally announced January 2023.

Spatially non-homogeneous metallization of VO2: a TDDFT+DMFT analysis

Authors: Jose Mario Galicia-Hernandez, Volodymyr Turkowski, Gregorio Hernandez-Cocoletzi, Talat S. Rahman

Abstract: We provide insights into the atomistic details of the ultrafast spatially-resolved breakdown of the insulating M1 phase in bulk VO2 employing an ab initio technique based on time-dependent density-functional theory and dynamical mean-field theory (DMFT-TDDFT). We find that the system is initially metallized preferentially along the vanadium-dimer chains (CR axis), with a subsequent growth of CR-el… ▽ More We provide insights into the atomistic details of the ultrafast spatially-resolved breakdown of the insulating M1 phase in bulk VO2 employing an ab initio technique based on time-dependent density-functional theory and dynamical mean-field theory (DMFT-TDDFT). We find that the system is initially metallized preferentially along the vanadium-dimer chains (CR axis), with a subsequent growth of CR-elongated metallic bubbles. Moreover, we trace the breakdown of the insulating phase to two types of oxygen atoms, resulting from vanadium dimerization, which produce an unusual charge-density modulation in the oxygen-atom chains with significant charge transfer to the inter-dimer distance. These results are in qualitative agreement with experimental data and shed light on the interplay between valence charge and lattice structure and its role in the ultrafast response of strongly correlated insulators. △ Less

Submitted 30 November, 2020; originally announced November 2020.

Ultrafast charge dynamics and photoluminescence in bilayer MoS2

Authors: Naseem Ud Din, Volodymyr Turkowski, Talat S Rahman

Abstract: Our examination of the interplay of ultrafast charge dynamics and electron-phonon interaction in bilayer MoS2 provides a microscopic basis for understanding the features (two peaks) in the emission spectrum. We demonstrate that while the initial accumulation of excited charge occurs at and near the Q point of the two-dimensional Brillioun zone, emission takes place predominantly through two pathwa… ▽ More Our examination of the interplay of ultrafast charge dynamics and electron-phonon interaction in bilayer MoS2 provides a microscopic basis for understanding the features (two peaks) in the emission spectrum. We demonstrate that while the initial accumulation of excited charge occurs at and near the Q point of the two-dimensional Brillioun zone, emission takes place predominantly through two pathways: direct charge recombination at the K point and indirect phonon-assisted recombination of electrons at the K valley and holes at Γがんま hill of the Brillouin zone. Analysis of the wave vector dependencies of the electron-phonon interaction traces the higher energy peak to phonon-assisted relaxation of the excited electrons from the Q to the K valley in the conduction band. Our results thus reveal the importance of ultrafast charge dynamics in understanding photoemissive properties of a few-layer transition-metal dichalcogenide. These calculations are based on time dependent density functional theory in the density matrix formulation. △ Less

Submitted 27 October, 2020; originally announced October 2020.

Electron Thermalization and Relaxation in Laser-Heated Nickel by Few-Femtosecond Core-Level Transient Absorption Spectroscopy

Authors: Hung-Tzu Chang, Alexander Guggenmos, Scott K. Cushing, Yang Cui, Naseem Ud Din, Shree Ram Acharya, Ilana J. Porter, Ulf Kleineberg, Volodymyr Turkowski, Talat S. Rahman, Daniel M. Neumark, Stephen R. Leone

Abstract: Direct measurements of photoexcited carrier dynamics in nickel are made using few-femtosecond extreme ultraviolet (XUV) transient absorption spectroscopy at the nickel M$_{2,3}$ edge. It is observed that the core-level absorption lineshape of photoexcited nickel can be described by a Gaussian broadening ($σしぐま$) and a red shift ($ωおめが_{s}$) of the ground state absorption spectrum. Theory predicts, and t… ▽ More Direct measurements of photoexcited carrier dynamics in nickel are made using few-femtosecond extreme ultraviolet (XUV) transient absorption spectroscopy at the nickel M$_{2,3}$ edge. It is observed that the core-level absorption lineshape of photoexcited nickel can be described by a Gaussian broadening ($σしぐま$) and a red shift ($ωおめが_{s}$) of the ground state absorption spectrum. Theory predicts, and the experimental results verify that after initial rapid carrier thermalization, the electron temperature increase ($ΔでるたT$) is linearly proportional to the Gaussian broadening factor $σしぐま$, providing quantitative real-time tracking of the relaxation of the electron temperature. Measurements reveal an electron cooling time for 50 nm thick polycrystalline nickel films of 640$\pm$80 fs. With hot thermalized carriers, the spectral red shift exhibits a power-law relationship with the change in electron temperature of $ωおめが_{s}\proptoΔでるたT^{1.5}$. Rapid electron thermalization via carrier-carrier scattering accompanies and follows the nominal 4 fs photoexcitation pulse until the carriers reach a quasi-thermal equilibrium. Entwined with a <6 fs instrument response function, carrier thermalization times ranging from 34 fs to 13 fs are estimated from experimental data acquired at different pump fluences and it is observed that the electron thermalization time decreases with increasing pump fluence. The study provides an initial example of measuring electron temperature and thermalization in metals in real time with XUV light, and it lays a foundation for further investigation of photoinduced phase transitions and carrier transport in metals with core-level absorption spectroscopy. △ Less

Submitted 5 January, 2021; v1 submitted 29 September, 2020; originally announced September 2020.

Comments: 20 pages, 8 figures

Journal ref: Phys. Rev. B 103, 064305 (2021)

Excited states in hydrogenated single-layer MoS$_2$

Authors: Naseem Ud Din, Volodymyr Turkowski, Talat S. Rahman

Abstract: We calculate the excitation spectrum of single-layer MoS$_2$ at several hydrogen coverages by using a method based on first-principles Density-Matrix Time-Dependent Density-Functional Theory (TDDFT). Our results show that the fully hydrogenated system is metallic, while in the low-coverage limit the spectrum of single-layer MoS$_2$ includes spin-polarized partially filled localized mid-gap states.… ▽ More We calculate the excitation spectrum of single-layer MoS$_2$ at several hydrogen coverages by using a method based on first-principles Density-Matrix Time-Dependent Density-Functional Theory (TDDFT). Our results show that the fully hydrogenated system is metallic, while in the low-coverage limit the spectrum of single-layer MoS$_2$ includes spin-polarized partially filled localized mid-gap states. These states arise from s-orbitals of H atoms which make a tilted bond with the surface S atoms. The calculated absorption spectrum of the system reveals standard excitonic peaks, which correspond to the bound valence-band hole and conduction-band electron, as well as excitonic peaks that involve the mid-gap charges. As in the case of pristine single-layer MoS$_2$, binding energies of the excitons of the hydrogenated system are found to be relatively large (few tens of meV), making their experimental detection facile and suggesting hydrogenation as a knob for tuning the optical properties of single-layer MoS$_2$. Importantly, we find hydrogenation to suppress visible light photoluminescence, in agreement with experimental observations. As an aside, we contrast the effects of hydrogen coverage to that of the next two elements in the same column of the periodic table (the lightest metals), Li and Na, on the spectral properties of single-layer MoS$_2$ which lead instead to the formation of n-doped non-magnetic semiconductors that do not allow excitonic states. △ Less

Submitted 27 August, 2020; originally announced August 2020.

Ultrafast demagnetization dynamics in Ni: role of electron correlations

Authors: Shree Ram Acharya, Volodymyr Turkowski, Guo-ping Zhang Talat S. Rahman

Abstract: Experimental observations of the ultrafast (less than 50 fs) demagnetization of Ni have so far defied theoretical explanations particularly since its spin-flipping time is much less than that resulting from spin-orbit and electron-lattice interactions. Through the application of an approach that benefits from spin-flip time-dependent density-functional theory and dynamical mean-field theory, we sh… ▽ More Experimental observations of the ultrafast (less than 50 fs) demagnetization of Ni have so far defied theoretical explanations particularly since its spin-flipping time is much less than that resulting from spin-orbit and electron-lattice interactions. Through the application of an approach that benefits from spin-flip time-dependent density-functional theory and dynamical mean-field theory, we show that proper inclusion of electron correlations and memory (time-dependence of electron-electron interaction) effects leads to demagnetization at the femtosecond scale, in good agreement with experimental observations. Furthermore, our calculations reveal that this ultrafast demagnetization results mainly from spin-flip transitions from occupied to unoccupied orbitals implying a dynamical reduction of exchange splitting. These conclusions are found to be valid for a wide range of laser pulse amplitudes. They also pave the way for ab initio investigations of ultrafast charge and spin dynamics in a variety of quantum materials in which electron correlations may play a definitive role. △ Less

Submitted 4 November, 2019; originally announced November 2019.

Journal ref: Phys. Rev. Lett. 125, 017202 (2020)

Plasmon excitations in mixed metallic nanoarrays

Authors: Kevin Conley, Neha Nayyar, Tuomas P. Rossi, Mikael Kuisma, Volodymyr Turkowski, Martti Puska, Talat S. Rahman

Abstract: We study the plasmonic properties of arrays of atomic chains which comprise noble (Cu, Ag, and Au) and transition (Pd, Pt) metal atoms using time-dependent density-functional theory. We show that the response to the electromagnetic radiation is related to both physics, the geometry-dependent confinement of sp-valence electrons, and chemistry, the energy position of d-electrons in the different ato… ▽ More We study the plasmonic properties of arrays of atomic chains which comprise noble (Cu, Ag, and Au) and transition (Pd, Pt) metal atoms using time-dependent density-functional theory. We show that the response to the electromagnetic radiation is related to both physics, the geometry-dependent confinement of sp-valence electrons, and chemistry, the energy position of d-electrons in the different atomic species and the hybridization between d and sp electrons. As a result it is possible to tune the position of the surface plasmon resonance, split it to several peaks, and eventually achieve broadband absorption of radiation. Mixing the arrays with transition metals can strongly attenuate the plasmonic behaviour. We analyze the origin of these phenomena and show that they arise from rich interactions between single-particle electron-hole and collective electron excitations. The tunability of the plasmonic response of arrays of atomic chains, which can be realized on solid surfaces, opens wide possibilities for their applications. In the present study we obtain guidelines how the desired properties can be achieved. △ Less

Submitted 30 September, 2018; originally announced October 2018.

Electron-phonon interaction, excitations and ultrafast photoemission from doped monolayer MoS2

Authors: Neha Nayyar, Duy Le, Volodymyr Turkowski, Talat S. Rahman

Abstract: We analyze the effect of electron-phonon coupling on photoemission properties and ultrafast response of doped monolayer MoS2. The analysis is based on combined DFT and many-body (Eliashberg theory) approaches. In particular, we have calculated the electronic and phonon spectra, the electron-phonon coupling and the electronic spectral function of the system at different values of doping. We have al… ▽ More We analyze the effect of electron-phonon coupling on photoemission properties and ultrafast response of doped monolayer MoS2. The analysis is based on combined DFT and many-body (Eliashberg theory) approaches. In particular, we have calculated the electronic and phonon spectra, the electron-phonon coupling and the electronic spectral function of the system at different values of doping. We have also analyzed the emissive properties and the response of the system to femtosecond (fs) laser pulses. It is shown that position of the emission peak of undoped system is in agreement with the experimental data if one takes into account the excitonic effects. The results for the self-energy and spectral functions of the doped systems suggest that one can expect ultrafast processes to be important in the system response , which makes the system attractive from the point of view of modern technological applications. Similar to graphene, the doped system demonstrates ultrafast (fs) relaxation of the electronic subsystem when excited by fs pulses, and a high ultrafast phonon relaxation-induced spectral fluence of visible light emission. Together with high carrier mobility, these features of monolayer MoS2 might be used in modern optoelectronic technologies. △ Less

Submitted 30 January, 2015; originally announced January 2015.

Nonadiabatic exchange-correlation kernel for strongly correlated materials

Authors: Volodymyr Turkowski, Talat S. Rahman

Abstract: We formulate a rigorous method for calculating a nonadiabatic (frequency-dependent) exchange-correlation (XC) kernel required for correct description of both equilibrium and nonequilibrium properties of strongly correlated systems within Time-Dependent Density Functional Theory (TDDFT). To do so we use the expression for charge susceptibility provided by Dynamical Mean Field Theory (DMFT) for the… ▽ More We formulate a rigorous method for calculating a nonadiabatic (frequency-dependent) exchange-correlation (XC) kernel required for correct description of both equilibrium and nonequilibrium properties of strongly correlated systems within Time-Dependent Density Functional Theory (TDDFT). To do so we use the expression for charge susceptibility provided by Dynamical Mean Field Theory (DMFT) for the effective multi-orbital Hubbard Model. We tested our formalism by applying it to the one-band Hubbard model: our nonadiabatic kernel leads to a significant modification of the excitation spectrum, shifting the peak that appears in adiabatic (simplified) solutions and disclosing a new one, in agreement with the DMFT solution. We also used our method to track the nonequilibrium charge-density response of a multi-orbital perovskite Mott insulator, YTiO3, to a perturbation by a femtosecond (fs) laser pulse. The results were quite different from those provided by the corresponding adiabatic formalism. These initial investigations indicate that electron-electron correlations and nonadiabatic features can significantly affect the spectrum and nonequilibrium properties of strongly correlated systems. △ Less

Submitted 10 October, 2014; v1 submitted 26 September, 2014; originally announced September 2014.

Magnetic anisotropy of FePt nanoparticles

Authors: Alamgir Kabir, Jun Hu, Volodymyr Turkowski, Ruqian Wu, Talat S. Rahman

Abstract: We carry out a systematic theoretical investigation of Magneto Crystalline Anisotropy (MCA) of L10 FePt clusters with alternating Fe and Pt planes along the (001) direction. We calculate the structural relaxation and magnetic moment of each cluster by using ab initio spin-polarized density functional theory (DFT), and the MCA with both spin-polarized DFT (including spin-orbit coupling self-consist… ▽ More We carry out a systematic theoretical investigation of Magneto Crystalline Anisotropy (MCA) of L10 FePt clusters with alternating Fe and Pt planes along the (001) direction. We calculate the structural relaxation and magnetic moment of each cluster by using ab initio spin-polarized density functional theory (DFT), and the MCA with both spin-polarized DFT (including spin-orbit coupling self-consistently) and the torque method. We find that the MCA of any composite structure of a given size is enhanced with respect to that of the same-sized pure Pt or pure Fe cluster as well as to that of any pair of Fe and Pt atoms in bulk L10 FePt. This enhancement results from the hybridization we observe between the 3d orbital of the Fe atoms and the 5d orbital of their Pt neighbors. This hybridization, however, affects the electronic properties of the component atoms in significantly different ways. While it somewhat increases the spin moment of the Fe atoms, it has little effect on their orbital moment; at the same time, it greatly increases both the spin and orbital moment of the Pt atoms. Given the fact that the spin-orbit coupling (SOC) constant of Pt is about 7 times greater than that of Fe, this Fe-induced jump in the orbital moment of the Pt atoms produces the increase in MCA of the composite structures over that of their pure counterparts. That any composite structure exhibits higher MCA than bulk L10 FePt results from the lower coordination of Pt atoms in the cluster, whether Fe or Pt predominates within it. We also find that bipyramidal clusters whose central layer is Pt have higher MCA than their same-sized counterparts whose central layer is Fe. This results from the fact that Pt atoms in such configurations are coordinated with more Fe atoms than in the latter. By thus participating in more instances of hybridization, they contribute higher orbital moments to the overall MCA of the unit. △ Less

Submitted 11 March, 2015; v1 submitted 15 September, 2014; originally announced September 2014.

Journal ref: Phys. Rev. B 92, 054424 (2015)

Tuning plasmon excitations in pure and transition metal-doped arrays of noble metal nanochains

Authors: Neha Nayyar, Volodymyr Turkowski, Talat S. Rahman

Abstract: We study the plasmonic properties of coupled noble-metal nanochains in the case of different number of coupled chains and doping by different transition-metal (TM) atoms within the time-dependent density-functional theory (TDDFT) approach. We find that as the number of chains in the array increases the plasmon peak shifts from the sub-eV towards the visible range. As doping with TM atoms increases… ▽ More We study the plasmonic properties of coupled noble-metal nanochains in the case of different number of coupled chains and doping by different transition-metal (TM) atoms within the time-dependent density-functional theory (TDDFT) approach. We find that as the number of chains in the array increases the plasmon peak shifts from the sub-eV towards the visible range. As doping with TM atoms increases, the visible absorption band broadens, owing to formation of additional plasmon peaks. The optical response is very sensitive to the type of doped atoms, their number and position; in particular, the additional peaks are most pronounced in the case of weak doping when they correspond to local plasmon oscillations around the impurity atom. These effects have a potential to be used in various modern technologies, from sensors to solar cells. Most of the studies of nano-plasmon effects have been focused on alkali- and noble-metal systems with extended s-electron states, while it was believed that doping with TM atoms with their more localized charge as a rule leads to an attenuation of the plasmon modes. We demonstrate that TM atoms can play a constructive role in plasmon generation in small chain systems, and that plasmonic modes can emerge even in some pure TM nanochains. △ Less

Submitted 11 September, 2014; originally announced September 2014.

Time-dependent density-matrix functional theory for trion excitations: application to monolayer MoS2

Authors: Alfredo Ramirez-Torres, Volodymyr Turkowski, Talat S. Rahman

Abstract: We study possible optically excited bound states in monolayer MoS2: excitons and trions. For this purpose we formulate and apply a generalized time-dependent density-matrix functional approach for bound states of multiple excitations. The approach was used in the cases of three different types of the exchange-correlation (XC) kernel: 1) two local kernels: a phenomenological contact and the adiabat… ▽ More We study possible optically excited bound states in monolayer MoS2: excitons and trions. For this purpose we formulate and apply a generalized time-dependent density-matrix functional approach for bound states of multiple excitations. The approach was used in the cases of three different types of the exchange-correlation (XC) kernel: 1) two local kernels: a phenomenological contact and the adiabatic local-density approximation (ALDA) (X and XC); 2) gradient-corrected X kernels: GEA, PW91 and PBE; and 3) two long-range (LR) kernels: a phenomenological (Coulomb) and Slater kernels. In the case of exciton, we find that LDA and its gradient-corrected kernels lead to too weak binding energy comparing to the experimental data, while the LR kernels are capable to reproduce the experimental results. Similarly, in the LR case (as well as in the case of local kernel), one can obtain the experimental value of the trion binding energy by taking into account the screening effects. Our results suggest that similar to the excitons, the LR structure of the XC kernel is necessary to describe the trion bound states. Our calculations for the first time confirm theoretically with time-dependent density-functional theory approach that in agreement with experimental data the exciton and trion binding energies are of order of hundreds (excitons) and tenth (trions) meVs, which can be used in different technological applications at the room temperature regime. The approach can be straightforwardly extended on the case of bound states and nonequilibrium response of systems with larger number of bound electrons and holes, including biexcitons. △ Less

Submitted 17 February, 2014; originally announced February 2014.

A DFT+Nonhomogeneous DMFT approach for Finite Systems

Authors: Alamgir Kabir, Volodymyr Turkowski, Talat S. Rahman

Abstract: For reliable and efficient inclusion of electron-electron correlation effects in nanosystems we propose a combined density-functional-theory/nonhomogeneous dynamical-mean-field-theory (DFT + DMFT) approach which employs an approximate Iterative Perturbative Theory (IPT) impurity solver. The validity of the method is demonstrated by successful examination of the size-dependent magnetic properties o… ▽ More For reliable and efficient inclusion of electron-electron correlation effects in nanosystems we propose a combined density-functional-theory/nonhomogeneous dynamical-mean-field-theory (DFT + DMFT) approach which employs an approximate Iterative Perturbative Theory (IPT) impurity solver. The validity of the method is demonstrated by successful examination of the size-dependent magnetic properties of iron nanoparticles containing 11-100 atoms. We show that the DFT+ DMFT solution is in very good agreement with experimental data, in particular it does not lead to the overestimation of magnetization that is found with the DFT and DFT+U techniques. More importantly, we demonstrate that DFT+DMFT approach can be used for accurate and realistic description of nanosystems containing about hundred atoms. △ Less

Submitted 13 November, 2013; originally announced November 2013.

Time-dependent density-functional theory of exciton-exciton correlations in the nonlinear optical response

Authors: Volodymyr Turkowski, Michael N. Leuenberger

Abstract: We analyze possible nonlinear exciton-exciton correlation effects in the optical response of semiconductors by using a time-dependent density-functional theory (TDDFT) approach. For this purpose, we derive the nonlinear (third-order) TDDFT equation for the excitonic polarization. In this equation, the nonlinear time-dependent effects are described by the time-dependent (non-adiabatic) part of the… ▽ More We analyze possible nonlinear exciton-exciton correlation effects in the optical response of semiconductors by using a time-dependent density-functional theory (TDDFT) approach. For this purpose, we derive the nonlinear (third-order) TDDFT equation for the excitonic polarization. In this equation, the nonlinear time-dependent effects are described by the time-dependent (non-adiabatic) part of the effective exciton-exciton interaction, which depends on the exchange-correlation (XC) kernel. We apply the approach to study the nonlinear optical response of a GaAs quantum well. In particular, we calculate the 2D Fourier spectra of the system and compare it with experimental data. We find that it is necessary to use a non-adiabatic XC kernel to describe excitonic bound states - biexcitons, which are formed due to the retarded TDDFT exciton-exciton interaction. △ Less

Submitted 19 December, 2013; v1 submitted 15 June, 2013; originally announced June 2013.

Nonadiabatic Time-Dependent Spin-Density Functional Theory for strongly correlated systems

Authors: Volodymyr Turkowski, Talat S. Rahman

Abstract: We propose a nonadiabatic time-dependent spin-density functional theory (TDSDFT) approach for studying the single-electron excited states and the ultrafast response of systems with strong electron correlations. The correlations are described by the correlation part of the nonadiabatic exchange-correlation (XC) kernel, which is constructed by using some exact results for the Hubbard model of strong… ▽ More We propose a nonadiabatic time-dependent spin-density functional theory (TDSDFT) approach for studying the single-electron excited states and the ultrafast response of systems with strong electron correlations. The correlations are described by the correlation part of the nonadiabatic exchange-correlation (XC) kernel, which is constructed by using some exact results for the Hubbard model of strongly correlated electrons. We demonstrate that the corresponding nonadiabatic XC kernel reproduces main features of the spectrum of the Hubbard dimer and infinite-dimensional Hubbard model, some of which are impossible to obtain within the adiabatic approach. The theory may be applied for DFT study of strongly correlated electron systems in- and out-of-equilibrium, including the important case of nanostructures, for which it leads to a dramatic reduction of necessary computational power. △ Less

Submitted 24 October, 2013; v1 submitted 23 October, 2012; originally announced October 2012.

Journal ref: J. Phys.: Condens. Matter 26, 022201 (2014) (Fast Track Communication)

Spectral moment sum rules for the retarded Green's function and self-energy of the inhomogeneous Bose-Hubbard model in equilibrium and nonequilibrium

Authors: J. K. Freericks, V. Turkowski, H. R. Krishnamurthy, M. Knap

Abstract: We derive expressions for the zeroth and the first three spectral moment sum rules for the retarded Green's function and for the zeroth and the first spectral moment sum rules for the retarded self-energy of the inhomogeneous Bose-Hubbard model in nonequilibrium, when the local on-site repulsion and the chemical potential are time-dependent, and in the presence of an external time-dependent electr… ▽ More We derive expressions for the zeroth and the first three spectral moment sum rules for the retarded Green's function and for the zeroth and the first spectral moment sum rules for the retarded self-energy of the inhomogeneous Bose-Hubbard model in nonequilibrium, when the local on-site repulsion and the chemical potential are time-dependent, and in the presence of an external time-dependent electromagnetic field. We also evaluate these moments for equilibrium, where all time dependence and external fields vanish. Unlike similar sum rules for the Fermi-Hubbard model, in the Bose-Hubbard model case, the sum rules often depend on expectation values that cannot be determined simply from parameters in the Hamiltonian like the interaction strength and chemical potential, but require knowledge of equal time many-body expectation values from some other source. We show how one can approximately evaluate these expectation values for the Mott-insulating phase in a systematic strong-coupling expansion in powers of the hopping divided by the interaction. We compare the exact moments to moments of spectral functions calculated from a variety of different approximations and use them to benchmark their accuracy. △ Less

Submitted 19 August, 2012; originally announced August 2012.

Comments: (24 pages, 5 figures, submitted to Phys. Rev. A)

Journal ref: Phys. Rev. A 87, 013628 (2013)

Transition metal impurity-induced generation of plasmonic collective modes in small gold clusters

Authors: N. Nayyar, A. Kabir, V. Turkowski, Talat S. Rahman

Abstract: We study the optical properties of small gold chains doped with different transition metal (TM) atoms (Ni,Rh,Fe) by using the time-dependent density-functional theory (TDDFT) approach. The optical absorption spectrum of such systems demonstrates a collective plasmon mode when the number of atoms is larger than approximately 10, and this mode splits into two peaks when the Au chain is doped with so… ▽ More We study the optical properties of small gold chains doped with different transition metal (TM) atoms (Ni,Rh,Fe) by using the time-dependent density-functional theory (TDDFT) approach. The optical absorption spectrum of such systems demonstrates a collective plasmon mode when the number of atoms is larger than approximately 10, and this mode splits into two peaks when the Au chain is doped with some of the TM (Ni,Fe) atoms. We associate the additional peak with a local plasmonic mode which corresponds to the charge oscillations around the potential created by the d-orbitals of the impurity atoms. The effect takes place when the number of TM atoms is much smaller than the number of Au atoms. This behavior is opposite to the case of larger bulk noble-metal-TM clusters (radius $>$1nm), where the doping with TM atoms does not lead to any generation of additional modes, and often leads to a suppression of the main plasmon peak. This phenomenon of tuning the optical properties of nanosystems with transition metal atoms can be used in many practical applications in the field of nanotechnology. △ Less

Submitted 5 September, 2011; originally announced September 2011.

Possible high-temperature superconductivity in multilayer graphane: can the cuprates be beaten?

Authors: V. M. Loktev, V. Turkowski

Abstract: We analyze a possible superconductivity in the hole-doped system of layered hydrogenized graphene by taking into account thermal fluctuations of the order parameter. In particular, we demonstrate that in the one-layer case the values of the high mean-field (MF) critical temperature $T_{c}^{MF}\sim 80-90K$, predicted recently by Savini et al \cite{Savini}, do not alter significantly due to the fluc… ▽ More We analyze a possible superconductivity in the hole-doped system of layered hydrogenized graphene by taking into account thermal fluctuations of the order parameter. In particular, we demonstrate that in the one-layer case the values of the high mean-field (MF) critical temperature $T_{c}^{MF}\sim 80-90K$, predicted recently by Savini et al \cite{Savini}, do not alter significantly due to the fluctuations, and the Berezinskii-Kosterlitz-Thouless (BKT) critical temperature of the vortex superconductivity is almost the same as the MF temperature at doping 0.01-0.1. We show that in the case of multilayer system, when the coupling between the layers stabilizes the superconducting phase in the form of fluxon superconductivity, the critical temperature $T_{c}$ can increase dramatically to the values $\sim 150K$, higher than the corresponding values in cuprates under ambient pressure. △ Less

Submitted 22 February, 2011; originally announced February 2011.

Time-dependent density-matrix functional theory for biexcitonic phenomena

Authors: Volodymyr Turkowski, Carsten A. Ullrich, Talat S. Rahman, Michael N. Leuenberger

Abstract: We formulate a time-dependent density-matrix functional theory (TDDMFT) approach for higher-order correlation effects like biexcitons in optical processes in solids based on the reduced two-particle density-matrix formalism within the normal orbital representation. A TDDMFT version of the Schrödinger equation for biexcitons in terms of one- and two-body reduced density matrices is derived, which l… ▽ More We formulate a time-dependent density-matrix functional theory (TDDMFT) approach for higher-order correlation effects like biexcitons in optical processes in solids based on the reduced two-particle density-matrix formalism within the normal orbital representation. A TDDMFT version of the Schrödinger equation for biexcitons in terms of one- and two-body reduced density matrices is derived, which leads to finite biexcitonic binding energies already with an adiabatic approximation. Biexcitonic binding energies for several bulk semiconductors are calculated using a contact biexciton model. △ Less

Submitted 9 August, 2010; originally announced August 2010.

Inhomogeneous spectral moment sum rules for the retarded Green function and self-energy of strongly correlated electrons or ultracold fermionic atoms in optical lattices

Authors: J. K. Freericks, V. Turkowski

Abstract: Spectral moment sum rules are presented for the inhomogeneous many-body problem described by the fermionic Falicov-Kimball or Hubbard models. These local sum rules allow for arbitrary hoppings, site energies, and interactions. They can be employed to quantify the accuracy of numerical solutions to the inhomogeneous many-body problem like strongly correlated multilayered devices, ultracold atoms… ▽ More Spectral moment sum rules are presented for the inhomogeneous many-body problem described by the fermionic Falicov-Kimball or Hubbard models. These local sum rules allow for arbitrary hoppings, site energies, and interactions. They can be employed to quantify the accuracy of numerical solutions to the inhomogeneous many-body problem like strongly correlated multilayered devices, ultracold atoms in an optical lattice with a trap potential, strongly correlated systems that are disordered, or systems with nontrivial spatial ordering like a charge density wave or a spin density wave. We also show how the spectral moment sum rules determine the asymptotic behavior of the Green function, self-energy, and dynamical mean field, when applied to the dynamical mean-field theory solution of the many body problem. In particular, we illustrate in detail how one can dramatically reduce the number of Matsubara frequencies needed to solve the Falicov-Kimball model, while still retaining high precision, and we sketch how one can incorporate these results into Hirsch-Fye quantum Monte Carlo solvers for the Hubbard (or more complicated) models. Since the solution of inhomogeneous problems is significantly more time consuming than periodic systems, efficient use of these sum rules can provide a dramatic speed up in the computational time required to solve the many-body problem. We also discuss how these sum rules behave in nonequilibrium situations as well, where the Hamiltonian has explicit time dependence due to a driving field or due to the time-dependent change of a parameter like the interaction strength or the origin of the trap potential. △ Less

Submitted 26 July, 2009; v1 submitted 7 July, 2009; originally announced July 2009.

Comments: (28 pages, 6 figures, ReVTeX) Paper updated to correct equations 11, 24, and 28

Journal ref: Phys. Rev. B 80, 115119 (2009)

Time-dependent density-functional approach for exciton binding energies

Authors: V. Turkowski, C. A. Ullrich

Abstract: Optical processes in insulators and semiconductors, including excitonic effects, can be described in principle exactly using time-dependent density-functional theory (TDDFT). Starting from a linearization of the TDDFT semiconductor Bloch equations in a two-band model, we derive a simple formalism for calculating excitonic binding energies. This formalism leads to a generalization of the standard… ▽ More Optical processes in insulators and semiconductors, including excitonic effects, can be described in principle exactly using time-dependent density-functional theory (TDDFT). Starting from a linearization of the TDDFT semiconductor Bloch equations in a two-band model, we derive a simple formalism for calculating excitonic binding energies. This formalism leads to a generalization of the standard Wannier equation for excitons, featuring a nonlocal effective electron-hole interaction determined by long-range and dynamical exchange-correlation (XC) effects. We calculate excitonic binding energies in several direct-gap semiconductors, using exchange-only and model XC kernels. △ Less

Submitted 27 February, 2009; originally announced March 2009.

Comments: 4 pages, 1 figure

Time-dependent density-functional theory for electronic excitations in materials: basics and perspectives

Authors: C. A. Ullrich, V. Turkowski

Abstract: Time-dependent density-functional theory (TDDFT) is widely used to describe electronic excitations in complex finite systems with large numbers of atoms, such as biomolecules and nanocrystals. The first part of this paper will give a simple and pedagogical explanation, using a two-level system, which shows how the basic TDDFT formalism for excitation energies works. There is currently an intense… ▽ More Time-dependent density-functional theory (TDDFT) is widely used to describe electronic excitations in complex finite systems with large numbers of atoms, such as biomolecules and nanocrystals. The first part of this paper will give a simple and pedagogical explanation, using a two-level system, which shows how the basic TDDFT formalism for excitation energies works. There is currently an intense effort underway to develop TDDFT methodologies for the charge and spin dynamics in extended systems, to calculate optical properties of bulk and nanostructured materials, and to study transport through molecular junctions. The second part of this paper highlights some challenges and recent advances of TDDFT in these areas. Two examples are discussed: excitonic effects in insulators and intersubband plasmon excitations in doped semiconductor quantum wells. △ Less

Submitted 14 August, 2008; originally announced August 2008.

Comments: 15 pages, 2 figures, International Conference on Materials Discovery and Databases: Materials Informatics and DFT

Spin fluctuations and the "strange metal" behavior of a weakly doped antiferromagnet

Authors: V. M. Loktev, V. Turkowski

Abstract: We analyze the spectral properties of a phenomenological model for a weakly doped two-dimensional antiferromagnet, in which the carriers move within one of the two sublattices where they were introduced. Such a constraint results in the free carrier spectra with the maxima at ${\bf k}=(\pm πぱい/2 ,\pm πぱい/2)$ observed in some cuprates. We consider the spectral properties of the model by taking into a… ▽ More We analyze the spectral properties of a phenomenological model for a weakly doped two-dimensional antiferromagnet, in which the carriers move within one of the two sublattices where they were introduced. Such a constraint results in the free carrier spectra with the maxima at ${\bf k}=(\pm πぱい/2 ,\pm πぱい/2)$ observed in some cuprates. We consider the spectral properties of the model by taking into account fluctuations of the spins in the antiferromagnetic background. We show that such fluctuations lead to a non-pole-like structure of the single-hole Green's function and these fluctuations can be responsible for some anomalous "strange metal" properties of underdoped cuprates in the nonsuperconducting regime. △ Less

Submitted 13 May, 2008; originally announced May 2008.

Comments: Submitted to Condensed Matter Physics

Nonequilibrium sum rules for the retarded self-energy of strongly correlated electrons

Authors: V. Turkowski, J. K. Freericks

Abstract: We derive the first two moment sum rules of the conduction electron retarded self-energy for both the Falicov-Kimball model and the Hubbard model coupled to an external spatially uniform and time-dependent electric field (this derivation also extends the known nonequilibrium moment sum rules for the Green's functions to the third moment). These sum rules are used to further test the accuracy of… ▽ More We derive the first two moment sum rules of the conduction electron retarded self-energy for both the Falicov-Kimball model and the Hubbard model coupled to an external spatially uniform and time-dependent electric field (this derivation also extends the known nonequilibrium moment sum rules for the Green's functions to the third moment). These sum rules are used to further test the accuracy of nonequilibrium solutions to the many-body problem; for example, we illustrate how well the self-energy sum rules are satisfied for the Falicov-Kimball model in infinite dimensions and placed in a uniform electric field turned on at time t=0. In general, the self-energy sum rules are satisfied to a significantly higher accuracy than the Green's functions sum rules. △ Less

Submitted 10 January, 2008; originally announced January 2008.

Journal ref: Phys. Rev. B 77, 205102-1--16 (2008)

Time-dependent density-functional theory for ultrafast interband excitations

Authors: V. Turkowski, C. A. Ullrich

Abstract: We formulate a time-dependent density functional theory (TDDFT) in terms of the density matrix to study ultrafast phenomena in semiconductor structures. A system of equations for the density matrix components, which is equivalent to the time-dependent Kohn-Sham equation, is derived. From this we obtain a TDDFT version of the semiconductor Bloch equations, where the electronic many-body effects a… ▽ More We formulate a time-dependent density functional theory (TDDFT) in terms of the density matrix to study ultrafast phenomena in semiconductor structures. A system of equations for the density matrix components, which is equivalent to the time-dependent Kohn-Sham equation, is derived. From this we obtain a TDDFT version of the semiconductor Bloch equations, where the electronic many-body effects are taken into account in principle exactly. As an example, we study the optical response of a three-dimensional two-band insulator to an external short-time pulsed laser field. We show that the optical absorption spectrum acquires excitonic features when the exchange-correlation potential contains a $1/q^{2}$ Coulomb singularity. A qualitative comparison of the TDDFT optical absorption spectra with the corresponding results obtained within the Hartree-Fock approximation is made. △ Less

Submitted 19 September, 2007; originally announced September 2007.

Nonequilibrium dynamical mean-field theory of strongly correlated electrons

Authors: V. Turkowski, J. K. Freericks

Abstract: We present a review of our recent work in extending the successful dynamical mean-field theory from the equilibrium case to nonequilibrium cases. In particular, we focus on the problem of turning on a spatially uniform, but possibly time varying, electric field (neglecting all magnetic field effects). We show how to work with a manifestly gauge-invariant formalism, and compare numerical calculat… ▽ More We present a review of our recent work in extending the successful dynamical mean-field theory from the equilibrium case to nonequilibrium cases. In particular, we focus on the problem of turning on a spatially uniform, but possibly time varying, electric field (neglecting all magnetic field effects). We show how to work with a manifestly gauge-invariant formalism, and compare numerical calculations from a transient-response formalism to different types of approximate treatments, including the semiclassical Boltzmann equation and perturbation theory in the interaction. In this review, we solve the nonequilibrium problem for the Falicov-Kimball model, which is the simplest many-body model and the easiest problem to illustrate the nonequilibrium behavior in both diffusive metals and Mott insulators. Due to space restrictions, we assume the reader already has some familiarity both with the Kadanoff-Baym-Keldysh nonequilibrium formalism and with equilibrium dynamical mean-field theory; we provide a guide to the literature where additional details can be found. △ Less

Submitted 18 December, 2006; originally announced December 2006.

Comments: To be published in "Strongly Correlated Systems, Coherence and Entanglement"

Journal ref: "Strongly Correlated Systems: Coherence and Entanglement", edited by J. M. P. Carmelo, J. M. B. Lopes dos Santos, V. Rocha Vieira, and P. D. Sacramento (World Scientific, Singapore, 2007), pp. 187--210.

Superconductivity and superconducting order parameter phase fluctuations in a weakly doped antiferromagnet

Authors: V. M. Loktev, V. M. Turkowski

Abstract: The superconducting properties of a recently proposed phenomenological model for a weakly doped antiferromagnet are analyzed, taking into account fluctuations of the phase of the order parameter. In this model, we assume that the doped charge carriers can't move out of the antiferromagnetic sublattice they were introduced. This case corresponds to the free carrier spectra with the maximum at… ▽ More The superconducting properties of a recently proposed phenomenological model for a weakly doped antiferromagnet are analyzed, taking into account fluctuations of the phase of the order parameter. In this model, we assume that the doped charge carriers can't move out of the antiferromagnetic sublattice they were introduced. This case corresponds to the free carrier spectra with the maximum at ${\bf k}=(\pm πぱい/2 ,\pm πぱい/2)$, as it was observed in ARPES experiments in some of the cuprates in the insulating state [1]. The doping dependence of the superconducting gap and the temperature-carrier density phase diagram of the model are studied in the case of the $d_{x^{2}-y^{2}}$ pairing symmetry and different values of the effective coupling. A possible relevance of the results to the experiments on high-temperature superconductors is discussed. △ Less

Submitted 30 October, 2006; originally announced October 2006.

Comments: 16 pages, 4 figures

Nonequilibrium perturbation theory of the spinless Falicov-Kimball model

Authors: V. Turkowski, J. K. Freericks

Abstract: We perform a perturbative analysis for the nonequilibrium Green functions of the spinless Falicov-Kimball model in the presence of an arbitrary external time-dependent but spatially uniform electric field. The conduction electron self-energy is found from a strictly truncated second-order perturbative expansion in the local electron-electron repulsion U. We examine the current at half-filling, a… ▽ More We perform a perturbative analysis for the nonequilibrium Green functions of the spinless Falicov-Kimball model in the presence of an arbitrary external time-dependent but spatially uniform electric field. The conduction electron self-energy is found from a strictly truncated second-order perturbative expansion in the local electron-electron repulsion U. We examine the current at half-filling, and compare to both the semiclassical Boltzmann equation and exact numerical solutions for the contour-ordered Green functions from a transient-response formalism (in infinite dimensions) on the Kadanoff-Baym-Keldysh contour. We find a strictly truncated perturbation theory in the two-time formalism cannot reach the long-time limit of the steady state; instead it illustrates pathological behavior for times larger than approximately 2/U. △ Less

Submitted 23 October, 2006; originally announced October 2006.

Journal ref: Phys. Rev. B.75, 125110-1--14 (2007)

Nonequilibrium dynamical mean-field theory

Authors: J. K. Freericks, V. M. Turkowski, V. Zlatic

Abstract: The many-body formalism for dynamical mean-field theory is extended to treat nonequilibrium problems. We illustrate how the formalism works by examining the transient decay of the oscillating current that is driven by a large electric field turned on at time t=0. We show how the Bloch oscillations are quenched by the electron-electron interactions, and how their character changes dramatically fo… ▽ More The many-body formalism for dynamical mean-field theory is extended to treat nonequilibrium problems. We illustrate how the formalism works by examining the transient decay of the oscillating current that is driven by a large electric field turned on at time t=0. We show how the Bloch oscillations are quenched by the electron-electron interactions, and how their character changes dramatically for a Mott insulator. △ Less

Submitted 3 July, 2006; originally announced July 2006.

Comments: 4 pages, 4 figures

Journal ref: Phys. Rev. Lett. 97, 266408-1--4 (2006)

Bosonic sector of the two-dimensional Hubbard model studied within a two-pole approximation

Authors: Adolfo Avella, Ferdinando Mancini, Volodymyr Turkowski

Abstract: The charge and spin dynamics of the two-dimensional Hubbard model in the paramagnetic phase is first studied by means of the two-pole approximation within the framework of the Composite Operator Method. The fully self-consistent scheme requires: no decoupling, the fulfillment of both Pauli principle and hydrodynamics constraints, the simultaneous solution of fermionic and bosonic sectors and a v… ▽ More The charge and spin dynamics of the two-dimensional Hubbard model in the paramagnetic phase is first studied by means of the two-pole approximation within the framework of the Composite Operator Method. The fully self-consistent scheme requires: no decoupling, the fulfillment of both Pauli principle and hydrodynamics constraints, the simultaneous solution of fermionic and bosonic sectors and a very rich momentum dependence of the response functions. The temperature and momentum dependencies, as well as the dependency on the Coulomb repulsion strength and the filling, of the calculated charge and spin susceptibilities and correlation functions are in very good agreement with the numerical calculations present in the literature. △ Less

Submitted 8 May, 2006; originally announced May 2006.

Journal ref: Phys. Rev. B 67, 115123 (2003)

Spectral moment sum rules for strongly correlated electrons in time-dependent electric fields

Authors: V. M. Turkowski, J. K. Freericks

Abstract: We derive exact operator average expressions for the first two spectral moments of nonequilibrium Green's functions for the Falicov-Kimball model and the Hubbard model in the presence of a spatially uniform, time-dependent electric field. The moments are similar to the well-known moments in equilibrium, but we extend those results to systems in arbitrary time-dependent electric fields. Moment su… ▽ More We derive exact operator average expressions for the first two spectral moments of nonequilibrium Green's functions for the Falicov-Kimball model and the Hubbard model in the presence of a spatially uniform, time-dependent electric field. The moments are similar to the well-known moments in equilibrium, but we extend those results to systems in arbitrary time-dependent electric fields. Moment sum rules can be employed to estimate the accuracy of numerical calculations; we compare our theoretical results to numerical calculations for the nonequilibrium dynamical mean-field theory solution of the Falicov-Kimball model at half-filling. △ Less

Submitted 2 November, 2005; originally announced November 2005.

Comments: (16 pages, submitted to Phys. Rev. B)

Journal ref: Phys. Rev. B 73 075108--1-15 (2006); Erratum, Phys. Rev. B 73, 209902(E)-1--1 (2006)

Steady-state nonequilibrium dynamical mean-field theory and the quantum Boltzman

Authors: J. K. Freericks, V. M. Turkowski

Abstract: We derive the formalism for steady state nonequilibrium dynamical mean-field theory in a real-time formalism along the Kadanoff-Baym contour. The resulting equations of motion are first transformed to Wigner coordinates (average and relative time), and then re-expressed in terms of differential operators. Finally, we perform a Fourier transform with respect to the relative time, and take the fir… ▽ More We derive the formalism for steady state nonequilibrium dynamical mean-field theory in a real-time formalism along the Kadanoff-Baym contour. The resulting equations of motion are first transformed to Wigner coordinates (average and relative time), and then re-expressed in terms of differential operators. Finally, we perform a Fourier transform with respect to the relative time, and take the first-order limit in the electric field to produce the quantum Boltzmann equation for dynamical mean-field theory. We next discuss the structure of the equations and their solutions, describing how these equations reduce to the Drude result in the limit of a constant relaxation time. We also explicitly demonstrate the equivalence between the Kubo and nonequilibrium approaches to linear response. There are a number of interesting modifications of the conventional quantum Boltzmann equation that arise due to the underlying bandstructure of the lattice. △ Less

Submitted 2 November, 2005; originally announced November 2005.

Comments: (14 pages, proceedings of the Workshop on Progress in Nonequilibrium Green's Functions III, Kiel Germany)

Journal ref: J. Phys.: Confer. Ser. 35, 39--52 (2006)

Real-time formalism for studying the nonlinear response of "smart" materials to an electric field

Authors: J. K. Freericks, V. Turkowski, and V. Zlatic

Abstract: The nonlinear response of a material to a large electric field (steady or pulsed) often determines the ultimate performance of the material for electronics applications. The formalism for understanding nonlinear effects in conventional semiconductors is well understood. The formalism is less well developed for so-called ``smart'' materials that are tuned to lie close to the metal-insulator trans… ▽ More The nonlinear response of a material to a large electric field (steady or pulsed) often determines the ultimate performance of the material for electronics applications. The formalism for understanding nonlinear effects in conventional semiconductors is well understood. The formalism is less well developed for so-called ``smart'' materials that are tuned to lie close to the metal-insulator transition. Here we show that the nonlinear response of a strongly correlated electronic material can be calculated with a massively parallel algorithm by discretizing a continuous matrix operator on the Kadanoff-Baym contour in real time. We benchmark the technique by examining the solutions when the field vanishes and comparing the results to exact results from an equilibrium formalism. We briefly discuss the numerical issues associated with the case of a large electric field and present results that show how the Bloch oscillations become damped as the scattering due to electron correlations increases. △ Less

Submitted 2 November, 2005; originally announced November 2005.

Comments: (8 pages, Proceedings of the 2005 DOD HPMC User's Group Conference)

Journal ref: in Proceedings of the HPCMP Users Group Conference 2005, Nashville, TN, June 28--30, 2005 edited by D. E. Post (IEEE Computer Society, Los Alamitos, CA, 2005), pp. 25--34

Nonlinear response of Bloch electrons in infinite dimensions

Authors: V. Turkowski, J. K. Freericks

Abstract: The exact nonlinear response of noninteracting (Bloch) electrons is examined within a nonequilibrium formalism on the infinite-dimensional hypercubic lattice. We examine the effects of a spatially uniform, but time-varying electric field (ignoring magnetic-field effects). The electronic Green's functions, Wigner density of states, and time-varying current are all determined and analyzed. We stud… ▽ More The exact nonlinear response of noninteracting (Bloch) electrons is examined within a nonequilibrium formalism on the infinite-dimensional hypercubic lattice. We examine the effects of a spatially uniform, but time-varying electric field (ignoring magnetic-field effects). The electronic Green's functions, Wigner density of states, and time-varying current are all determined and analyzed. We study both constant and pulsed electric fields, focusing on the transient response region. These noninteracting Green's functions are an important input into nonequilibrium dynamical mean field theory for the nonlinear response of strongly correlated electrons. △ Less

Submitted 20 August, 2004; originally announced August 2004.

Comments: 11 pages 6 figures

Journal ref: Phys. Rev. B 71, 085104--1-11 (2005)

Parallelizing the Keldysh formalism for strongly correlated electrons

Authors: J. K. Freericks, V. Turkowski, V. Zlatic

Abstract: Nonequilibrium quantum mechanics can be solved with the Keldysh formalism, which evolves the quantum mechanical states forward in time in the presence of a time-dependent field, and then evolves them backward in time, undoing the effect of the time-dependent field. The Feynman path integral over the Keldysh contour is employed to calculate the strongly correlated Green's function. We examine the… ▽ More Nonequilibrium quantum mechanics can be solved with the Keldysh formalism, which evolves the quantum mechanical states forward in time in the presence of a time-dependent field, and then evolves them backward in time, undoing the effect of the time-dependent field. The Feynman path integral over the Keldysh contour is employed to calculate the strongly correlated Green's function. We examine the accuracy of this procedure for the simplest problem that requires a nonequilibrium formulation: the f-electron spectral function of the spinless Falicov-Kimball model. △ Less

Submitted 15 July, 2004; originally announced July 2004.

Comments: (8 pages, 9 figures, conference proceeding for the 2004 HPMC User's group conference)

Journal ref: in Proceedings of the HPCMP Users Group Conference 2004, Williamsburg, VA, June 7--11, 2004 edited by R. E. Peterkin (IEEE Computer Society, Los Alamitos, CA, 2004), p. 7--16

F-electron spectral function near a quantum critical point

Authors: J. K. Freericks, V. Turkowski, V. Zlatic

Abstract: We calculate the f-electron spectral function using a Keldysh formalism for the Falicov-Kimball model in infinite dimensions. We study the region close to the quantum critical point on both the hypercubic and Bethe lattices. We calculate the f-electron spectral function using a Keldysh formalism for the Falicov-Kimball model in infinite dimensions. We study the region close to the quantum critical point on both the hypercubic and Bethe lattices. △ Less

Submitted 15 July, 2004; originally announced July 2004.

Comments: (2 pages, 3 figures, submitted to SCES'04)

Journal ref: Physica 359--361C, 684--686 (2005)

F-electron spectral function of the Falicov-Kimball model in infinite dimensions: the half-filled case

Authors: J. K. Freericks, V. Turkowski, V. Zlatic

Abstract: The f-electron spectral function of the Falicov-Kimball model is calculated via a Keldysh-based many-body formalism originally developed by Brandt and Urbanek. We provide results for both the Bethe lattice and the hypercubic lattice at half filling. Since the numerical computations are quite sensitive to the discretization along the Kadanoff-Baym contour and to the maximum cutoff in time that is… ▽ More The f-electron spectral function of the Falicov-Kimball model is calculated via a Keldysh-based many-body formalism originally developed by Brandt and Urbanek. We provide results for both the Bethe lattice and the hypercubic lattice at half filling. Since the numerical computations are quite sensitive to the discretization along the Kadanoff-Baym contour and to the maximum cutoff in time that is employed, we analyze the accuracy of the results using a variety of different moment sum-rules and spectral formulas. We find that the f-electron spectral function has interesting temperature dependence becoming a narrow single-peaked function for small U and developing a gap, with two broader peaks for large U. △ Less

Submitted 15 July, 2004; originally announced July 2004.

Comments: (13 pages, 11 figures, typeset in RevTex 4)

Journal ref: Phys. Rev. B 71, 115111--1-12 (2005)

On the theory of pseudogap anisotropy in the cuprate superconductors

Authors: V. M. Loktev, V. Turkowski

Abstract: We show by means of the theory of the order parameter phase fluctuations that the temperature of "closing" (or "opening") of the gap (and pseudogap) in the electron spectra of superconductors with anisotropic order parameter takes place within a finite temperature range. Every Fourier-component of the order parameter has its own critical temperature. We show by means of the theory of the order parameter phase fluctuations that the temperature of "closing" (or "opening") of the gap (and pseudogap) in the electron spectra of superconductors with anisotropic order parameter takes place within a finite temperature range. Every Fourier-component of the order parameter has its own critical temperature. △ Less

Submitted 15 April, 2004; v1 submitted 14 April, 2004; originally announced April 2004.

Superconducting properties of the 2D models with different types of inter-particle coupling

Authors: V. M. Loktev, V. Turkowski

Abstract: The superconducting properties of the 2D fermion system with local and different types of the indirect boson-exchange attractions in the cases of s-wave and d-wave pairing are reviewed and analysed at T=0. In particular, the possibility of the crossover from the Bose-Einstein condensation regime to the BCS-like superconductivity with the carrier density and coupling changing in the case of diffe… ▽ More The superconducting properties of the 2D fermion system with local and different types of the indirect boson-exchange attractions in the cases of s-wave and d-wave pairing are reviewed and analysed at T=0. In particular, the possibility of the crossover from the Bose-Einstein condensation regime to the BCS-like superconductivity with the carrier density and coupling changing in the case of different pairing channels is discussed. Gaussian fluctuations of the order parameter are taken into account and the carrier density dependence of the gap is studied in this case. The role of the interaction form in the physical behavior of the systems is also discussed. △ Less

Submitted 8 May, 2003; originally announced May 2003.

Comments: To be published in 'Horizons in Superconductivity Research', Nova Science

Temperature-doping phase diagram of layered superconductors

Authors: V. M. Loktev, V. Turkowski

Abstract: The superconducting properties of a layered system are analyzed for the cases of zero- and non-zero angular momentum of the pairs. The effective thermodynamic potential for the quasi-2D XY-model for the gradients of the phase of the order parameter is derived from the microscopic superconducting Hamiltonian. The dependence of the superconducting critical temperature T_c on doping, or carrier den… ▽ More The superconducting properties of a layered system are analyzed for the cases of zero- and non-zero angular momentum of the pairs. The effective thermodynamic potential for the quasi-2D XY-model for the gradients of the phase of the order parameter is derived from the microscopic superconducting Hamiltonian. The dependence of the superconducting critical temperature T_c on doping, or carrier density, is studied at different values of coupling and inter-layer hopping. It is shown that the critical temperature T_c of the layered system can be lower than the critical temperature of the two-dimensional Berezinskii-Kosterlitz-Thouless transition T_BKT at some values of the model parameters, contrary to the case when the parameters of the XY-model do not depend on the microscopic Hamiltonian parameters. △ Less

Submitted 8 May, 2003; originally announced May 2003.

Comments: To be published in Phys. Rev. B

Journal ref: Phys. Rev. B 67, 214510 (2003)

Non-equilibrium properties of the S=1/2 Heisenberg model in a time-dependent magnetic field

Authors: V. Turkowski, V. R. Vieira, P. D. Sacramento

Abstract: The time-dependent behavior of the Heisenberg model in contact with a phonon heat bath and in an external time-dependent magnetic field is studied by means of a path integral approach. The action of the phonon heat bath is taken into account up to the second order in the coupling to the heath bath. It is shown that there is a minimal value of the magnetic field below which the average magnetizat… ▽ More The time-dependent behavior of the Heisenberg model in contact with a phonon heat bath and in an external time-dependent magnetic field is studied by means of a path integral approach. The action of the phonon heat bath is taken into account up to the second order in the coupling to the heath bath. It is shown that there is a minimal value of the magnetic field below which the average magnetization of the system does not relax to equilibrium when the external magnetic field is flipped. This result is in qualitative agreement with the mean field results obtained within $φふぁい^{4}$-theory. △ Less

Submitted 7 May, 2003; originally announced May 2003.

Comments: To be published in Physica A

Fluctuating order parameter in doped cuprate superconductors

Authors: V. M. Loktev, Yu. G. Pogorelov, V. M. Turkowski

Abstract: We discuss static fluctuations of the d-wave superconducting order parameter $Δでるた$ in CuO$_2$ planes, due to quasiparticle scattering by charged dopants. The analysis of two-particle anomalous Green functions at $T = 0$ permits to estimate the mean-square fluctuation $δでるた^2 = <Δでるた^2> - <Δでるた>^2$, averaged in random dopant configurations, to the lowest order in doping level $c$. Since $Δでるた$ is found to satu… ▽ More We discuss static fluctuations of the d-wave superconducting order parameter $Δでるた$ in CuO$_2$ planes, due to quasiparticle scattering by charged dopants. The analysis of two-particle anomalous Green functions at $T = 0$ permits to estimate the mean-square fluctuation $δでるた^2 = <Δでるた^2> - <Δでるた>^2$, averaged in random dopant configurations, to the lowest order in doping level $c$. Since $Δでるた$ is found to saturate with growing doping level while $δでるた$ remains to grow, this can explain the collapse of $T_c$ at overdoping. Also we consider the spatial correlations $<Δでるた(0)Δでるた({\bf R})>$ for order parameter in different points of the plane. △ Less

Submitted 14 February, 2003; originally announced February 2003.

Comments: RevTex4, 3 pages, to be published in Proceedings of New3SC-4 International Conference, San Diego, California, January 15-21, 2003

Crossover from superfluidity to superconductivity in a system with doping dependent attraction

Authors: V. M. Loktev, V. Turkowski

Abstract: Zero temperature crossover from superfluidity to superconductivity with carrier density increasing is studied for a two-dimensional system in the s-wave and d-wave pairing channels. It was assumed that the particle attraction correlation length depends on carrier density $n_f$ as $r_0\sim1/\sqrt{n_f}$. Such a dependence was found experimentally for the radius of magnetic correlations in… ▽ More Zero temperature crossover from superfluidity to superconductivity with carrier density increasing is studied for a two-dimensional system in the s-wave and d-wave pairing channels. It was assumed that the particle attraction correlation length depends on carrier density $n_f$ as $r_0\sim1/\sqrt{n_f}$. Such a dependence was found experimentally for the radius of magnetic correlations in $La_{2-x}Sr_xCuO_4$. The short range Coulomb repulsion was also taken into account. It is shown that the behaviour of the system with doping is fundamentally different from the case with $r_0(n_f)=const$. In particular, similarly to the d-wave case, the crossover in the s-channel takes place only if the coupling is larger of some minimal value, otherwise the Cooper pairing scenario takes place at any small carrier density. The relevance of the model to high-temperature superconductors is discussed. △ Less

Submitted 22 March, 2002; originally announced March 2002.

Comments: Accepted for publication in Physica C

Journal ref: Physica C 383, 256 (2002)

Green function of fermions in 2D superconducting Frohlich model with inhomogeneous order parameter

Authors: V. M. Loktev, V. M. Turkowski

Abstract: The fermion Green function and spectral characteristics for the 2D Frohlich model of superconductivity at static fluctuations in the phase of the order parameter are calculated. The results demonstrate strongly non-Fermi-liquid properties of the system at finite temperatures and relate with the pseudogap behavior of high-$T_{c}$ superconductors at relatively small charge carrier densities. The fermion Green function and spectral characteristics for the 2D Frohlich model of superconductivity at static fluctuations in the phase of the order parameter are calculated. The results demonstrate strongly non-Fermi-liquid properties of the system at finite temperatures and relate with the pseudogap behavior of high-$T_{c}$ superconductors at relatively small charge carrier densities. △ Less

Submitted 24 July, 2000; originally announced July 2000.

Journal ref: Low. Temp. Phys. 26, 97 (2000) (Fiz.Nizk.Temp. 26, 135 (2000))

Spin magnetic susceptibility in the two-layer Hubbard model

Authors: Ferdinando Mancini, Volodymyr Turkowski

Abstract: The two-layer Hubbard model is studied by using the composite operator method. Magnetic properties of the model in the normal state are analyzed by studying the uniform static spin magnetic susceptibility as a function of doping and temperature. The relevance of the model to describe the physical behavior of two-layered cuprate high-$T_{c}$ superconductors is discussed. The two-layer Hubbard model is studied by using the composite operator method. Magnetic properties of the model in the normal state are analyzed by studying the uniform static spin magnetic susceptibility as a function of doping and temperature. The relevance of the model to describe the physical behavior of two-layered cuprate high-$T_{c}$ superconductors is discussed. △ Less

Submitted 24 July, 2000; originally announced July 2000.

Journal ref: Physica B 284-288, 1575 (2000)

On temperature versus doping phase diagram of high critiical temperature superconductors

Authors: V. M. Loktev, V. M. Turkowski

Abstract: The attempt to describe the bell-shape dependence of the critical temperature of high-$T_{c}$ superconductors on charge carriers density is made. Its linear increase in the region of small densities (underdoped regime) is proposed to explain by the role of the order parameter phase 2D fluctuations which become less at this density growth. The critical temperature suppression in the region of lar… ▽ More The attempt to describe the bell-shape dependence of the critical temperature of high-$T_{c}$ superconductors on charge carriers density is made. Its linear increase in the region of small densities (underdoped regime) is proposed to explain by the role of the order parameter phase 2D fluctuations which become less at this density growth. The critical temperature suppression in the region of large carrier densities (overdoped regime) is connected with the appearance (because of doping) of the essential damping of long-wave bosons which in the frame of the model proposed define the mechanism of indirect inter-fermion attraction. △ Less

Submitted 2 March, 1998; originally announced March 1998.

Comments: 15 pages, 3 figures, EMTEX

On the Theory of the Pseudogap Formation in 2D Attracting Fermion Systems

Authors: V. M. Loktev, V. M. Turkowski

Abstract: Two-dimensional system of the fermions with the indirect Einstein phonon-exchange attraction and added local four-fermion interaction is considered. It is shown that in such a system at resulting attraction between particles a new nonsuperconducting phase arises along with the normal and superconducting phases. In this, called "abnormal normal", or pseudogap, phase the absolute value of the orde… ▽ More Two-dimensional system of the fermions with the indirect Einstein phonon-exchange attraction and added local four-fermion interaction is considered. It is shown that in such a system at resulting attraction between particles a new nonsuperconducting phase arises along with the normal and superconducting phases. In this, called "abnormal normal", or pseudogap, phase the absolute value of the order parameter is finite but its phase is a random quantity. It is important that the new phase really exists at low carrier density only, i.e. it shrinks with doping increasing in the case of phonon attraction. The relevance of the results for high-temperature superconductors is speculated. Key words: 2D metal, arbitrary carrier density, normal phase, abnormal normal phase, pseudogap, suderconducting phase, Berezinskii-Kosterlitz-Thouless phase, electron-electron and electron-hole pairing △ Less

Submitted 18 July, 1997; originally announced July 1997.

Comments: 19 pages, 2 figures (emtex)

Phase diagram in 2D Fröhlich model of metal at arbitrary carrier density: pseudogap versus doping

Authors: V. M. Loktev, S. G. Sharapov, V. M. Turkowski

Abstract: It is shown that in 2D system of the fermions with simplest indirect boson-induced attraction (through the Einstein phonon exchange as an example) along with the normal and superconducting phases there arises a new (called "abnormal normal" or pseudogap) one where the absolute value of the order parameter is finite but its phase is a random quantity. It is important that this new phase really ex… ▽ More It is shown that in 2D system of the fermions with simplest indirect boson-induced attraction (through the Einstein phonon exchange as an example) along with the normal and superconducting phases there arises a new (called "abnormal normal" or pseudogap) one where the absolute value of the order parameter is finite but its phase is a random quantity. It is important that this new phase really exists at low carriers density only, i.e. it shrinks when doping increases. The relevance of the results obtained with dependence of pseudogap on doping in high-temperature superconductors is speculated. △ Less

Submitted 13 March, 1997; v1 submitted 7 March, 1997; originally announced March 1997.

Comments: 13 pages, EMTeX, with 1 figure added

Journal ref: Physica C. 296 (1998) 84