Chemical Physics (physics.chem-ph)

  • PDF
    Spectral density functions are central to the simulation of complex many body systems. Their determination requires to make approximations not only to the dynamics but also to the underlying electronic structure theory. Here, blending different methods bears the danger of an inconsistent description. To solve this issue we propose an all-DFTB approach to determine spectral densities for the description of Frenkel excitons in molecular assemblies. The protocol is illustrated for a model of a PTCDI crystal, which involves the calculation of monomeric excitation energies, Coulomb couplings between monomer transitions, as well as their spectral distributions due to thermal fluctuations of the nuclei. Using dynamically defined normal modes, a mapping onto the standard harmonic oscillator spectral densities is achieved.
  • PDF
    FeMoco, MoFe$_7$S$_9$C, has been shown to be the active catalytic site for the reduction of nitrogen to ammonia in the nitrogenase protein. An understanding of its electronic structure including strong electron correlation is key to designing mimic catalysts capable of ambient nitrogen fixation. Active spaces ranging from [54, 54] to [65, 57] have been predicted for a quantitative description of FeMoco's electronic structure. However, a wavefunction approach for a singlet state using a [54,54] active space would require 10$^{29}$ variables. In this work, we systematically explore the active-space size necessary to qualitatively capture strong correlation in FeMoco and two related moieties, MoFe$_3$S$_7$ and Fe$_4$S$_7$. Using CASSCF and 2-RDM methods, we consider active-space sizes up to [14,14] and [30,30], respectively, with STO-3G, 3-21G, and DZP basis sets and use fractional natural-orbital occupation numbers to assess the level of multireference electron correlation, an examination of which reveals a competition between single-reference and multi-reference solutions to the electronic Schrödinger equation for smaller active spaces and a consistent multi-reference solution for larger active spaces.
  • PDF
    An alternative formulation of the fundamental gap problem in density-functional theory (DFT) is presented. In the new approach, referred to as N-centered ensemble DFT, the infamous derivative discontinuity correction to the gap appears automatically through the weight dependence of the ensemble exchange-correlation energy, like in the Gross-Oliveira-Kohn DFT formulation of the optical gap problem. A two-weight extension is proposed for a separate extraction of the ionization potential and electron affinity. The exact theory is then applied to the asymmetric Hubbard dimer.
  • PDF
    The structure of entangled polymers is perturbed by strong shear, resulting in a non-linear response. However, practical limitations make it difficult to probe the molecular structure far from equilibrium, and its understanding is currently incomplete. Here we report \emphin situ measurements of the form factor of a semi-dilute polymer solution under a steady shear flow, as revealed by small angle neutron scattering. A simple analytical formula is derived with only four parameters, giving an excellent fit to the anisotropic 2D scattering pattern in the flow-vorticity plane. There emerges a length scale below which the polymer is not perturbed by shear, while above it there is elongation along the flow and contraction perpendicular to it. Similar patterns have earlier been reported for polymer melts, which now amounts to two greatly different dynamical systems producing the same time-averaged form factor.
  • PDF
    The quadratic Jahn-Teller effect of C$_{60}^{n-}$ ($n=$ 1-5) is investigated from the first principles. Employing the density functional theory calculations with hybrid functional, the quadratic vibronic coupling constants of C$_{60}^-$ were derived. The warping of the adiabatic potential energy surface of C$_{60}^-$ by the quadratic vibronic coupling is estimated about 2 meV, which is much smaller than the Jahn-Teller stabilization energy ($\approx$ 50 meV). Because of the selection rule and the vibronic reduction, the quadratic coupling slightly modifies the vibronic states of C$_{60}$ anions. Particularly, in the case of C$_{60}^{3-}$, parity and symmetry selection rule significantly reduces the effect of quadratic coupling on vibronic states. The present results confirm that the low-energy vibronic dynamics of C$_{60}^{n-}$ is of pseudorotational type.

Recent comments

Dan Elton Mar 16 2018 04:36 UTC

Comments are appreciated. Message me here or on twitter @moreisdifferent

Code is open source and available at :


Jarrod McClean Apr 02 2014 18:33 UTC

Ryan is exactly correct. The method would work with any of the clock constructions, however we decided that the machinery they developed to make sure it was implementable in qubits was unnecessary overhead for a classical implementation, where having a qudit with large d does not present a great ch

Ryan Babbush Mar 23 2014 19:18 UTC

I'll let Jarrod correct me if he disagrees but I think the point Jarrod is making is that the original clocks of Feynman and Kitaev were designed with physical implementation in mind whereas this proposal is for a classical algorithm and as such, should not be measured by the same considerations (su

Omar Shehab Mar 22 2014 04:06 UTC


Would there be any issue if you had used the original Feynman's clock or Kitaev's clock for simulating quantum systems classically?

Omar Shehab Mar 22 2014 00:28 UTC

Thank you for your quick reply. In that case, as long as we are using qubits for your scheme we should be good. Having said that I am afraid it may jeopardize the local structure recommended by Feynman to some extent.

Jarrod McClean Mar 21 2014 14:16 UTC

Thank you for your interest in our paper! The reason for the apparent discrepancy is that our paper is designed for classical simulation of quantum systems using ideas (namely Feynman's Clock) from quantum computation. When performing the simulation on a classical computer, the underlying qubit st

Omar Shehab Mar 21 2014 00:14 UTC

In equation 7 of the paper, the clock register is different from Feynman's original proposal. According to this paper, for a four clock steps quantum circuit, the sequence of the state of the clock register will be $|00\rangle \to |01\rangle \to |10\rangle \to |11\rangle$. It infers that we will nee

Sergey Filippov Aug 02 2013 12:17 UTC

Physicists may be also interested in a realistic proposal of DFS computation with solid-state qubits: Phys. Lett. A 374, 3285 (2010), arXiv:0903.1056 [quant-ph]