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

Ilya G. Ryabinkin

Ilya G. Ryabinkin 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

14works
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

14 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

On the Importance of Well-Defined Thermal Correlation Functions in Simulating Vibronic Spectra

Two difficulties associated with the computations of thermal vibrational correlation functions are discussed. The first one is the lack of a well-behaved expression that is valid at both high-temperature and $T \to 0$ K limits. Specifically, if the partition function and the propagator are considered separately, then thermal vibrational correlation functions may have an indeterminate form 0/0 in the limit $T \to 0$ K. This difficulty is resolved when the partition function and the propagator are jointly considered in the harmonic approximation, which allows a problematic term that emanates from the zero-point energy to be cancelled out thereby producing a thermal correlation function with a determinate form in $T \to 0$ K limit. The second difficulty is related to the multivaluedness of the vibrational correlation function. We show numerically that an improper selection of branch leads to discontinuities in the computed correlation function and an incorrect vibronic spectra. We propose a phase tracking procedure that ensures continuity of both real and imaginary parts of the correlation function to recover the correct spectra. We support our findings by simulating the UV-vis absorption spectra of pentacene at 4 K and benzene at 298 K. Both are found to be in good agreement with their experimental counterparts.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Estimating Phosphorescent Emission Energies in Ir(III) Complexes using Large-Scale Quantum Computing Simulations

Quantum chemistry simulations that accurately predict the properties of materials are among the most highly anticipated applications of quantum computing. It is widely believed that simulations running on quantum computers will allow for higher accuracy, but there has not yet been a convincing demonstration that quantum methods are competitive with existing classical methods at scale. Here we apply the iterative qubit coupled cluster (iQCC) method on classical hardware to the calculation of the $T_1 \to S_0$ transition energies in nine phosphorescent iridium complexes, to determine if quantum simulations have any advantage over traditional computing methods. Phosphorescent iridium complexes are integral to the widespread commercialization of organic light-emitting diode (OLED) technology, yet accurate computational prediction of their emission energies remains a challenge. Our simulations would require a gate-based quantum computer with a minimum of 72 fully-connected and error-corrected logical qubits. Since such devices do not yet exist, we demonstrate the iQCC quantum method using a special purpose quantum simulator on classical hardware. The results are compared to a selection of common density-functional theory (DFT) functionals (B3LYP, CAM-B3LYP, LC-wHPBE), ab initio methods (HF and MP2), and experimental data. The iQCC quantum method is found to match the accuracy of the fine-tuned DFT functionals, has a better Pearson correlation coefficient, and still has considerable potential for systematic improvement. Based on these results, we anticipate that the iQCC quantum method will have the required accuracy to design organometallic complexes when deployed on emerging quantum hardware.

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

Reduction of Electronic Wavefunctions to Kohn-Sham Effective Potentials

A method for calculating the Kohn--Sham exchange-correlation potential, $v_\text{XC}(\mathbf{r})$, from a given electronic wavefunction is devised and implemented. It requires on input one- and two-electron density matrices and involves construction of the generalized Fock matrix. The method is free from numerical limitations and basis-set artifacts of conventional schemes for constructing $v_\text{XC}(\mathbf{r})$ in which the potential is recovered from a given electron density, and is simpler than various many-body techniques. The chief significance of this development is that it allows one to directly probe the functional derivative of the true exchange-correlation energy functional and to rigorously test and improve various density-functional approximations.

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

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

Accurate and efficient approximation to the optimized effective potential for exchange

We devise an efficient practical method for computing the Kohn-Sham exchange-correlation potential corresponding to a Hartree-Fock electron density. This potential is almost indistinguishable from the exact-exchange optimized effective potential (OEP) and, when used as an approximation to the OEP, is vastly better than all existing models. Using our method one can obtain unambiguous, nearly exact OEPs for any finite one-electron basis set at the same low cost as the Krieger-Li-Iafrate and Becke-Johnson potentials. For all practical purposes, this solves the long-standing problem of black-box construction of OEPs in exact-exchange calculations.

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