Chemical Physics (physics.chem-ph)

• In the framework of the Gibbs statistical theory, the question of the size of the particles forming the statistical system is investigated. This task is relevant for a wide variety of applications. The distribution for particle sizes and the moments of this quantity are determined from probabilistic considerations. The results are compared with different models and approximations. The particle size depends on the interactions in the system, on the compressibility factor, on the number of interacting particles. The cases of the absence and presence of a rigid incompressible particle core are considered.
• We introduce a simple ansatz for the wavefunction of a many-body system based on coupled forward and backward-propagating semiclassical trajectories. This method is primarily aimed at, but not limited to, treating nonequilibrium dynamics in electron-phonon systems. The time-evolution of the system is obtained from the Euler-Lagrange variational principle, and we show that this ansatz yields Ehrenfest mean field theory in the limit that the forward and backward trajectories are orthogonal, and in the limit that they coalesce. We investigate accuracy and performance of this method by simulating electronic relaxation in the spin-boson model and the Holstein model. Although this method involves only pairs of semiclassical trajectories, it shows a substantial improvement over mean field theory, capturing quantum coherence of nuclear dynamics as well as electron-nuclear correlations. This improvement is particularly evident in nonadiabatic systems, where the accuracy of this coupled trajectory method extends well beyond the perturbative electron-phonon coupling regime. This approach thus provides an attractive route forward to the ab-initio description of relaxation processes, such as thermalization, in condensed phase systems.
• Automatic Differentiation (AD) is a powerful tool that allows calculating derivatives of implemented algorithms with respect to all of their parameters up to machine precision, without the need to explicitly add any additional functions. Thus, AD has great potential in quantum chemistry, where gradients are omnipresent but also difficult to obtain, and researchers typically spend a considerable amount of time finding suitable analytical forms when implementing derivatives. Here, we demonstrate that automatic differentiation can be used to compute gradients with respect to any parameter throughout a complete quantum chemistry method. We implement DiffiQult, a fully autodifferentiable Hartree-Fock (HF) algorithm, which serves as a proof-of-concept that illustrates the capabilities of AD for quantum chemistry. We leverage the obtained gradients to optimize the parameters of one-particle basis sets in the context of the floating Gaussian framework.
• Ordinary differential equation models have become a standard tool for the mechanistic description of biochemical processes. If parameters are inferred from experimental data, such mechanistic models can provide accurate predictions about the behavior of latent variables or the process under new experimental conditions. Complementarily, inference of model structure can be used to identify the most plausible model structure from a set of candidates, and thus gain novel biological insight. Several toolboxes can infer model parameters and structure for small- to medium-scale mechanistic models out of the box. However, models for highly multiplexed datasets can require hundreds to thousands of state variables and parameters. For the analysis of such large-scale models, most algorithms require intractably high computation times. This chapter provides an overview of state-of-the-art methods for parameter and model inference, with an emphasis on scalability.
• The magnetic tuning of the low rotational levels in the Xsigma+(0,0,0), Api(0,0,0), and Bsigma+ (0,0,0) electronic states of strontium hydroxide, SrOH, have been experimentally investigated using high resolution optical Zeeman spectroscopy of a cold molecular beam sample. The observed Zeeman shifts and splittings are successfully modeled using a traditional effective Hamiltonian approach to account for the interaction between the Api and Bsigma+ states. The determined magnetic g-factors for the Xsigma+, Api , and Bsigma+ states are compared to those predicted by perturbation theory. The dispersed fluorescence resulting from laser excitation of rotationally resolved branch features of the Bsigma+ (0,0,0)<Xsigma+(0,0,0) Api(0,0,0)< Xsigma+(0,0,0) transitions have been recorded and analyzed. The measured fluorescence branching ratios are compared with Franck-Condon calculations. The required bending motion wave functions are derived using a discrete variable representation (DVR) method. Implications for laser slowing and magneto-optical trapping experiments for SrOH are described.

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]