Researcher profile

B. J. Powell

B. J. Powell contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
25works
0followers
13topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2022arXiv

$C_3$ symmetry breaking metal-insulator transitions in a flat band in the half-filled Hubbard model on the decorated honeycomb lattice

We study the single-orbital Hubbard model on the half-filled decorated honeycomb lattice. In the non-interacting theory at half-filling the Fermi energy lies within a flat band where strong correlations are enhanced. The lattice is highly frustrated. We find a correlation driven first-order metal-insulator transition to two different insulating ground states - a dimer valence bond solid Mott insulator when inter-triangle correlations dominate, and a broken $\mathcal{C}_3$-symmetry antiferromagnet that arises from frustration when intra-triangle correlations dominate. The metal-insulator transitions into these two phases have very different characters. The metal-broken $\mathcal{C}_3$ antiferromagnetic transition is driven by spontaneous $\mathcal{C}_3$ symmetry breaking that lifts the topologically required degeneracy at the Fermi energy and opens an energy gap in the quasiparticle spectrum. The metal-dimer valence bond solid transition breaks no symmetries of the Hamiltonian. It is caused by strong correlations renormalizing the electronic structure into a phase that is adiabatically connected to both the trivial band insulator and the ground state of the spin-1/2 Heisenberg model in the relevant parameter regime. Therefore, neither of these metal-insulator transitions can be understood in either the Brinkmann-Rice or Slater paradigms.

preprint2020arXiv

Frustration, ring exchange, and the absence of long-range order in EtMe$_3$Sb[Pd(dmit)$_2$]$_2$: from first principles to many-body theory

We parameterize Hubbard and spin models for EtMe$_3$Sb[Pd(dmit)$_2$]$_2$ from broken symmetry density functional calculations. This gives a scalene triangular model where the largest net exchange interaction is three times larger than the mean interchain coupling. The chain random phase approximation shows that the difference in the interchain couplings is equivalent to a bipartite interchain coupling, favoring long-range magnetic order. This competes with ring exchange, which favors quantum disorder. Ring exchange wins. We predict that the thermal conductivity, $κ$, along the chain direction is much larger than that along the crystallographic axes and that $κ/T\rightarrow0$ as $T\rightarrow0$ along the crystallographic axes, but that $κ/T\rightarrow\textrm{a constant}>0$ as $T\rightarrow0$ along the chain direction.

preprint2019arXiv

Interplay of dipoles and spins in $κ$-(BEDT-TTF)$_2X$, where $X=$ Hg(SCN)$_2$Cl, Hg(SCN)$_2$Br, Cu[N(CN)$_2$]Cl, Cu[N(CN)$_2$]Br, and Ag$_2$(CN)$_3$

We combine first principles density functional calculations with empirical relationships for the Coulomb interactions in the `monomer' model of $κ$-(BEDT-TTF)$_2X$. This enables us to calculate the parameters for the model of coupled dipolar and spin degrees of freedom proposed by Hotta [Phys. Rev. B $\textbf{82}$, 241104 (2010)], and Naka and Ishihara [J. Phys. Soc. Japan $\textbf{79}$, 063707 (2010)]. In all materials studied, retaining only the largest interactions leads to a transverse field Ising model of the dipoles. This quantifies, justifies and confirms recent claims that the dipoles are of crucial importance for understanding these materials. We show that two effects are responsible for a range of behaviors found in the dipoles in different $κ$-(BEDT-TTF)$_2X$ salts. (i) The inter-dimer hopping, $t_{b1}$, which gives rise to the "transverse field" in the Ising model for the dipoles ($H^T=2t_{b1}$), is between a third and a tenth smaller in the mercuric materials than for the mercury-free salts. (ii) The Ising model of dipoles is in the quasi-one-dimensional limit for the mercuric salts, but quasi-two-dimensional (between the square and isotropic triangular limits) for the mercury-free materials. Thus, the dimensionless critical field is much smaller in the mercuric compounds. Effect (ii) is much larger than effect (i). Simple explanations of both effects based on the band structures of the different salts are given. We show that dipolar order and even short-range dipolar correlations have a profound impact on the nature of the interdimer magnetic (superexchange) interactions. For example, dipole crystallization drives the materials towards quasi-one-dimensional magnetic interactions, which could be important for understanding the spin liquids found in some of these materials.

preprint2019arXiv

Spin-state ice in geometrically frustrated spin-crossover materials

Spin crossover materials contain metal ions that can access two spin-states: one low-spin (LS), the other high-spin (HS). We propose that frustrated elastic interactions can give rise to spin-state ices -- phases of matter without long-range order, characterized by a local constraint or `ice rule'. The low-energy physics of spin-state ices is described by an emergent divergence-less gauge field with a gap to topological excitations that are deconfined quasi-particles with spin fractionalized midway between the spins of the LS and HS states.

preprint2018arXiv

A theory of the quantum spin liquid in the hyper-honeycomb metal-organic framework [(C$_2$H$_5$)$_3$NH]$_2$Cu$_2$(C$_2$O$_4$)$_3$ from first principles

We construct a tight-binding model of [(C$_2$H$_5$)$_3$NH]$_2$Cu$_2$(C$_2$O$_4$)$_3$ from Wannier orbital overlaps. Including interactions within the Jahn-Teller distorted Cu-centered $e_g$ Wannier orbitals leads to an effective Heisenberg model. The hyper-honeycomb lattice contains two symmetry distinct sublattices of Cu atoms arranged in coupled chains. One sublattice is strongly dimerized, the other forms isotropic antiferromagnetic chains. Integrating out the strongest (intradimer) exchange interactions leaves extremely weakly coupled Heisenberg chains, consistent with the observed low temperature physics.

preprint2015arXiv

Conservation laws, radiative decay rates, and excited state localization in organometallic complexes with strong spin-orbit coupling

There is longstanding fundamental interest in 6-fold coordinated $d^6$ ($t_{2g}^6$) transition metal complexes such as [Ru(bpy)$_3$]$^{2+}$ and Ir(ppy)$_3$, particularly their phosphorescence. This interest has increased with the growing realisation that many of these complexes have potential uses in applications including photovoltaics, imaging, sensing, and light-emitting diodes. In order to design new complexes with properties tailored for specific applications a detailed understanding of the low-energy excited states, particularly the lowest energy triplet state, $T_1$, is required. Here we describe a model of pseudo-octahedral complexes based on a pseudo-angular momentum representation and show that the predictions of this model are in excellent agreement with experiment - even when the deviations from octahedral symmetry are large. This model gives a natural explanation of zero-field splitting of $T_1$ and of the relative radiative rates of the three sublevels in terms of the conservation of time-reversal parity and total angular momentum modulo two. We show that the broad parameter regime consistent with the experimental data implies significant localization of the excited state.

preprint2015arXiv

Theories of phosphorescence in organo-transition metal complexes - from relativistic effects to simple models and design principles for organic light-emitting diodes

We review theories of phosphorescence in cyclometalated complexes. We focus primarily on pseudooctahedrally coordinated $t_{2g}^6$ metals (e.g., [Os(II)(bpy)$_3$]$^{2+}$, Ir(III)(ppy)$_3$ and Ir(III)(ptz)$_3$) as, for reasons that are explored in detail, these show particularly strong phosphorescence. We discuss both first principles approaches and semi-empirical models, e.g., ligand field theory. We show that together these provide a clear understanding of the photophysics and in particular the lowest energy triplet excitation, T$_1$. In order to build a good model relativistic effects need to be included. The role of spin-orbit coupling is well-known, but scalar relativistic effects are also large - and are therefore also introduced and discussed. No expertise in special relativity or relativistic quantum mechanics is assumed and a pedagogical introduction to these subjects is given. It is shown that, once both scalar relativistic effects and spin-orbit coupling are included, time dependent density functional theory (TDDFT) provides quantitatively accurate predictions of the radiative decay rates of the substates of T$_1$ in phosphorescent organotranstion-metal complexes. We describe the pseudo-angular momentum model, and show that it reproduces the key experimental findings. For example, this model provides a simple explanation of the relative radiative rates of the substates of T$_1$, which differ by orders of magnitude. Special emphasis is placed on materials with potential applications as active materials in organic light-emitting diodes (OLEDs) and principles for the design of new complexes are identified on the basis of the insights provided by the theories reviewed. We discuss the remaining theoretical challenges, which include deepening our understanding of solvent effects and, vitally, understanding and predicting non-radiative decay rates.

preprint2015arXiv

When is the Kadowaki-Woods ratio universal?

We calculate the Kadowaki-Woods ratio (KWR) in Fermi liquids with arbitrary band structures. We find that, contrary to the single band case, the ratio is not generally independent of the effects of electronic correlations (universal). This is very surprising given the experimental findings of a near universal KWR in many multiband metals. We identify a limit where the universality of the ratio, which has been observed experimentally in many strongly correlated electron systems, is recovered. We discuss the KWR in Dirac semimetals in two and three dimensions. In the two-dimensional case we also generalize the KWR to account for the logarithmic factor in the self-energy. In both cases we find that the KWR is independent of correlations, but strongly dependent on the doping of the system: for massless fermions the KWR is proportional to the inverse square of the carrier density, whereas the KWR for systems with massive quasiparticles is proportional to the inverse of the carrier density.

preprint2014arXiv

Topological spin liquid phase in a low-dimensional organic molecular compound

We report the discovery of a correlated insulator with a bulk gap at two-thirds filling in a geometrically frustrated Hubbard model that describes the low-energy physics of Mo$_3$S$_7$(dmit)$_3$. This is very different from the Mott insulator expected at half-filling. We show that the insulating phase, which persists even for very weak electron-electron interactions ($U$), is adiabatically connected to the Haldane phase and is consistent with experiments on Mo$_3$S$_7$(dmit)$_3$

preprint2013arXiv

In-plane superfluid density and microwave conductivity of the organic superconductor κ-(BEDT-TTF)2Cu[N(CN)2]Br: evidence for d-wave pairing and resilient quasiparticles

We report the in-plane microwave surface impedance of a high quality single crystal of κ-(BEDT-TTF)2Cu[N(CN)2]Br. In the superconducting state, we find three independent signatures of d-wave pairing: i) a strong, linear temperature dependence of superfluid density; ii) deep in the superconducting state the quasiparticle scattering rate Γ~ T^3; and iii) no BCS coherence peak is observed in the quasiparticle conductivity. Above T_c, the Kadowaki-Woods ratio and the temperature dependence of the in-plane conductivity show that the normal state is a Fermi liquid below ~ 23 K, yet resilient quasiparticles dominate the transport up to ~ 50 K.

preprint2012arXiv

Geometrical frustration in the spin liquid beta'-Me3EtSb[Pd(dmit)2]2 and the valence bond solid Me3EtP[Pd(dmit)2]2

We show that the electronic structures of the title compounds predicted by density functional theory (DFT) are well described by tight binding models. We determine the frustration ratio, J&#39;/J, of the Heisenberg model on the anisotropic triangular lattice, which describes the spin degrees of freedom in the Mott insulating phase for a range of Pd(dmit)2 salts. All of the antiferromagnetic materials studied have J&#39;/J < 0.5 or J&#39;/J > 0.9, consistent with predictions for the Heisenberg model. All salts with 0.5 < J&#39;/J < 0.9, where many-body theories find a number of competing ground states, are known, experimentally, to be charge ordered, valence bond solids or spin liquids.

preprint2011arXiv

Electronic correlations in organometallic complexes

We investigate an effective model for organometallic complexes (with potential uses in optoelectronic devices) via both exact diagonalisation and the configuration interaction singles (CIS) approximation. This model captures a number of important features of organometallic complexes, notably the sensitivity of the radiative decay rate to small chemical changes. We find that for large parameter ranges the CIS approximation accurately reproduces the low energy excitations and hence the photophysical properties of the exact solution. This suggests that electronic correlations do \emph{not} play an important role in these complexes. This explains why time-dependent density functional theory works surprisingly well in these complexes.

preprint2011arXiv

Quantum frustration in organic Mott insulators: from spin liquids to unconventional superconductors

We review the interplay of frustration and strong electronic correlations in quasi-two-dimensional organic charge transfer salts, such as k-(BEDT-TTF)_2X and Et_nMe_{4-n}Pn[Pd(dmit)2]2. These two forces drive a range of exotic phases including spin liquids, valence bond crystals, pseudogapped metals, and unconventional superconductivity. Of particular interest is that in several materials there is a direct transition as a function of pressure from a spin liquid Mott insulating state to a superconducting state. Experiments on these materials raise a number of profound questions about the quantum behaviour of frustrated systems, particularly the intimate connection between spin liquids and superconductivity. Insights into these questions have come from a wide range of theoretical techniques including first principles electronic structure, quantum many-body theory and quantum field theory. In this review we introduce the basic ideas of the field by discussing a simple frustrated Heisenberg model with four spins. We then describe the key experimental results, emphasizing that for two materials, k-(BEDT-TTF)_2Cu_2(CN)_3 and EtMe_3Sb[Pd(dmit)_2]_2, there is strong evidence for a spin liquid ground state, and for EtMe_3P[Pd(dmit)_2]_2, a valence bond solid ground state. We review theoretical attempts to explain these phenomena, arguing that this can be captured by a Hubbard model on the anisotropic triangular lattice at half filling, and that resonating valence bond wavefunctions can capture most of the essential physics. We review evidence that this model can have a spin liquid ground state for a range of parameters that are realistic for the relevant materials. We conclude by summarising the progress made thus far and identifying some of the key questions still to be answered.

preprint2010arXiv

An introduction to effective low-energy Hamiltonians in condensed matter physics and chemistry

These lecture notes introduce some simple effective Hamiltonians (also known as semi-empirical models) that have widespread applications to solid state and molecular systems. They are aimed as an introduction to a beginning graduate student. I also hope that it may help to break down the divide between the physics and chemistry literatures. After a brief introduction to second quantisation notation, which is used extensively, I focus of the &#34;four H&#39;s&#34;: the Huckel (or tight binding), Hubbard, Heisenberg and Holstein models. Some other related models, such as the Pariser-Parr-Pople model, the extended Hubbard model, multi-orbital models and the ionic Hubbard model, are also discussed. Finally, I discuss the epistemological basis of effective Hamiltonians and compare and contrast this with that of ab initio methods as well as discussing the problem of parametrising effective Hamiltonians.

preprint2010arXiv

Competition between Superconductivity and Weak Localization in Metal-Mixed Ion-Implanted Polymers

We study the effects of varying the pre-implant film thickness and implant temperature on the electrical and superconducting properties of metal-mixed ion-implanted polymers. We show that it is possible to drive a superconductor-insulator transition in these materials via control of the fabrication parameters. We observe peaks in the magnetoresistance and demonstrate that these are caused by the interplay between superconductivity and weak localization in these films, which occurs due to their granular structure. We compare these magnetoresistance peaks with those seen in unimplanted films and other organic superconductors, and show that they are distinctly different.

preprint2010arXiv

Sensitivity of the photo-physical properties of organometallic complexes to small chemical changes

We investigate an effective model Hamiltonian for organometallic complexes that are widely used in optoelectronic devices. The two most important parameters in the model are $J$, the effective exchange interaction between the $π$ and $π^*$ orbitals of the ligands, and $ε^*$, the renormalized energy gap between the highest occupied orbitals on the metal and on the ligand. We find that the degree of metal-to-ligand charge transfer (MLCT) character of the lowest triplet state is strongly dependent on the ratio $ε^*/J$. $ε^*$ is purely a property of the complex and can be changed significantly by even small variations in the complex&#39;s chemistry, such as replacing substituents on the ligands. We find that that small changes in $ε^*/J$ can cause large changes in the properties of the complex, including the lifetime of the triplet state and the probability of injected charges (electrons and holes) forming triplet excitations. These results give some insight into the observed large changes in the photophysical properties of organometallic complexes caused by small changes in the ligands.

preprint2009arXiv

Effective Coulomb interactions within BEDT-TTF dimers

We calculate the effective Coulomb interactions between holes in dimers of the organic molecule BEDT-TTF in vacuo. We use density functional theory (DFT) to parameterise Hubbard models for beta and kappa phase organic charge transfer salts. We focus on the intra-dimer Coulomb repulsion, U_d, and the inter-monomer Coulomb repulsion, V_m. We find that U_d = 3.22 \pm 0.09 eV and V_m = 2.71 \pm 0.10 eV for 23 experimental geometries taken from a range of materials in both polymorphs. The quoted error is one standard deviation over the set of conformations studied. We conclude that U_d and V_m are essentially the same for an isolated dimer with the geometries present in all of the compounds studied. We highlight the disagreement between our parameterisation of the Hubbard model and previous results from both DFT and Huckel methods and show that this is caused by the failure of an assumption made in previous calculations (that U_m >> V_m, where U_m is the effective intra-monomer Coulomb repulsion). We discuss the implications of our calculations for theories of the BEDT-TTF salts based on the Hubbard model on the 2D anisotropic triangular lattice and explore the role of conformational disorder in these materials.

preprint2009arXiv

Electronic and magnetic properties of the ionic Hubbard model on the striped triangular lattice at 3/4 filling

We report a detailed study of a model Hamiltonian which exhibits a rich interplay of geometrical spin frustration, strong electronic correlations, and charge ordering. The character of the insulating phase depends on the magnitude of Delta/|t| and on the sign of t. We find a Mott insulator for Delta >> U >> |t|; a charge transfer insulator for U >> Δ>> |t|; and a correlated covalent insulator for U >> Δ~ |t|. The charge transfer insulating state is investigated using a strong coupling expansion. The frustration of the triangular lattice can lead to antiferromagnetism or ferromagnetism depending on the sign of the hopping matrix element, t. We identify the &#34;ring&#34; exchange process around a triangular plaquette which determines the sign of the magnetic interactions. Exact diagonalization calculations are performed on the model for a wide range of parameters and compared to the strong coupling expansion. The regime U >> Δ~ |t| and t<0 is relevant to Na05CoO2. The calculated optical conductivity and the spectral density are discussed in the light of recent experiments on Na05CoO2.

preprint2009arXiv

Interplay of frustration, magnetism, charge ordering, and covalency in a model of Na0.5CoO2

We investigate an effective Hamiltonian for Na0.5CoO2 that includes the electrostatic potential due to the ordered Na ions and strong electronic correlations. This model displays a subtle interplay between metallic and insulating phases and between charge and magnetic order. For realistic parameters, the model predicts an insulating phase with similarities to a covalent insulator. We show that this interpretation gives a consistent explanation of experiments on Na0.5CoO2, including the small degree of charge ordering, the small charge gap, the large moment, and the optical conductivity.

preprint2009arXiv

Ionic Hubbard model on a triangular lattice for Na_0.5CoO_2, Rb_0.5CoO_2 and K_0.5CoO_2: Mean-field slave boson theory

We introduce a strongly correlated mean-field theory of the ionic Hubbard model on the triangular lattice with alternating stripes of site energy using Barnes-Coleman slave bosons. We study the paramagnetic phases of this theory at three quarters filling, where it is a model of Na_0.5CoO_2, Rb_0.5CoO_2, and K_0.5CoO_2. This theory has two bands of fermionic quasi-particles: one of which is filled or nearly filled and hence weakly correlated; the other is half-filled or nearly half-filled and hence strongly correlated. Further results depend strongly on the sign of the hopping integral, t. The light band is always filled for t>0, but only becomes filled for |Delta/t|>=1.5 for t<0, where Delta is the difference in the site energies of the two sublattices. A metal--charge transfer insulator transition occurs at |Delta/t|=5.0 for t>0 and |Delta/t|=8.0 for t<0. In the charge transfer insulator complete charge disproportionation occurs: one sublattice is filled and the other in half filled. We compare our results with exact diagonalisation calculations and experiments on Na_0.5CoO_2, and discuss the relevance of our results to Rb_0.5CoO_2 and K_0.5CoO_2. In particular we propose a resolution of seemingly contradictory experimental results on Na_0.5CoO_2. Many experiments suggest that there is a charge gap, yet quantum oscillations are observed suggesting the existence of quasiparticle states at arbitrarily low excitation energies. We argue that the heavy band is gapped while the light band, which contains less than 1 charge carrier per 100 unit cells, remains ungapped.

preprint2009arXiv

Spin Fluctuations and the Pseudogap in Organic Superconductors

We show that there are strong similarities in the spin lattice relaxation of non-magnetic organic charge transfer salts, and that these similarities can be understood in terms of spin fluctuations. Further, we show that, in all of the kappa-phase organic superconductors for which there is nuclear magnetic resonance data, the energy scale for the spin fluctuations coincides with the energy scale for the pseudogap. This suggests that the pseudogap is caused by short-range spin correlations. In the weakly frustrated metals k-(BEDT-TTF)_2Cu[N(CN)_2]Br, k-(BEDT-TTF)_2Cu(NCS)_2, and k-(BEDT-TTF)_2Cu[N(CN)_2]Cl (under pressure) the pseudogap opens at the same temperature as coherence emerges in the (intralayer) transport. We argue that this is because the spin correlations are cut off by the loss of intralayer coherence at high temperatures. We discuss what might happen to these two energy scales at high pressures, where the electronic correlations are weaker. In these weakly frustrated materials the data is well described by the chemical pressure hypothesis (that anion substitution is equivalent to hydrostatic pressure). However, we find important differences in the metallic state of k-(BEDT-TTF)_2Cu_2(CN)_3, which is highly frustrated and displays a spin liquid insulating phase. We also show that the characteristic temperature scale of the spin fluctuations in (TMTSF)_2ClO_4 is the same as superconducting critical temperature, which may be evidence that spin fluctuations mediate the superconductivity in the Bechgaard salts.

preprint2007arXiv

Antiferromagnetic Spin Fluctuations and the Pseudogap in the Paramagnetic Phases of Quasi-Two-Dimensional Organic Superconductors

We give a quantitative analysis of the published results of nuclear magnetic resonance (NMR) experiments for k-(ET)2X family of organic charge transfer salts by using the phenomenological spin fluctuation model of Moriya, and Millis, Monien and Pines (M-MMP). For temperatures above T_NMR ~ 50 K, the model gives a good quantitative description of the data for the paramagnetic metallic phase of several k-(ET)2X materials, with an antiferromagnetic correlation length which increases with decreasing temperature; growing to several lattice constants by T_NMR. It is shown that the fact that the dimensionless Korringa ratio is much larger than unity is inconsistent with a broad class of theoretical models (such as dynamical mean-field theory) which neglect spatial correlations and/or vertex corrections. For materials close to the Mott insulating phase, 1/T1, Ks, and K all decrease significantly with decreasing temperature below T_NMR. This cannot be described by the M-MMP model and the most natural explanation is that a pseudogap opens up in the density of states below T_NMR, as in, for example, the cuprates. We show that the NMR measurements reported for k-(ET)2Cu2(CN)3 are qualitatively inconsistent with this material having a ground state with long range magnetic order. A pseudogap in the metallic state of organic superconductors is an important prediction of the resonating valence bond theory of superconductivity. We propose specific new experiments on organic superconductors to elucidate the nature, origin, and momentum dependence of the pseudogap and its relationship to superconductivity.

preprint2005arXiv

Electron-phonon coupling as an order-one problem

The coupling between electrons and phonons plays important roles in physics, chemistry and biology. However, the accurate calculation of the electron-phonon coupling constants is computationally expensive as it involves solving the Schrodinger equation for O(N) nuclear configurations, where N is the number of nuclei. Herein we show that by considering the forces on the nuclei caused by the addition or subtraction of an arbitrarily small electronic charge one may calculate the electron-phonon coupling constants from O(1) solutions of the Schrodinger equation. We show that Janak&#39;s theorem means that this procedure is exact within the density functional formalism. We demonstrate that the O(1) approach produces numerically accurate results by calculating the electron-phonon coupling constants for a series of molecules ranging in size from H_2 to C_60. We use our approach to introduce a computationally fast approximation for the adiabatic ionisation potentials and electron affinities which is shown to be accurate for large molecules. We also show that our approach allows for the calculation isotope effects in O(0) time and discuss the deuteration driven Mott transition in kappa-(BEDT-TTF)_2Cu[N(CN)_2]Br.

preprint2004arXiv

5,6-dihydroxyindole-2-carboxylic acid (DHICA): a First Principles Density-Functional Study

We report first principles density functional calculations for 5,6-dihydroxyindole-2-carboxylic acid (DHICA) and several reduced forms. DHICA and 5,6-dihydroxyindole (DHI) are believed to be the basic building blocks of the eumelanins. Our results show that carboxylation has a significant effect on the physical properties of the molecules. In particular, the relative stabilities and the HOMO-LUMO gaps (calculated with the $Δ$SCF method) of the various redox forms are strongly affected. We predict that, in contrast to DHI, the density of unpaired electrons, and hence the ESR signal, in DHICA is negligibly small.

preprint2004arXiv

A First Principles Density-Functional Calculation of the Electronic and Vibrational Structure of the Key Melanin Monomers

We report first principles density functional calculations for hydroquinone (HQ), indolequinone (IQ) and semiquinone (SQ). These molecules are believed to be the basic building blocks of the eumelanins, a class of bio-macromolecules with important biological functions (including photoprotection) and with potential for certain bioengineering applications. We have used the DeltaSCF (difference of self consistent fields) method to study the energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), Delta_HL. We show that Delta_HL is similar in IQ and SQ but approximately twice as large in HQ. This may have important implications for our understanding of the observed broad band optical absorption of the eumelanins. The possibility of using this difference in Delta_HL to molecularly engineer the electronic properties of eumelanins is discussed. We calculate the infrared and Raman spectra of the three redox forms from first principles. Each of the molecules have significantly different infrared and Raman signatures, and so these spectra could be used in situ to non-destructively identify the monomeric content of macromolecules. It is hoped that this may be a helpful analytical tool in determining the structure of eumelanin macromolecules and hence in helping to determine the structure-property-function relationships that control the behaviour of the eumelanins.