BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Artur F. Izmaylov

Artur F. Izmaylov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.chem-phquant-phphysics.comp-phcond-mat.mes-hallcond-mat.mtrl-sci

Catalog footprint

What is connected

29works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2026arXiv

Multistate iterative qubit coupled cluster (MS-iQCC): a quantum-inspired, state-averaged approach to ground- and excited-state energies

We introduce the multistate iterative qubit coupled cluster (MS-iQCC) method, a quantum-inspired algorithm that runs efficiently on classical hardware and is designed to predict both ground and excited electronic states of molecules. Accurate excited-state energetics are essential for interpreting spectroscopy and chemical reactivity, but standard electronic structure methods are either too computationally expensive for larger systems or lose reliability in the presence of strong electron correlation. MS-iQCC addresses this challenge by simultaneously optimizing multiple electronic states in a single, state-averaged procedure that treats ground and excited states on equal footing. This removes the energetic bias that is introduced when excited states are computed one at a time and constrained to remain orthogonal to previously optimized states. The approach supports multireference electronic structure by allowing multideterminantal initial guesses and by adaptively building a compact exponential ansatz from a pool of qubit excitation generators. We apply MS-iQCC to H$_4$, H$_2$O, N$_2$, and C$_2$, including strongly correlated geometries, and observe robust convergence of all targeted state energies to chemically meaningful accuracy across their potential energy surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Controlling energy conservation in quantum dynamics with independently moving basis functions: Application to Multi-Configuration Ehrenfest

Application of the time-dependent variational principle to a linear combination of frozen-width Gaussians describing the nuclear wavefunction provides a formalism where the total energy is conserved. The computational downside of this formalism is that trajectories of individual Gaussians are solutions of a coupled system of differential equations, limiting implementation to serial propagation algorithms. To allow for parallelization and acceleration of the computation, independent trajectories based on simplified equations of motion were suggested. Unfortunately, within practical realizations involving finite Gaussian bases, this simplification leads to breaking the energy conservation. We offer a solution for this problem by using Lagrange multipliers to ensure the energy and norm conservation regardless of basis function trajectories or basis completeness. We illustrate our approach within the Multi-Configuration Ehrenfest method considering a linear vibronic coupling model.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Deterministic improvements of quantum measurements with grouping of compatible operators, non-local transformations, and covariance estimates

Obtaining the expectation value of an observable on a quantum computer is a crucial step in the variational quantum algorithms. For complicated observables such as molecular electronic Hamiltonians, a common strategy is to present the observable as a linear combination of measurable fragments. The main problem of this approach is a large number of measurements required for accurate sampling of the observable's expectation value. We consider several partitioning schemes based on grouping of commuting multi-qubit Pauli products with the goal of minimizing the number of measurements. Three main directions are explored: 1) grouping commuting operators using the greedy approach, 2) involving non-local unitary transformations for measuring, and 3) taking advantage of compatibility of some Pauli products with several measurable groups. The last direction gives rise to a general framework that not only provides improvements over previous methods but also connects measurement grouping approaches with recent advances in techniques of shadow tomography. Following this direction, we develop two new measurement schemes that achieve a severalfold reduction in the number of measurements for a set of model molecules compared to previous state-of-the-art methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fidelity overhead for non-local measurements in variational quantum algorithms

Measuring quantum observables by grouping terms that can be rotated to sums of only products of Pauli $\hat z$ operators (Ising form) is proven to be efficient in near term quantum computing algorithms. This approach requires extra unitary transformations to rotate the state of interest so that the measurement of a fragment's Ising form would be equivalent to measurement of the fragment for the unrotated state. These extra rotations allow one to perform a fewer number of measurements by grouping more terms into the measurable fragments with a lower overall estimator variance. However, previous estimations of the number of measurements did not take into account non-unit fidelity of quantum gates implementing the additional transformations. Through a circuit fidelity reduction, additional transformations introduce extra uncertainty and increase the needed number of measurements. Here we consider a simple model for errors introduced by additional gates needed in schemes involving grouping of commuting Pauli products. For a set of molecular electronic Hamiltonians, we confirm that the numbers of measurements in schemes using non-local qubit rotations are still lower than those in their local qubit rotation counterparts, even after accounting for uncertainties introduced by additional gates.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fluid fermionic fragments for optimizing quantum measurements of electronic Hamiltonians in the variational quantum eigensolver

Measuring the expectation value of the molecular electronic Hamiltonian is one of the challenging parts of the variational quantum eigensolver. A widely used strategy is to express the Hamiltonian as a sum of measurable fragments using fermionic operator algebra. Such fragments have an advantage of conserving molecular symmetries that can be used for error mitigation. The number of measurements required to obtain the Hamiltonian expectation value is proportional to a sum of fragment variances. Here, we introduce a new method for lowering the fragments' variances by exploiting flexibility in the fragments' form. Due to idempotency of the occupation number operators, some parts of two-electron fragments can be turned into one-electron fragments, which then can be partially collected in a purely one-electron fragment. This repartitioning does not affect the expectation value of the Hamiltonian but has non-vanishing contributions to the variance of each fragment. The proposed method finds the optimal repartitioning by employing variances estimated using a classically efficient proxy for the quantum wavefunction. Numerical tests on several molecules show that repartitioning of one-electron terms lowers the number of measurements by more than an order of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Improving quantum measurements by introducing "ghost" Pauli products

Reducing the number of measurements required to estimate the expectation value of an observable is crucial for the variational quantum eigensolver to become competitive with state-of-the-art classical algorithms. To measure complicated observables such as a molecular electronic Hamiltonian, one of the common strategies is to partition the observable into linear combinations (fragments) of mutually commutative Pauli products. The total number of measurements for obtaining the expectation value is then proportional to the sum of variances of individual fragments. We propose a method that lowers individual fragment variances by modifying the fragments without changing the total observable expectation value. Our approach is based on adding Pauli products ("ghosts") that are compatible with members of multiple fragments. The total expectation value does not change because a sum of coefficients for each "ghost" Pauli product introduced to several fragments is zero. Yet, these additions change individual fragment variances because of the non-vanishing contributions of "ghost" Pauli products within each fragment. The proposed algorithm minimizes individual fragment variances using a classically efficient approximation of the quantum wavefunction for variance estimations. Numerical tests on a few molecular electronic Hamiltonian expectation values show several-fold reductions in the number of measurements in the "ghost" Pauli algorithm compared to those in the other recently developed techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum Kinetic Rates within the Nonequilibrium Steady State

The nonequilibrium steady state (NESS) of a quantum network is central to a host of physical and biological scenarios. Examples include natural processes such as vision and photosynthesis, as well as technical devices such as photocells, both activated by incoherent light (e.g. sunlight) and leading to quantum transport. Here, a completely general approach to defining components of a quantum network in the NESS, and obtaining rates of processes between these components is provided. Quantum effects are explicitly included throughout, both in (a) defining network components via projection operators, and (b) in determining the role of coherences in rate processes. As examples, the methodology is applied to model cases, two versions of the V-level system, and to the spin-boson model, wherein the role of the environment and of internal system properties in determining the rates is examined. In addition, the role of Markovian vs. non-Markovian contributions is quantified, exposing conditions under which NESS rates can be obtained by perturbing the nonequilibrium steady state.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Analytic gradients in variational quantum algorithms: Algebraic extensions of the parameter-shift rule to general unitary transformations

Optimization of unitary transformations in Variational Quantum Algorithms benefits highly from efficient evaluation of cost function gradients with respect to amplitudes of unitary generators. We propose several extensions of the parametric-shift-rule to formulating these gradients as linear combinations of expectation values for generators with general eigen-spectrum (i.e. with more than two eigenvalues). Our approaches are exact and do not use any auxiliary qubits, instead they rely on a generator eigen-spectrum analysis. Two main directions in the parametric-shift-rule extensions are 1) polynomial expansion of the exponential unitary operator based on a limited number of different eigenvalues in the generator and 2) decomposition of the generator as a linear combination of low-eigenvalue operators (e.g. operators with only 2 or 3 eigenvalues). These techniques have a range of scalings for the number of needed expectation values with the number of generator eigenvalues from quadratic (for polynomial expansion) to linear and even $\log_2$ (for generator decompositions). This allowed us to propose efficient differentiation schemes superior to previous approaches for commonly used 2-qubit transformations (e.g. match-gates, transmon and fSim gates) and $\hat S^2$-conserving fermionic operators for the variational quantum eigensolver.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Tequila: A platform for rapid development of quantum algorithms

Variational quantum algorithms are currently the most promising class of algorithms for deployment on near-term quantum computers. In contrast to classical algorithms, there are almost no standardized methods in quantum algorithmic development yet, and the field continues to evolve rapidly. As in classical computing, heuristics play a crucial role in the development of new quantum algorithms, resulting in high demand for flexible and reliable ways to implement, test, and share new ideas. Inspired by this demand, we introduce tequila, a development package for quantum algorithms in python, designed for fast and flexible implementation, prototyping, and deployment of novel quantum algorithms in electronic structure and other fields. Tequila operates with abstract expectation values which can be combined, transformed, differentiated, and optimized. On evaluation, the abstract data structures are compiled to run on state-of-the-art quantum simulators or interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measurement Optimization in the Variational Quantum Eigensolver Using a Minimum Clique Cover

Solving the electronic structure problem using the Variational Quantum Eigensolver (VQE) technique involves measurement of the Hamiltonian expectation value. Current hardware can perform only projective single-qubit measurements, and thus, the Hamiltonian expectation value is obtained by measuring parts of the Hamiltonian rather than the full Hamiltonian. This restriction makes the measurement process inefficient because the number of terms in the Hamiltonian grows as $O(N^4)$ with the size of the system, $N$. To optimize VQE measurement one can try to group as many Hamiltonian terms as possible for their simultaneous measurement. Single-qubit measurements allow one to group only the terms that commute within corresponding single-qubit subspaces or qubit-wise commuting. We found that qubit-wise commutativity between the Hamiltonian terms can be expressed as a graph and the problem of the optimal grouping is equivalent of finding a minimum clique cover (MCC) for the Hamiltonian graph. The MCC problem is NP-hard but there exist several polynomial heuristic algorithms to solve it approximately. Several of these heuristics were tested in this work for a set of molecular electronic Hamiltonians. On average, grouping qubit-wise commuting terms reduced the number of operators to measure three times compared to the total number of terms in the considered Hamiltonians.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring all compatible operators in one series of a single-qubit measurements using unitary transformations

The Variational Quantum Eigensolver approach to the electronic structure problem on a quantum computer involves measurement of the Hamiltonian expectation value. Formally, quantum mechanics allows one to measure all mutually commuting or compatible operators simultaneously. Unfortunately, the current hardware permits measuring only a much more limited subset of operators that share a common tensor product eigen-basis. We introduce unitary transformations that transform any fully commuting group of operators to a group that can be measured on current hardware. These unitary operations can be encoded as a sequence of Clifford gates and let us not only measure much larger groups of terms but also to obtain these groups efficiently on a classical computer. The problem of finding the minimum number of fully commuting groups of terms covering the whole Hamiltonian is found to be equivalent to the minimum clique cover problem for a graph representing Hamiltonian terms as vertices and commutativity between them as edges. Tested on few molecular electronic Hamiltonians with up to 50 thousand terms, the introduced technique allows for the reduction of the number of separately measurable operator groups down to few hundreds, thus achieving up to 2 orders of magnitude reduction. It was also shown that this gain scales at least linearly with the number of qubits.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Exact and approximate symmetry projectors for the electronic structure problem on a quantum computer

Solving the electronic structure problem on a universal-gate quantum computer within the variational quantum eigensolver (VQE) methodology requires constraining the search procedure to a subspace defined by relevant physical symmetries. Ignoring symmetries results in convergence to the lowest eigenstate of the Fock space for the second quantized electronic Hamiltonian. Moreover, this eigenstate can be symmetry broken due to limitations of the wavefunction ansatz. To address this VQE problem, we introduce and assess methods of exact and approximate projection operators to irreducible eigen-subspaces of available physical symmetries. Feasibility of symmetry projection operators in the VQE framework is discussed, and their efficiency is compared with symmetry constraint optimization procedures. Generally, projectors introduce higher numbers of terms for VQE measurement compared to the constraint approach. On the other hand, the projection formalism improves accuracy of the variational wavefunction ansatz without introducing additional unitary transformations, which is beneficial for reducing depths of quantum circuits.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A diabatic definition of geometric phase effects

Electronic wave-functions in the adiabatic representation acquire nontrivial geometric phases (GPs) when corresponding potential energy surfaces undergo conical intersection (CI). These GPs have profound effects on the nuclear quantum dynamics and cannot be eliminated in the adiabatic representation without changing the physics of the system. To define dynamical effects arising from the GP presence the nuclear quantum dynamics of the CI containing system is compared with that of the system with artificially removed GP. We explore a new construction of the system with removed GP via a modification of the diabatic representation for the original CI containing system. Using an absolute value function of diabatic couplings we remove the GP while preserving adiabatic potential energy surfaces and CI. We assess GP effects in dynamics of a two-dimensional linear vibronic coupling model both for ground and excited state dynamics. Results are compared with those obtained with a conventional removal of the GP by ignoring double-valued boundary conditions of the real electronic wave-functions. Interestingly, GP effects appear similar in two approaches only for the low energy dynamics. In contrast with the conventional approach, a new approach does not have substantial GP effects in the ultra-fast excited state dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Entanglement in the Born-Oppenheimer Approximation

The role of electron-nuclear entanglement on the validity of the Born-Oppenheimer (BO) approximation is investigated. While nonadiabatic couplings generally lead to entanglement and to a failure of the BO approximation, surprisingly the degree of electron-nuclear entanglement is found to be uncorrelated with the degree of validity of the BO approximation. This is because while the degree of entanglement of BO states is determined by their deviation from the corresponding states in the crude BO approximation, the accuracy of the BO approximation is dictated, instead, by the deviation of the BO states from the exact electron-nuclear states. In fact, in the context of a minimal avoided crossing model, extreme cases are identified where an adequate BO state is seen to be maximally entangled, and where the BO approximation fails but the associated BO state remains approximately unentangled. Further, the BO states are found to not preserve the entanglement properties of the exact electron-nuclear eigenstates, and to be completely unentangled only in the limit in which the BO approximation becomes exact.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nuclear dynamics investigation of the initial electron transfer in the cyclobutane pyrimidine dimer lesion repair process by photolyases

Photolyases are proteins capable of harvesting the sunlight to repair DNA damages caused by UV light. In this work we focus on the first step in the repair process of the cyclobutane pyrimidine dimer photoproduct (CPD) lesion, which is an electron transfer (ET) from a flavine cofactor to CPD, and study the role of various nuclear degrees of freedom (DOF) in this step. The ET step has been experimentally studied using transient spectroscopy and the corresponding data provide excellent basis for testing the quality of quantum dynamical models. Based on previous theoretical studies of electronic structure and conformations of the protein active site, we present a procedure to build a diabatic Hamiltonian for simulating the ET reaction in a molecular complex mimicking the enzyme's active site. We generate a reduced nuclear dimensional model that provides a first non-empirical quantum dynamical description of the structural features influencing the ET rate. By varying the nuclear DOF parametrization in the model to assess the role of different nuclear motions, we demonstrate that the low frequency flavin butterfly bending mode slows ET by reducing Franck-Condon overlaps between donor and acceptor states and also induces decoherence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the inclusion of the diagonal Born-Oppenheimer correction in surface hopping methods

The diagonal Born-Oppenheimer correction (DBOC) stems from the diagonal second derivative coupling term in the adiabatic representation, and it can have an arbitrary large magnitude when a gap between neighbouring Born-Oppenheimer (BO) potential energy surfaces (PESs) is closing. Nevertheless, DBOC is typically neglected in mixed quantum-classical methods of simulating nonadiabatic dynamics (e.g., fewest-switch surface hopping (FSSH) method). A straightforward addition of DBOC to BO PESs in the FSSH method, FSSH+D, has been shown to lead to numerically much inferior results for models containing conical intersections. More sophisticated variation of the DBOC inclusion, phase-space surface-hopping (PSSH) was more successful than FSSH+D but on model problems without conical intersections. This work comprehensively assesses the role of DBOC in nonadiabatic dynamics of two electronic state problems and the performance of FSSH, FSSH+D, and PSSH methods in variety of one- and two-dimensional models. Our results show that the inclusion of DBOC can enhance the accuracy of surface hopping simulations when two conditions are simultaneously satisfied: 1) nuclei have kinetic energy lower than DBOC and 2) PESs are not strongly nonadiabatically coupled. The inclusion of DBOC is detrimental in situations where its energy scale becomes very high or even diverges, because in these regions PESs are also very strongly coupled. In this case, the true quantum formalism heavily relies on an interplay between diagonal and off-diagonal nonadiabatic couplings while surface hopping approaches treat diagonal terms as PESs and off-diagonal ones stochastically.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Topologically correct quantum nonadiabatic formalism for on-the-fly dynamics

On-the-fly quantum nonadiabatic dynamics for large systems greatly benefits from the adiabatic representation readily available from the electronic structure programs. However, frequently occurring in this representation conical intersections introduce non-trivial geometric or Berry phases which require a special treatment for adequate modelling of the nuclear dynamics. We analyze two approaches for nonadiabatic dynamics using the time-dependent variational principle and the adiabatic representation. The first approach employs adiabatic electronic functions with global parametric dependence on the nuclear coordinates. The second approach uses adiabatic electronic functions obtained only at the centres of moving localized nuclear basis functions (e.g. frozen-width Gaussians). Unless a gauge transformation is used to enforce single-valued boundary conditions, the first approach fails to capture the geometric phase. In contrast, the second approach accounts for the geometric phase naturally because of the absence of the global nuclear coordinate dependence in the electronic functions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast numerical evaluation of time-derivative nonadiabatic couplings for mixed quantum-classical methods

We have developed a numerical differentiation scheme which eliminates evaluation of overlap determinants in calculating the time-derivative non-adiabatic couplings (TDNACs). Evaluation of these determinants was the bottleneck in previous implementations of mixed quantum-classical methods using numerical differentiation of electronic wave functions in the Slater-determinant representation. The central idea of our approach is, first, to reduce the analytic time derivatives of Slater determinants to time derivatives of molecular orbitals, and then to apply a finite-difference formula. Benchmark calculations prove the efficiency of the proposed scheme showing impressive several-order-of-magnitude speedups of the TDNAC calculation step for midsize molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Problem-free time-dependent variational principle for open quantum systems

Methods of quantum nuclear wave-function dynamics have become very efficient in simulating large isolated systems using the time-dependent variational principle (TDVP). However, a straightforward extension of the TDVP to the density matrix framework gives rise to methods that do not conserve the energy in the isolated system limit and the total system population for open systems where only energy exchange with the environment is allowed. These problems arise when the system density is in a mixed state and is simulated using an incomplete basis. Thus, the basis set incompleteness, which is inevitable in practical calculations, creates artificial channels for energy and population dissipation. To overcome this unphysical behavior, we have introduced a constrained Lagrangian formulation of TDVP applied to the non-stochastic open system Schrodinger equation (NOSSE) [L. Joubert-Doriol, I. G. Ryabinkin, and A. F. Izmaylov, J. Chem. Phys. 141, 234112 (2014)]. While our formulation can be applied to any variational ansatz for the system density matrix, derivation of working equations and numerical assessment are done within the variational multiconfiguration Gaussian approach for a two-dimensional linear vibronic coupling model system interacting with a harmonic bath.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Why do mixed quantum-classical methods describe short-time dynamics through conical intersections so well? Analysis of geometric phase effects

Adequate simulation of non-adiabatic dynamics through conical intersection requires account for a non-trivial geometric phase (GP) emerging in electronic and nuclear wave-functions in the adiabatic representation. Popular mixed quantum-classical (MQC) methods, surface hopping and Ehrenfest, do not carry a nuclear wave-function to be able to incorporate the GP into nuclear dynamics. Surprisingly, the MQC methods reproduce ultra-fast interstate crossing dynamics generated with the exact quantum propagation so well as if they contained information about the GP. Using two-dimensional linear vibronic coupling models we unravel how the MQC methods can effectively mimic the most significant dynamical GP effects: 1) compensation for repulsive diagonal second order non-adiabatic couplings and 2) transfer enhancement for a fully cylindrically symmetric component of a nuclear distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A perturbative formalism for electronic transitions through conical intersections in a fully quadratic vibronic model

We consider a fully quadratic vibronic model Hamiltonian for studying photoinduced electronic transitions through conical intersections. Using a second order perturbative approximation for diabatic couplings we derive an analytical expression for the time evolution of electronic populations at a given temperature. This formalism extends upon a previously developed perturbative technique for a linear vibronic coupling Hamiltonian. The advantage of the quadratic model Hamiltonian is that it allows one to use separate quadratic representations for potential energy surfaces of different electronic states and a more flexible representation of interstate couplings. We explore features introduced by the quadratic Hamiltonian in a series of 2D models, and then apply our formalism to the 2,6-bis(methylene) adamantyl cation, and its dimethyl derivative. The Hamiltonian parameters for the molecular systems have been obtained from electronic structure calculations followed by a diabatization procedure. The evolution of electronic populations in the molecular systems using the perturbative formalism shows a good agreement with that from variational quantum dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An efficient implementation of the localized operator partitioning method for electronic energy transfer

The localized operator partitioning method [Y. Khan and P. Brumer, J. Chem. Phys. 137, 194112 (2012)] rigorously defines the electronic energy on any subsystem within a molecule and gives a precise meaning to the subsystem ground and excited electronic energies, which is crucial for investigating electronic energy transfer from first principles. However, an efficient implementation of this approach has been hindered by complicated one- and two-electron integrals arising in its formulation. Using a resolution of the identity in the definition of partitioning we reformulate the method in a computationally efficient manner that involves standard one- and two-electron integrals. We apply the developed algorithm to the 9-((1-naphthyl)-methyl)-anthracene (A1N) molecule by partitioning A1N into anthracenyl and CH2-naphthyl groups as subsystems, and examine their electronic energies and populations for several excited states using Configuration Interaction Singles method. The implemented approach shows a wide variety of different behaviors amongst the excited electronic states.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Local operator partitioning of electronic energy for electronic energy transfer: An efficient algorithm

An efficient computational algorithm to implement a local operator approach to partitioning electronic energy in general molecular systems is presented. This approach, which rigorously defines the electronic energy on any subsystem within a molecule, gives a precise meaning to the subsystem ground and excited electronic energies, which is crucial for investigating electronic energy transfer from first principles. We apply the technique to the $9-$(($1-$naphthyl)$-$methyl)-anthracene (A1N) molecule by partitioning A1N into anthracenyl and CH$_2-$naphthyl groups as subsystems, and examine their electronic energies and populations for several excited states using Configuration Interaction Singles method. The implemented approach shows a wide variety of different behaviors amongst these excited electronic states.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Non-stochastic matrix Schrödinger equation for open systems

We propose an extension of the Schrödinger equation for a quantum system interacting with environment. This equation describes dynamics of auxiliary wave-functions $\mathbf{m}$, from which the system density matrix can be reconstructed as $\hatρ = \mathbf{m} \mathbf{m}^\dagger$. We formulate a compatibility condition, which ensures that the reconstructed density satisfies a given quantum master equation for the system density. The resulting non-stochastic evolution equation preserves positive-definiteness of the system density and is applicable to both Markovian and non-Markovian system-bath treatments. Our formalism also resolves a long-standing problem of energy non-conservation in the time-dependent variational principle applied to mixed states of closed systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

When do we need to account for the geometric phase in excited state dynamics ?

We investigate the role of the geometric phase (GP) in an internal conversion process when the system changes its electronic state by passing through a conical intersection (CI). Local analysis of a two-dimensional linear vibronic coupling (LVC) model Hamiltonian near the CI shows that the role of the GP is twofold. First, it compensates for a repulsion created by the so-called diagonal Born-Oppenheimer correction (DBOC). Second, the GP enhances the non-adiabatic transition probability for a wave-packet part that experiences a central collision with the CI. To assess the significance of both GP contributions we propose two indicators that can be computed from parameters of electronic surfaces and initial conditions. To generalize our analysis to N-dimensional systems we introduce a reduction of a general N-dimensional LVC model to an effective 2D LVC model using a mode transformation that preserves short-time dynamics of the original N-dimensional model. Using examples of the bis(methylene) adamantyl and butatriene cations, and the pyrazine molecule we have demonstrated that their effective 2D models reproduce the short-time dynamics of the corresponding full dimensional models, and the introduced indicators are very reliable in assessing GP effects.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Analysis of geometric phase effects in the quantum-classical Liouville formalism

We analyze two approaches to the quantum-classical Liouville (QCL) formalism that differ in the order of two operations: Wigner transformation and projection onto adiabatic electronic states. The analysis is carried out on a two-dimensional linear vibronic model where geometric phase (GP) effects arising from a conical intersection profoundly affect nuclear dynamics. We find that the Wigner-then-Adiabatic (WA) QCL approach captures GP effects, whereas the Adiabatic-then-Wigner (AW) QCL approach does not. Moreover, the Wigner transform in AW-QCL leads to an ill-defined Fourier transform of double-valued functions. The double-valued character of these functions stems from the nontrivial GP of adiabatic electronic states in the presence of a conical intersection. In contrast, WA-QCL avoids this issue by starting with the Wigner transform of single-valued quantities of the full problem. Since the WA-QCL approach uses solely the adiabatic potentials and non-adiabatic derivative couplings as an input, our results indicate that WA-QCL can capture GP effects in general two-state crossing problems using first-principles electronic structure calculations without prior diabatization or introduction of explicit phase factors.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Geometric phase effects in dynamics near conical intersections: Symmetry breaking and spatial localization

We show that finite systems with conical intersections can exhibit spontaneous symmetry breaking which manifests itself in spatial localization of eigenstates. This localization has a geometric phase origin and is robust against variation of model parameters. The transition between localized and delocalized eigenstate regimes resembles a continuous phase transition. The localization slows down the low-energy quantum nuclear dynamics at zero and low temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Geometric phase effects in low-energy dynamics near conical intersections: A study of the multidimensional linear vibronic coupling model

In molecular systems containing conical intersections (CIs), a nontrivial geometric phase (GP) appears in the nuclear and electronic wave-functions in the adiabatic representation. We study GP effects in nuclear dynamics of an N-dimensional linear vibronic coupling (LVC) model. The main impact of GP on low-energy nuclear dynamics is reduction of population transfer between the local minima of the LVC lower energy surface. For the LVC model, we proposed an isometric coordinate transformation that confines non-adiabatic effects within a two-dimensional subsystem interacting with an N-2 dimensional environment. Since environmental modes do not couple electronic states, all GP effects originate from nuclear dynamics within the subsystem. We explored when the GP affects nuclear dynamics of the isolated subsystem, and how the subsystem-environment interaction can interfere with GP effects. Comparing quantum dynamics with and without GP allowed us to devise simple rules to determine significance of the GP for nuclear dynamics in this model.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Analytical infrared intensities for periodic systems with local basis sets

We report a method for the efficient evaluation of analytic infrared (IR) intensities within generalized Kohn-Sham density functional theory using Gaussian orbitals and periodic boundary conditions. A discretized form of the Berry phase is used to evaluate a periodic dipole moment and its derivatives with respect to in-phase nuclear coordinate displacements. Benchmark calculations are presented for one-dimensional chains of water molecules and poly(paraphenylenevinylene).

0 citations0 reviewsNo code linkedOpen access