Researcher profile

Jiri Vanicek

Jiri Vanicek contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

15 published item(s)

preprint2020arXiv

A combined on-the-fly/interpolation procedure for evaluating energy values needed in molecular simulations

We propose an algorithm for molecular dynamics or Monte Carlo simulations that uses an interpolation procedure to estimate potential energy values from energies and gradients evaluated previously at points of a simplicial mesh. We chose an interpolation procedure which is exact for harmonic systems and considered two possible mesh types: Delaunay triangulation and an alternative anisotropic triangulation designed to improve performance in anharmonic systems. The mesh is generated and updated on the fly during the simulation. The procedure is tested on two-dimensional quartic oscillators and on the path integral Monte Carlo evaluation of HCN/DCN equilibrium isotope effect.

preprint2016arXiv

Path integral approach to the quantum fidelity amplitude

The Loschmidt echo is a measure of quantum irreversibility and is determined by the fidelity amplitude of an imperfect time-reversal protocol. Fidelity amplitude plays an important role both in the foundations of quantum mechanics and its applications, such as time-resolved electronic spectroscopy. We derive an exact path integral formula for the fidelity amplitude and use it to obtain a series of increasingly accurate semiclassical approximations by truncating an exact expansion of the path integral exponent. While the zeroth-order expansion results in a remarkably simple, yet nontrivial approximation for the fidelity amplitude, the first-order expansion yields an alternative derivation of the so-called "dephasing representation", circumventing the use of semiclassical propagator as in the original derivation. We also obtain approximate expression for fidelity based on the second-order expansion, which resolves several shortcomings of the dephasing representation. The rigorous derivation from the path integral permits the identification of sufficient conditions under which various approximations obtained become exact.

preprint2015arXiv

Spectral properties of electrons in fractal nanowires

In view of promising applications of fractal nanostructures, we analyze the spectra of quantum particles in the Sierpinski carpet and study the non-correlated electron gas in this geometry. We show that the spectrum exhibits scale invariance with almost arbitrary spacing between energy levels, including large energy gaps at high energies. These features disappear in the analogous random fractal---where Anderson localization dominates---and in the regular lattice of equally sized holes---where only two length scales are present. The fractal structure amplifies microscopic effects, resulting in the presence of quantum behavior of the electron gas even at high temperatures. Our results demonstrate the potential of fractal nanostructures to improve the light-matter interaction at any frequency, with possible applications, e.g., in the development of solar cells with a wide absorption spectrum, artificial photosynthesis, or nanometamaterials with tailored Fermi levels and band gaps, operating in a wide range of frequencies, and with extended operating temperature range.

preprint2014arXiv

On-the-fly ab initio semiclassical dynamics: Identifying degrees of freedom essential for emission spectra of oligothiophenes

Vibrationally resolved spectra provide a stringent test of the accuracy of theoretical calculations. We combine the thawed Gaussian approximation (TGA) with an on-the-fly ab initio (OTF-AI) scheme to calculate the vibrationally resolved emission spectra of oligothiophenes with up to five rings. The efficiency of the OTF-AI-TGA permits treating all vibrational degrees of freedom on an equal footing even in pentathiophene with 105 vibrational degrees of freedom, thus obviating the need for the global harmonic approximation, popular for large systems. Besides reproducing almost perfectly the experimental emission spectra, in order to provide a deeper insight into the associated physical and chemical processes, we also develop a novel systematic approach to assess the importance and coupling between individual vibrational degrees of freedom during the dynamics. This allows us to explain how the vibrational line shapes of the oligothiophenes change with increasing number of rings. Furthermore, we observe the dynamical interplay between the quinoid and aromatic characters of individual rings in the oligothiophene chain during the dynamics and confirm that the quinoid character prevails in the center of the chain.

preprint2013arXiv

Improving the accuracy and efficiency of time-resolved electronic spectra calculations: Cellular dephasing representation with a prefactor

Time-resolved electronic spectra can be obtained as the Fourier transform of a special type of time correlation function known as fidelity amplitude, which, in turn, can be evaluated approximately and efficiently with the dephasing representation. Here we improve both the accuracy of this approximation---with an amplitude correction derived from the phase-space propagator---and its efficiency---with an improved cellular scheme employing inverse Weierstrass transform and optimal scaling of the cell size. We demonstrate the advantages of the new methodology by computing dispersed time-resolved stimulated emission spectra in the harmonic potential, pyrazine, and the NCO molecule. In contrast, we show that in strongly chaotic systems such as the quartic oscillator the original dephasing representation is more appropriate than either the cellular or prefactor-corrected methods.

preprint2012arXiv

Evaluation of the importance of spin-orbit couplings in the nonadiabatic quantum dynamics with quantum fidelity and with its efficient "on-the-fly" ab initio semiclassical approximation

We propose to measure the importance of spin-orbit couplings (SOCs) in the nonadiabatic molecular quantum dynamics rigorously with quantum fidelity. To make the criterion practical, quantum fidelity is estimated efficiently with the multiple-surface dephasing representation (MSDR). The MSDR is a semiclassical method that includes nuclear quantum effects through interference of mixed quantum-classical trajectories without the need for the Hessian of potential energy surfaces. Two variants of the MSDR are studied, in which the nuclei are propagated either with the fewest-switches surface hopping or with the locally mean field dynamics. The fidelity criterion and MSDR are first tested on one-dimensional model systems amenable to numerically exact quantum dynamics. Then, the MSDR is combined with "on-the-fly" computed electronic structure to measure the importance of SOCs and nonadiabatic couplings (NACs) in the photoisomerization dynamics of CH2NH2+ considering 20 electronic states and in the collision of F + H2 considering six electronic states.

preprint2012arXiv

Role of the sampling weight in evaluating classical time autocorrelation functions

We analyze how the choice of the sampling weight affects the efficiency of the Monte Carlo evaluation of classical time autocorrelation functions. Assuming uncorrelated sampling or sampling with constant correlation length, we propose a sampling weight for which the number of trajectories needed for convergence is independent of the correlated quantity, dimensionality, dynamics, and phase-space density. In contrast, it is shown that the computational cost of the "standard" intuitive algorithm which samples directly from the phase-space density may scale exponentially with the number of degrees of freedom. Yet, for the stationary Gaussian distribution of harmonic systems and for the autocorrelation function of a linear function of phase-space coordinates, the computational cost of this standard algorithm is also independent of dimensionality.

preprint2011arXiv

Beating the efficiency of both quantum and classical simulations with semiclassics

While rigorous quantum dynamical simulations of many-body systems are extremely difficult (or impossible) due to the exponential scaling with dimensionality, corresponding classical simulations completely ignore quantum effects. Semiclassical methods are generally more efficient but less accurate than quantum methods, and more accurate but less efficient than classical methods. We find a remarkable exception to this rule by showing that a semiclassical method can be both more accurate and faster than a classical simulation. Specifically, we prove that for the semiclassical dephasing representation the number of trajectories needed to simulate quantum fidelity is independent of dimensionality and also that this semiclassical method is even faster than the most efficient corresponding classical algorithm. Analytical results are confirmed with simulations of quantum fidelity in up to 100 dimensions with 2^1700-dimensional Hilbert space.

preprint2011arXiv

Breaking the exponential wall in classical simulations of fidelity

We analyze the efficiency of available algorithms for the simulation of classical fidelity and show that their computational costs increase exponentially with the number of degrees of freedom for almost all initial states. Then we present an algorithm whose cost is independent of the system's dimensionality and show that, within a continuous family of algorithms, our algorithm is the only one with this property. Simultaneously we propose a general analytical approach to estimate efficiency of trajectory-based methods.

preprint2011arXiv

Measuring nonadiabaticity of molecular quantum dynamics with quantum fidelity and with its efficient semiclassical approximation

We propose to measure nonadiabaticity of molecular quantum dynamics rigorously with the quantum fidelity between the Born-Oppenheimer and fully nonadiabatic dynamics. It is shown that this measure of nonadiabaticity applies in situations where other criteria, such as the energy gap criterion or the extent of population transfer, fail. We further propose to estimate this quantum fidelity efficiently with a generalization of the dephasing representation to multiple surfaces. Two variants of the multiple-surface dephasing representation (MSDR) are introduced, in which the nuclei are propagated either with the fewest-switches surface hopping (FSSH) or with the locally mean field dynamics (LMFD). The LMFD can be interpreted as the Ehrenfest dynamics of an ensemble of nuclear trajectories, and has been used previously in the nonadiabatic semiclassical initial value representation. In addition to propagating an ensemble of classical trajectories, the MSDR requires evaluating nonadiabatic couplings and solving the Schrödinger (or more generally, the quantum Liouville-von Neumann) equation for a single discrete degree of freedom. The MSDR can be also used to measure the importance of other terms present in the molecular Hamiltonian, such as diabatic couplings, spin-orbit couplings, or couplings to external fields, and to evaluate the accuracy of quantum dynamics with an approximate nonadiabatic Hamiltonian. The method is tested on three model problems introduced by Tully, on a two-surface model of dissociation of NaI, and a three-surface model including spin-orbit interactions. An example is presented that demonstrates the importance of often-neglected second-order nonadiabatic couplings.

preprint2010arXiv

Direct evaluation of the temperature dependence of the rate constant based on the quantum instanton approximation

A general method for the direct evaluation of the temperature dependence of the quantum-mechanical reaction rate constant in many-dimensional systems is described. The method is based on the quantum instanton approximation for the rate constant, thermodynamic integration with respect to the inverse temperature, and the path integral Monte Carlo evaluation. It can describe deviations from the Arrhenius law due to the coupling of rotations and vibrations, zero-point energy, tunneling, corner-cutting, and other nuclear quantum effects. The method is tested on the Eckart barrier and the full-dimensional H + H_2 -> H_2 + H reaction. In the temperature range from 300K to 1500K, the error of the present method remains within 13% despite the very large deviations from the Arrhenius law. The direct approach makes the calculations much more efficient, and the efficiency is increased even further (by up to two orders of magnitude in the studied reactions) by using optimal estimators for reactant and transition state thermal energies. Which of the estimators is optimal, however, depends on the system and the strength of constraint in a constrained simulation.

preprint2010arXiv

Evaluation of the nondiabaticity of quantum molecular dynamics with the dephasing representation of quantum fidelity

We propose an approximate method for evaluating the importance of non-Born-Oppenheimer effects on the quantum dynamics of nuclei. The method uses a generalization of the dephasing representation (DR) of quantum fidelity to several diabatic potential energy surfaces and its computational cost is the cost of dynamics of a classical phase space distribution. It can be implemented easily into any molecular dynamics program and also can utilize on-the-fly ab initio electronic structure information. We test the methodology on three model problems introduced by Tully and on the photodissociation of NaI. The results show that for dynamics close to the diabatic limit the decay of fidelity due to nondiabatic effects is described accurately by the DR. In this regime, unlike the mixed quantum-classical methods such as surface hopping or Ehrenfest dynamics, the DR can capture more subtle quantum effects than the population transfer between potential energy surfaces. Hence we propose using the DR to estimate the dynamical importance of diabatic, spin-orbit, or other couplings between potential energy surfaces. The acquired information can help reduce the complexity of a studied system without affecting the accuracy of the quantum simulation.

preprint2010arXiv

Three applications of path integrals: equilibrium and kinetic isotope effects, and the temperature dependence of the rate constant of the [1,5] sigmatropic hydrogen shift in (Z)-1,3-pentadiene

Recent experiments have confirmed the importance of nuclear quantum effects even in large biomolecules at physiological temperature. Here we describe how the path integral formalism can be used to describe rigorously the nuclear quantum effects on equilibrium and kinetic properties of molecules. Specifically, we explain how path integrals can be employed to evaluate the equilibrium (EIE) and kinetic (KIE) isotope effects, and the temperature dependence of the rate constant. The methodology is applied to the [1,5] sigmatropic hydrogen shift in pentadiene. Both the KIE and the temperature dependence of the rate constant confirm the importance of tunneling and other nuclear quantum effects as well as of the anharmonicity of the potential energy surface. Moreover, previous results on the KIE were improved by using a combination of a high level electronic structure calculation within the harmonic approximation with a path integral anharmonicity correction using a lower level method.

preprint2009arXiv

Efficient evaluation of accuracy of molecular quantum dynamics using dephasing representation

Ab initio methods for electronic structure of molecules have reached a satisfactory accuracy for calculation of static properties, but remain too expensive for quantum dynamical calculations. We propose an efficient semiclassical method for evaluating the accuracy of a lower level quantum dynamics, as compared to a higher level quantum dynamics, without having to perform any quantum dynamics. The method is based on dephasing representation of quantum fidelity and its feasibility is demonstrated on the photodissociation dynamics of CO2. We suggest how to implement the method in existing molecular dynamics codes and describe a simple test of its applicability.

preprint2009arXiv

Path integral evaluation of equilibrium isotope effects

A general and rigorous methodology to compute the quantum equilibrium isotope effect is described. Unlike standard approaches, ours does not assume separability of rotational and vibrational motions and does not make the harmonic approximation for vibrations or rigid rotor approximation for the rotations. In particular, zero point energy and anharmonicity effects are described correctly quantum mechanically. The approach is based on the thermodynamic integration with respect to the mass of isotopes and on the Feynman path integral representation of the partition function. An efficient estimator for the derivative of free energy is used whose statistical error is independent of the number of imaginary time slices in the path integral, speeding up calculations by a factor of 60 at 500 K. We describe the implementation of the methodology in the molecular dynamics package Amber 10. The method is tested on three [1,5] sigmatropic hydrogen shift reactions. Because of the computational expense, we use ab initio potentials to evaluate the equilibrium isotope effects within the harmonic approximation, and then the path integral method together with semiempirical potentials to evaluate the anharmonicity corrections. Our calculations show that the anharmonicity effects amount up to 30% of the symmetry reduced reaction free energy. The numerical results are compared with recent experiments of Doering and coworkers, confirming the accuracy of the most recent measurement on 2,4,6,7,9-pentamethyl-5-(5,5-$^2$H$_2$)methylene-11,11a-dihydro-12H-naphthacene as well as concerns about compromised accuracy, due to side reactions, of another measurement on 2-methyl-10-(10,10-$^2$H$_2$)methylenebicyclo[4.4.0]dec-1-ene.