Showing 1–50 of 140 results for author: McKenzie, R H

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

Authors: H. L. Nourse, Ross H. McKenzie, B. J. Powell

Abstract: 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 insul… ▽ More 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. △ Less

Submitted 17 May, 2022; v1 submitted 28 February, 2022; originally announced March 2022.

Journal ref: Phys. Rev. B 105, 205119 (2022)

Spin-$0$ Mott insulator to metal to spin-$1$ Mott insulator transition in the single-orbital Hubbard model on the decorated honeycomb lattice

Authors: H. L. Nourse, Ross H. McKenzie, B. J. Powell

Abstract: We study the interplay of strong electron correlations and intra-triangle spin exchange at two-thirds filling of the single-orbital Hubbard model on the decorated honeycomb lattice using rotationally invariant slave bosons. We find that the spin exchange tunes between a spin-$1$ Mott insulator, a metal, and a spin-$0$ Mott insulator when the exchange is antiferromagnetic. The Mott insulators occur… ▽ More We study the interplay of strong electron correlations and intra-triangle spin exchange at two-thirds filling of the single-orbital Hubbard model on the decorated honeycomb lattice using rotationally invariant slave bosons. We find that the spin exchange tunes between a spin-$1$ Mott insulator, a metal, and a spin-$0$ Mott insulator when the exchange is antiferromagnetic. The Mott insulators occur from effective intra-triangle multi-orbital interactions and are adiabatically connected to the ground state of an isolated triangle. An antiferromagnetic spin exchange, as determined by the Goodenough-Kanamori rules, may occur in coordination polymers from kinetic exchange via the ligands. We characterize the magnetism in the regime where spin-triplets dominate. For small $U$ a spin-$1$ Slater insulator occurs with antiferromagnetic order between triangles. Magnetism in the spin-$1$ Mott insulator is described by a spin-$1$ Heisenberg model on a honeycomb lattice, whose ground state is Néel ordered. △ Less

Submitted 4 August, 2021; v1 submitted 1 April, 2021; originally announced April 2021.

Journal ref: Phys. Rev. B 104, 075104 (2021)

Hidden Devil's staircase in a two-dimensional elastic model of spin crossover materials

Authors: Gian Ruzzi, Jace Cruddas, Ross H. McKenzie, Ben J. Powell

Abstract: Spin crossover (SCO) materials are reversible molecular switches found in a wide range of transition metal complexes and metal organic frameworks (MOFs). They exhibit diverse spin state orderings and transitions between them. We present an exact mapping from an elastic lattice mismatch model to a long-range Ising model, with an inverse square decay of the interaction strengths at large distances (… ▽ More Spin crossover (SCO) materials are reversible molecular switches found in a wide range of transition metal complexes and metal organic frameworks (MOFs). They exhibit diverse spin state orderings and transitions between them. We present an exact mapping from an elastic lattice mismatch model to a long-range Ising model, with an inverse square decay of the interaction strengths at large distances (on the square lattice). This provides a microscopic justification for an Ising model description, which has previously only been justified on phenomenological grounds. Elastic frustration is required for non-zero Ising interactions, but whether or not the short-range interactions in the Ising model are geometrically frustrated depends on the ratio of the bulk and shear moduli or equivalently Poisson's ratio. We show that, for a simple square lattice model with realistic parameters, sweeping the enthalpy difference between the two spin-states at zero temperature leads to a large (probably infinite) number of spin-state orderings and corresponding steps in the fraction of high-spin ions, consistent with a Devil's staircase. The staircase can also be climbed by varying the temperature, but then some of the steps are hidden and only a finite number remain, consistent with experiments on relevant framework materials, such as {(Fe[Hg(SCN)$_3$]$_2$(4,4'-bipy)$_2$)}$_n$. Our results are also relevant to other binary systems with lattice mismatch, e.g., heterogeneous solids. △ Less

Submitted 26 May, 2021; v1 submitted 19 August, 2020; originally announced August 2020.

Comments: Added discussion about a devil's staircase found in model Updated sections 3 and 4 Expanded section 5 (Collective behaviour) Added appendix B about isotope effect

Multiple insulating phases due to the interplay of strong correlations and lattice geometry in a single-orbital Hubbard model

Authors: H. L. Nourse, Ross H. McKenzie, B. J. Powell

Abstract: We report ten ground states arising from strong correlations in the single-orbital Hubbard model on the decorated honeycomb lattice; including Dirac metals, flat-band ferromagnets, real-space Mott insulators, dimer and trimer Mott insulators, and a spin-$1$ Mott insulator. The rich phase diagram arises from structures within the unit cell. Hence, such states are absent on simpler lattices. We argu… ▽ More We report ten ground states arising from strong correlations in the single-orbital Hubbard model on the decorated honeycomb lattice; including Dirac metals, flat-band ferromagnets, real-space Mott insulators, dimer and trimer Mott insulators, and a spin-$1$ Mott insulator. The rich phase diagram arises from structures within the unit cell. Hence, such states are absent on simpler lattices. We argue that such insulating phases are prevalent on decorated lattices. These are found in many materials and common in coordination polymers, providing a playground to explore this physics. △ Less

Submitted 1 March, 2021; v1 submitted 10 March, 2020; originally announced March 2020.

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

Orbital-selective bad metals due to Hund's rule and orbital anisotropy: a finite-temperature slave-spin treatment of the two-band Hubbard model

Authors: Alejandro Mezio, Ross H. McKenzie

Abstract: We study the finite-temperature properties of the half-filled two-band Hubbard model in the presence of Hund's rule coupling and orbital anisotropy. We use the mean-field treatment of the $Z_2$ slave-spin theory with a finite-temperature extension of the zero-temperature gauge variable previously developed by Hassan and de' Medici http://dx.doi.org/10.1103/PhysRevB.81.035106}{Phys. Rev. B {\bf 81}… ▽ More We study the finite-temperature properties of the half-filled two-band Hubbard model in the presence of Hund's rule coupling and orbital anisotropy. We use the mean-field treatment of the $Z_2$ slave-spin theory with a finite-temperature extension of the zero-temperature gauge variable previously developed by Hassan and de' Medici http://dx.doi.org/10.1103/PhysRevB.81.035106}{Phys. Rev. B {\bf 81}, 035106 (2010)}]. We consider the instability of the Fermi liquid phases and how it is enhanced by the Hund's rule. We identify paramagnetic solutions that have zero quasi-particle weight with bad metallic phases, and the first-order transition temperature between it and the Fermi liquid phase as a coherence temperature that signals the crossover to the bad metallic state. When orbital anisotropy is present, we found an intermediate transition to an orbital-selective bad metal (OSBM), where the narrow band becomes a bad metal while the wide band remains a renormalised Fermi liquid. The temperatures $T_\textrm{coh}$ and $T_\textrm{OSBM}$ at which the system transitions to the bad metal phases can be orders of magnitude less than the Fermi temperature associated with the non-interacting band. The parameter dependence of the temperature at which the OSBM is destroyed can be understood in terms of a ferromagnetic Kondo-Hubbard lattice model. In general, Hund's rule coupling enhances the bad metallic phases, reduce interorbital charge fluctuations and increase spin fluctuations. The qualitative difference found in the ground state whether the Hund's rule is present or not, related to the degeneracy of the low energy manifold, is also maintained for finite temperatures. △ Less

Submitted 2 April, 2019; v1 submitted 17 March, 2019; originally announced March 2019.

Comments: 16 pages, 14 figures, 1 table

Journal ref: Phys. Rev. B 100, 205134 (2019)

Low quasiparticle coherence temperature in the one band-Hubbard model: A slave-boson approach

Authors: Alejandro Mezio, Ross H. McKenzie

Abstract: We use the Kotliar-Ruckenstein slave-boson formalism to study the temperature dependence of paramagnetic phases of the one-band Hubbard model for a variety of band structures. We calculate the Fermi liquid quasiparticle spectral weight $Z$ and identify the temperature at which it decreases significantly to a crossover to a bad metal region. Near the Mott metal-insulator transition, this coherence… ▽ More We use the Kotliar-Ruckenstein slave-boson formalism to study the temperature dependence of paramagnetic phases of the one-band Hubbard model for a variety of band structures. We calculate the Fermi liquid quasiparticle spectral weight $Z$ and identify the temperature at which it decreases significantly to a crossover to a bad metal region. Near the Mott metal-insulator transition, this coherence temperature $T_\textrm{coh}$ is much lower than the Fermi temperature of the uncorrelated Fermi gas, as is observed in a broad range of strongly correlated electron materials. After a proper rescaling of temperature and interaction, we find a universal behavior that is independent of the band structure of the system. We obtain the temperature-interaction phase diagram as a function of doping, and we compare the temperature dependence of the double occupancy, entropy, and charge compressibility with previous results obtained with Dynamical Mean-Field Theory. We analyse the stability of the method by calculating the charge compressibility. △ Less

Submitted 13 July, 2017; v1 submitted 9 February, 2017; originally announced February 2017.

Comments: 12 figures, 11 pages

Journal ref: Phys. Rev. B 96, 035121 (2017)

Quantum critical local spin dynamics near the Mott metal-insulator transition in infinite dimensions

Authors: Nagamalleswararao Dasari, N. S. Vidhyadhiraja, Mark Jarrell, Ross H. McKenzie

Abstract: Finding microscopic models for metallic states that exhibit quantum critical properties such as $ω/T$ scaling is a major theoretical challenge. We calculate the local dynamical spin susceptibility $χ(T,ω)$ for a Hubbard model at half filling using Dynamical Mean-Field Theory, which is exact in infinite dimensions. Qualitatively distinct behavior is found in the different regions of the phase diagr… ▽ More Finding microscopic models for metallic states that exhibit quantum critical properties such as $ω/T$ scaling is a major theoretical challenge. We calculate the local dynamical spin susceptibility $χ(T,ω)$ for a Hubbard model at half filling using Dynamical Mean-Field Theory, which is exact in infinite dimensions. Qualitatively distinct behavior is found in the different regions of the phase diagram: Mott insulator, Fermi liquid metal, bad metal, and a quantum critical region above the finite temperature critical point. The signature of the latter is $ω/T$ scaling where $T$ is the temperature. Our results are consistent with previous results showing scaling of the dc electrical conductivity and are relevant to experiments on organic charge transfer salts. △ Less

Submitted 8 March, 2017; v1 submitted 6 November, 2016; originally announced November 2016.

Comments: 9 pages, 9 figures

Journal ref: Phys. Rev. B 95, 165105 (2017)

Effect of hydrogen bonding on infrared absorption intensity

Authors: Bijyalaxmi Athokpam, Sai G. Ramesh, Ross H. McKenzie

Abstract: We consider how the infrared intensity of an O-H stretch in a hydrogen bonded complex varies as the strength of the H-bond varies from weak to strong. We obtain trends for the fundamental and overtone transitions as a function of donor-acceptor distance R, which is a common measure of H-bond strength. Our calculations use a simple two-diabatic state model that permits symmetric and asymmetric bond… ▽ More We consider how the infrared intensity of an O-H stretch in a hydrogen bonded complex varies as the strength of the H-bond varies from weak to strong. We obtain trends for the fundamental and overtone transitions as a function of donor-acceptor distance R, which is a common measure of H-bond strength. Our calculations use a simple two-diabatic state model that permits symmetric and asymmetric bonds, i.e. where the proton affinity of the donor and acceptor are equal and unequal, respectively. The dipole moment function uses a Mecke form for the free OH dipole moment, associated with the diabatic states. The transition dipole moment is calculated using one-dimensional vibrational eigenstates associated with the H-atom transfer coordinate on the ground state adiabatic surface of our model. Over 20-fold intensity enhancements for the fundamental are found for strong H-bonds, where there are significant non-Condon effects. The isotope effect on the intensity yields a non-monotonic H/D intensity ratio as a function of R, and is enhanced by the secondary geometric isotope effect. The first overtone intensity is found to vary non-monotonically with H-bond strength; strong enhancements are possible for strong H-bonds. Modifying the dipole moment through the Mecke parameters is found to have a stronger effect on the overtone than the fundamental. We compare our findings with those for specific molecular systems analysed through experiments and theory in earlier works. Our model results compare favourably for strong and medium strength symmetric H-bonds. However, for weak asymmetric bonds we find much smaller effects than in earlier work. △ Less

Submitted 14 June, 2016; originally announced June 2016.

Comments: 12 pages, 10 figures

Journal ref: Chemical Physics 488 & 489 (2017) Page 43

Spin-orbit coupling and odd-parity superconductivity in the quasi-one-dimensional compound Li$_{0.9}$Mo$_6$O$_{17}$

Authors: Christian Platt, Weejee Cho, Ross H. McKenzie, Ronny Thomale, Sri Raghu

Abstract: Previous theoretical studies [W. Cho, C. Platt, R. H. McKenzie, and S. Raghu, Phys. Rev. B 92, 134514 (2015); N. Lera and J. V. Alvarez, Phys. Rev. B 92, 174523 (2015)] have suggested that Li$_{0.9}$Mo$_6$O$_{17}$, a quasi-one dimensional "purple bronze" compound, exhibits spin-triplet superconductivity and that the gap function changes sign across the two nearly degenerate Fermi surface sheets. W… ▽ More Previous theoretical studies [W. Cho, C. Platt, R. H. McKenzie, and S. Raghu, Phys. Rev. B 92, 134514 (2015); N. Lera and J. V. Alvarez, Phys. Rev. B 92, 174523 (2015)] have suggested that Li$_{0.9}$Mo$_6$O$_{17}$, a quasi-one dimensional "purple bronze" compound, exhibits spin-triplet superconductivity and that the gap function changes sign across the two nearly degenerate Fermi surface sheets. We investigate the role of spin-orbit coupling (SOC) in determining the symmetry and orientation of the $d$-vector associated with the superconducting order parameter. We propose that the lack of local inversion symmetry within the four-atom unit cell leads to a staggered spin-orbit coupling analogous to that proposed for graphene, MoS$_2$, or SrPtAs. In addition, from a weak-coupling renormalization group treatment of an effective model Hamiltonian, we find that SOC favors the odd parity $A_{1u}$ state with $S_z = \pm 1$ over the $B$ states with $S_z=0$, where $z$ denotes the least-conducting direction. We discuss possible definitive experimental signatures of this superconducting state. △ Less

Submitted 18 January, 2016; originally announced January 2016.

Journal ref: Phys. Rev. B 93, 214515 (2016)

Spin-triplet superconductivity in a weak-coupling Hubbard model for the quasi-one-dimensional compound Li$_{0.9}$Mo$_6$O$_{17}$

Authors: Weejee Cho, Christian Platt, Ross H. McKenzie, S. Raghu

Abstract: The purple bronze Li$_{0.9}$Mo$_6$O$_{17}$ is of interest due to its quasi-one-dimensional electronic structure and the possible Luttinger liquid behavior resulting from it. For sufficiently low temperatures, it is a superconductor with a pairing symmetry that is still to be determined. To shed light on this issue, we analyze a minimal Hubbard model for this material involving four molybdenum orbi… ▽ More The purple bronze Li$_{0.9}$Mo$_6$O$_{17}$ is of interest due to its quasi-one-dimensional electronic structure and the possible Luttinger liquid behavior resulting from it. For sufficiently low temperatures, it is a superconductor with a pairing symmetry that is still to be determined. To shed light on this issue, we analyze a minimal Hubbard model for this material involving four molybdenum orbitals per unit cell near quarter filling, using asymptotically exact perturbative renormalization group methods. We find that spin triplet odd-parity superconductivity is the dominant instability. Approximate nesting properties of the two quasi-one-dimensional Fermi surfaces enhance certain second-order processes, which play crucial roles in determining the structure of the pairing gap. Notably, we find that the gap has accidental nodes, i.e. it has more sign changes than required by the point-group symmetry. △ Less

Submitted 9 November, 2015; v1 submitted 22 July, 2015; originally announced July 2015.

Comments: Updated

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

Holon-Doublon Binding as the Mechanism for the Mott transition

Authors: Peter Prelovšek, Jure Kokalj, Zala Lenarčič, Ross H. McKenzie

Abstract: We study the binding of a holon to a doublon in a half-filled Hubbard model as the mechanism of the zero-temperature metal-insulator transition. In a spin polarized system and a non-bipartite lattice a single holon-doublon (HD) pair exhibits a binding transition (e.g., on a face-centred cubic lattice), or a sharp crossover (e.g., on a triangular lattice) corresponding well to the standard Mott tra… ▽ More We study the binding of a holon to a doublon in a half-filled Hubbard model as the mechanism of the zero-temperature metal-insulator transition. In a spin polarized system and a non-bipartite lattice a single holon-doublon (HD) pair exhibits a binding transition (e.g., on a face-centred cubic lattice), or a sharp crossover (e.g., on a triangular lattice) corresponding well to the standard Mott transition in unpolarized systems. We extend the HD-pair study towards non-polarized systems by considering more general spin background and by treating the finite HD density within a BCS-type approximation. Both approaches lead to a discontinuous transition away from the fully polarized system and give density correlations consistent with numerical results on a triangular lattice. △ Less

Submitted 22 December, 2015; v1 submitted 6 May, 2015; originally announced May 2015.

Comments: 6 pages, 4 figures

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

Shear viscosity of strongly interacting fermionic quantum fluids

Authors: Nandan Pakhira, Ross H. McKenzie

Abstract: Eighty years ago Eyring proposed that the shear viscosity of a liquid, $η$, has a quantum limit $η\gtrsim n\hbar$ where $n$ is the density of the fluid. Using holographic duality and the AdS/CFT correspondence in string theory Kovtun, Son, and Starinets (KSS) conjectured a universal bound $\fracη{s}\geq \frac{\hbar}{4πk_{B}}$ for the ratio between the shear viscosity and the entropy density, $s$.… ▽ More Eighty years ago Eyring proposed that the shear viscosity of a liquid, $η$, has a quantum limit $η\gtrsim n\hbar$ where $n$ is the density of the fluid. Using holographic duality and the AdS/CFT correspondence in string theory Kovtun, Son, and Starinets (KSS) conjectured a universal bound $\fracη{s}\geq \frac{\hbar}{4πk_{B}}$ for the ratio between the shear viscosity and the entropy density, $s$. Using Dynamical Mean-Field Theory (DMFT) we calculate the shear viscosity and entropy density for a fermionic fluid described by a single band Hubbard model at half filling. Our calculated shear viscosity as a function of temperature is compared with experimental data for liquid $^{3}$He. At low temperature the shear viscosity is found to be well above the quantum limit and is proportional to the characteristic Fermi liquid $1/T^{2}$ dependence, where $T$ is the temperature. With increasing temperature and interaction strength $U$ there is significant deviation from the Fermi liquid form. Also, the shear viscosity violates the quantum limit near the crossover from coherent quasi-particle based transport to incoherent transport (the bad metal regime). Finally, the ratio of the shear viscosity to the entropy density is found to be comparable to the KSS bound for parameters appropriate to liquid $^{3}$He. However, this bound is found to be strongly violated in the bad metal regime for parameters appropriate to lattice electronic systems such as organic charge transfer salts. △ Less

Submitted 17 August, 2015; v1 submitted 27 April, 2015; originally announced April 2015.

Comments: Revised manuscript with added references, 14 pages 5 figures

Isotopic fractionation in proteins as a measure of hydrogen bond length

Authors: Ross H. McKenzie, Bijyalaxmi Athokpam, Sai Ramesh

Abstract: If a deuterated molecule containing strong intramolecular hydrogen bonds is placed in a hydrogenated solvent it may preferentially exchange deuterium for hydrogen. This preference is due to the difference between the vibrational zero-point energy for hydrogen and deuterium. It is found that the associated fractionation factor $Φ$ is correlated with the strength of the intramolecular hydrogen bonds… ▽ More If a deuterated molecule containing strong intramolecular hydrogen bonds is placed in a hydrogenated solvent it may preferentially exchange deuterium for hydrogen. This preference is due to the difference between the vibrational zero-point energy for hydrogen and deuterium. It is found that the associated fractionation factor $Φ$ is correlated with the strength of the intramolecular hydrogen bonds. This correlation has been used to determine the length of the H-bonds (donor-acceptor separation) in a diverse range of enzymes and has been argued to support the existence of short low-barrier H-bonds. Starting with a potential energy surface based on a simple diabatic state model for H-bonds we calculate $Φ$ as a function of the proton donor-acceptor distance $R$. For numerical results, we use a parameterization of the model for symmetric O-H.... O bonds. We consider the relative contributions of the O-H stretch vibration, O-H bend vibrations (both in plane and out of plane), tunnelling splitting effects at finite temperature, and the secondary geometric isotope effect. We compare our total $Φ$ as a function of $R$ with NMR experimental results for enzymes, and in particular with an empirical parametrisation $Φ(R)$, used previously to determine bond lengths. △ Less

Submitted 17 July, 2015; v1 submitted 20 April, 2015; originally announced April 2015.

Comments: Final version, accepted for publication in Journal of Chemical Physics. Minor changes, including more extensive discussion of relevant of model to proteins

Journal ref: J. Chem. Phys. 143, 044309 (2015)

Enhancement of the thermal expansion of organic charge transfer salts by strong electronic correlations

Authors: J. Kokalj, Ross H. McKenzie

Abstract: Organic charge transfer salts exhibit thermal expansion anomalies similar to those found in other strongly correlated electron systems. The thermal expansion can be anisotropic and have a non-monotonic temperature dependence. We show how these anomalies can arise from electronic effects and be significantly enhanced, particularly at temperatures below 100 K, by strong electronic correlations. For… ▽ More Organic charge transfer salts exhibit thermal expansion anomalies similar to those found in other strongly correlated electron systems. The thermal expansion can be anisotropic and have a non-monotonic temperature dependence. We show how these anomalies can arise from electronic effects and be significantly enhanced, particularly at temperatures below 100 K, by strong electronic correlations. For the relevant Hubbard model the thermal expansion is related to the dependence of the entropy on the parameters ($t$, $t'$, and $U$) in the Hamiltonian or the temperature dependence of bond orders and double occupancy. The latter are calculated on finite lattices with the Finite Temperature Lanczos Method. Although many features seen in experimental data, in both the metallic and Mott insulating phase, are described qualitatively, the calculated magnitude of the thermal expansion is smaller than that observed experimentally. △ Less

Submitted 4 November, 2014; originally announced November 2014.

Comments: 9 pages, 6 figures

Journal ref: Phys. Rev. B 91, 205121 (2015)

Valence-bond Non-equilibrium Solvation Model for a Twisting Monomethine Cyanine

Authors: Sean McConnell, Ross H. McKenzie, Seth Olsen

Abstract: We propose and analyze a two-state valence-bond model of non-equilibrium solvation effects on the excited-state twisting reaction of monomethine cyanines. Suppression of this reaction is thought responsible for environment-dependent fluorescence yield enhancement in these dyes. Fluorescence is quenched because twisting is accompanied via the formation of dark twisted intramolecular charge-transfer… ▽ More We propose and analyze a two-state valence-bond model of non-equilibrium solvation effects on the excited-state twisting reaction of monomethine cyanines. Suppression of this reaction is thought responsible for environment-dependent fluorescence yield enhancement in these dyes. Fluorescence is quenched because twisting is accompanied via the formation of dark twisted intramolecular charge-transfer (TICT) states. For monomethine cyanines, where the ground state is a superposition of structures with different bond and charge localization, there are two possible twisting pathways with different charge localization in the excited state. For parameters corresponding to symmetric monomethines, the model predicts two low-energy twisting channels on the excited-state surface that lead to a manifold of twisted intramolecular charge-transfer (TICT) states. For typical monomethines, twisting on the excited state surface will occur with a small barrier or no barrier. Changes in the solvation configuration can differentially stabilize TICT states in channels corresponding to different bonds, and that the position of a conical intersection between adiabatic states moves in response to solvation to stabilize either one channel or the other. There is a conical intersection seam that grows along the bottom of the excited-state potential with increasing solvent polarity. For monomethine cyanines with modest-sized terminal groups in moderately polar solution, the bottom of the excited-state potential surface is completely spanned by a conical intersection seam. △ Less

Submitted 28 January, 2015; v1 submitted 28 October, 2014; originally announced October 2014.

Comments: To appear in the Journal of Chemical Physics

Journal ref: J. Chem. Phys. 142, 084502 (2015)

Enhancement of the thermoelectric power by electronic correlations in bad metals: a study of the Kelvin formula

Authors: J. Kokalj, Ross. H. McKenzie

Abstract: In many strongly correlated electron metals the thermoelectric power has a non-monotonic temperature dependence and values that are orders of magnitude larger than for elemental metals. Kelvin proposed a particularly simple expression for the thermopower in terms of the temperature dependence of the chemical potential. We consider a Hubbard model on an anisotropic triangular lattice at half fillin… ▽ More In many strongly correlated electron metals the thermoelectric power has a non-monotonic temperature dependence and values that are orders of magnitude larger than for elemental metals. Kelvin proposed a particularly simple expression for the thermopower in terms of the temperature dependence of the chemical potential. We consider a Hubbard model on an anisotropic triangular lattice at half filling, a minimal effective Hamiltonian for several classes of organic charge transfer salts. The finite temperature Lanczos method is used to calculate the temperature dependence of the thermopower using the Kelvin formula. We find that electronic correlations significantly enhance the magnitude of the thermopower and lead to a non-monotonic temperature dependence. The latter reflects a crossover with increasing temperature from a Fermi liquid to a bad metal. Although, the Kelvin formula gives a semi-quantitative description of some experimental results it cannot describe the directional dependence of the sign of the thermopower in some materials. △ Less

Submitted 3 October, 2014; originally announced October 2014.

Comments: 10 pages, 3 figures

Journal ref: Phys. Rev. B 91, 125143 (2015)

Absence of a quantum limit to charge diffusion in bad metals

Authors: Nandan Pakhira, Ross H. McKenzie

Abstract: Good metals are characterised by diffusive transport of coherent quasi-particle states and the resistivity is much less than the Mott-Ioffe-Regel (MIR) limit, $\frac{ha}{e^{2}}$, where $a$ is the lattice constant. In bad metals, such as many strongly correlated electron materials, the resistivity exceeds the Mott-Ioffe-Regel limit and the transport is incoherent in nature. Hartnoll, loosely motiva… ▽ More Good metals are characterised by diffusive transport of coherent quasi-particle states and the resistivity is much less than the Mott-Ioffe-Regel (MIR) limit, $\frac{ha}{e^{2}}$, where $a$ is the lattice constant. In bad metals, such as many strongly correlated electron materials, the resistivity exceeds the Mott-Ioffe-Regel limit and the transport is incoherent in nature. Hartnoll, loosely motivated by holographic duality (AdS/CFT correspondence) in string theory, recently proposed a lower bound to the charge diffusion constant, $D \gtrsim \hbar v_{F}^{2}/(k_{B}T)$, in the incoherent regime of transport, where $v_F$ is the Fermi velocity and $T$ the temperature. Using dynamical mean field theory (DMFT) we calculate the charge diffusion constant in a single band Hubbard model at half filling. We show that in the strongly correlated regime the Hartnoll's bound is violated in the crossover region between the coherent Fermi liquid region and the incoherent (bad metal) local moment region. The violation occurs even when the bare Fermi velocity $v_F$ is replaced by its low temperature renormalised value, $v_F^*$.The bound is satisfied at all temperatures in the weakly and moderately correlated systems as well as in strongly correlated systems in the high temperature region where the resistivity is close to linear in temperature. Our calculated charge diffusion constant, in the incoherent regime of transport, also strongly violates a proposed quantum limit of spin diffusion, $D_{s} \sim 1.3 \hbar/m$, where $m$ is the fermion mass, experimentally observed and theoretically calculated in a cold degenerate Fermi gas in the unitary limit of scattering. △ Less

Submitted 15 February, 2015; v1 submitted 19 September, 2014; originally announced September 2014.

Comments: Final version, 11 pages, 6 figures, accepted for publication in Physical Review B

A diabatic state model for double proton transfer in hydrogen bonded complexes

Authors: Ross H. McKenzie

Abstract: Four diabatic states are used to construct a simple model for double proton transfer in hydrogen bonded complexes. Key parameters in the model are the proton donor-acceptor separation R and the ratio, D_1/D_2, between the proton affinity of a donor with one and two protons. Depending on the values of these two parameters the model describes four qualitatively different ground state potential energ… ▽ More Four diabatic states are used to construct a simple model for double proton transfer in hydrogen bonded complexes. Key parameters in the model are the proton donor-acceptor separation R and the ratio, D_1/D_2, between the proton affinity of a donor with one and two protons. Depending on the values of these two parameters the model describes four qualitatively different ground state potential energy surfaces, having zero, one, two, or four saddle points. In the limit D_2=D_1 the model reduces to two decoupled hydrogen bonds. As R decreases a transition can occur from a concerted to a sequential mechanism for double proton transfer. △ Less

Submitted 14 July, 2014; originally announced July 2014.

Comments: 13 pages, 4 figures

Journal ref: J. Chem. Phys. 141, 104314 (2014)

Effect of quantum nuclear motion on hydrogen bonding

Authors: Ross H. McKenzie, Christiaan Bekker, Bijyalaxmi Athokpam, Sai G. Ramesh

Abstract: This work considers how the properties of hydrogen bonded complexes, D-H....A, are modified by the quantum motion of the shared proton. Using a simple two-diabatic state model Hamiltonian, the analysis of the symmetric case, where the donor (D) and acceptor (A) have the same proton affinity, is carried out. For quantitative comparisons, a parametrization specific to the O-H....O complexes is used.… ▽ More This work considers how the properties of hydrogen bonded complexes, D-H....A, are modified by the quantum motion of the shared proton. Using a simple two-diabatic state model Hamiltonian, the analysis of the symmetric case, where the donor (D) and acceptor (A) have the same proton affinity, is carried out. For quantitative comparisons, a parametrization specific to the O-H....O complexes is used. The vibrational energy levels of the one-dimensional ground state adiabatic potential of the model are used to make quantitative comparisons with a vast body of condensed phase data, spanning a donor-acceptor separation (R) range of about 2.4-3.0 A, i.e., from strong to weak bonds. The position of the proton and its longitudinal vibrational frequency, along with the isotope effects in both are discussed. An analysis of the secondary geometric isotope effects, using a simple extension of the two-state model, yields an improved agreement of the predicted variation with R of frequency isotope effects. The role of the bending modes in also considered: their quantum effects compete with those of the stretching mode for certain ranges of H-bond strengths. In spite of the economy in the parametrization of the model used, it offers key insights into the defining features of H-bonds, and semi-quantitatively captures several experimental trends. △ Less

Submitted 21 April, 2014; v1 submitted 15 January, 2014; originally announced January 2014.

Comments: 12 pages, 8 figures. Notation clarified. Revised figure including the effect of bending vibrations on secondary geometric isotope effect. Final version, accepted for publication in Journal of Chemical Physics

Journal ref: J. Chem. Phys. 140, 174508 (2014)

Collective Modes and Nonlinear Acoustics in Superfluid 3He-B

Authors: Ross H. McKenzie, J. A. Sauls

Abstract: We discuss the relationship of collisionless sound propagation and attenuation to the order parameter collective modes of superfluid 3He-B. These modes, which reflect the symmetries of the normal state as well as the broken gauge and relative rotational symmetries of the condensate, play a central role in the high-frequency acoustic response of the superfluid. We review the theory of acoustic spec… ▽ More We discuss the relationship of collisionless sound propagation and attenuation to the order parameter collective modes of superfluid 3He-B. These modes, which reflect the symmetries of the normal state as well as the broken gauge and relative rotational symmetries of the condensate, play a central role in the high-frequency acoustic response of the superfluid. We review the theory of acoustic spectroscopy based on the linear coupling of the J=2+/- modes. We develop the theory of the nonlinear acoustic response of superfluid 3He-B for the case of a three-wave resonance between two sound waves of different frequencies and the J=2+/- modes, as well as the generation of higher harmonic sound waves. The J=2+ modes couple nonlinearly to two zero-sound waves. The predictions for the nonlinear absorption and velocity anomalies resulting from the J=2+ modes depend on relatively well understood material properties of 3He. △ Less

Submitted 23 September, 2013; originally announced September 2013.

Comments: 62 pages, 7 figures, with typographical corrections to the published version

Journal ref: Volume 26, Chapter 5, pp. 255-311 of Modern Problems in Condensed Matter Physics, Ed. by W. P. Halperin and L. P. Pitaevskii, North-Holland, 1990

Signatures of the Berry curvature in the frequency dependent interlayer magnetoresistance in tilted magnetic fields

Authors: Anthony R. Wright, Ross H. McKenzie

Abstract: We show that in a layered metal, the angle dependent, finite frequency, interlayer magnetoresistance is altered due to the presence of a non-zero Berry curvature at the Fermi surface. At zero frequency, we find a conservation law which demands that the `magic angle' condition for interlayer magnetoresistance extrema as a function of magnetic field tilt angle is essentially both field and Berry cur… ▽ More We show that in a layered metal, the angle dependent, finite frequency, interlayer magnetoresistance is altered due to the presence of a non-zero Berry curvature at the Fermi surface. At zero frequency, we find a conservation law which demands that the `magic angle' condition for interlayer magnetoresistance extrema as a function of magnetic field tilt angle is essentially both field and Berry curvature independent. In the finite frequency case, however, we find that surprisingly large signatures of a finite Berry curvature occur in the periodic orbit resonances. We outline a method whereby the presence and magnitude of the Berry curvature at the Fermi surface can be extracted. △ Less

Submitted 30 September, 2013; v1 submitted 19 June, 2013; originally announced June 2013.

Comments: We corrected an error in the original preprint. A term in the equations of motion was overlooked which leads to a cancellation of terms involving the Berry curvature in the magic angle condition. We show this surprising result follows from a strict conservation law. However, the finite frequency interlayer conductivity is altered by the Berry curvature due to its effect on the cyclotron frequency

Thermodynamics of a bad metal-Mott insulator transition in the presence of frustration

Authors: J. Kokalj, Ross H. McKenzie

Abstract: We study a range of thermodynamic properties (charge susceptibility, specific heat, entropy and spin susceptibility) of the Hubbard model on the anisotropic triangular lattice at half filling by means of the numerical finite-temperature Lanczos method. We observe clear signatures of a metal-Mott insulator transition in the charge susceptibility, and show that the metallic phase is characterized by… ▽ More We study a range of thermodynamic properties (charge susceptibility, specific heat, entropy and spin susceptibility) of the Hubbard model on the anisotropic triangular lattice at half filling by means of the numerical finite-temperature Lanczos method. We observe clear signatures of a metal-Mott insulator transition in the charge susceptibility, and show that the metallic phase is characterized by a small charge susceptibility, large entropy, low coherence temperature, large renormalized quasiparticle mass, and large spin susceptibility. We show that the local magnetic moment in the metallic phase is large and comparable to the local moment in the insulating phase. These are all characteristics of bad metallic behavior. In addition, we show how frustration increases the density of low-lying spin excitations in the Mott insulating phase and decreases longer range spin correlations. △ Less

Submitted 26 April, 2013; v1 submitted 6 November, 2012; originally announced November 2012.

Comments: 5+9 pages, 17 figures, updated version

Journal ref: Phys. Rev. Lett. 110, 206402 (2013)

Quantum oscillations and Berry's phase in topological insulator surface states with broken particle-hole symmetry

Authors: Anthony R. Wright, Ross H. McKenzie

Abstract: Quantum oscillations can be used to determine properties of the Fermi surface of metals by varying the magnitude and orientation of an external magnetic field. Topological insulator surface states are an unusual mix of normal and Dirac fermions. Unlike in graphene and simple metals, Berry's geometric phase in topological insulator surface states is not necessarily quantised. We show that reliably… ▽ More Quantum oscillations can be used to determine properties of the Fermi surface of metals by varying the magnitude and orientation of an external magnetic field. Topological insulator surface states are an unusual mix of normal and Dirac fermions. Unlike in graphene and simple metals, Berry's geometric phase in topological insulator surface states is not necessarily quantised. We show that reliably extracting this geometric phase from the phase offset associated with the quantum oscillations is subtle. This is especially so in the presence of a Dirac gap such as that associated with the Zeeman splitting or interlayer tunneling. We develop a semi-classical theory for general mixed normal-Dirac systems in the presence of a gap, and in doing so clarify the role of topology and broken particle-hole symmetry. We propose a systematic procedure of fitting Landau level index plots at large filling factors to reliably extract the phase offset associated with Berry's phase. △ Less

Submitted 29 October, 2012; v1 submitted 15 October, 2012; originally announced October 2012.

Comments: 12 pages, 6 figures. Included effect of bulk Fermi surface

Journal ref: Phys. Rev. B 87, 085411 (2013)

A two-state model of twisted intramolecular chargetransfer in monomethine dyes

Authors: Seth Olsen, Ross H. McKenzie

Abstract: A two-state model Hamiltonian is proposed to model the coupling of twisting displacements to charge-transfer behavior in the ground and excited states of a general monomethine dye molecule. This coupling may be relevant to the molecular mechanism of environment-dependent fluorescence yield enhancement. The model is parameterized against quantum chemical calculations on different protonation states… ▽ More A two-state model Hamiltonian is proposed to model the coupling of twisting displacements to charge-transfer behavior in the ground and excited states of a general monomethine dye molecule. This coupling may be relevant to the molecular mechanism of environment-dependent fluorescence yield enhancement. The model is parameterized against quantum chemical calculations on different protonation states of the green fluorescent protein chromophore (GFP), which are chosen to sample different regimes of detuning from the cyanine (resonant) limit. The model provides a simple yet realistic description of the charge transfer character along two possible excited state twisting channels associated with the methine bridge. It describes qualitatively different behavior in three regions that can be classified by their relationship to the resonant (cyanine) limit. The regimes differ by the presence or absence of twist-dependent polarization reversal and the occurrence of conical intersections. We find that selective biasing of one twisting channel over another by an applied diabatic biasing potential can only be achieved in a finite range of parameters near the cyanine limit. △ Less

Submitted 14 October, 2012; v1 submitted 12 August, 2012; originally announced August 2012.

Comments: 45 pages, 9 Figures (incl. 2 chemical schemes). Accepted for publication by the Journal of Chemical Physics. Changes include 2 additional figures to and expanded discussion of key points felt to be important, and condensed discussion of some points felt to be less important

Effective Hamiltonian for the electronic properties of the quasi-one-dimensional material Li0.9Mo6O17

Authors: Jaime Merino, Ross H. McKenzie

Abstract: The title material has a quasi-one-dimensional electronic structure and is of considerable interest because it has a metallic phase with properties different from a simple Fermi liquid, a poorly understood "insulating" phase, and a superconducting phase which may involve spin triplet Cooper pairs. Using the Slater-Koster approach and comparison with published band structure calculations we present… ▽ More The title material has a quasi-one-dimensional electronic structure and is of considerable interest because it has a metallic phase with properties different from a simple Fermi liquid, a poorly understood "insulating" phase, and a superconducting phase which may involve spin triplet Cooper pairs. Using the Slater-Koster approach and comparison with published band structure calculations we present the simplest possible tight-binding model for the electronic band structure near the Fermi energy. This describes a set of ladders with weak (and frustrated) inter-ladder hopping. In the corresponding lattice model the system is actually close to one-quarter filling (i.e., one electron per pair of sites) rather than half-filling, as has often been claimed. We consider the simplest possible effective Hamiltonian that may capture the subtle competition between unconventional superconducting, charge ordered, and non-Fermi liquid metal phases. We argue that this is an extended Hubbard model with long-range Coulomb interactions. Estimates of the relevant values of the parameters in the Hamiltonian are given. NMR relaxation rate experiments should be performed to clarify the role of charge fluctuations in Li0.9Mo6O17 associated with the proximity to a Coulomb driven charge ordering transition. △ Less

Submitted 25 April, 2012; originally announced April 2012.

Comments: 15 pages, 7 figures, submitted to PRB

Journal ref: Phys. Rev. B 85, 235128 (2012)

An Effective Hamiltonian for Symmetric Diarylmethanes from a Series of Analogous Quantum Chemical Models

Authors: Seth Olsen, Ross H. McKenzie

Abstract: We propose a single effective Hamiltonian to describe the low-energy electronic structure of a series of symmetric cationic diarylmethanes, which are all bridge-substituted derivatives of Michler's Hydrol Blue. Three-state diabatic Hamiltonians for the dyes are calculated using four-electron three-orbital state-averaged complete active space self-consistent field and multi-state multi-reference pe… ▽ More We propose a single effective Hamiltonian to describe the low-energy electronic structure of a series of symmetric cationic diarylmethanes, which are all bridge-substituted derivatives of Michler's Hydrol Blue. Three-state diabatic Hamiltonians for the dyes are calculated using four-electron three-orbital state-averaged complete active space self-consistent field and multi-state multi-reference perturbation theory models. The approach takes advantage of an isolobal analogy that can be established between the orbitals spanning the active spaces of the different substituted dyes. The solutions of the chemical problem are expressed in a diabatic Hilbert space that is analogous to classical resonance models. The effective Hamiltonians for all dyes can be fit to a single functional form that depends on the mixing angle between a bridge-charged diabatic state and a superposition representing the canonical resonance. We find that the structure of the bridge-charged state changes in a regular fashion across the series. The change is consistent with an inversion of the sign of the charge carrier on the bridge, which changes from an electron pair to a hole as the series is traversed. △ Less

Submitted 26 March, 2012; originally announced March 2012.

Comments: 24 pages, 8 figures

Transport properties of the metallic state of overdoped cuprate superconductors from an anisotropic marginal Fermi liquid model

Authors: Jure Kokalj, Nigel E. Hussey, Ross H. McKenzie

Abstract: We consider the implications of a phenomenological model self-energy for the charge transport properties of the metallic phase of the overdoped cuprate superconductors. The self-energy is the sum of two terms with characteristic dependencies on temperature, frequency, location on the Fermi surface, and doping. The first term is isotropic over the Fermi surface, independent of doping, and has the f… ▽ More We consider the implications of a phenomenological model self-energy for the charge transport properties of the metallic phase of the overdoped cuprate superconductors. The self-energy is the sum of two terms with characteristic dependencies on temperature, frequency, location on the Fermi surface, and doping. The first term is isotropic over the Fermi surface, independent of doping, and has the frequency and temperature dependence characteristic of a Fermi liquid. The second term is anisotropic over the Fermi surface (vanishing at the same points as the superconducting energy gap), strongly varies with doping (scaling roughly with $T_c$, the superconducting transition temperature), and has the frequency and temperature dependence characteristic of a marginal Fermi liquid. Previously it has been shown this self-energy can describe a range of experimental data including angle-dependent magnetoresistance (ADMR) and quasi-particle renormalisations determined from specific heat, quantum oscillations, and angle-resolved photo-emission spectroscopy (ARPES). Without introducing new parameters and neglecting vertex corrections we show that this model self-energy can give a quantitative description of the temperature and doping dependence of a range of reported transport properties of Tl2201 samples. These include the intra-layer resistivity, the frequency dependent optical conductivity, the intra-layer magnetoresistance, and the Hall coefficient. The temperature dependence of the latter two are particularly sensitive to the anisotropy of the scattering rate and to the shape of the Fermi surface. In contrast, the temperature dependence of the Hall angle is dominated by the Fermi liquid contribution to the self-energy that determines the scattering rate in the nodal regions of the Fermi surface. △ Less

Submitted 10 July, 2012; v1 submitted 22 February, 2012; originally announced February 2012.

Comments: 17 pages, 16 figures

Journal ref: Phys. Rev. B 86, 045132 (2012)

Unified description of hydrogen bonding and of proton transfer by a two-state effective Hamiltonian

Authors: Ross H. McKenzie

Abstract: An effective Hamiltonian is considered for hydrogen bonding between two molecules due to the quantum mechanical interaction between the orbitals of the H-atom and the donor and acceptor atoms in the molecules. The Hamiltonian acts on two diabatic states and has a simple chemically motivated form for its matrix elements. The model gives insight into the "H-bond puzzle", describes different classes… ▽ More An effective Hamiltonian is considered for hydrogen bonding between two molecules due to the quantum mechanical interaction between the orbitals of the H-atom and the donor and acceptor atoms in the molecules. The Hamiltonian acts on two diabatic states and has a simple chemically motivated form for its matrix elements. The model gives insight into the "H-bond puzzle", describes different classes of bonds, and empirical correlations between the donor-acceptor distance R and binding energies, bond lengths, and the softening of vibrational frequencies. A key prediction is the UV photo-dissociation of H-bonded complexes via an excited electronic state with an exalted vibrational frequency. △ Less

Submitted 21 September, 2011; v1 submitted 9 September, 2011; originally announced September 2011.

Comments: 5 pages, 7 figures, minor changes

Quantum entanglement between electronic and vibrational degrees of freedom in molecules

Authors: Laura K. McKemmish, Ross H. McKenzie, Noel S. Hush, Jeffrey R. Reimers

Abstract: We consider the quantum entanglement of the electronic and vibrational degrees of freedom in molecules with a tendency towards double welled potentials using model coupled harmonic diabatic potential-energy surfaces. The von Neumann entropy of the reduced density matrix is used to quantify the electron-vibration entanglement for the lowest two vibronic wavefunctions in such a bipartite system. Sig… ▽ More We consider the quantum entanglement of the electronic and vibrational degrees of freedom in molecules with a tendency towards double welled potentials using model coupled harmonic diabatic potential-energy surfaces. The von Neumann entropy of the reduced density matrix is used to quantify the electron-vibration entanglement for the lowest two vibronic wavefunctions in such a bipartite system. Significant entanglement is found only in the region in which the ground vibronic state contains a density profile that is bimodal (i.e., contains two separate local minima). However, in this region two distinct types of entanglement are found: (1) entanglement that arises purely from the degeneracy of energy levels in the two potential wells and which is destroyed by slight asymmetry, and (2) entanglement that involves strongly interacting states in each well that is relatively insensitive to asymmetry. These two distinct regions are termed fragile degeneracy-induced entanglement and persistent entanglement, respectively. Six classic molecular systems describable by two diabatic states are considered: ammonia, benzene, semibullvalene, pyridine excited triplet states, the Creutz-Taube ion, and the radical cation of the "special pair" of chlorophylls involved in photosynthesis. These chemically diverse systems are all treated using the same general formalism and the nature of the entanglement that they embody is elucidated. △ Less

Submitted 6 July, 2011; originally announced July 2011.

Journal ref: Journal of Chemical Physics 135, 244110 (2011)

Consistent description of the metallic phase of overdoped cuprate superconductors as an anisotropic marginal Fermi liquid

Authors: J. Kokalj, Ross H. McKenzie

Abstract: We consider a model self-energy consisting of an isotropic Fermi liquid term and a Marginal Fermi liquid term which is anisotropic over the Fermi surface, vanishing in the same directions as the superconducting gap and the pseudogap. This model self-energy gives a consistent description of experimental results from Angle-Dependent Magneto-Resistance (ADMR), specific heat, de Haas-van Alphen, and m… ▽ More We consider a model self-energy consisting of an isotropic Fermi liquid term and a Marginal Fermi liquid term which is anisotropic over the Fermi surface, vanishing in the same directions as the superconducting gap and the pseudogap. This model self-energy gives a consistent description of experimental results from Angle-Dependent Magneto-Resistance (ADMR), specific heat, de Haas-van Alphen, and measurements of the quasi-particle dispersion near the Fermi surface from photoemission. In particular, we reconcile the strongly doping dependent anomalous scattering rate observed in ADMR with the almost doping independent specific heat. △ Less

Submitted 18 July, 2011; v1 submitted 11 May, 2011; originally announced May 2011.

Comments: 5+5 pages, 6 figures, revised version

Journal ref: Phys. Rev. Lett. 107, 147001 (2011)

Quantum entanglement between a nonlinear nanomechanical resonator and a microwave field

Authors: Charles P. Meaney, Ross H. McKenzie, G. J. Milburn

Abstract: We consider a theoretical model for a nonlinear nanomechanical resonator coupled to a superconducting microwave resonator. The nanomechanical resonator is driven parametrically at twice its resonance frequency, while the superconducting microwave resonator is driven with two tones that differ in frequency by an amount equal to the parametric driving frequency. We show that the semi-classical appro… ▽ More We consider a theoretical model for a nonlinear nanomechanical resonator coupled to a superconducting microwave resonator. The nanomechanical resonator is driven parametrically at twice its resonance frequency, while the superconducting microwave resonator is driven with two tones that differ in frequency by an amount equal to the parametric driving frequency. We show that the semi-classical approximation of this system has an interesting fixed point bifurcation structure. In the semi-classical dynamics a transition from stable fixed points to limit cycles is observed as one moves from positive to negative detuning. We show that signatures of this bifurcation structure are also present in the full dissipative quantum system and further show that it leads to mixed state entanglement between the nanomechanical resonator and the microwave cavity in the dissipative quantum system that is a maximum close to the semi-classical bifurcation. Quantum signatures of the semi-classical limit-cycles are presented. △ Less

Submitted 29 November, 2010; v1 submitted 19 October, 2010; originally announced October 2010.

Comments: 36 pages, 18 figures

Journal ref: Phys. Rev. E 83, 056202 (2011)

Proposed measurements of the interlayer magnetoresistance of underdoped cuprate superconductors can distinguish closed pockets from open arcs in the Fermi surface

Authors: M. F. Smith, Ross H. McKenzie

Abstract: An outstanding question concerning the underdoped cuprate concerns the true nature of their Fermi surface which appears as a set of disconnected arcs. Theoretical models have proposed two distinct possibilities: (1) each arc is the observable part of a partially-hidden closed pocket, and (2) each arc is open, truncated at its apparent ends. We show that measurements of the variation of the interla… ▽ More An outstanding question concerning the underdoped cuprate concerns the true nature of their Fermi surface which appears as a set of disconnected arcs. Theoretical models have proposed two distinct possibilities: (1) each arc is the observable part of a partially-hidden closed pocket, and (2) each arc is open, truncated at its apparent ends. We show that measurements of the variation of the interlayer resistance with the direction of a magnetic field parallel to the layers can qualitatively distinguish closed pockets from open arcs. This is possible because the field can be oriented such that all electrons on arcs encounter a large Lorentz force and resulting magnetoresistance whereas some electrons on pockets escape the effect by moving parallel to the field. △ Less

Submitted 8 September, 2010; v1 submitted 6 September, 2010; originally announced September 2010.

Comments: 10 pages, 2 figures

Bond Alternation, Polarizability and Resonance Detuning in Methine Dyes

Authors: Seth Olsen, Ross H. McKenzie

Abstract: We derive structure-property relationships for methine ("Brooker") dyes relating the color of the dye and its symmetric parents to its bond alternation in the ground state and also to the dipole properties associated with its low-lying charge-resonance (or charge-transfer) transition. We calibrate and test these relationships on an array of different protonation states of the green fluorescent pro… ▽ More We derive structure-property relationships for methine ("Brooker") dyes relating the color of the dye and its symmetric parents to its bond alternation in the ground state and also to the dipole properties associated with its low-lying charge-resonance (or charge-transfer) transition. We calibrate and test these relationships on an array of different protonation states of the green fluorescent protein chromophore motif (an asymmetric halochromic methine dye) and its symmetric parent dyes. The relationships rely on the assumption that the diabatic states that define the Platt model for methine dye color [J.R. Platt, J. Chem. Phys. 25 80 (1956)] can also be distinguished by their single-double bond alternation and by their charge localization character. These assumptions are independent of the primary constraint that defines the diabatic states in the Platt model - specifically, the Brooker deviation rule for methine dyes [L.G.S. Brooker, Rev. Mod. Phys. 14 275 (1942)]. Taking these assumptions, we show that the Platt model offers an alternate route to known structure-property relationships between the bond length alternation and the quadratic nonlinear polarizability β. We show also that the Platt model can be parameterized without the need for synthesis of the symmetric parents of a given dye, using dipole data obtained through spectroscopic measurements. This suggests that the Platt model parameters may be used as independent variables in free-energy relationships for chromophores whose symmetric parents cannot be synthesized or chromophores strongly bound to biomolecular environments. The latter category includes several recently characterized biomolecular probe constructs. We illustrate these concepts by an analysis of previously reported electroabsorption and second-harmonic generation experiments on green fluorescent proteins. △ Less

Submitted 21 February, 2011; v1 submitted 29 August, 2010; originally announced August 2010.

Comments: Keywords: Structure-property relationships, Brooker dyes, methine dyes, push-pull polyenes, nonlinear optical response, bond length alternation, bond order - bond length relationships, fluorescent proteins

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

Authors: B. J. Powell, Ross H. McKenzie

Abstract: 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… ▽ More 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. △ Less

Submitted 13 January, 2011; v1 submitted 30 July, 2010; originally announced July 2010.

Comments: Major rewrite. New material added and many typos corrected. 67 pages, 41 figures. Thanks to those who commented on the previous version

Journal ref: Rep. Prog. Phys. 74, 056501 (2011)

Models of organometallic complexes for optoelectronic applications

Authors: A. C. Jacko, Ross H. McKenzie, B. J. Powell

Abstract: Organometallic complexes have potential applications as the optically active components of organic light emitting diodes (OLEDs) and organic photovoltaics (OPV). Development of more effective complexes may be aided by understanding their excited state properties. Here we discuss two key theoretical approaches to investigate these complexes: first principles atomistic models and effective Hamiltoni… ▽ More Organometallic complexes have potential applications as the optically active components of organic light emitting diodes (OLEDs) and organic photovoltaics (OPV). Development of more effective complexes may be aided by understanding their excited state properties. Here we discuss two key theoretical approaches to investigate these complexes: first principles atomistic models and effective Hamiltonian models. We review applications of these methods, such as, determining the nature of the emitting state, predicting the fraction of injected charges that form triplet excitations, and explaining the sensitivity of device performance to small changes in the molecular structure of the organometallic complexes. △ Less

Submitted 2 July, 2010; originally announced July 2010.

Comments: To appear in themed issue of J. Mat. Chem. on the modelling of materials

Journal ref: J. Mater. Chem. 20, 10301 (2010)

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

Authors: A. C. Jacko, B. J. Powell, Ross H. McKenzie

Abstract: 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… ▽ More 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'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. △ Less

Submitted 25 July, 2010; v1 submitted 21 May, 2010; originally announced May 2010.

Comments: Accepted for publication in J. Chem. Phys. 14 pages, 9 figures, Supplementary Info: 15 pages, 17 figures

Journal ref: J. Chem. Phys. 133, 124314 (2010)

Thermal and electrical conductivity of Fermi pocket models of underdoped cuprates

Authors: M. F. Smith, Ross H. McKenzie

Abstract: Several models of the electronic spectrum in the pseudogap state of underdoped cuprates have been proposed to explain ARPES and STM measurements, which reveal only truncated Fermi pockets instead of a full metallic Fermi surface. We consider the transport properties expected of four physically distinct models, and calculate the thermal and electrical conductivity of the electronic quasiparticles.… ▽ More Several models of the electronic spectrum in the pseudogap state of underdoped cuprates have been proposed to explain ARPES and STM measurements, which reveal only truncated Fermi pockets instead of a full metallic Fermi surface. We consider the transport properties expected of four physically distinct models, and calculate the thermal and electrical conductivity of the electronic quasiparticles. By proposing transport currents that reflect the close correspondence between quasiparticles on the Fermi pockets in the pseudogap and those near nodes in the superconducting state, we show that measurable transport coefficients provide stringent tests of pseudogap models. △ Less

Submitted 18 May, 2010; v1 submitted 22 April, 2010; originally announced April 2010.

Comments: 5 pages, 3 figures

A Dark Excited State of Fluorescent Protein Chromophores, Considered as Brooker Dyes

Authors: Seth Olsen, Ross H. McKenzie

Abstract: The green fluorescent protein (GFP) chromophore is an asymmetric monomethine dye system. In the resonance color theory of dyes, a strong optical excitation arises from interactions of two valence-bond structures with a third, higher structure. We use correlated quantum chemistry to show that the anionic chromophore is a resonant Brooker dye, and that the third structure corresponds to a higher s… ▽ More The green fluorescent protein (GFP) chromophore is an asymmetric monomethine dye system. In the resonance color theory of dyes, a strong optical excitation arises from interactions of two valence-bond structures with a third, higher structure. We use correlated quantum chemistry to show that the anionic chromophore is a resonant Brooker dye, and that the third structure corresponds to a higher stationary electronic state of this species. The excitation energy of this state should be just below the first excitation energy of the neutral form. This has implications for excited state mechanism in GFPs, which we discuss. △ Less

Submitted 14 January, 2010; originally announced January 2010.

Comments: This paper has been submitted for publication in Chemical Physics Letters

Quantum Transition State Theory for proton transfer reactions in enzymes

Authors: Jacques P. Bothma, Joel Gilmore, Ross H. McKenzie

Abstract: We consider the role of quantum effects in the transfer of hyrogen-like species in enzyme-catalysed reactions. This study is stimulated by claims that the observed magnitude and temperature dependence of kinetic isotope effects imply that quantum tunneling below the energy barrier associated with the transition state significantly enhances the reaction rate in many enzymes. We use a path integra… ▽ More We consider the role of quantum effects in the transfer of hyrogen-like species in enzyme-catalysed reactions. This study is stimulated by claims that the observed magnitude and temperature dependence of kinetic isotope effects imply that quantum tunneling below the energy barrier associated with the transition state significantly enhances the reaction rate in many enzymes. We use a path integral approach which provides a general framework to understand tunneling in a quantum system which interacts with an environment at non-zero temperature. Here the quantum system is the active site of the enzyme and the environment is the surrounding protein and water. Tunneling well below the barrier only occurs for temperatures less than a temperature $T_0$ which is determined by the curvature of potential energy surface near the top of the barrier. We argue that for most enzymes this temperature is less than room temperature. For physically reasonable parameters quantum transition state theory gives a quantitative description of the temperature dependence and magnitude of kinetic isotope effects for two classes of enzymes which have been claimed to exhibit signatures of quantum tunneling. The only quantum effects are those associated with the transition state, both reflection at the barrier top and tunneling just below the barrier. We establish that the friction due to the environment is weak and only slightly modifies the reaction rate. Furthermore, at room temperature and for typical energy barriers environmental degrees of freedom with frequencies much less than 1000 cm$^{-1}$ do not have a significant effect on quantum corrections to the reaction rate. △ Less

Submitted 6 October, 2009; originally announced October 2009.

Comments: Aspects of the article are discussed at condensedconcepts.blogspot.com

Journal ref: New Journal of Physics 12, 055002 (2010)

A Consistent Picture of Charge Distribution in Reduced Ceria Phases

Authors: E. Shoko, M. F. Smith, Ross H. McKenzie

Abstract: We consider the implications of the bond valence model (BVM) description of charge distribution in reduced ceria phases (CeO$_{2-x}$)\cite{Shoko2009,Shoko2009b} to the models used to describe electronic and ionic conductivity in these phases. We conclude that the BVM is consistent with both the small polaron model (SPM) and the atomistic models which describe the electronic and ionic conductivit… ▽ More We consider the implications of the bond valence model (BVM) description of charge distribution in reduced ceria phases (CeO$_{2-x}$)\cite{Shoko2009,Shoko2009b} to the models used to describe electronic and ionic conductivity in these phases. We conclude that the BVM is consistent with both the small polaron model (SPM) and the atomistic models which describe the electronic and ionic conductivities respectively. For intermediate phases, i.e., $x \sim 0.3$, we suggest the possibility of low temperature metallic conductivity. This has not yet been experimentally observed. We contrast the BVM results and the conventional description of charge distribution in reduced ceria phases. △ Less

Submitted 22 December, 2009; v1 submitted 5 October, 2009; originally announced October 2009.

Comments: 5 figures Changes in V2: 1. New title. 2. Added a theory section. 3. Additional references. 4. Additional calculations

Conical Intersections, Charge Transfer and Photoisomerization Pathway Selection in a Minimal Model of a Degenerate Monomethine Dye

Authors: Seth Olsen, Ross H. McKenzie

Abstract: We propose a minimal model Hamiltonian for the electronic structure of a monomethine dye, in order to describe the photoisomerization of such dyes. The model describes interactions between three diabatic electronic states, each of which can be associated with a valence bond structure. Monomethine dyes are characterized by a charge-transfer resonance; the indeterminacy of the single-double bondin… ▽ More We propose a minimal model Hamiltonian for the electronic structure of a monomethine dye, in order to describe the photoisomerization of such dyes. The model describes interactions between three diabatic electronic states, each of which can be associated with a valence bond structure. Monomethine dyes are characterized by a charge-transfer resonance; the indeterminacy of the single-double bonding structure dictated by the resonance is reflected in a duality of photoisomerization pathways corresponding to the different methine bonds. The possible multiplicity of decay channels complicates mechanistic models of the effect of the environment on fluorescent quantum yields, as well as coherent control strategies. We examine the extent and topology of intersection seams between the electronic states of the dye, and how they relate to charge localization and selection between different decay pathways. We find that intersections between the S1 and S0 surfaces only occur for large twist angles. In contrast, S2/S1 intersections can occur near the Franck-Condon region. When the molecule has left-right symmetry, all intersections are associated with con- or dis-rotations and never with single bond twists. For asymmetric molecules (i.e. where the bridge couples more strongly to one end) then the S2 and S1 surfaces bias torsion about different bonds. Charge localization and torsion pathway biasing are correlated. We relate our observations to several recent experimental and theoretical results, which have been obtained for dyes with similar structure. △ Less

Submitted 5 November, 2009; v1 submitted 28 September, 2009; originally announced September 2009.

Comments: This manuscript is a significant rewrite of an earlier manuscript (arXiv:0906.3914). The current version expands significantly on the previous model. It includes several new figures and the text has been completely rewritten. A shortened version of this manuscript has been accepted for publication in The Journal of Chemical Physics

Charge Distribution Near Oxygen Vacancies in Reduced Ceria

Authors: E. Shoko, M. F. Smith, Ross H. McKenzie

Abstract: Understanding the electronic charge distribution around oxygen vacancies in transition metal and rare earth oxides is a scientific challenge of considerable technological importance. We show how significant information about the charge distribution around vacancies in cerium oxide can be gained from a study of high resolution crystal structures of higher order oxides which exhibit ordering of ox… ▽ More Understanding the electronic charge distribution around oxygen vacancies in transition metal and rare earth oxides is a scientific challenge of considerable technological importance. We show how significant information about the charge distribution around vacancies in cerium oxide can be gained from a study of high resolution crystal structures of higher order oxides which exhibit ordering of oxygen vacancies. Specifically, we consider the implications of a bond valence sum analysis of Ce$_{7}$O$_{12}$ and Ce$_{11}$O$_{20}$. To illuminate our analysis we show alternative representations of the crystal structures in terms of orderly arrays of co-ordination defects and in terms of flourite-type modules. We found that in Ce$_{7}$O$_{12}$, the excess charge resulting from removal of an oxygen atom delocalizes among all three triclinic Ce sites closest to the O vacancy. In Ce$_{11}$O$_{20}$, the charge localizes on the next nearest neighbour Ce atoms. Our main result is that the charge prefers to distribute itself so that it is farthest away from the O vacancies. This contradicts \emph{the standard picture of charge localisation} which assumes that each of the two excess electrons localises on one of the cerium ions nearest to the vacancy. This standard picture is assumed in most calculations based on density functional theory (DFT). Based on the known crystal structure of Pr$_{6}$O$_{11}$, we also predict that the charge in Ce$_{6}$O$_{11}$ will be found in the second coordination shell of the O vacancy. Although this review focuses on bulk cerium oxides our approach to characterising electronic properties of oxygen vacancies and the physical insights gained should also be relevant to surface defects and to other rare earth and transition metal oxides. △ Less

Submitted 28 December, 2009; v1 submitted 7 July, 2009; originally announced July 2009.

Comments: 20 pages, 23 figures. The replacement file has a new format for the figures are the document layout but no change in content. v3 has the following main changes: 1. The abstract and introduction were extensively revised. 2. Sec. IV was removed. 3. The Conclusion was rewritten

Journal ref: J. Phys.: Condens. Matter 22 223201 (2010)

Conical Intersections and Charge Localization in a Minimal Model of a Monomethine Dye System: Degenerate Case

Authors: Seth Olsen, Ross H. McKenzie

Abstract: We investigate a minimal valence-bond-based model Hamiltonian for photoisomerization of monomethine dyes. These include a broad class of large conjugated organic molecules including diarylmethane and cyanine dyes, and the chromophores of fluorescent proteins. The model fulfills the most basic expectations for a double bond photoisomerization model: 1) it describes favorable twisting in the excit… ▽ More We investigate a minimal valence-bond-based model Hamiltonian for photoisomerization of monomethine dyes. These include a broad class of large conjugated organic molecules including diarylmethane and cyanine dyes, and the chromophores of fluorescent proteins. The model fulfills the most basic expectations for a double bond photoisomerization model: 1) it describes favorable twisting in the excited state and 2) it describes conical intersections which can mediate radidationless decay. We examine the topology, topography and extent of the conical intersection seams. Important properties of the model include 1) two regimes of the coupling ratio where the first excited state intersects the ground and a higher excited state along orthogonal twisting coordinates or along the same coordinate 2) first and higher excited state intersections occurring in regions of little or no twisting and 3) charge localization connected with excited-state twisting which changes discontinuously across a conical intersection. The model is promising as a simple extension of empirical force fields in order to describe photoisomerization processes in methine-bridged systems and as a starting point for more elaborate models which may include antisymmetric internal couplings or which may be interfaced with environmental or quantum dynamical models. △ Less

Submitted 21 June, 2009; originally announced June 2009.

Interlayer magnetoresistance in an anisotropic pseudogap state

Authors: M. F. Smith, Ross H. McKenzie

Abstract: The interlayer magnetoresistance of a quasi-two-dimensional layered metal with a d-wave pseudogap is calculated semiclassically. An expression for the interlayer resistivity as a function of the strength and direction of the magnetic field, the magnitude of the pseudogap, temperature, and scattering rate is obtained. We find that the pseudogap, by introducing low-energy nodal quasiparticle conto… ▽ More The interlayer magnetoresistance of a quasi-two-dimensional layered metal with a d-wave pseudogap is calculated semiclassically. An expression for the interlayer resistivity as a function of the strength and direction of the magnetic field, the magnitude of the pseudogap, temperature, and scattering rate is obtained. We find that the pseudogap, by introducing low-energy nodal quasiparticle contours, smooths the dependence on field direction in a manner characteristic of its anisotropy. We thus propose that interlayer resistance measurements under a strong field of variable orientation can be used to fully characterize an anisotropic pseudogap. The general result is applied to the case of a magnetic field parallel to the conducting layers using a model band structure appropriate for overdoped T$\ell$2201. △ Less

Submitted 4 June, 2009; originally announced June 2009.

Comments: 18 pages, 3 figures

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

Authors: J. Merino, R. H. McKenzie, B. J. Powell

Abstract: 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 fo… ▽ More 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 "ring" 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. △ Less

Submitted 22 April, 2009; originally announced April 2009.

Comments: 15 pages, 15 figures

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

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

Authors: B. J. Powell, J. Merino, Ross H. McKenzie

Abstract: 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 fille… ▽ More 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. △ Less

Submitted 13 August, 2009; v1 submitted 3 April, 2009; originally announced April 2009.

Comments: To appear in PRB. 19 pages, 19 figures, 6 movies. References added

Spin Fluctuations and the Pseudogap in Organic Superconductors

Authors: B. J. Powell, Eddy Yusuf, Ross H. McKenzie

Abstract: 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… ▽ More 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. △ Less

Submitted 9 July, 2009; v1 submitted 16 March, 2009; originally announced March 2009.

Comments: 7 pages, 4 figures; to appear in PRB

The Jahn-Teller instability in dissipative quantum electromechanical systems

Authors: Charles P. Meaney, Tim Duty, Ross H. McKenzie, Gerard J. Milburn

Abstract: We consider the steady states of a harmonic oscillator coupled so strongly to a two-level system (a qubit) that the rotating wave approximation cannot be made. The Hamiltonian version of this model is known as the $E\otimesβ$ Jahn-Teller model. The semiclassical version of this system exhibits a fixed point bifurcation, which in the quantum model leads to a ground state with substantial entangle… ▽ More We consider the steady states of a harmonic oscillator coupled so strongly to a two-level system (a qubit) that the rotating wave approximation cannot be made. The Hamiltonian version of this model is known as the $E\otimesβ$ Jahn-Teller model. The semiclassical version of this system exhibits a fixed point bifurcation, which in the quantum model leads to a ground state with substantial entanglement between the oscillator and the qubit. We show that the dynamical bifurcation survives in a dissipative quantum description of the system, amidst an even richer bifurcation structure. We propose two experimental implementations of this model based on superconducting cavities: a parametrically driven nonlinear nanomechanical resonator coupled capacitively to a coplanar microwave cavity and a superconducting junction in the central conductor of a coplanar waveguide. △ Less

Submitted 15 March, 2009; originally announced March 2009.

Comments: 24 pages, 13 figures

Journal ref: Phys. Rev. A 81, 043805 (2010)

A Diabatic Three-State Representation of Photoisomerization in the Green Fluorescent Protein Chromophore

Authors: Seth Olsen, Ross H. McKenzie

Abstract: We give a quantum chemical description of bridge photoisomerization reaction of green fluorescent protein (GFP) chromophores using a representation over three diabatic states. Bridge photoisomerization leads to non-radiative decay, and competes with fluorescence in these systems. In the protein, this pathway is suppressed, leading to fluorescence. Understanding the electronic structure of the ph… ▽ More We give a quantum chemical description of bridge photoisomerization reaction of green fluorescent protein (GFP) chromophores using a representation over three diabatic states. Bridge photoisomerization leads to non-radiative decay, and competes with fluorescence in these systems. In the protein, this pathway is suppressed, leading to fluorescence. Understanding the electronic structure of the photoisomerization is a prerequisite to understanding how the protein suppresses this pathway and preserves the emitting state of the chromophore. We present a solution to the state-averaged complete active space problem, which is spanned at convergence by three fragment-localized orbitals. We generate the diabatic-state representation by applying a block diagonalization transformation to the Hamiltonian calculated for the anionic chromophore model HBDI with multi-reference, multi-state perturbation theory. The diabatic states that emerge are charge-localized structures with a natural valence-bond interpretation. At planar geometries, the diabatic picture recaptures the charge transfer resonance of the anion. The strong S0-S1 excitation at these geometries is reasonably described within a two-state model, but extension to a three-state model is necessary to describe decay via two possible pathways associated with photoisomerization of the (methine) bridge. Parametric Hamiltonians based on the three-state ansatz can be fit directly to data generated using the underlying active space. We provide an illustrative example of such a parametric Hamiltonian. △ Less

Submitted 14 May, 2009; v1 submitted 22 February, 2009; originally announced February 2009.

Journal ref: S. Olsen and R.H. McKenzie, The Journal of Chemical Physics, 130, 184302 (2009)

Mixed valency in cerium oxide crystallographic phases: Determination of valence of the different cerium sites by the bond valence method

Authors: E. Shoko, M. F. Smith, Ross H. McKenzie

Abstract: We have applied the bond valence method to cerium oxides to determine the oxidation states of the Ce ion at the various site symmetries of the crystals. The crystals studied include cerium dioxide and the two sesquioxides along with some selected intermediate phases which are crystallographically well characterized. Our results indicate that cerium dioxide has a mixed-valence ground state with a… ▽ More We have applied the bond valence method to cerium oxides to determine the oxidation states of the Ce ion at the various site symmetries of the crystals. The crystals studied include cerium dioxide and the two sesquioxides along with some selected intermediate phases which are crystallographically well characterized. Our results indicate that cerium dioxide has a mixed-valence ground state with an f-electron population on the Ce site of 0.27 while both the A- and C-sesquioxides have a nearly pure f^1 configuration. The Ce sites in most of the intermediate oxides have non-integral valences. Furthermore, many of these valences are different from the values predicted from a naive consideration of the stoichiometric valence of the compound. △ Less

Submitted 27 November, 2008; originally announced November 2008.

Journal ref: Phys. Rev. B 79 134108 (2009)