Showing 1–50 of 55 results for author: Whaley, K B

Performance of wave function and Green's functions based methods for non equilibrium many-body dynamics

Authors: Cian C. Reeves, Gaurav Harsha, Avijit Shee, Yuanran Zhu, Chao Yang, K Birgitta Whaley, Dominika Zgid, Vojtech Vlcek

Abstract: Theoretical descriptions of non equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cas… ▽ More Theoretical descriptions of non equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cases of each of the methods and address how these complementary approaches capture the dynamics driven by intense field perturbations to non equilibrium states. Under strong driving, the systems are characterized by strong entanglement of the single particle density matrix and natural populations approaching those of a strongly interacting equilibrium system. We generate a representative set of results that are numerically exact and form a basis for critical comparison of the distinct families of methods. We demonstrate that the compressed formulation based on similarity transformed Hamiltonians (coupled cluster approach) is practically exact in weak fields and, hence, weakly or moderately correlated systems. Coupled cluster, however, struggles for strong driving fields, under which the system exhibits strongly correlated behavior, as measured by the von Neumann entropy of the single particle density matrix. The dynamics predicted by Green's functions in the (widely popular) GW approximation are less accurate by improve significantly upon the mean-field results in the strongly driven regime. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 11 pages, 2 figures in main text. 11 page supplemental information with 6 figures

Comparing self-consistent GW and vertex corrected G0W0 (G0W0Γがんま) accuracy for molecular ionization potentials

Authors: Ming Wen, Vibin Abraham, Gaurav Harsha, Avijit Shee, K. Birgitta Whaley, Dominika Zgid

Abstract: We test the performance of self-consistent GW and several representative implementations of vertex corrected G0W0 (G0W0Γがんま). These approaches are tested on benchmark data sets covering full valence spectra (first ionization potentials and some inner valence shell excitations). For small molecules, when comparing against state of the art wave function techniques, our results show that performing full… ▽ More We test the performance of self-consistent GW and several representative implementations of vertex corrected G0W0 (G0W0Γがんま). These approaches are tested on benchmark data sets covering full valence spectra (first ionization potentials and some inner valence shell excitations). For small molecules, when comparing against state of the art wave function techniques, our results show that performing full self-consistency in the GW scheme either systematically outperforms vertex corrected G0W0 or gives results of at least the same quality. Moreover, the G0W0Γがんま results in additional computational cost when compared to G0W0 or self-consistent GW and the G0W0Γがんま dependency on the starting mean-filed solution is frequently larger than the magnitude of the vertex correction. Consequently, for molecular systems self-consistent GW performed on imaginary axis and then followed by modern analytical continuation techniques offers a more reliable approach to make predictions of IP spectra. △ Less

Submitted 5 April, 2024; v1 submitted 20 November, 2023; originally announced November 2023.

Enhancement of vibrationally assisted energy transfer by proximity to exceptional points, probed by fluorescence-detected vibrational spectroscopy

Authors: Zeng-Zhao Li, K. Birgitta Whaley

Abstract: Emulation of energy transfer processes in natural systems on quantum platforms can further our understanding of complex dynamics in nature. One notable example is the demonstration of vibrationally assisted energy transfer (VAET) on a trapped-ion quantum emulator, which offers insights for the energetics of light harvesting. In this work, we expand the study of VAET simulation with trapped ions to… ▽ More Emulation of energy transfer processes in natural systems on quantum platforms can further our understanding of complex dynamics in nature. One notable example is the demonstration of vibrationally assisted energy transfer (VAET) on a trapped-ion quantum emulator, which offers insights for the energetics of light harvesting. In this work, we expand the study of VAET simulation with trapped ions to a non-Hermitian quantum system comprising a $\mathscr{PT}$-symmetric chromophore dimer weakly coupled to a vibrational mode. We first characterize exceptional points (EPs) and non- Hermitian features of the excitation energy transfer processes in the absence of the vibration, finding a degenerate pair of second-order EPs. Exploring the non-Hermitian dynamics of the whole system including vibrations, we find that energy transfer accompanied by absorption of phonons from a vibrational mode can be significantly enhanced near such a degenerate EP. Our calculations reveal a unique spectral feature accompanying the coalescing of eigenstates and eigenenergies that provides a novel approach to probe the degenerate EP by fluorescence-detected vibrational spectroscopy. Enhancement of the VAET process near the EP is found to be due to maximal favorability of phonon absorption at the degenerate EP, enabling multiple simultaneous excitations. Our work on improving VAET processes in non-Hermitian quantum systems paves the way for leveraging non-Hermiticity in quantum dynamics related to excitation energy transfer. △ Less

Submitted 23 April, 2024; v1 submitted 6 September, 2023; originally announced September 2023.

Comments: 17 pages, 14 figures, 1 table

Journal ref: Phys. Rev. Research 6, 023149 (2024)

Error Mitigated Metasurface-Based Randomized Measurement Schemes

Authors: Hang Ren, Yipei Zhang, Ze Zheng, Cuifeng Ying, Lei Xu, Mohsen Rahmani, K. Birgitta Whaley

Abstract: Estimating properties of quantum states via randomized measurements has become a significant part of quantum information science. In this paper, we design an innovative approach leveraging metasurfaces to perform randomized measurements on photonic qubits, together with error mitigation techniques that suppress realistic metasurface measurement noise. Through fidelity and purity estimation, we con… ▽ More Estimating properties of quantum states via randomized measurements has become a significant part of quantum information science. In this paper, we design an innovative approach leveraging metasurfaces to perform randomized measurements on photonic qubits, together with error mitigation techniques that suppress realistic metasurface measurement noise. Through fidelity and purity estimation, we confirm the capability of metasurfaces to implement randomized measurements and the unbiased nature of our error-mitigated estimator. Our findings show the potential of metasurface-based randomized measurement schemes in achieving robust and resource-efficient estimation of quantum state properties. △ Less

Submitted 19 April, 2024; v1 submitted 16 August, 2023; originally announced August 2023.

Performing quantum entangled biphoton spectroscopy using classical light pulses

Authors: Liwen Ko, Robert L. Cook, K. Birgitta Whaley

Abstract: We show that for a class of quantum light spectroscopy (QLS) experiments using n = 0,1,2,$\cdots$ classical light pulses and an entangled photon pair (a biphoton state) where one photon acts as a reference without interacting with the matter sample, identical signals can be obtained by replacing the biphotons with classical-like coherent states of light, where these are defined explicitly in terms… ▽ More We show that for a class of quantum light spectroscopy (QLS) experiments using n = 0,1,2,$\cdots$ classical light pulses and an entangled photon pair (a biphoton state) where one photon acts as a reference without interacting with the matter sample, identical signals can be obtained by replacing the biphotons with classical-like coherent states of light, where these are defined explicitly in terms of the parameters of the biphoton states. An input-output formulation of quantum nonlinear spectroscopy is used to prove this equivalence. We demonstrate the equivalence numerically by comparing a classical pump - quantum probe experiment with the corresponding classical pump - classical probe experiment. This analysis shows that understanding the equivalence between entangled biphoton probes and carefully designed classical-like coherent state probes leads to quantum-inspired classical experiments and provides insights for future design of QLS experiments that could provide a true quantum advantage. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Say NO to Optimization: A Non-Orthogonal Quantum Eigensolver

Authors: Unpil Baek, Diptarka Hait, James Shee, Oskar Leimkuhler, William J. Huggins, Torin F. Stetina, Martin Head-Gordon, K. Birgitta Whaley

Abstract: A balanced description of both static and dynamic correlations in electronic systems with nearly degenerate low-lying states presents a challenge for multi-configurational methods on classical computers. We present here a quantum algorithm utilizing the action of correlating cluster operators to provide high-quality wavefunction ansätze employing a non-orthogonal multireference basis that captures… ▽ More A balanced description of both static and dynamic correlations in electronic systems with nearly degenerate low-lying states presents a challenge for multi-configurational methods on classical computers. We present here a quantum algorithm utilizing the action of correlating cluster operators to provide high-quality wavefunction ansätze employing a non-orthogonal multireference basis that captures a significant portion of the exact wavefunction in a highly compact manner, and that allows computation of the resulting energies and wavefunctions at polynomial cost with a quantum computer. This enables a significant improvement over the corresponding classical non-orthogonal solver, which incurs an exponential cost when evaluating off-diagonal matrix elements between the ansatz states, and is therefore intractable. We implement the non-orthogonal quantum eigensolver (NOQE) here with an efficient ansatz parameterization inspired by classical quantum chemistry methods that succeed in capturing significant amounts of electronic correlation accurately. By taking advantage of classical methods for chemistry, NOQE provides a flexible, compact, and rigorous description of both static and dynamic electronic correlation, making it an attractive method for the calculation of electronic states of a wide range of molecular systems. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: 27+11 pages, 14+4 figures

Dynamics of photosynthetic light harvesting systems interacting with N-photon Fock states

Authors: Liwen Ko, Robert L. Cook, K. Birgitta Whaley

Abstract: We develop a method to simulate the excitonic dynamics of realistic photosynthetic light harvesting systems including non-Markovian coupling to phonon degrees of freedom, under excitation by N-photon Fock state pulses. This method combines the input-output formalism and the hierarchical equations of motion (HEOM) formalism into a double hierarchy of coupled linear equations in density matrices. We… ▽ More We develop a method to simulate the excitonic dynamics of realistic photosynthetic light harvesting systems including non-Markovian coupling to phonon degrees of freedom, under excitation by N-photon Fock state pulses. This method combines the input-output formalism and the hierarchical equations of motion (HEOM) formalism into a double hierarchy of coupled linear equations in density matrices. We show analytically that, under weak field excitation relevant to natural photosynthesis conditions, an N-photon Fock state input and a corresponding coherent state input give rise to equal density matrices in the excited manifold. However, an important difference is that an N-photon Fock state input has no off-diagonal coherence between the ground and excited subspaces, in contrast with the coherences created by a coherent state input. We derive expressions for the probability to absorb a single Fock state photon, with or without the influence of phonons. For short pulses (or equivalently, wide bandwidth pulses), we show that the absorption probability has a universal behavior that depends only upon a system-dependent effective energy spread parameter Δでるた and an exciton-light coupling constant Γがんま. This holds for a broad range of chromophore systems and for a variety of pulse shapes. We also analyse the absorption probability in the opposite long pulses (narrow bandwidth) regime. We also derive an expression for the long time emission rate in the presence of phonons and use it to study the difference between collective versus independent emission. Finally, we present a numerical simulation for the LHCII monomer (14-mer) system under single photon excitation that illustrates the use of the double hierarchy for calculation of Fock state excitation of a light harvesting system including chromophore coupling to a non-Markovian phonon bath. △ Less

Submitted 24 November, 2021; v1 submitted 12 November, 2021; originally announced November 2021.

Comments: 54 pages, 17 figures; revised abstract, added acknowledgements

Journal ref: J. Chem. Phys. 156, 244108 (2022)

A Quantum trajectory picture of single photon absorption and energy transport in photosystem II

Authors: Robert L. Cook, Liwen Ko, K. Birgitta Whaley

Abstract: In this work we study the first step in photosynthesis for the limiting case of a single photon interacting with photosystem II (PSII). We model our system using quantum trajectory theory, which allows us to consider not only the average evolution, but also the conditional evolution of the system given individual realizations of idealized measurements of photons that have been absorbed and subsequ… ▽ More In this work we study the first step in photosynthesis for the limiting case of a single photon interacting with photosystem II (PSII). We model our system using quantum trajectory theory, which allows us to consider not only the average evolution, but also the conditional evolution of the system given individual realizations of idealized measurements of photons that have been absorbed and subsequently emitted as fluorescence. The quantum nature of the single photon input requires a fully quantum model of both the input and output light fields. We show that PSII coupled to the field via three collective ``bright states'', whose orientation and distribution correlate strongly with its natural geometry. Measurements of the transmitted beam strongly affects the system state, since a (null) detection of the outgoing photon confirms that the system must be in the electronic (excited) ground state. Using numerical and analytical calculations we show that observing the null result transforms a state with a low excited state population $O( 10^{-5} )$ to a state with nearly all population contained in the excited states. This is solely a property of the single photon input, as we confirm by comparing this behavior with that for excitation by a coherent state possessing an average of one photon, using a smaller five site ``pentamer'' system. We also examine the effect of a dissipative phononic environment on the conditional excited state dynamics. We show that the environment has a strong effect on the observed rates of fluorescence, which could act as a new photon-counting witness of excitonic coherence. The long time evolution of the phononic model predicts an experimentally consistent quantum efficiency of 92%. △ Less

Submitted 8 August, 2022; v1 submitted 24 October, 2021; originally announced October 2021.

Comments: 22 pages, 8 figures

Interplay of vibration- and environment-assisted energy transfer

Authors: Zeng-Zhao Li, Liwen Ko, Zhibo Yang, Mohan Sarovar, K. Birgitta Whaley

Abstract: We study the interplay between two environmental influences on excited state energy transfer in photosynthetic light harvesting complexes, namely, vibrationally assisted energy transfer (VAET) and environment-assisted quantum transport (ENAQT), considering a dimeric chromophore donor-acceptor model as a prototype for larger systems. We demonstrate how the basic features of the excitonic energy tra… ▽ More We study the interplay between two environmental influences on excited state energy transfer in photosynthetic light harvesting complexes, namely, vibrationally assisted energy transfer (VAET) and environment-assisted quantum transport (ENAQT), considering a dimeric chromophore donor-acceptor model as a prototype for larger systems. We demonstrate how the basic features of the excitonic energy transfer are influenced by these two environments, both separately and together, with the environment being fully quantum in the case of VAET and treated in the Haken-Strobl-Reineker classical limit in the case of ENAQT. Our results reveal that in the weak noise regime, the presence of a classical noise source is detrimental to the energy transfer that is resonantly assisted by the exciton-vibration interactions intrinsic to VAET. In the strong noise regime we reproduce all the features of ENAQT including the turnover into a Zeno regime where energy transfer is suppressed, and VAET is insignificant. △ Less

Submitted 25 March, 2022; v1 submitted 16 October, 2021; originally announced October 2021.

Comments: 6 pages + Appendices

Journal ref: New J. Phys. 24 033032 (2022)

Unraveling excitation energy transfer assisted by collective behaviors of vibrations

Authors: Zeng-Zhao Li, Liwen Ko, Zhibo Yang, Mohan Sarovar, K. Birgitta Whaley

Abstract: We investigate how collective behaviors of vibrations such as cooperativity and interference can enhance energy transfer in a nontrivial way, focusing on an example of a donor-bridge-acceptor trimeric chromophore system coupled to two vibrational degrees of freedom. Employing parameters selected to provide an overall uphill energy transfer from donor to acceptor, we use numerical calculations of d… ▽ More We investigate how collective behaviors of vibrations such as cooperativity and interference can enhance energy transfer in a nontrivial way, focusing on an example of a donor-bridge-acceptor trimeric chromophore system coupled to two vibrational degrees of freedom. Employing parameters selected to provide an overall uphill energy transfer from donor to acceptor, we use numerical calculations of dynamics in a coupled exciton-vibration basis, together with perturbation-based analytics and calculation of vibronic spectra, we identify clear spectral features of single- and multi-phonon vibrationally-assisted energy transfer (VAET) dynamics, where the latter include up to six-phonon contributions. We identify signatures of vibrational cooperation and interference that provide enhancement of energy transfer relative to that obtained from VAET with a single vibrational mode. We observe a phononic analogue of two-photon absorption, as well as a novel heteroexcitation mechanism in which a single phonon gives rise to simultaneous excitation of both the trimeric system and the vibrational degrees of freedom. The impact of vibrations and of the one- and two-phonon VAET processes on the energy transfer are seen to be quite different in the weak and strong site-vibration coupling regimes. In the weak coupling regime, two-phonon processes dominate, whereas in the strong coupling regime up to six-phonon VAET processes can be induced. The VAET features are seen to be enhanced with increasing temperature and site-vibration coupling strength, and are reduced in the presence of dissipation. We analyze the dependence of these phenomena on the explicit form of the chromophore-vibration couplings, with comparison of VAET spectra for local and non-local couplings. △ Less

Submitted 8 October, 2020; originally announced October 2020.

Journal ref: New J. Phys. 23, 073012 (2021)

Are multi-quasiparticle interactions important in molecular ionization?

Authors: Carlos Mejuto-Zaera, Guorong Weng, Mariya Romanova, Stephen J. Cotton, K. Birgitta Whaley, Norm M. Tubman, Vojtěch Vlček

Abstract: Photo-emission spectroscopy directly probes individual electronic states, ranging from single excitations to high-energy satellites, which simultaneously represent multiple quasiparticles (QPs) and encode information about electronic correlation. First-principles description of the spectra requires an efficient and accurate treatment of all many-body effects. This is especially challenging for inn… ▽ More Photo-emission spectroscopy directly probes individual electronic states, ranging from single excitations to high-energy satellites, which simultaneously represent multiple quasiparticles (QPs) and encode information about electronic correlation. First-principles description of the spectra requires an efficient and accurate treatment of all many-body effects. This is especially challenging for inner valence excitations where the single QP picture breaks down. Here, we provide the full valence spectra of small closed-shell molecules, exploring the independent and interacting quasiparticle regimes, computed with the fully-correlated adaptive sampling configuration interaction (ASCI) method. We critically compare these results to calculations with the many-body perturbation theory, based on the $GW$ and vertex corrected $GWΓがんま$ approaches. The latter explicitly accounts for two-QP quantum interactions, which have been often neglected. We demonstrate that for molecular systems, the vertex correction universally improves the theoretical spectra, and it is crucial for accurate prediction of QPs as well as capturing the rich satellite structures of high-energy excitations. $GWΓがんま$ offers a unified description across all relevant energy scales. Our results suggest that the multi-QP regime corresponds to dynamical correlations, which can be described via perturbation theory. △ Less

Submitted 1 December, 2020; v1 submitted 4 September, 2020; originally announced September 2020.

Comments: 14 pages, 3 figures, plus SI

Journal ref: J. Chem. Phys. 154, 121101 (2021)

The Ground State Electronic Energy of Benzene

Authors: Janus J. Eriksen, Tyler A. Anderson, J. Emiliano Deustua, Khaldoon Ghanem, Diptarka Hait, Mark R. Hoffmann, Seunghoon Lee, Daniel S. Levine, Ilias Magoulas, Jun Shen, Norman M. Tubman, K. Birgitta Whaley, Enhua Xu, Yuan Yao, Ning Zhang, Ali Alavi, Garnet Kin-Lic Chan, Martin Head-Gordon, Wenjian Liu, Piotr Piecuch, Sandeep Sharma, Seiichiro L. Ten-no, C. J. Umrigar, Jürgen Gauss

Abstract: We report on the findings of a blind challenge devoted to determining the frozen-core, full configuration interaction (FCI) ground state energy of the benzene molecule in a standard correlation-consistent basis set of double-$ζぜーた$ quality. As a broad international endeavour, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire of… ▽ More We report on the findings of a blind challenge devoted to determining the frozen-core, full configuration interaction (FCI) ground state energy of the benzene molecule in a standard correlation-consistent basis set of double-$ζぜーた$ quality. As a broad international endeavour, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire offered by modern electronic structure theory. In our assessment, the evaluated high-level methods are all found to qualitatively agree on a final correlation energy, with most methods yielding an estimate of the FCI value around $-863$ m$E_{\text{H}}$. However, we find the root-mean-square deviation of the energies from the studied methods to be considerable (1.3 m$E_{\text{H}}$), which in light of the acclaimed performance of each of the methods for smaller molecular systems clearly displays the challenges faced in extending reliable, near-exact correlation methods to larger systems. While the discrepancies exposed by our study thus emphasize the fact that the current state-of-the-art approaches leave room for improvement, we still expect the present assessment to provide a valuable community resource for benchmark and calibration purposes going forward. △ Less

Submitted 7 October, 2020; v1 submitted 6 August, 2020; originally announced August 2020.

Comments: 29 pages, 1 figure, 2 tables. SI as an ancillary file

Journal ref: J. Phys. Chem. Lett., 11, 8922 (2020)

Efficient phase-factor evaluation in quantum signal processing

Authors: Yulong Dong, Xiang Meng, K. Birgitta Whaley, Lin Lin

Abstract: Quantum signal processing (QSP) is a powerful quantum algorithm to exactly implement matrix polynomials on quantum computers. Asymptotic analysis of quantum algorithms based on QSP has shown that asymptotically optimal results can in principle be obtained for a range of tasks, such as Hamiltonian simulation and the quantum linear system problem. A further benefit of QSP is that it uses a minimal n… ▽ More Quantum signal processing (QSP) is a powerful quantum algorithm to exactly implement matrix polynomials on quantum computers. Asymptotic analysis of quantum algorithms based on QSP has shown that asymptotically optimal results can in principle be obtained for a range of tasks, such as Hamiltonian simulation and the quantum linear system problem. A further benefit of QSP is that it uses a minimal number of ancilla qubits, which facilitates its implementation on near-to-intermediate term quantum architectures. However, there is so far no classically stable algorithm allowing computation of the phase factors that are needed to build QSP circuits. Existing methods require the usage of variable precision arithmetic and can only be applied to polynomials of relatively low degree. We present here an optimization based method that can accurately compute the phase factors using standard double precision arithmetic operations. We demonstrate the performance of this approach with applications to Hamiltonian simulation, eigenvalue filtering, and the quantum linear system problems. Our numerical results show that the optimization algorithm can find phase factors to accurately approximate polynomials of degree larger than $10,000$ with error below $10^{-12}$. △ Less

Submitted 10 July, 2021; v1 submitted 26 February, 2020; originally announced February 2020.

Comments: 26 pages, 17 figures

Journal ref: Phys. Rev. A 103, 042419 (2021)

CASSCF with Extremely Large Active Spaces using the Adaptive Sampling Configuration Interaction Method

Authors: Daniel S. Levine, Diptarka Hait, Norm M. Tubman, Susi Lehtola, K. Birgitta Whaley, Martin Head-Gordon

Abstract: The complete active space self-consistent field (CASSCF) method is the principal approach employed for studying strongly correlated systems. However, exact CASSCF can only be performed on small active spaces of ~20 electrons in ~20 orbitals due to exponential growth in the computational cost. We show that employing the Adaptive Sampling Configuration Interaction (ASCI) method as an approximate Ful… ▽ More The complete active space self-consistent field (CASSCF) method is the principal approach employed for studying strongly correlated systems. However, exact CASSCF can only be performed on small active spaces of ~20 electrons in ~20 orbitals due to exponential growth in the computational cost. We show that employing the Adaptive Sampling Configuration Interaction (ASCI) method as an approximate Full CI solver in the active space allows CASSCF-like calculations within chemical accuracy (<1 kcal/mol for relative energies) in active spaces with more than ~50 active electrons in ~50 active orbitals, significantly increasing the sizes of systems amenable to accurate multiconfigurational treatment. The main challenge with using any selected CI-based approximate CASSCF is the orbital optimization problem; they tend to exhibit large numbers of local minima in orbital space due to their lack of invariance to active-active rotations (in addition to the local minima that exist in exact CASSCF). We highlight methods that can avoid spurious local extrema as a practical solution to the orbital optimization problem. We employ ASCI-SCF to demonstrate lack of polyradical character in moderately sized periacenes with up to 52 correlated electrons and compare against heat-bath CI on an iron porphyrin system with more than 40 correlated electrons. △ Less

Submitted 4 February, 2020; v1 submitted 18 December, 2019; originally announced December 2019.

Comments: 44 pages, 5 figures

Journal ref: J. Chem. Theory Comput. 16, 2340 (2020)

A Non-Orthogonal Variational Quantum Eigensolver

Authors: William J. Huggins, Joonho Lee, Unpil Baek, Bryan O'Gorman, K. Birgitta Whaley

Abstract: Variational algorithms for strongly correlated chemical and materials systems are one of the most promising applications of near-term quantum computers. We present an extension to the variational quantum eigensolver that approximates the ground state of a system by solving a generalized eigenvalue problem in a subspace spanned by a collection of parametrized quantum states. This allows for the sys… ▽ More Variational algorithms for strongly correlated chemical and materials systems are one of the most promising applications of near-term quantum computers. We present an extension to the variational quantum eigensolver that approximates the ground state of a system by solving a generalized eigenvalue problem in a subspace spanned by a collection of parametrized quantum states. This allows for the systematic improvement of a logical wavefunction ansatz without a significant increase in circuit complexity. To minimize the circuit complexity of this approach, we propose a strategy for efficiently measuring the Hamiltonian and overlap matrix elements between states parametrized by circuits that commute with the total particle number operator. We also propose a classical Monte Carlo scheme to estimate the uncertainty in the ground state energy caused by a finite number of measurements of the matrix elements. We explain how this Monte Carlo procedure can be extended to adaptively schedule the required measurements, reducing the number of circuit executions necessary for a given accuracy. We apply these ideas to two model strongly correlated systems, a square configuration of H$_4$ and the $πぱい$-system of Hexatriene (C$_6$H$_8$). △ Less

Submitted 12 February, 2020; v1 submitted 19 September, 2019; originally announced September 2019.

Efficient and Noise Resilient Measurements for Quantum Chemistry on Near-Term Quantum Computers

Authors: William J. Huggins, Jarrod McClean, Nicholas Rubin, Zhang Jiang, Nathan Wiebe, K. Birgitta Whaley, Ryan Babbush

Abstract: Variational algorithms are a promising paradigm for utilizing near-term quantum devices for modeling electronic states of molecular systems. However, previous bounds on the measurement time required have suggested that the application of these techniques to larger molecules might be infeasible. We present a measurement strategy based on a low rank factorization of the two-electron integral tensor.… ▽ More Variational algorithms are a promising paradigm for utilizing near-term quantum devices for modeling electronic states of molecular systems. However, previous bounds on the measurement time required have suggested that the application of these techniques to larger molecules might be infeasible. We present a measurement strategy based on a low rank factorization of the two-electron integral tensor. Our approach provides a cubic reduction in term groupings over prior state-of-the-art and enables measurement times three orders of magnitude smaller than those suggested by commonly referenced bounds for the largest systems we consider. Although our technique requires execution of a linear-depth circuit prior to measurement, this is compensated for by eliminating challenges associated with sampling non-local Jordan-Wigner transformed operators in the presence of measurement error, while enabling a powerful form of error mitigation based on efficient postselection. We numerically characterize these benefits with noisy quantum circuit simulations for ground state energies of strongly correlated electronic systems. △ Less

Submitted 25 June, 2020; v1 submitted 30 July, 2019; originally announced July 2019.

Journal ref: npj Quantum Inf 7, 23 (2021)

Efficient hybridization fitting for dynamical mean-field theory via semi-definite relaxation

Authors: Carlos Mejuto-Zaera, Leonardo Zepeda-Núñez, Michael Lindsey, Norm Tubman, K. Birgitta Whaley, Lin Lin

Abstract: We introduce a nested optimization procedure using semi-definite relaxation for the fitting step in Hamiltonian-based cluster dynamical mean-field theory (DMFT) methodologies. We show that the proposed method is more efficient and flexible than state-of-the-art fitting schemes, which allows us to treat as large a number of bath sites as the impurity solver at hand allows. We characterize its robus… ▽ More We introduce a nested optimization procedure using semi-definite relaxation for the fitting step in Hamiltonian-based cluster dynamical mean-field theory (DMFT) methodologies. We show that the proposed method is more efficient and flexible than state-of-the-art fitting schemes, which allows us to treat as large a number of bath sites as the impurity solver at hand allows. We characterize its robustness to initial conditions and symmetry constraints, thus providing conclusive evidence that in the presence of a large bath, our semi-definite relaxation approach can find the correct set of bath parameters without needing to include \emph{a priori} knowledge of the properties that are to be described. We believe this method will be of great use for Hamiltonian-based calculations, simplifying and improving one of the key steps in cluster dynamical mean-field theory calculations. △ Less

Submitted 27 January, 2020; v1 submitted 16 July, 2019; originally announced July 2019.

Comments: 23 pages, 12 figures, 5 tables. SI: 10 pages, 1 figure, 16 tables

Journal ref: Phys. Rev. B 101, 035143 (2020)

Generalized swap networks for near-term quantum computing

Authors: Bryan O'Gorman, William J. Huggins, Eleanor G. Rieffel, K. Birgitta Whaley

Abstract: The practical use of many types of near-term quantum computers requires accounting for their limited connectivity. One way of overcoming limited connectivity is to insert swaps in the circuit so that logical operations can be performed on physically adjacent qubits, which we refer to as solving the `routing via matchings' problem. We address the routing problem for families of quantum circuits def… ▽ More The practical use of many types of near-term quantum computers requires accounting for their limited connectivity. One way of overcoming limited connectivity is to insert swaps in the circuit so that logical operations can be performed on physically adjacent qubits, which we refer to as solving the `routing via matchings' problem. We address the routing problem for families of quantum circuits defined by a hypergraph wherein each hyperedge corresponds to a potential gate. Our main result is that any unordered set of $k$-qubit gates on distinct $k$-qubit subsets of $n$ logical qubits can be ordered and parallelized in $O(n^{k-1})$ depth using a linear arrangement of $n$ physical qubits; the construction is completely general and achieves optimal scaling in the case where gates acting on all $\binom{n}{k}$ sets of $k$ qubits are desired. We highlight two classes of problems for which our method is particularly useful. First, it applies to sets of mutually commuting gates, as in the (diagonal) phase separators of Quantum Alternating Operator Ansatz (Quantum Approximate Optimization Algorithm) circuits. For example, a single level of a QAOA circuit for Maximum Cut can be implemented in linear depth, and a single level for $3$-SAT in quadratic depth. Second, it applies to sets of gates that do not commute but for which compilation efficiency is the dominant criterion in their ordering. In particular, it can be adapted to Trotterized time-evolution of fermionic Hamiltonians under the Jordan-Wigner transformation, and also to non-standard mixers in QAOA. Using our method, a single Trotter step of the electronic structure Hamiltonian in an arbitrary basis of $n$ orbitals can be done in $O(n^3)$ depth while a Trotter step of the unitary coupled cluster singles and doubles method can be implemented in $O(n^2 ηいーた)$ depth, where $ηいーた$ is the number of electrons. △ Less

Submitted 13 May, 2019; originally announced May 2019.

Generalized Unitary Coupled Cluster Wavefunctions for Quantum Computation

Authors: Joonho Lee, William J. Huggins, Martin Head-Gordon, K. Birgitta Whaley

Abstract: We introduce a unitary coupled-cluster (UCC) ansatz termed $k$-UpCCGSD that is based on a family of sparse generalized doubles (D) operators which provides an affordable and systematically improvable unitary coupled-cluster wavefunction suitable for implementation on a near-term quantum computer. $k$-UpCCGSD employs $k$ products of the exponential of pair coupled-cluster double excitation operator… ▽ More We introduce a unitary coupled-cluster (UCC) ansatz termed $k$-UpCCGSD that is based on a family of sparse generalized doubles (D) operators which provides an affordable and systematically improvable unitary coupled-cluster wavefunction suitable for implementation on a near-term quantum computer. $k$-UpCCGSD employs $k$ products of the exponential of pair coupled-cluster double excitation operators (pCCD), together with generalized single (S) excitation operators. We compare its performance in both efficiency of implementation and accuracy with that of the generalized UCC ansatz employing the full generalized SD excitation operators (UCCGSD), as well as with the standard ansatz employing only SD excitations (UCCSD). $k$-UpCCGSD is found to show the best scaling for quantum computing applications, requiring a circuit depth of $\mathcal O(kN)$, compared with $\mathcal O(N^3)$ for UCCGSD and $\mathcal O((N-ηいーた)^2 ηいーた)$ for UCCSD where $N$ is the number of spin orbitals and $ηいーた$ is the number of electrons. We analyzed the accuracy of these three ansätze by making classical benchmark calculations on the ground state and the first excited state of H$_4$ (STO-3G, 6-31G), H$_2$O (STO-3G), and N$_2$ (STO-3G), making additional comparisons to conventional coupled cluster methods. The results for ground states show that $k$-UpCCGSD offers a good tradeoff between accuracy and cost, achieving chemical accuracy for lower cost of implementation on quantum computers than both UCCGSD and UCCSD. Excited states are calculated with an orthogonally constrained variational quantum eigensolver approach. This is seen to generally yield less accurate energies than for the corresponding ground states. We demonstrate that using a specialized multi-determinantal reference state constructed from classical linear response calculations allows these excited state energetics to be improved. △ Less

Submitted 3 January, 2019; v1 submitted 4 October, 2018; originally announced October 2018.

Postponing the orthogonality catastrophe: efficient state preparation for electronic structure simulations on quantum devices

Authors: Norm M. Tubman, Carlos Mejuto-Zaera, Jeffrey M. Epstein, Diptarka Hait, Daniel S. Levine, William Huggins, Zhang Jiang, Jarrod R. McClean, Ryan Babbush, Martin Head-Gordon, K. Birgitta Whaley

Abstract: Despite significant work on resource estimation for quantum simulation of electronic systems, the challenge of preparing states with sufficient ground state support has so far been largely neglected. In this work we investigate this issue in several systems of interest, including organic molecules, transition metal complexes, the uniform electron gas, Hubbard models, and quantum impurity models ar… ▽ More Despite significant work on resource estimation for quantum simulation of electronic systems, the challenge of preparing states with sufficient ground state support has so far been largely neglected. In this work we investigate this issue in several systems of interest, including organic molecules, transition metal complexes, the uniform electron gas, Hubbard models, and quantum impurity models arising from embedding formalisms such as dynamical mean-field theory. Our approach uses a state-of-the-art classical technique for high-fidelity ground state approximation. We find that easy-to-prepare single Slater determinants such as the Hartree-Fock state often have surprisingly robust support on the ground state for many applications of interest. For the most difficult systems, single-determinant reference states may be insufficient, but low-complexity reference states may suffice. For this we introduce a method for preparation of multi-determinant states on quantum computers. △ Less

Submitted 14 September, 2018; originally announced September 2018.

Comments: 8 pages + SI, 5 figures

Molecular Mechanics Simulations and Improved Tight-binding Hamiltonians for Artificial Light Harvesting Systems: Predicting Geometric Distributions, Disorder, and Spectroscopy of Chromophores in a Protein Environment

Authors: Joonho Lee, Donghyun Lee, Aleksey Kocherzhenko, Loren Greenman, Daniel T. Finley, Matthew B. Francis, K. Birgitta Whaley

Abstract: We present molecular mechanics {and spectroscopic} calculations on prototype artificial light harvesting systems consisting of chromophores attached to a tobacco mosaic virus (TMV) protein scaffold. These systems have been synthesized and characterized spectroscopically, but information about the microscopic configurations and geometry of these TMV-templated chromophore assemblies is largely unkno… ▽ More We present molecular mechanics {and spectroscopic} calculations on prototype artificial light harvesting systems consisting of chromophores attached to a tobacco mosaic virus (TMV) protein scaffold. These systems have been synthesized and characterized spectroscopically, but information about the microscopic configurations and geometry of these TMV-templated chromophore assemblies is largely unknown. We use a Monte Carlo conformational search algorithm to determine the preferred positions and orientations of two chromophores, Coumarin 343 together with its linker, and Oregon Green 488, when these are attached at two different sites (104 and 123) on the TMV protein. The resulting geometric information shows that the extent of disorder and aggregation properties, and therefore the optical properties of the TMV-templated chromophore assembly, are highly dependent on the choice of chromophores and protein site to which they are bound. We used the results of the conformational search as geometric parameters together with an improved tight-binding Hamiltonian to simulate the linear absorption spectra and compare with experimental spectral measurements. The ideal dipole approximation to the Hamiltonian is not valid since the distance between chromophores can be very small. We found that using the geometries from the conformational search is necessary to reproduce the features of the experimental spectral peaks. △ Less

Submitted 3 January, 2019; v1 submitted 7 September, 2018; originally announced September 2018.

An efficient deterministic perturbation theory for selected configuration interaction methods

Authors: Norm M. Tubman, Daniel S. Levine, Diptarka Hait, Martin Head-Gordon, K. Birgitta Whaley

Abstract: The interplay between advances in stochastic and deterministic algorithms has recently led to development of interesting new selected configuration interaction (SCI) methods for solving the many body Schrödinger equation. The performance of these SCI methods can be greatly improved with a second order perturbation theory (PT2) correction, which is often evaluated in a stochastic or hybrid-stochast… ▽ More The interplay between advances in stochastic and deterministic algorithms has recently led to development of interesting new selected configuration interaction (SCI) methods for solving the many body Schrödinger equation. The performance of these SCI methods can be greatly improved with a second order perturbation theory (PT2) correction, which is often evaluated in a stochastic or hybrid-stochastic manner. In this work, we present a highly efficient, fully deterministic PT2 algorithm for SCI methods and demonstrate that our approach is orders of magnitude faster than recent proposals for stochastic SCI+PT2. We also show that it is important to have a compact reference SCI wave function, in order to obtain optimal SCI+PT2 energies. This indicates that it advantageous to use accurate search algorithms such as 'ASCI search' rather than more approximate approaches. Our deterministic PT2 algorithm is based on sorting techniques that have been developed for modern computing architectures and is inherently straightforward to use on parallel computing architectures. Related architectures such as GPU implementations can be also used to further increase the efficiency. Overall, we demonstrate that the algorithms presented in this work allow for efficient evaluation of trillions of PT2 contributions with modest computing resources. △ Less

Submitted 6 August, 2018; originally announced August 2018.

Comments: 13 pages, 2 figures

Modern Approaches to Exact Diagonalization and Selected Configuration Interaction with the Adaptive Sampling CI Method

Authors: Norm M. Tubman, C. Daniel Freeman, Daniel S. Levine, Diptarka Hait, Martin Head-Gordon, K. Birgitta Whaley

Abstract: Recent advances in selected CI, including the adaptive sampling configuration interaction (ASCI) algorithm and its heat bath extension, have made the ASCI approach competitive with the most accurate techniques available, and hence an increasingly powerful tool in solving quantum Hamiltonians. In this work, we show that a useful paradigm for generating efficient selected CI/exact diagonalization al… ▽ More Recent advances in selected CI, including the adaptive sampling configuration interaction (ASCI) algorithm and its heat bath extension, have made the ASCI approach competitive with the most accurate techniques available, and hence an increasingly powerful tool in solving quantum Hamiltonians. In this work, we show that a useful paradigm for generating efficient selected CI/exact diagonalization algorithms is driven by fast sorting algorithms, much in the same way iterative diagonalization is based on the paradigm of matrix vector multiplication. We present several new algorithms for all parts of performing a selected CI, which includes new ASCI search, dynamic bit masking, fast orbital rotations, fast diagonal matrix elements, and residue arrays. The algorithms presented here are fast and scalable, and we find that because they are built on fast sorting algorithms they are more efficient than all other approaches we considered. After introducing these techniques we present ASCI results applied to a large range of systems and basis sets in order to demonstrate the types of simulations that can be practically treated at the full-CI level with modern methods and hardware, presenting double- and triple-zeta benchmark data for the G1 dataset. The largest of these calculations is Si$_{2}$H$_{6}$ which is a simulation of 34 electrons in 152 orbitals. We also present some preliminary results for fast deterministic perturbation theory simulations that use hash functions to maintain high efficiency for treating large basis sets. △ Less

Submitted 28 December, 2019; v1 submitted 2 July, 2018; originally announced July 2018.

Comments: 22 pages,8 figures, 15 tables (added supplemental information on Cr2 in the svp basis)

Journal ref: J. Chem. Theory Comput. 2020, 16, 4, 2139-2159

Quantum phases of dipolar rotors on two-dimensional lattices

Authors: Brendan P. Abolins, Robert E. Zillich, K. Birgitta Whaley

Abstract: The quantum phase transitions of dipoles confined to the vertices of two dimensional (2D) lattices of square and triangular geometry is studied using path integral ground state quantum Monte Carlo (PIGS). We analyze the phase diagram as a function of the strength of both the dipolar interaction and a transverse electric field. The study reveals the existence of a class of orientational phases of q… ▽ More The quantum phase transitions of dipoles confined to the vertices of two dimensional (2D) lattices of square and triangular geometry is studied using path integral ground state quantum Monte Carlo (PIGS). We analyze the phase diagram as a function of the strength of both the dipolar interaction and a transverse electric field. The study reveals the existence of a class of orientational phases of quantum dipolar rotors whose properties are determined by the ratios between the strength anisotropic dipole-dipole interaction, the strength of the applied transverse field, and the rotational constant. For the triangular lattice, the generic orientationally disordered phase found at zero and weak values of both dipolar interaction strength and applied field, is found to show a transition to a phase characterized by net polarization in the lattice plane as the strength of the dipole-dipole interaction is increased, independent of the strength of the applied transverse field, in addition to the expected transition to a transverse polarized phase as the electric field strength increases. The square lattice is also found to exhibit a transition from a disordered phase to an ordered phase as the dipole-dipole interaction strength is increased, as well as the expected transition to a transverse polarized phase as the electric field strength increases. In contrast to the situation with a triangular lattice, on square lattices the ordered phase at high dipole-dipole interaction strength possesses a striped ordering. The properties of these quantum dipolar rotor phases are dominated by the anisotropy of the interaction and provide useful models for developing quantum phases beyond the well-known paradigms of spin Hamiltonian models, realizing in particular a novel physical realization of a quantum rotor-like Hamiltonian that possesses an anisotropic long range interaction. △ Less

Submitted 8 March, 2018; v1 submitted 6 March, 2018; originally announced March 2018.

Comments: Updated credit line and changed line spacing

Journal ref: J. Chem. Phys. 148, 102338 (2018)

Single-photon absorption by single photosynthetic light-harvesting complexes

Authors: Herman C. H. Chan, Omar E. Gamel, Graham R. Fleming, K. Birgitta Whaley

Abstract: We provide a unified theoretical approach to the quantum dynamics of absorption of single photons and subsequent excitonic energy transfer in photosynthetic light-harvesting complexes. Our analysis combines a continuous mode <n>-photon quantum optical master equation for the chromophoric system with the hierarchy of equations of motion describing excitonic dynamics in presence of non-Markovian cou… ▽ More We provide a unified theoretical approach to the quantum dynamics of absorption of single photons and subsequent excitonic energy transfer in photosynthetic light-harvesting complexes. Our analysis combines a continuous mode <n>-photon quantum optical master equation for the chromophoric system with the hierarchy of equations of motion describing excitonic dynamics in presence of non-Markovian coupling to vibrations of the chromophores and surrounding protein. We apply the approach to simulation of absorption of single-photon coherent states by pigment-protein complexes containing between one and seven chromophores, and compare with results obtained by excitation using a thermal radiation field. We show that the values of excitation probability obtained under single-photon absorption conditions can be consistently related to bulk absorption cross-sections. Analysis of the timescale and efficiency of single-photon absorption by light-harvesting systems within this full quantum description of pigment-protein dynamics coupled to a quantum radiation field reveals a non-trivial dependence of the excitation probability and the excited state dynamics induced by exciton-phonon coupling during and subsequent to the pulse, on the bandwidth of the incident photon pulse. For bandwidths equal to the spectral bandwidth of Chlorophyll a, our results yield an estimation of an average time of ~0.09 s for a single chlorophyll chromophore to absorb the energy equivalent of one (single-polarization) photon under irradiation by single-photon states at the intensity of sunlight. △ Less

Submitted 13 January, 2018; originally announced January 2018.

Comments: 24 pages, 12 figures

Journal ref: J. Phys. B: At. Mol. Opt. Phys. 51, 054002 (2018)

Cluster decomposition of full configuration interaction wave functions: a tool for chemical interpretation of systems with strong correlation

Authors: Susi Lehtola, Norm M. Tubman, K. Birgitta Whaley, Martin Head-Gordon

Abstract: Approximate full configuration interaction (FCI) calculations have recently become tractable for systems of unforeseen size thanks to stochastic and adaptive approximations to the exponentially scaling FCI problem. The result of an FCI calculation is a weighted set of electronic configurations, which can also be expressed in terms of excitations from a reference configuration. The excitation ampli… ▽ More Approximate full configuration interaction (FCI) calculations have recently become tractable for systems of unforeseen size thanks to stochastic and adaptive approximations to the exponentially scaling FCI problem. The result of an FCI calculation is a weighted set of electronic configurations, which can also be expressed in terms of excitations from a reference configuration. The excitation amplitudes contain information on the complexity of the electronic wave function, but this information is contaminated by contributions from disconnected excitations, i.e. those excitations that are just products of independent lower-level excitations. The unwanted contributions can be removed via a cluster decomposition procedure, making it possible to examine the importance of connected excitations in complicated multireference molecules which are outside the reach of conventional algorithms. We present an implementation of the cluster decomposition analysis and apply it to both true FCI wave functions, as well as wave functions generated from the adaptive sampling CI (ASCI) algorithm. The cluster decomposition is useful for interpreting calculations in chemical studies, as a diagnostic for the convergence of various excitation manifolds, as well as as a guidepost for polynomially scaling electronic structure models. Applications are presented for (i) the double dissociation of water, (ii) the carbon dimer, (iii) the πぱい space of polyacenes, as well as (iv) the chromium dimer. While the cluster amplitudes exhibit rapid decay with increasing rank for the first three systems, even connected octuple excitations still appear important in Cr$_2$, suggesting that spin-restricted single-reference coupled-cluster approaches may not be tractable for some problems in transition metal chemistry. △ Less

Submitted 24 October, 2017; v1 submitted 13 July, 2017; originally announced July 2017.

Comments: 15 pages, 5 figures

Journal ref: J. Chem. Phys. 147, 154105 (2017)

Optimized pulses for Raman excitation through the continuum: verification using multi-configurational time-dependent Hartree-Fock

Authors: Loren Greenman, K. Birgitta Whaley, Daniel J. Haxton, C. William McCurdy

Abstract: We have verified a mechanism for Raman excitation of atoms through continuum levels previously obtained by quantum optimal control using the multi-configurational time-dependent Hartree-Fock (MCTDHF) method. For the optimal control, which requires running multiple propagations to determine the optimal pulse sequence, we used the computationally inexpensive time-dependent configuration interaction… ▽ More We have verified a mechanism for Raman excitation of atoms through continuum levels previously obtained by quantum optimal control using the multi-configurational time-dependent Hartree-Fock (MCTDHF) method. For the optimal control, which requires running multiple propagations to determine the optimal pulse sequence, we used the computationally inexpensive time-dependent configuration interaction singles (TDCIS) method. TDCIS captures all of the necessary correlation of the desired processes but assumes that ionization pathways reached via double excitations are not present. MCTDHF includes these pathways and all multiparticle correlations in a set of time-dependent orbitals. The mechanism that was determined to be optimal in the Raman excitation of the Ne $1s^22s^22p^53p^1$ valence state via the metastable $1s^22s^12p^63p^1$ resonance state involves a sequential resonance-valence excitation. First, a long pump pulse excites the core-hole state, and then a shorter Stokes pulse transfers the population to the valence state. This process represents the first step in a multidimensional x-ray spectroscopy scheme that will provide a local probe of valence electronic correlations. Although at the optimal pulse intensities at the TDCIS level of theory the MCTDHF method predicts multiple ionization of the atom, at slightly lower intensities (reduced by a factor of about 4) the TDCIS mechanism is shown to hold qualitatively. Quantitatively, the MCTDHF populations are reduced from the TDCIS calculations by a factor of 4. △ Less

Submitted 28 December, 2016; v1 submitted 15 September, 2016; originally announced September 2016.

Journal ref: Phys. Rev. A 96, 013411 (2017)

A deterministic alternative to the full configuration interaction quantum Monte Carlo method

Authors: Norm M. Tubman, Joonho Lee, Tyler Y. Takeshita, Martin Head-Gordon, K. Birgitta Whaley

Abstract: Development of exponentially scaling methods has seen great progress in tackling larger systems than previously thought possible. One such technique, full configuration interaction quantum Monte Carlo, is a useful algorithm that allows exact diagonalization through stochastically sampling determinants. The method derives its utility from the information in the matrix elements of the Hamiltonian, a… ▽ More Development of exponentially scaling methods has seen great progress in tackling larger systems than previously thought possible. One such technique, full configuration interaction quantum Monte Carlo, is a useful algorithm that allows exact diagonalization through stochastically sampling determinants. The method derives its utility from the information in the matrix elements of the Hamiltonian, along with a stochastic projected wave function, to find the important parts of Hilbert space. However, the stochastic representation of the wave function is not required to search Hilbert space efficiently, and here we describe a highly efficient deterministic method to achieve chemical accuracy for a wide range of systems, including the difficult Cr$_{2}$ dimer. In addition our method also allows efficient calculation of excited state energies, for which we illustrate with benchmark results for the excited states of C$_{2}$. △ Less

Submitted 8 March, 2016; originally announced March 2016.

Comments: 4 pages, 2 figures

Journal ref: J. Chem. Phys. 145, 044112 (2016)

Long-range coherent energy transport in Photosystem II

Authors: Jan J. J. Roden, Doran I. G. Bennett, K. Birgitta Whaley

Abstract: We simulate the long-range inter-complex electronic energy transfer in Photosystem II -- from the antenna complex, via a core complex, to the reaction center -- using a non-Markovian (ZOFE) quantum master equation description that allows us to quantify the electronic coherence involved in the energy transfer. We identify the pathways of the energy transfer in the network of coupled chromophores, u… ▽ More We simulate the long-range inter-complex electronic energy transfer in Photosystem II -- from the antenna complex, via a core complex, to the reaction center -- using a non-Markovian (ZOFE) quantum master equation description that allows us to quantify the electronic coherence involved in the energy transfer. We identify the pathways of the energy transfer in the network of coupled chromophores, using a description based on excitation probability currents. We investigate how the energy transfer depends on the initial excitation -- localized, coherent initial excitation versus delocalized, incoherent initial excitation -- and find that the energy transfer is remarkably robust with respect to such strong variations of the initial condition. To explore the importance of vibrationally enhanced transfer and to address the question of optimization in the system parameters, we vary the strength of the coupling between the electronic and the vibrational degrees of freedom. We find that the original parameters lie in a (broad) region that enables optimal transfer efficiency, and that the energy transfer appears to be very robust with respect to variations in the vibronic coupling. Nevertheless, vibrationally enhanced transfer appears to be crucial to obtain a high transfer efficiency. We compare our quantum simulation to a "classical" rate equation based on a modified-Redfield/generalized-Förster description that was previously used to simulate energy transfer dynamics in the entire Photosystem II complex, and find very good agreement between quantum and rate-equation simulation of the overall energy transfer dynamics. △ Less

Submitted 27 January, 2015; originally announced January 2015.

Laser pulses for coherent xuv Raman excitation

Authors: Loren Greenman, Christiane P. Koch, K. Birgitta Whaley

Abstract: We combine multi-channel electronic structure theory with quantum optimal control to derive Raman pulse sequences that coherently populate a valence excited state. For a neon atom, Raman target populations of up to 13% are obtained. Superpositions of the ground and valence Raman states with a controllable relative phase are found to be reachable with up to 4.5% population and phase control facilit… ▽ More We combine multi-channel electronic structure theory with quantum optimal control to derive Raman pulse sequences that coherently populate a valence excited state. For a neon atom, Raman target populations of up to 13% are obtained. Superpositions of the ground and valence Raman states with a controllable relative phase are found to be reachable with up to 4.5% population and phase control facilitated by the pump pulse carrier envelope phase. Our results open a route to creating core-hole excitations in molecules and aggregates that locally address specific atoms and represent the first step towards realization of multidimensional spectroscopy in the xuv and x-ray regimes. △ Less

Submitted 27 September, 2014; originally announced September 2014.

Journal ref: Phys. Rev. A 92, 013407 (2015)

Photo-activated biological processes as quantum measurements

Authors: Atac Imamoglu, K. Birgitta Whaley

Abstract: We outline a framework for describing photo-activated biological reactions as generalized quantum measurements of external fields, for which the biological system takes on the role of a quantum meter. By using general arguments regarding the Hamiltonian that describes the measurement interaction, we identify the cases where it is essential for a complex chemical or biological system to exhibit non… ▽ More We outline a framework for describing photo-activated biological reactions as generalized quantum measurements of external fields, for which the biological system takes on the role of a quantum meter. By using general arguments regarding the Hamiltonian that describes the measurement interaction, we identify the cases where it is essential for a complex chemical or biological system to exhibit non-equilibrium quantum coherent dynamics in order to achieve the requisite functionality. We illustrate the analysis by considering measurement of the solar radiation field in photosynthesis and measurement of the earth's magnetic field in avian magnetoreception. △ Less

Submitted 21 August, 2014; originally announced August 2014.

Comments: 9 pages, 2 figures

Infrared-dressed entanglement of cold open-shell polar molecules for universal matchgate quantum computing

Authors: Felipe Herrera, Yudong Cao, Sabre Kais, K. Birgitta Whaley

Abstract: Implementing a scalable quantum information processor using polar molecules in optical lattices requires precise control over the long-range dipole-dipole interaction between molecules in selected lattice sites. We present here a scheme using trapped open-shell $^2Σしぐま$ polar molecules that allows dipolar exchange processes between nearest and next-nearest neighbors to be controlled to construct a ge… ▽ More Implementing a scalable quantum information processor using polar molecules in optical lattices requires precise control over the long-range dipole-dipole interaction between molecules in selected lattice sites. We present here a scheme using trapped open-shell $^2Σしぐま$ polar molecules that allows dipolar exchange processes between nearest and next-nearest neighbors to be controlled to construct a generalized transverse Ising spin Hamiltonian with tunable $XX$, $YY$ and $XY$ couplings in the rotating frame of the driving lasers. The scheme requires a moderately strong bias magnetic field with near-infrared light to provide local tuning of the qubit energy gap, and mid-infrared pulses to perform rotational state transfer via stimulated Raman adiabatic passage. No interaction between qubits is present in the absence of the infrared driving. We analyze the fidelity of the resulting two-qubit matchgate, and demonstrate its robustness as a function of the driving parameters. We discuss a realistic application of the system for universal matchgate quantum computing in optical lattices. △ Less

Submitted 3 February, 2014; originally announced February 2014.

Comments: 18 pages, 4 figures

Journal ref: New J. Phys. 16, 075001 (2014)

Generalized master equation for modular exciton density transfer

Authors: Seogjoo Jang, Stephan Hoyer, Graham Fleming, K. Birgitta Whaley

Abstract: A generalized master equation (GME) governing quantum evolution of modular exciton density (MED) is derived for large scale light harvesting systems composed of weakly interacting modules of multiple chromophores. The GME-MED offers a practical framework to incorporate real time coherent quantum dynamics calculations at small length scales into dynamics over large length scales, without assumption… ▽ More A generalized master equation (GME) governing quantum evolution of modular exciton density (MED) is derived for large scale light harvesting systems composed of weakly interacting modules of multiple chromophores. The GME-MED offers a practical framework to incorporate real time coherent quantum dynamics calculations at small length scales into dynamics over large length scales, without assumptions of time scale separation or specific forms of intra-module quantum dynamics. A test of the GME-MED for four sites of the Fenna-Matthews-Olson complex demonstrates how coherent dynamics of excitonic populations over many coupled chromophores can be accurately described by transitions between subgroups (modules) of delocalized excitons. △ Less

Submitted 8 November, 2013; originally announced November 2013.

Comments: 5 pages, 3 figures

Measurement- and comparison-based sizes of Schrödinger cat states of light

Authors: T. J. Volkoff, K. B. Whaley

Abstract: We extend several measurement-based definitions of effective "cat-size" to coherent state superpositions with branches composed of either single coherent states or tensor products of coherent states. These effective cat-size measures depend on determining the maximal quantum distinguishability of certain states associated with the superposition state: e.g., in one measure, the maximal distinguisha… ▽ More We extend several measurement-based definitions of effective "cat-size" to coherent state superpositions with branches composed of either single coherent states or tensor products of coherent states. These effective cat-size measures depend on determining the maximal quantum distinguishability of certain states associated with the superposition state: e.g., in one measure, the maximal distinguishability of the branches of the superposition is considered as in quantum binary decision theory; in another measure, the maximal distinguishability of the initial superposition and its image after a one-parameter evolution generated by a local Hermitian operator is of interest. The cat-size scaling with the number of modes and mode intensity (i.e., photon number) is compared to the scaling derived directly from the Wigner function of the superposition and to that estimated experimentally from decoherence. We also apply earlier comparison-based methods for determining macroscopic superposition size that require a reference GHZ state. The case of a hierarchical Schrödinger cat state with branches composed of smaller superpositions is also analyzed from a measurement-based perspective. △ Less

Submitted 22 October, 2013; originally announced October 2013.

Journal ref: Phys. Rev. A 89, 012122 (2014)

Ab-Initio Calculation of Molecular Aggregation Effects: a Coumarin-343 Case Study

Authors: Donghyun Lee, Loren Greenman, Mohan Sarovar, K. Birgitta Whaley

Abstract: We present time-dependent density functional theory (TDDFT) calculations for single and dimerized Coumarin-343 molecules in order to investigate the quantum mechanical effects of chromophore aggregation in extended systems designed to function as a new generation of sensors and light-harvesting devices. Using the single-chromophore results, we describe the construction of effective Hamiltonians to… ▽ More We present time-dependent density functional theory (TDDFT) calculations for single and dimerized Coumarin-343 molecules in order to investigate the quantum mechanical effects of chromophore aggregation in extended systems designed to function as a new generation of sensors and light-harvesting devices. Using the single-chromophore results, we describe the construction of effective Hamiltonians to predict the excitonic properties of aggregate systems. We compare the electronic coupling properties predicted by such effective Hamiltonians to those obtained from TDDFT calculations of dimers, and to the coupling predicted by the transition density cube (TDC) method. We determine the accuracy of the dipole-dipole approximation and TDC with respect to the separation distance and orientation of the dimers. In particular, we investigate the effects of including Coulomb coupling terms ignored in the typical tight-binding effective Hamiltonian. We also examine effects of orbital relaxation which cannot be captured by either of these models. △ Less

Submitted 20 September, 2013; originally announced September 2013.

Realistic and verifiable coherent control of excitonic states in a light harvesting complex

Authors: Stephan Hoyer, Filippo Caruso, Simone Montangero, Mohan Sarovar, Tommaso Calarco, Martin B. Plenio, K. Birgitta Whaley

Abstract: We explore the feasibility of coherent control of excitonic dynamics in light harvesting complexes, analyzing the limits imposed by the open nature of these quantum systems. We establish feasible targets for phase and phase/amplitude control of the electronically excited state populations in the Fenna-Mathews-Olson (FMO) complex and analyze the robustness of this control with respect to orientatio… ▽ More We explore the feasibility of coherent control of excitonic dynamics in light harvesting complexes, analyzing the limits imposed by the open nature of these quantum systems. We establish feasible targets for phase and phase/amplitude control of the electronically excited state populations in the Fenna-Mathews-Olson (FMO) complex and analyze the robustness of this control with respect to orientational and energetic disorder, as well as decoherence arising from coupling to the protein environment. We further present two possible routes to verification of the control target, with simulations for the FMO complex showing that steering of the excited state is experimentally verifiable either by extending excitonic coherence or by producing novel states in a pump-probe setup. Our results provide a first step toward coherent control of these complex biological quantum systems in an ultrafast spectroscopy setup. △ Less

Submitted 20 January, 2014; v1 submitted 17 July, 2013; originally announced July 2013.

Comments: 12 pages, 8 figures

Journal ref: New J. Phys. 16 045007 (2014)

Collective Effects in Linear Spectroscopy of Dipole-Coupled Molecular Arrays

Authors: A. A. Kocherzhenko, J. Dawlaty, B. P. Abolins, F. Herrera, D. B. Abraham, K. B. Whaley

Abstract: We present a consistent analysis of linear spectroscopy for arrays of nearest neighbor dipole-coupled two-level molecules that reveals distinct signatures of weak and strong coupling regimes separated for infinite size arrays by a quantum critical point. In the weak coupling regime, the ground state of the molecular array is disordered, but in the strong coupling regime it has (anti)ferroelectric… ▽ More We present a consistent analysis of linear spectroscopy for arrays of nearest neighbor dipole-coupled two-level molecules that reveals distinct signatures of weak and strong coupling regimes separated for infinite size arrays by a quantum critical point. In the weak coupling regime, the ground state of the molecular array is disordered, but in the strong coupling regime it has (anti)ferroelectric ordering. We show that multiple molecular excitations (odd/even in weak/strong coupling regime) can be accessed directly from the ground state. We analyze the scaling of absorption and emission with system size and find that the oscillator strengths show enhanced superradiant behavior in both ordered and disordered phases. As the coupling increases, the single excitation oscillator strength rapidly exceeds the well known Heitler-London value. In the strong coupling regime we show the existence of a unique spectral transition with excitation energy that can be tuned by varying the system size and that asymptotically approaches zero for large systems. The oscillator strength for this transition scales quadratically with system size, showing an anomalous one-photon superradiance. For systems of infinite size, we find a novel, singular spectroscopic signature of the quantum phase transition between disordered and ordered ground states. We outline how arrays of ultra cold dipolar molecules trapped in an optical lattice can be used to access the strong coupling regime and observe the anomalous superradiant effects associated with this regime. △ Less

Submitted 16 October, 2014; v1 submitted 10 July, 2013; originally announced July 2013.

Comments: 12 pages, 7 figures main text

Entanglement creation in cold molecular gases using strong laser pulses

Authors: Felipe Herrera, Sabre Kais, K. Birgitta Whaley

Abstract: While many-particle entanglement can be found in natural solids and strongly interacting atomic and molecular gases, generating highly entangled states between weakly interacting particles in a controlled and scalable way presents a significant challenge. We describe here a one-step method to generate entanglement in a dilute gas of cold polar molecules. For molecules in optical traps separated by… ▽ More While many-particle entanglement can be found in natural solids and strongly interacting atomic and molecular gases, generating highly entangled states between weakly interacting particles in a controlled and scalable way presents a significant challenge. We describe here a one-step method to generate entanglement in a dilute gas of cold polar molecules. For molecules in optical traps separated by a few micrometers, we show that maximally entangled states can be created using the strong off-resonant pulses that are routinely used in molecular alignment experiments. We show that the resulting alignment-mediated entanglement can be detected by measuring laser-induced fluorescence with single-site resolution and that signatures of this molecular entanglement also appear in the microwave absorption spectra of the molecular ensemble. We analyze the robustness of these entangled molecular states with respect to intensity fluctuations of the trapping laser and discuss possible applications of the system for quantum information processing. △ Less

Submitted 26 February, 2013; originally announced February 2013.

Comments: 29 pages, 7 figures

Design principles and fundamental trade-offs in biomimetic light harvesting

Authors: Mohan Sarovar, K. Birgitta Whaley

Abstract: Recent developments in synthetic and supramolecular chemistry have created opportunities to design organic systems with tailored nanoscale structure for various technological applications. A key application area is the capture of light energy and its conversion into electrochemical or chemical forms for photovoltaic or sensing applications. In this work we consider cylindrical assemblies of chromo… ▽ More Recent developments in synthetic and supramolecular chemistry have created opportunities to design organic systems with tailored nanoscale structure for various technological applications. A key application area is the capture of light energy and its conversion into electrochemical or chemical forms for photovoltaic or sensing applications. In this work we consider cylindrical assemblies of chromophores that model structures produced by several supramolecular techniques. Our study is especially guided by the versatile structures produced by virus-templated assembly. We use a multi-objective optimization framework to determine design principles and limitations in light harvesting performance for such assemblies, both in the presence and absence of disorder. We identify a fundamental trade-off in cylindrical assemblies that is encountered when attempting to maximize both efficiency of energy transfer and absorption bandwidth. We also rationalize the optimal design strategies and provide explanations for why various structures provide optimal performance. Most importantly, we find that the optimal design strategies depend on the amount of energetic and structural disorder in the system. The aim of these studies is to develop a program of quantum-informed rational design for construction of organic assemblies that have the same degree of tailored nanoscale structure as biological photosynthetic light harvesting complexes, and also have the potential to reproduce their remarkable light harvesting performance. △ Less

Submitted 25 March, 2013; v1 submitted 7 October, 2012; originally announced October 2012.

Comments: 35 pages, 8 figures. Published version

Journal ref: New J. Phys. 15, 013030 (2013)

Inverting pump-probe spectroscopy for state tomography of excitonic systems

Authors: Stephan Hoyer, K. Birgitta Whaley

Abstract: We propose a two-step protocol for inverting ultrafast spectroscopy experiments on a molecular aggregate to extract the time-evolution of the excited state density matrix. The first step is a deconvolution of the experimental signal to determine a pump-dependent response function. The second step inverts the quantum state of the system from this response function, given a model for how the system… ▽ More We propose a two-step protocol for inverting ultrafast spectroscopy experiments on a molecular aggregate to extract the time-evolution of the excited state density matrix. The first step is a deconvolution of the experimental signal to determine a pump-dependent response function. The second step inverts the quantum state of the system from this response function, given a model for how the system evolves following the probe interaction. We demonstrate this inversion analytically and numerically for a dimer model system, and evaluate the feasibility of scaling it to larger molecular aggregates such as photosynthetic protein-pigment complexes. Our scheme provides a direct alternative to the approach of determining all Hamiltonian parameters and then simulating excited state dynamics. △ Less

Submitted 2 April, 2013; v1 submitted 28 September, 2012; originally announced September 2012.

Comments: 15 pages, 6 figures. This is the accepted version for publication in the Journal of Chemical Physics

Journal ref: J. Chem. Phys. 138, 164102 (2013)

Spatial propagation of excitonic coherence enables ratcheted energy transfer

Authors: Stephan Hoyer, Akihito Ishizaki, K. Birgitta Whaley

Abstract: Experimental evidence shows that a variety of photosynthetic systems can preserve quantum beats in the process of electronic energy transfer, even at room temperature. However, whether this quantum coherence arises in vivo and whether it has any biological function have remained unclear. Here we present a theoretical model that suggests that the creation and recreation of coherence under natural c… ▽ More Experimental evidence shows that a variety of photosynthetic systems can preserve quantum beats in the process of electronic energy transfer, even at room temperature. However, whether this quantum coherence arises in vivo and whether it has any biological function have remained unclear. Here we present a theoretical model that suggests that the creation and recreation of coherence under natural conditions is ubiquitous. Our model allows us to theoretically demonstrate a mechanism for a ratchet effect enabled by quantum coherence, in a design inspired by an energy transfer pathway in the Fenna-Matthews-Olson complex of the green sulfur bacteria. This suggests a possible biological role for coherent oscillations in spatially directing energy transfer. Our results emphasize the importance of analyzing long-range energy transfer in terms of transfer between inter-complex coupling (ICC) states rather than between site or exciton states. △ Less

Submitted 17 September, 2012; v1 submitted 15 June, 2011; originally announced June 2011.

Comments: Accepted version for Phys. Rev. E. 14 pages, 7 figures

Journal ref: Phys. Rev. E 86, 041911 (2012)

Quantum entanglement phenomena in photosynthetic light harvesting complexes

Authors: K. Birgitta Whaley, Mohan Sarovar, Akihito Ishizaki

Abstract: We review recent theoretical calculations of quantum entanglement in photosynthetic light harvesting complexes. These works establish, for the first time, a manifestation of this characteristically quantum mechanical phenomenon in biologically functional structures. We begin by summarizing calculations on model biomolecular systems that aim to reveal non-trivial characteristics of quantum entangle… ▽ More We review recent theoretical calculations of quantum entanglement in photosynthetic light harvesting complexes. These works establish, for the first time, a manifestation of this characteristically quantum mechanical phenomenon in biologically functional structures. We begin by summarizing calculations on model biomolecular systems that aim to reveal non-trivial characteristics of quantum entanglement in non-equilibrium biological environments. We then discuss and compare several calculations performed recently of excitonic dynamics in the Fenna-Matthews-Olson light harvesting complex and of the entanglement present in this widely studied pigment-protein structure. We point out the commonalities between the derived results and also identify and explain the differences. We also discuss recent work that examines entanglement in the structurally more intricate light harvesting complex II (LHCII). During this overview, we take the opportunity to clarify several subtle issues relating to entanglement in such biomolecular systems, including the role of entanglement in biological function, the complexity of dynamical modeling that is required to capture the salient features of entanglement in such biomolecular systems, and the relationship between entanglement and other quantum mechanical features that are observed and predicted in light harvesting complexes. Finally, we suggest possible extensions of the current work and also review the options for experimental confirmation of the predicted entanglement phenomena in light harvesting complexes. △ Less

Submitted 18 December, 2010; originally announced December 2010.

Comments: 11 pages, 3 figures. Review article. To appear in the proceedings of the 22nd Solvay Conference in Chemistry on "Quantum Effects in Chemistry and Biology"

Environmental correlation effects on excitation energy transfer in photosynthetic light harvesting

Authors: Mohan Sarovar, Yuan-Chung Cheng, K. Birgitta Whaley

Abstract: Several recent studies of energy transfer in photosynthetic light harvesting complexes have revealed a subtle interplay between coherent and decoherent dynamic contributions to the overall transfer efficiency in these open quantum systems. In this work we systematically investigate the impact of temporal and spatial correlations in environmental fluctuations on excitation transport in the Fenna-Ma… ▽ More Several recent studies of energy transfer in photosynthetic light harvesting complexes have revealed a subtle interplay between coherent and decoherent dynamic contributions to the overall transfer efficiency in these open quantum systems. In this work we systematically investigate the impact of temporal and spatial correlations in environmental fluctuations on excitation transport in the Fenna-Matthews-Olson photosynthetic complex. We demonstrate that the exact nature of the correlations can have a large impact on the efficiency of light harvesting. In particular, we find that (i) spatial correlations can enhance coherences in the site basis while at the same time slowing transport, and (ii) the overall efficiency of transport is optimized at a finite temporal correlation that produces maximum overlap between the environmental power spectrum and the excitonic energy differences, which in turn results in enhanced driving of transitions between excitonic states. △ Less

Submitted 20 February, 2011; v1 submitted 28 November, 2009; originally announced November 2009.

Comments: 9 pages, 2 figures. New results in this version. Published version

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

Limits of quantum speedup in photosynthetic light harvesting

Authors: Stephan Hoyer, Mohan Sarovar, K. Birgitta Whaley

Abstract: It has been suggested that excitation transport in photosynthetic light harvesting complexes features speedups analogous to those found in quantum algorithms. Here we compare the dynamics in these light harvesting systems to the dynamics of quantum walks, in order to elucidate the limits of such quantum speedups. For the Fenna-Matthews-Olson (FMO) complex of green sulfur bacteria, we show that whi… ▽ More It has been suggested that excitation transport in photosynthetic light harvesting complexes features speedups analogous to those found in quantum algorithms. Here we compare the dynamics in these light harvesting systems to the dynamics of quantum walks, in order to elucidate the limits of such quantum speedups. For the Fenna-Matthews-Olson (FMO) complex of green sulfur bacteria, we show that while there is indeed speedup at short times, this is short lived (70 fs) despite longer lived (ps) quantum coherence. Remarkably, this time scale is independent of the details of the decoherence model. More generally, we show that the distinguishing features of light-harvesting complexes not only limit the extent of quantum speedup but also reduce rates of diffusive transport. These results suggest that quantum coherent effects in biological systems are optimized for efficiency or robustness rather than the more elusive goal of quantum speedup. △ Less

Submitted 7 May, 2010; v1 submitted 12 October, 2009; originally announced October 2009.

Comments: 9 pages, 6 figures. To appear in New Journal Physics, special issue on "Quantum Effects and Noise in Biomolecules." Updated to accepted version

Journal ref: New J. Phys. 12, 065041 (2010)

Quantum entanglement in photosynthetic light harvesting complexes

Authors: Mohan Sarovar, Akihito Ishizaki, Graham R. Fleming, K. Birgitta Whaley

Abstract: Light harvesting components of photosynthetic organisms are complex, coupled, many-body quantum systems, in which electronic coherence has recently been shown to survive for relatively long time scales despite the decohering effects of their environments. Within this context, we analyze entanglement in multi-chromophoric light harvesting complexes, and establish methods for quantification of entan… ▽ More Light harvesting components of photosynthetic organisms are complex, coupled, many-body quantum systems, in which electronic coherence has recently been shown to survive for relatively long time scales despite the decohering effects of their environments. Within this context, we analyze entanglement in multi-chromophoric light harvesting complexes, and establish methods for quantification of entanglement by presenting necessary and sufficient conditions for entanglement and by deriving a measure of global entanglement. These methods are then applied to the Fenna-Matthews-Olson (FMO) protein to extract the initial state and temperature dependencies of entanglement. We show that while FMO in natural conditions largely contains bipartite entanglement between dimerized chromophores, a small amount of long-range and multipartite entanglement exists even at physiological temperatures. This constitutes the first rigorous quantification of entanglement in a biological system. Finally, we discuss the practical utilization of entanglement in densely packed molecular aggregates such as light harvesting complexes. △ Less

Submitted 7 June, 2010; v1 submitted 23 May, 2009; originally announced May 2009.

Comments: 14 pages, 7 figures. Improved presentation, published version

Journal ref: Nature Physics, 6, 462 (2010)

Nanoscale Quantum Solvation of para-H$_2$ around the Linear OCS Molecule inside $^4$He Droplets

Authors: Yongkyung Kwon, K. Birgitta Whaley

Abstract: We present a microscopic analysis of the quantum solvation structures of para-H$_2$ around the OCS molecule when embedded in low temperature $^4$He droplets. The structures of clusters containing M=5 and 6 para-H$_2$ molecules are compared with corresponding structures for M=1 (OCS-H$_2$ complex) and M=17 (a full solvation shell), as well as with the clusters in the absence of helium. We find th… ▽ More We present a microscopic analysis of the quantum solvation structures of para-H$_2$ around the OCS molecule when embedded in low temperature $^4$He droplets. The structures of clusters containing M=5 and 6 para-H$_2$ molecules are compared with corresponding structures for M=1 (OCS-H$_2$ complex) and M=17 (a full solvation shell), as well as with the clusters in the absence of helium. We find that the helium has negligible effect on the structures for the small and large OCS(H$_2$)$_M$ clusters, but that it modifies the cluster structure for M=6. We discuss implications of these results for the onset of superfluidity in the solvating hydrogen shell and for spectroscopic measurements. △ Less

Submitted 12 December, 2003; originally announced December 2003.

Comments: 4 pages, 2 figures, accepted for publication in J. Low Temp. Phys

Structure and energetics of helium adsorption on nanosurfaces

Authors: Patrick Huang, Heather D. Whitley, K. Birgitta Whaley

Abstract: The ground and excited state properties of small helium clusters, 4He_N, containing nanoscale (~3-10 Angstroms) planar aromatic molecules have been studied with quantum Monte Carlo methods. Ground state structures and energies are obtained from importance-sampled, rigid-body diffusion Monte Carlo. Excited state energies due to helium vibrational motion are evaluated using the projection operator… ▽ More The ground and excited state properties of small helium clusters, 4He_N, containing nanoscale (~3-10 Angstroms) planar aromatic molecules have been studied with quantum Monte Carlo methods. Ground state structures and energies are obtained from importance-sampled, rigid-body diffusion Monte Carlo. Excited state energies due to helium vibrational motion are evaluated using the projection operator, imaginary time spectral evolution technique. We examine the adsorption of N helium atoms (N less than or equal to 24) on a series of planar aromatic molecules (benzene, naphthalene, anthracene, tetracene, phthalocyanine). The first layer of helium atoms is well-localized on the molecule surface, and we find well-defined localized excitations due to in-plane vibrational motion of helium on the molecule surface. We discuss the implications of these confined excitations for the molecule spectroscopy. △ Less

Submitted 17 November, 2003; originally announced November 2003.

Comments: 6 pages, 2 figures, QFS 2003 Symposium, submitted to J. Low Temp. Phys

Comparison of Rotational Energies and Rigidity of OCS-paraH_2 and OCS-4He complexes

Authors: R. E. Zillich, K. B. Whaley

Abstract: We analyze the nature of the rotational energy level structure of the OCS-He and OCS-H_2 complexes with a comparison of exact calculations to several differentdynamical approximations. We compare with the clamped coordinate quasiadiabatic approximation that introduces an effective potential for each asymmetric rotor level, with an effective rotation Hamiltonian constructed from ground state aver… ▽ More We analyze the nature of the rotational energy level structure of the OCS-He and OCS-H_2 complexes with a comparison of exact calculations to several differentdynamical approximations. We compare with the clamped coordinate quasiadiabatic approximation that introduces an effective potential for each asymmetric rotor level, with an effective rotation Hamiltonian constructed from ground state averages of the inverse of the inertial matrix, and investigate the usefulness of the Eckart condition to decouple rotations and vibrations of these weakly bound complexes between linear OCS and 4He or H_2. Comparison with exact results allows an assessment of the accuracies of the different approximate methods and indicates which approaches are suitable for larger clusters of OCS with 4He and with H_2. We find the OCS-H_2 complex is considerably more rigid than the OCS-4He complex, suggesting that semi-rigid models are useful for analysis of larger clusters of H_2 with OCS. △ Less

Submitted 17 September, 2003; originally announced September 2003.

Comments: accepted by Chem. Phys., 2003

OCS in small para-hydrogen clusters: energetics and structure with N=1-8 complexed hydrogen molecules

Authors: F. Paesani, R. E. Zillich, K. B. Whaley

Abstract: We determine the structure and energetics of complexes of the linear OCS molecule with small numbers of para-hydrogen molecules, N=1-8, using zero temperature quantum Monte Carlo methods. Ground state calculations are carried out with importance-sampled rigid body diffusion Monte Carlo (IS-RBDMC) and excited state calculations with the projection operator imaginary time spectral evolution (POITS… ▽ More We determine the structure and energetics of complexes of the linear OCS molecule with small numbers of para-hydrogen molecules, N=1-8, using zero temperature quantum Monte Carlo methods. Ground state calculations are carried out with importance-sampled rigid body diffusion Monte Carlo (IS-RBDMC) and excited state calculations with the projection operator imaginary time spectral evolution (POITSE) methodology. The ground states are found to be highly structured, with a gradual build up of two axial rings as N increases to 8. Analysis of the azimuthal density correlations around the OCS molecule shows that these rings are quite delocalized for small N values, but become strongly localized for N \geq 5 . Excited state calculations are made for a range of total cluster angular momentum values and the rotational energy levels fitted to obtain effective rotational and distortion constants of the complexed OCS molecule as a function of cluster size N. Detailed analysis of these spectroscopic constants indicates that the complexes of OCS with para-hydrogen have an unusually rich variation in dynamical behavior, with sizes N=1-2 showing near rigid behavior, sizes N=3-4 showing extremely floppy behavior, and the larger sizes N=5-8 showing more rigid behavior again. The large values of the distortion constant D obtained for N=3-4 are rationalized in terms of the coupling between the OCS rotations and the "breathing" mode of the first, partially filled ring of para-hydrogen molecules. △ Less

Submitted 16 September, 2003; originally announced September 2003.

Comments: 26 pages, 11 figures. accepted for publication in the Journal of Chemical Physics

Localized helium excitations in 4He_N-benzene clusters

Authors: Patrick Huang, K. Birgitta Whaley

Abstract: We compute ground and excited state properties of small helium clusters 4He_N containing a single benzene impurity molecule. Ground-state structures and energies are obtained for N=1,2,3,14 from importance-sampled, rigid-body diffusion Monte Carlo (DMC). Excited state energies due to helium vibrational motion near the molecule surface are evaluated using the projection operator, imaginary time s… ▽ More We compute ground and excited state properties of small helium clusters 4He_N containing a single benzene impurity molecule. Ground-state structures and energies are obtained for N=1,2,3,14 from importance-sampled, rigid-body diffusion Monte Carlo (DMC). Excited state energies due to helium vibrational motion near the molecule surface are evaluated using the projection operator, imaginary time spectral evolution (POITSE) method. We find excitation energies of up to ~23 K above the ground state. These states all possess vibrational character of helium atoms in a highly anisotropic potential due to the aromatic molecule, and can be categorized in terms of localized and collective vibrational modes. These results appear to provide precursors for a transition from localized to collective helium excitations at molecular nanosubstrates of increasing size. We discuss the implications of these results for analysis of anomalous spectral features in recent spectroscopic studies of large aromatic molecules in helium clusters. △ Less

Submitted 27 January, 2003; originally announced January 2003.

Comments: 15 pages, 5 figures, submitted to Phys. Rev. B