# Chemical Physics (physics.chem-ph)

• The reliability of popular density functionals was studied for the description of torsional profiles of 36 molecules: glyoxal, oxalyl halides and their thiocarbonyl derivatives. Sixteen functionals of varying complexity, from local density to range-separated hybrid approximations, have been considered and benchmarked against CCSD(T)-level rotational profiles. For molecules containing heavy halogens, all functionals except M05-2X and M06-2X fail to reproduce barrier heights accurately and a number of functionals introduce spurious minima. Dispersion corrections show no improvement. Calibrated torsion-corrected atom-centered potentials rectify the shortcomings of PBE and also improve on $\sigma$-hole based intermolecular binding in dimers and crystals.
• The analytic energy gradients in the atomic orbital representation have recently been published (J. Chem. Phys. 146, 014102, 2017) within the framework of the natural orbital functional theory (NOFT). We provide here an alternative expression for them in terms of natural orbitals, and use it to derive the analytic second-order energy derivatives with respect to nuclear displacements in the NOFT. The computational burden is shifted to the calculation of perturbed natural orbitals and occupancies, since a set of linear coupled-perturbed equations obtained from the variational Euler equations must be solved to attain the analytic Hessian at the perturbed geometry. The linear response of both natural orbitals and occupation numbers to nuclear geometry displacements need only specify the reconstruction of the second-order reduced density matrix in terms of occupation numbers.
• Experiments on single crystal Li$_X$FePO$_4$ (LFP) nanoparticles indicate rich nonequilibrium phase behavior, such as suppression of phase separation at high lithiation rates, striped patterns of coherent phase boundaries, nucleation by binarysolid surface wetting and intercalation waves. These observations have been successfully predicted (prior to the experiments) by 1D depth-averaged phase-field models, which neglect any subsurface phase separation. In this paper, using an electro-chemo-mechanical phase-field model, we investigate the coherent non-equilibrium subsurface phase morphologies that develop in the $ab$- plane of platelet-like single-crystal platelet-like Li$_X$FePO$_4$ nanoparticles. Finite element simulations are performed for 2D plane-stress conditions in the $ab$- plane, and validated by 3D simulations, showing similar results. We show that the anisotropy of the interfacial tension tensor, coupled with electroautocatalytic surface intercalation reactions, plays a crucial role in determining the subsurface phase morphology. With isotropic interfacial tension, subsurface phase separation is observed, independent of the reaction kinetics, but for strong anisotropy, phase separation is controlled by surface reactions, as assumed in 1D models. Moreover, the driven intercalation reaction suppresses phase separation during lithiation, while enhancing it during delithiation, by electro-autocatalysis, in quantitative agreement with \it in operando imaging experiments in single-crystalline nanoparticles, given measured reaction rate constants.
• The paper has two goals: It presents basic ideas, notions, and methods for reduction of reaction kinetics models: quasi-steady-state, quasi-equilibrium, slow invariant manifolds, and limiting steps. It describes briefly the current state of the art and some latest achievements in the broad area of model reduction in chemical and biochemical kinetics, including new results in methods of invariant manifolds, computation singular perturbation, bottleneck methods, asymptotology, tropical equilibration, and reaction mechanism skeletonisation.
• We give an explicit and general description of the energy, linear momentum, angular momentum and boost momentum of a molecule to order $1/c^2$, where it necessary to take account of kinetic contributions made by the electrons and nuclei as well as electromagnetic contributions made by the intramolecular field. A wealth of interesting subtleties are encountered that are not seen at order $1/c^0$, including relativistic Hall shifts, anomalous velocities and hidden momenta. Some of these have well known analogues in solid state physics.

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

...(continued)
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

...(continued)
Omar Shehab Mar 22 2014 04:06 UTC

Jarrod,

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

...(continued)
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

...(continued)
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]