Semi‐classical formulation of time‐dependent discrete variable representation method

In this article, we apply a novel time‐dependent discrete variable representation (TDDVR) method proposed by Barkakaty and Adhikari to investigate tunneling through an Eckart barrier. This semi‐classical method is theoretically rigorous and straightforward to implement. Among the TDDVR formulations, this report presents the first derivation of a rigorous form of quantum force (QF) for the present perspective. The validity of this semi‐classical approach is demanded based on the excellent agreement of the tunneling probability with the corresponding quantum results. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Nonsingular constraints in time‐dependent variational principle for parametrized wave functions

In this work, we consider two conditions required for the nonsingularity of constraints in the time‐dependent variational principle (TDVP) for parametrized wave functions. One is the regularity condition which assures the static nonsingularity of the constraint surface. The other condition is the second‐class condition of constraints which assures the dynamic nonsingularity of the constraint surface with a symplectic metric. For analytic wave functions for complex TDVP‐parameters, the regularity and the second‐class conditions become equivalent. The second‐class condition for expectation values is reduced to the noncommutability of the corresponding quantum operators. The symplectic singularity of the equation of motion of TDVP is also shown to be a local breakdown of the second‐class condition in an extended canonical phase‐space. © 2012 Wiley Periodicals, Inc.     

A global potential energy surface and time‐dependent quantum wave packet calculation of Au + H2 reaction

A global potential energy surface (PES) corresponding to the ground state of AuH₂ system has been constructed based on 22 853 ab initio energies calculated by the multireference configuration interaction method with a Davidson correction. The neural network method is used to fit the PES, and the root mean square error is only 1.87 meV. The topographical features of the novel global PES are compared with previous PES which is constructed by Zanchet et al. (Zanchet PES). The global minimum energy reaction paths on the two PESs both have a well and a barrier. Relative to the Au + H₂ reactants, the energy of well is 0.316 eV on the new PES, which is 0.421 eV deeper than Zanchet PES. The calculation of Au(²S) + H₂(X¹Σ_g⁺) → AuH(X¹Σ⁺) + H(²S) dynamical reaction is carried out on new PES, by the time‐dependent quantum wave packet method (TDWP) with second order split operator. The reaction probabilities, integral cross‐sections (ICSs) and differential cross‐sections are obtained from the dynamics calculation. The threshold in the reaction is about 1.46 eV, which is 0.07 eV smaller than Zanchet PES due to the different endothermic energies on the two PESs. At low collision energy (<2.3 eV), the total ICS is larger than the result obtained on Zanchet PES, which can be attributed to the difference of the wells and endothermic energies.     

Theory of time‐dependent transport in quantum dot systems

The formulation of the time‐dependent tunnel current in quantum dot (QD) systems in the (generalized) transfer Hamiltonian formalism is reconsidered, taking into account the nonorthogonality between the subsystems. Exploiting the fact that only the total charge in the system is conserved, in general, gives rise to a new formulation of the transport theory. As a result, it is shown necessary to depart from the orthodox picture, in which the current is treated as a local property of the system. A general formula for the current is derived. By expressing properties of the quantum dot in the Hubbard I approximation with the loop correction, thereby including strong correlations of the QD states, numerical results are given for various time‐dependent source‐drain voltages. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Efficient calculation of open quantum system dynamics and time‐resolved spectroscopy with distributed memory HEOM (DM‐HEOM)

Time‐ and frequency‐resolved optical signals provide insights into the properties of light‐harvesting molecular complexes, including excitation energies, dipole strengths and orientations, as well as in the exciton energy flow through the complex. The hierarchical equations of motion (HEOM) provide a unifying theory, which allows one to study the combined effects of system‐environment dissipation and non‐Markovian memory without making restrictive assumptions about weak or strong couplings or separability of vibrational and electronic degrees of freedom. With increasing system size the exact solution of the open quantum system dynamics requires memory and compute resources beyond a single compute node. To overcome this barrier, we developed a scalable variant of HEOM. Our distributed memory HEOM, DM‐HEOM, is a universal tool for open quantum system dynamics. It is used to accurately compute all experimentally accessible time‐ and frequency‐resolved processes in light‐harvesting molecular complexes with arbitrary system‐environment couplings for a wide range of temperatures and complex sizes. © 2018 Wiley Periodicals, Inc.     

Electronuclear multiconfiguration time‐dependent hartree calculations on the confined H atom with mobile electron and nucleus

A multiconfiguration time‐dependent Hartree method oriented toward calculations of a non‐Born‐Oppenheimer nature has been applied to the calculation of the dynamical properties of a confined H atom. The calculation is fully six‐dimensional and does not take into account constraints arising from linear or angular momentum conservation. The orbital evolution is monitored and the energy level spectrum of the system, as well as the dependence of the results on the decomposition of the Hamiltonian and on the correlation between radial degrees of freedom, is determined. © 2012 Wiley Periodicals, Inc.     

Dielectric relaxation in proteins: a continuum electrostatics model incorporating dielectric heterogeneity of the protein and time‐dependent charges

A boundary element formulation of continuum electrostatics is used to examine time‐independent dielectric relaxation and screening in two proteins, and time‐dependent relaxation in two simpler solutes. Cytochrome c oxidation is modeled by inserting partial charges on the heme, using one to three dielectric regions in the protein. It was suggested recently that for charge insertion on a protein‐bound ligand, all or part of the ligand should be treated as a cavity within the protein medium. Here, the effect of an internal cavity surrounding the central heme atoms is examined, considering separately the static and relaxation (or reorganization) free energies. The former is the free energy to remove the redox electron while maintaining the rest of the structure and charge distribution fixed; the latter is the free energy associated with the relaxation into the product state after the corresponding constraints are released. The effect of the cavity is found to be small for the static free energy, while for the relaxation free energy it is large, as polarization of groups immediately around the heme dominates the relaxation. If the protein surface groups are treated as a distinct medium with a dielectric of 25 (as suggested by recent molecular dynamics simulations), the relaxation free energy decreases significantly (from −37.0 to −43.9 kcal/mol), compared to a model where the whole protein has a dielectric constant of two. Therefore, with this model, although polarization of groups immediately around the heme still dominates the relaxation, polar groups near the protein surface also contribute significantly, and solvent negligibly. The screening of an applied field within myoglobin is calculated, with the protein surrounded by either a low‐dielectric or a high‐dielectric glass. In the vicinity of the CO ligand, the screening is approximately isotropic with a low‐dielectric glass. It is anisotropic with a high‐dielectric glass, but the applied and local fields are still approximately parallel. This has implications for experiments that probe dielectric screening in proteins with the newly developed technique of vibrational Stark spectroscopy: with a high‐dielectric glass, a single, rotationally averaged screening factor can be used, the local field being about 1.65 times the applied field. Finally, we calculate the time‐dependent relaxation in response to instantaneous charge insertion within a spherical cavity in a Debye solvent, and to photoexcitation of a tryptophan solute, illustrating the extension of the boundary element formulation to time‐dependent problems. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 290–305, 2001     

Development of spin‐dependent relativistic open‐shell Hartree–Fock theory with time‐reversal symmetry (I): The unrestricted approach

An open‐shell Hartree–Fock (HF) theory for spin‐dependent, two‐component relativistic calculations, termed the Kramers‐unrestricted HF (KUHF) method, is developed. The present KUHF method, which is formulated as a relativistic counterpart of nonrelativistic UHF, is based on quaternion algebra and partly uses time‐reversal symmetry. The fundamental characteristics of KUHF are discussed in this study. From numerical assessments, it was revealed that KUHF gives a corresponding solution to nonrelativistic UHF; furthermore, KUHF properly describes spin‐orbit interactions. In addition, KUHF can improve the self‐consistent field convergence behavior in spin‐dependent calculations, for example, for f‐block elements.     

Highly efficient implementation of pseudospectral time‐dependent density‐functional theory for the calculation of excitation energies of large molecules

We have developed and implemented pseudospectral time‐dependent density‐functional theory (TDDFT) in the quantum mechanics package Jaguar to calculate restricted singlet and restricted triplet, as well as unrestricted excitation energies with either full linear response (FLR) or the Tamm–Dancoff approximation (TDA) with the pseudospectral length scales, pseudospectral atomic corrections, and pseudospectral multigrid strategy included in the implementations to improve the chemical accuracy and to speed the pseudospectral calculations. The calculations based on pseudospectral time‐dependent density‐functional theory with full linear response (PS‐FLR‐TDDFT) and within the Tamm–Dancoff approximation (PS‐TDA‐TDDFT) for G2 set molecules using B3LYP/6‐31G*^* show mean and maximum absolute deviations of 0.0015 eV and 0.0081 eV, 0.0007 eV and 0.0064 eV, 0.0004 eV and 0.0022 eV for restricted singlet excitation energies, restricted triplet excitation energies, and unrestricted excitation energies, respectively; compared with the results calculated from the conventional spectral method. The application of PS‐FLR‐TDDFT to OLED molecules and organic dyes, as well as the comparisons for results calculated from PS‐FLR‐TDDFT and best estimations demonstrate that the accuracy of both PS‐FLR‐TDDFT and PS‐TDA‐TDDFT. Calculations for a set of medium‐sized molecules, including C_n fullerenes and nanotubes, using the B3LYP functional and 6‐31G^** basis set show PS‐TDA‐TDDFT provides 19‐ to 34‐fold speedups for C_n fullerenes with 450–1470 basis functions, 11‐ to 32‐fold speedups for nanotubes with 660–3180 basis functions, and 9‐ to 16‐fold speedups for organic molecules with 540–1340 basis functions compared to fully analytic calculations without sacrificing chemical accuracy. The calculations on a set of larger molecules, including the antibiotic drug Ramoplanin, the 46‐residue crambin protein, fullerenes up to C₅₄₀ and nanotubes up to 14×(6,6), using the B3LYP functional and 6‐31G^** basis set with up to 8100 basis functions show that PS‐FLR‐TDDFT CPU time scales as N^2.05 with the number of basis functions. © 2016 Wiley Periodicals, Inc.     

Development of spin‐dependent relativistic open‐shell Hartree–Fock theory with time‐reversal symmetry (II): The restricted open‐shell approach

An open‐shell Hartree–Fock (HF) theory for spin‐dependent two‐component relativistic calculations, termed the Kramers‐restricted open‐shell HF (KROHF) method, is developed. The present KROHF method is defined as a relativistic analogue of ROHF using time‐reversal symmetry and quaternion algebra, based on the Kramers‐unrestricted HF (KUHF) theory reported in our previous study (Int. J. Quantum Chem., doi: 10.1002/qua.25356 ). As seen in the nonrelativistic ROHF theory, the ambiguity of the KROHF Fock operator gives physically meaningless spinor energies. To avoid this problem, the canonical parametrization of KROHF to satisfy Koopmans' theorem is also discussed based on the procedure proposed by Plakhutin et al. (J. Chem. Phys. 2006 , 125, 204110). Numerical assessments confirmed that KROHF using Plakhutin's canonicalization procedure correctly gives physical spinor energies within the frozen‐orbital approximation under spin–orbit interactions.     

Nonstationarity and related measures for time‐dependent hartree–fock and multiconfigurational models

Based on an earlier article (Eberly and Singh, Phys. Rev. D 1973 , 7, 359) and related works on short‐time evolution, this article proposes a many‐electron formulation for the nonstationarity degree which can be assigned to quantum system at each time point. The key measure introduced, , is a nonstationarity index that can be thought of as an inverse nominal lifetime at each instance of time. The index is directly computed from the time derivative of one‐electron density matrix and is a size‐consistent quantity. In this article, the approach is developed for the time‐dependent Hartree–Fock (TDHF), single‐excitation (TDCIS), and time‐dependent full configuration interaction (TDFCI) models. As a rule, nonstationarity effects are more pronounced in correlated electron systems, and a joint analysis of and the multiconfigurational character of wave functions apparently provide a deeper insight into dynamical molecular processes. The performed calculations on small molecules in laser fields show a preference for the TDCIS model when comparing TDCIS and TDHF with the “exact” TDFCI model. © 2013 Wiley Periodicals, Inc.     

Applications of automatic differentiation to the time‐dependent Schrödinger equation

A new approach based upon the Taylor series method is proposed for propagating solutions of the time‐dependent Schrödinger equation. Replacing the spatial derivative of the wave function with finite difference formulas, we derive a recursive formula for the evaluation of Taylor coefficients. The automatic differentiation technique is used to recursively calculate the required Taylor coefficients. We also develop an implicit scheme for the recursive evaluation of these coefficients. We then advance the solution in time using a Taylor series expansion. Excellent computational results are obtained when this method is applied to a one‐dimensional reflectionless potential and a two‐dimensional barrier transmission problem. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

A quantitative analysis of light‐driven charge transfer processes using voronoi partitioning of time dependent DFT‐derived electron densities

An analytical method is presented that provides quantitative insight into light‐driven electron density rearrangement using the output of standard time‐dependent density functional theory (TD‐DFT) computations on molecular compounds. Using final and initial electron densities for photochemical processes, the subtraction of summed electron density in each atom‐centered Voronoi polyhedron yields the electronic charge difference, Q ^VECD. This subtractive method can also be used with Bader, Mulliken and Hirshfeld charges. A validation study shows Q ^VECD to have the most consistent performance across basis sets and good conservation of charge between electronic states. Besides vertical transitions, relaxation processes can be investigated as well. Significant electron transfer is computed for isomerization on the excited state energy surface of azobenzene. A number of linear anilinepyridinium donor‐bridge‐acceptor chromophores was examined using Q ^VECD to unravel the influence of its pi‐conjugated bridge on charge separation. Finally, the usefulness of the presented method as a tool in optimizing charge transfer is shown for a homologous series of organometallic pigments. The presented work allows facile calculation of a novel, relevant quantity describing charge transfer processes at the atomic level. © 2017 The Authors Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Excited states of OsO4: A comprehensive time‐dependent relativistic density functional theory study

A large number of scalar as well as spinor excited states of OsO₄, in the experimentally accessible energy range of 3–11 eV, have been captured by time‐dependent relativistic density functional linear response theory based on an exact two‐component Hamiltonian resulting from the symmetrized elimination of the small component. The results are grossly in good agreement with those by the singles and doubles coupled‐cluster linear response theory in conjunction with relativistic effective core potentials. The simulated‐excitation spectrum is also in line with the available experiment. Furthermore, combined with detailed analysis of the excited states, the nature of the observed optical transitions is clearly elucidated. It is found that a few scalar states of ³T₁ and ³T₂ symmetries are split significantly by the spin‐orbit coupling. The possible source for the substantial spin‐orbit splittings of ligand molecular orbitals is carefully examined, leading to a new interpretation on the primary valence photoelectron ionization spectrum of OsO₄. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Efficient blue‐emitting Ir(III) complexes with phenyl‐methyl‐benzimidazolyl and picolinate ligands: A DFT and time‐dependent DFT study

The series of heteroleptic cyclometalated Ir(III) complexes for organic light‐emitting display application have been investigated theoretically to explore their electronic structures and spectroscopic properties. The geometries, electronic structures, and the lowest‐lying singlet absorptions and triplet emissions of Ir‐(pmb)₃ and theoretically designed models Ir‐(Rpmb)₂pic were investigated with density functional theory (DFT)‐based approaches, where pmb = phenyl‐methyl‐benzimidazolyl, pic = picolinate, and R = H/F. Their structures in the ground and excited states have been optimized at the DFT/B3LYP/LANL2DZ and TDDFT/B3LYP/LANL2DZ levels, and the lowest absorptions and emissions were evaluated at B3LYP and M062X level of theory, respectively. The mobility of holes and electrons were studied computationally based on the Marcus theory. Calculations of ionization potentials were used to evaluate the injection abilities of holes into these complexes. The reasons for the higher electroluminescence efficiency and phosphorescence quantum yields in Ir‐(Rpmb)₂pic than in Ir‐(pmb)₃ have been investigated. The designed moleculars are expected to be highly emissive in pure‐blue region. © 2013 Wiley Periodicals, Inc.