Impurity model for non-equilibrium steady states

We propose an out-of-equilibrium impurity model for the dynamical mean-field description of the Hubbard model driven by a finite electric field. The out-of-equilibrium impurity environment is represented by a collection of equilibrium reservoirs at different chemical potentials. We discuss the validity of the impurity model and propose a non-perturbative method, based on a quantum Monte Carlo solver, which provides the steady-state solutions of the impurity and original lattice problems. We discuss the relevance of this approach to other non-equilibrium steady-state contexts.Comment: 5 pages, 4 figure

Emergence of novel magnetic order at finite temperature in overdoped pnictides

We examine the temperature dependence of the magnetic ordering in the frustrated Heisenberg $J_1-J_2$ model in presence of two different kind of dopants: vacancies or magnetic impurities. We demonstrate that, irrespective to their magnetic ratio, the introduction of impurities quenches the order by disorder selection mechanism associated with an Ising-like phase transition at low temperatures and gives way to a $90^\circ$ (anticollinear) order . The presence of dopants triggers a non trivial competition between entropically selected states (collinear) and energetically favoured ones (anticollinear) in dependence of both dilution and temperature. While in case of magnetic impurity, the interesting magnetic phases are observed for full range of temperature and doping, in case of nonmagnetic impurities every magnetic order is destroyed at all temperatures above $12\%$ dilution. At fixed low temperature and tuning the doping we show a first order phase transition leading to the re-entrance of the Ising-like order with percolation of islands of $90^\circ$ order. At fixed doping and varying the temperature we observe a transition from the anticollinear to the collinear phase assisted by a new emerging magnetic phase in the presence of magnetic impurities, whilst in case of vacancies this transition is characterised by a coexistent region of both. Furthermore, tuning the magnetic moment of the impurities, a complete collapse of the Ising-like order is attained. This is in agreement with observations of Ir dopant atoms in superconducting Ba(Fe$_{1-x}$Ir$_x$)$_2$As$_2$ with $x<0.047$

Renormalization of myoglobin-ligand binding energetics by quantum many-body effects

We carry out a first-principles atomistic study of the electronic mechanisms of ligand binding and discrimination in the myoglobin protein. Electronic correlation effects are taken into account using one of the most advanced methods currently available, namely a linear-scaling density functional theory (DFT) approach wherein the treatment of localized iron 3d electrons is further refined using dynamical mean-field theory (DMFT). This combination of methods explicitly accounts for dynamical and multi-reference quantum physics, such as valence and spin fluctuations, of the 3d electrons, whilst treating a significant proportion of the protein (more than 1000 atoms) with density functional theory. The computed electronic structure of the myoglobin complexes and the nature of the Fe-O2 bonding are validated against experimental spectroscopic observables. We elucidate and solve a long standing problem related to the quantum-mechanical description of the respiration process, namely that DFT calculations predict a strong imbalance between O2 and CO binding, favoring the latter to an unphysically large extent. We show that the explicit inclusion of many body-effects induced by the Hund's coupling mechanism results in the correct prediction of similar binding energies for oxy- and carbonmonoxymyoglobin.Comment: 7 pages, 5 figures. Accepted for publication in the Proceedings of the National Academy of Sciences of the United States of America (2014). For the published article see http://www.pnas.org/content/early/2014/04/09/1322966111.abstrac

Metal-insulator transition in copper oxides induced by apex displacements

High temperature superconductivity has been found in many kinds of compounds built from planes of Cu and O, separated by spacer layers. Understanding why critical temperatures are so high has been the subject of numerous investigations and extensive controversy. To realize high temperature superconductivity, parent compounds are either hole-doped, such as {La$_{2}$CuO$_4$} (LCO) with Sr (LSCO), or electron doped, such as {Nd$_{2}$CuO$_4$} (NCO) with Ce (NCCO). In the electron doped cuprates, the antiferromagnetic phase is much more robust than the superconducting phase. However, it was recently found that the reduction of residual out-of-plane apical oxygens dramatically affects the phase diagram, driving those compounds to a superconducting phase. Here we use a recently developed first principles method to explore how displacement of the apical oxygen (A-O) in LCO affects the optical gap, spin and charge susceptibilities, and superconducting order parameter. By combining quasiparticle self-consistent GW (QS\emph{GW}) and dynamical mean field theory (DMFT), that LCO is a Mott insulator; but small displacements of the apical oxygens drive the compound to a metallic state through a localization/delocalization transition, with a concomitant maximum $d$-wave order parameter at the transition. We address the question whether NCO can be seen as the limit of LCO with large apical displacements, and elucidate the deep physical reasons why the behaviour of NCO is so different than the hole doped materials. We shed new light on the recent correlation observed between T$_c$ and the charge transfer gap, while also providing a guide towards the design of optimized high-Tc superconductors. Further our results suggest that strong correlation, enough to induce Mott gap, may not be a prerequisite for high-Tc superconductivity

Calculating DMFT forces in ab-initio ultrasoft pseudopotential formalism

In this paper, we show how to calculate analytical atomic forces within self-consistent density functional theory + dynamical mean-field theory (DFT+DMFT) approach in the case when ultra-soft or norm-conserving pseudopotentials are used. We show how to treat the non-local projection terms arising within the pseudopotential formalism and circumvent the problem of non-orthogonality of the Kohn-Sham eigenvectors. Our approach is, in principle, independent of the DMFT solver employed, and here was tested with the Hubbard I solver. We benchmark our formalism by comparing against the forces calculated in Ce$_{2}$O$_{3}$ and PrO$_2$ by numerical differentiation of the total free energy, as well as by comparing the energy profiles against the numerically integrated analytical forces.Comment: 12 pages, 3 figure

Local selection rules that can determine specific pathways of DNA unknotting by type II DNA topoisomerases

We performed numerical simulations of DNA chains to understand how local geometry of juxtaposed segments in knotted DNA molecules can guide type II DNA topoisomerases to perform very efficient relaxation of DNA knots. We investigated how the various parameters defining the geometry of inter-segmental juxtapositions at sites of inter-segmental passage reactions mediated by type II DNA topoisomerases can affect the topological consequences of these reactions. We confirmed the hypothesis that by recognizing specific geometry of juxtaposed DNA segments in knotted DNA molecules, type II DNA topoisomerases can maintain the steady-state knotting level below the topological equilibrium. In addition, we revealed that a preference for a particular geometry of juxtaposed segments as sites of strand-passage reaction enables type II DNA topoisomerases to select the most efficient pathway of relaxation of complex DNA knots. The analysis of the best selection criteria for efficient relaxation of complex knots revealed that local structures in random configurations of a given knot type statistically behave as analogous local structures in ideal geometric configurations of the corresponding knot typ

Local selection rules that can determine specific pathways of DNA unknotting by type II DNA topoisomerases

We performed numerical simulations of DNA chains to understand how local geometry of juxtaposed segments in knotted DNA molecules can guide type II DNA topoisomerases to perform very efficient relaxation of DNA knots. We investigated how the various parameters defining the geometry of inter-segmental juxtapositions at sites of inter-segmental passage reactions mediated by type II DNA topoisomerases can affect the topological consequences of these reactions. We confirmed the hypothesis that by recognizing specific geometry of juxtaposed DNA segments in knotted DNA molecules, type II DNA topoisomerases can maintain the steady-state knotting level below the topological equilibrium. In addition, we revealed that a preference for a particular geometry of juxtaposed segments as sites of strand-passage reaction enables type II DNA topoisomerases to select the most efficient pathway of relaxation of complex DNA knots. The analysis of the best selection criteria for efficient relaxation of complex knots revealed that local structures in random configurations of a given knot type statistically behave as analogous local structures in ideal geometric configurations of the corresponding knot type