Chemical Physics (physics.chem-ph)

  • PDF
    A central difficulty of state-specific Multi-Reference Coupled Cluster (MR-CC) formalisms concerns the definition of the amplitudes of the single and double excitation operators appearing in the exponential wave operator. If the reference space is a complete active space (CAS) the number of these amplitudes is larger than the number of singly and doubly excited determinants on which one may project the eigenequation, and one must impose additional conditions. The present work first defines a state-specific reference-independent operator $\hat{\tilde{T}}^m$ which acting on the CAS component of the wave function $|\Psi_0^m \rangle$ maximizes the overlap between $(1+\hat{\tilde{T}}^m)|\Psi_0^m \rangle$ and the eigenvector of the CAS-SD CI matrix $|\Psi_{\rm CAS-SD}^m \rangle$. This operator may be used to generate approximate coefficients of the Triples and Quadruples, and a dressing of the CAS-SD CI matrix, according to the intermediate Hamiltonian formalism. The process may be iterated to convergence. As a refinement towards a strict Coupled Cluster formalism, one may exploit reference-independent amplitudes provided by $(1+\hat{\tilde{T}}^m)|\Psi_0^m \rangle$ to define a reference-dependent operator $\hat{T}^m$ by fitting the eigenvector of the (dressed) CAS-SD CI matrix. The two variants, which are internally uncontracted, give rather similar results. The new MR-CC version has been tested on the ground state potential energy curves of 6 molecules (up to triple-bond breaking) and a two excited states. The non-parallelism error with respect to the Full-CI curves is of the order of 1 m$E_{\rm h}$.
  • PDF
    Metal atoms and small clusters introduced into superfluid helium (He II) concentrate there in quantized vortices to form (by further coagulation) the thin nanowires. The nanowires' thickness and structure are well predicted by a double-staged mechanism. On the first stage the coagulation of cold particles in the vortex cores leads to melting of their fusion product, which acquires a spherical shape due to surface tension. Then (second stage) when these particles reach a certain size they do not possess sufficient energy to melt and eventually coalesce into the nanowires. Nevertheless the assumption of melting for such refractory metal as tungsten, especially in He II, which possesses an extremely high thermal conductivity, induces natural skepticism. That is why we decided to register directly the visible thermal emission accompanying metals coagulation in He II. The brightness temperatures of this radiation for the tungsten, molybdenum, and platinum coagulation were found to be noticeably higher than even the metals' melting temperatures. The region of He II that contained suspended metal particles expanded with the velocity of 50 m/s, being close to the Landau velocity, but coagulation took place even more quickly, so that the whole process of nanowire growth is completed at distances about 1.5 mm from the place of metal injection into He II. High rate of coagulation of guest metal particles as well as huge local overheating are associated with them concentrating in quantized vortex cores. The same process should take place not only for metals but for any atoms, molecules and small clusters embedded into He II.
  • PDF
    Wavefunctions restricted to electron-pair states are promising models to describe static/nondynamic electron correlation effects encountered, for instance, in bond-dissociation processes and transition-metal and actinide chemistry. To reach spectroscopic accuracy, however, the missing dynamic electron correlation effects that cannot be described by electron-pair states need to be included \textita posteriori. In this article, we extend the previously presented perturbation theory models with an Antisymmetric Product of 1-reference orbital Geminal (AP1roG) reference function that allow us to describe both static/nondynamic and dynamic electron correlation effects. Specifically, our perturbation theory models combine a diagonal and off-diagonal zero-order Hamiltonian, a single-reference and multi-reference dual state, and different excitation operators used to construct the projection manifold. We benchmark all proposed models as well as an \textita posteriori linearized coupled cluster correction on top of AP1roG against CR-CCSD(T) reference data for reaction energies of several closed-shell molecules that are extrapolated to the basis set limit. Moreover, we test the performance of our new methods for multiple bond breaking processes in the N$_2$, C$_2$, and BN dimers against MRCI-SD and MRCI-SD+Q reference data. Our numerical results indicate that the best performance is obtained from a linearized coupled cluster correction as well as second-order perturbation theory corrections employing a diagonal and off-diagonal zero-order Hamiltonian and a single-determinant dual state. These dynamic corrections on top of AP1roG allow us to reliably model molecular systems dominated by static/nondynamic as well as dynamic electron correlation.
  • PDF
    First-principles calculations combining density-functional theory and continuum solvation models enable realistic theoretical modeling and design of electrochemical systems. When a reaction proceeds in such systems, the number of electrons in the portion of the system treated quantum mechanically changes continuously, with a balancing charge appearing in the continuum electrolyte. A grand-canonical ensemble of electrons at a chemical potential set by the electrode potential is therefore the ideal description of such systems that directly mimics the experimental condition. We present two distinct algorithms, a self-consistent field method (GC-SCF) and a direct variational free energy minimization method using auxiliary Hamiltonians (GC-AuxH), to solve the Kohn-Sham equations of electronic density-functional theory directly in the grand canonical ensemble at fixed potential. Both methods substantially improve performance compared to a sequence of conventional fixed-number calculations targeting the desired potential, with the GC-AuxH method additionally exhibiting reliable and smooth exponential convergence of the grand free energy. Finally, we apply grand-canonical DFT to the under-potential deposition of copper on platinum from chloride-containing electrolytes and show that chloride desorption, not partial copper monolayer formation, is responsible for the second voltammetric peak.
  • PDF
    Coherent control of chaotic molecular systems, using laser-assisted alignment of sulphur dioxide (SO$_2$) molecules in the presence of a static electric field as an example, is considered. Conditions for which the classical version of this system is chaotic are established, and the quantum and classical analogs are shown to be in very good correspondence. It is found that the chaos present in the classical system does not impede the alignment, neither in the classical nor the quantum system. Using the results of numerical calculations, we suggest that laser-assisted alignment is stable against rotational chaos for all asymmetric top molecules.

Recent comments

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]