Restrictions on second order response properties of quantum states

We distinguish two extreme classes of perturbation problems depending on the signs of second-order response properties. The first class refers to a positive value of the same for any state, and is overwhelmingly more probable. The other category offers all-but-one negative values, or at least some negative values for highly excited states. The classes are seen to differ in reproducing results of finite-dimensional matrix Hamiltonian perturbations, allowing the emergence of a type of sum rule. A few analytical findings are employed for direct demonstration. The outcomes provide notable restrictions on second order response properties of quantum states.     

Full optimization of Jastrow-Slater wave functions with application to the first-row atoms and homonuclear diatomic molecules

We pursue the development and application of the recently introduced linear optimization method for determining the optimal linear and nonlinear parameters of Jastrow-Slater wave functions in a variational Monte Carlo framework. In this approach, the optimal parameters are found iteratively by diagonalizing the Hamiltonian matrix in the space spanned by the wave function and its first-order derivatives, making use of a strong zero-variance principle. We extend the method to optimize the exponents of the basis functions, simultaneously with all the other parameters, namely, the Jastrow, configuration state function, and orbital parameters. We show that the linear optimization method can be thought of as a so-called augmented Hessian approach, which helps explain the robustness of the method and permits us to extend it to minimize a linear combination of the energy and the energy variance. We apply the linear optimization method to obtain the complete ground-state potential energy curve of the C(2) molecule up to the dissociation limit and discuss size consistency and broken spin-symmetry issues in quantum Monte Carlo calculations. We perform calculations for the first-row atoms and homonuclear diatomic molecules with fully optimized Jastrow-Slater wave functions, and we demonstrate that molecular well depths can be obtained with near chemical accuracy quite systematically at the diffusion Monte Carlo level for these systems.     

Ab initio density functional theory applied to quasidegenerate problems

Ab initio density functional theory (DFT), previously applied primarily at the second-order many-body perturbation theory (MBPT) level, is generalized to selected infinite-order effects by using a new coupled-cluster perturbation theory (CCPT). This is accomplished by redefining the unperturbed Hamiltonian in ab initio DFT to correspond to the CCPT2 orbital dependent functional. These methods are applied to the Be-isoelectronic systems as an example of a quasidegenerate system. The CCPT2 variant shows better convergence to the exact quantum Monte Carlo correlation potential for Be than any prior attempt. When using MBPT2, the semicanonical choice of unperturbed Hamiltonian, plays a critical role in determining the quality of the obtained correlation potentials and obtaining convergence, while the usual Kohn-Sham choice invariably diverges. However, without the additional infinite-order effects, introduced by CCPT2, the final potentials and energies are not sufficiently accurate. The issue of the effects of the single excitations on the divergence in ordinary OEP2 is addressed, and it is shown that, whereas their individual values are small, their infinite-order summation is essential to the good convergence of ab initio DFT.     

Algorithmic dimensionality reduction for molecular structure analysis

Dimensionality reduction approaches have been used to exploit the redundancy in a Cartesian coordinate representation of molecular motion by producing low-dimensional representations of molecular motion. This has been used to help visualize complex energy landscapes, to extend the time scales of simulation, and to improve the efficiency of optimization. Until recently, linear approaches for dimensionality reduction have been employed. Here, we investigate the efficacy of several automated algorithms for nonlinear dimensionality reduction for representation of trans, trans-1,2,4-trifluorocyclo-octane conformation--a molecule whose structure can be described on a 2-manifold in a Cartesian coordinate phase space. We describe an efficient approach for a deterministic enumeration of ring conformations. We demonstrate a drastic improvement in dimensionality reduction with the use of nonlinear methods. We discuss the use of dimensionality reduction algorithms for estimating intrinsic dimensionality and the relationship to the Whitney embedding theorem. Additionally, we investigate the influence of the choice of high-dimensional encoding on the reduction. We show for the case studied that, in terms of reconstruction error root mean square deviation, Cartesian coordinate representations and encodings based on interatom distances provide better performance than encodings based on a dihedral angle representation.     

Quantum solution of coupled harmonic oscillator systems beyond normal coordinates

Normal coordinates can be defined as orthogonal linear combinations of coordinates that remove the second order couplings in coupled harmonic oscillator systems. In this paper we go further and explore the possibility of using linear although non-orthogonal coordinate transformations to get the quantum solution of coupled systems. The idea is to use as non-orthogonal linear coordinates those which allow us to express the second-order Hamiltonian matrix in a block diagonal form. To illustrate the viability of this treatment, we first apply it to a system of two bilinearly coupled harmonic oscillators which admits analytical exact solutions. The method provides in this case, as an extra mathematical result, the analytical expressions for the eigenvalues of a certain type of symmetrical tridiagonal matrices. Second, we carry out a numerical application to the Barbanis coupled oscillators system, which contains a third order coupling term and cannot be solved in closed form. We demonstrate that the non-orthogonal coordinates used, named oblique coordinates, are much more efficient than normal coordinates to determine the energy levels and eigenfunctions of this system variationally.     

Generalized electronic diabatic scheme: Diagonalizing the electronic Hamiltonian for artificial molecular systems. How do molecular meccanos move?

A quantum–classic model is presented and used to describe systems ranging from normal molecules up to electronic systems sensed in real space. The quantum system is a set of n‐electrons; a positive background in real space completes the model. A generalized electronic diabatic (GED) theory is introduced. The diabatic functions diagonalize the electronic Hamiltonian for any arrangement of the positive background. Physical quantum states are represented as linear superpositions in the diabatic basis; this latter is always fixed. For systems sensed in real space, the coefficients of the linear superposition are functions of the real space configuration coordinates. Physical changes are produced by interactions with external sources/sinks of energy. An interaction couples different diabatic states; diagonalizing the electronic Hamiltonian plus the couplings leads to new coefficients describing physical states. Among other things, these couplings can be used to simulate the effects produced by scanning tunneling microscopy, atomic force, and transmission electron microscopy on substrates located in real space. The important thing is that time–evolution in electronic Hilbert space can be related to actual motion in real space. The experiment of lateral hopping of a substrate on a metallic surface induced by vibration excitation and followed with scanning tunneling microscope is discussed. A result of the present work is that motion of molecular meccanos reflects then time–evolution in electronic Hilbert space. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Solving the electron and electron-nuclear Schrodinger equations for the excited states of helium atom with the free iterative-complement-interaction method

Very accurate variational calculations with the free iterative-complement-interaction (ICI) method for solving the Schrodinger equation were performed for the 1sNs singlet and triplet excited states of helium atom up to N=24. This is the first extensive applications of the free ICI method to the calculations of excited states to very high levels. We performed the calculations with the fixed-nucleus Hamiltonian and moving-nucleus Hamiltonian. The latter case is the Schrodinger equation for the electron-nuclear Hamiltonian and includes the quantum effect of nuclear motion. This solution corresponds to the nonrelativistic limit and reproduced the experimental values up to five decimal figures. The small differences from the experimental values are not at all the theoretical errors but represent the physical effects that are not included in the present calculations, such as relativistic effect, quantum electrodynamic effect, and even the experimental errors. The present calculations constitute a small step toward the accurately predictive quantum chemistry.     

Does ℏ play a role in multidimensional spectroscopy? Reduced hierarchy equations of motion approach to molecular vibrations

To investigate the role of quantum effects in vibrational spectroscopies, we have carried out numerically exact calculations of linear and nonlinear response functions for an anharmonic potential system nonlinearly coupled to a harmonic oscillator bath. Although one cannot carry out the quantum calculations of the response functions with full molecular dynamics (MD) simulations for a realistic system which consists of many molecules, it is possible to grasp the essence of the quantum effects on the vibrational spectra by employing a model Hamiltonian that describes an intra- or intermolecular vibrational motion in a condensed phase. The present model fully includes vibrational relaxation, while the stochastic model often used to simulate infrared spectra does not. We have employed the reduced quantum hierarchy equations of motion approach in the Wigner space representation to deal with nonperturbative, non-Markovian, and nonsecular system-bath interactions. Taking the classical limit of the hierarchy equations of motion, we have obtained the classical equations of motion that describe the classical dynamics under the same physical conditions as in the quantum case. By comparing the classical and quantum mechanically calculated linear and multidimensional spectra, we found that the profiles of spectra for a fast modulation case were similar, but different for a slow modulation case. In both the classical and quantum cases, we identified the resonant oscillation peak in the spectra, but the quantum peak shifted to the red compared with the classical one if the potential is anharmonic. The prominent quantum effect is the 1-2 transition peak, which appears only in the quantum mechanically calculated spectra as a result of anharmonicity in the potential or nonlinearity of the system-bath coupling. While the contribution of the 1-2 transition is negligible in the fast modulation case, it becomes important in the slow modulation case as long as the amplitude of the frequency fluctuation is small. Thus, we observed a distinct difference between the classical and quantum mechanically calculated multidimensional spectra in the slow modulation case where spectral diffusion plays a role. This fact indicates that one may not reproduce the experimentally obtained multidimensional spectrum for high-frequency vibrational modes based on classical molecular dynamics simulations if the modulation that arises from surrounding molecules is weak and slow. A practical way to overcome the difference between the classical and quantum simulations was discussed.     

On exact calculation of response properties of oscillators in static electric field: A Fourier grid Hamiltonian approach. I. One-dimensional systems

The Fourier grid Hamiltonian method is used to calculate the response properties of different types of 1-d (one-dimensional) quantum oscillators in a uniform static electric field. The calculations are potentially exact. Excepting the harmonic oscillator, the other model oscillators studied are seen to possess nonlinear polarizabilities. In general, the polarizabilities are not monotonic functions of appropriate vibrational quantum numbers. The exact nature of this vibrational-state dependence of polarizabilities is shown to depend on the type of mechanical anharmonicity in which the nuclei move and the nature of electrical anharmonicity characterizing the field–oscillator coupling. The large vibrational contribution to nonlinear polarizabilities often predicted for real diatomics could therefore originate from the mechanical and electrical anharmonicities of the potential in which the nuclei move when placed in a static electric field. © 1994 John Wiley & Sons, Inc.     

New algorithms for iterative matrix‐free eigensolvers in quantum chemistry

New algorithms for iterative diagonalization procedures that solve for a small set of eigen‐states of a large matrix are described. The performance of the algorithms is illustrated by calculations of low and high‐lying ionized and electronically excited states using equation‐of‐motion coupled‐cluster methods with single and double substitutions (EOM‐IP‐CCSD and EOM‐EE‐CCSD). We present two algorithms suitable for calculating excited states that are close to a specified energy shift (interior eigenvalues). One solver is based on the Davidson algorithm, a diagonalization procedure commonly used in quantum‐chemical calculations. The second is a recently developed solver, called the “Generalized Preconditioned Locally Harmonic Residual (GPLHR) method.” We also present a modification of the Davidson procedure that allows one to solve for a specific transition. The details of the algorithms, their computational scaling, and memory requirements are described. The new algorithms are implemented within the EOM‐CC suite of methods in the Q‐Chem electronic structure program. © 2014 Wiley Periodicals, Inc.     

Second-harmonic generation and effect of adsorbates for free-electron metal and semiconductor surfaces

We discuss aspects of a developing microscopic theory of SHG from simple metal and semiconductor surfaces. For semiconductors calculations of the dynamical nonlinear susceptibility on the basis of realistic tight-binding parametrizations of the electronic Hamiltonian provide a practical scheme. In the resulting spectra the effect of the dangling bonds on SHG is clearly seen together with a strong decrease upon saturation with H atoms. In the metal case the adsorbate induced changes of the static nonlinear electron density can be calculated self-consistently by applying density functional theory to the jellium model. The second-order dipole moment determines the effect of adsorbates on the SHG intensity in the adiabatic limit. Quite general a correlation with the nature of the adsorbate expressed by its electronegativity and the characteristic charge transfer, adsorption dipole and polarizabilities in first and second order is found.     

Assignment and extraction of dynamics of a small molecule with a complex vibrational spectrum: thiophosgene

The dispersed fluorescence spectrum of the ground electronic state of thiophosgene, SCCl2, is analyzed in a very complex region of vibrational excitation, 7000-9000 cm(-1). The final result is that most of the inferred excited vibrational levels are assigned in terms of approximate constants of the motion. Furthermore, each level is associated with a rung on a ladder of quantum states on the basis of common reduced dimension fundamental motions. The resulting ladders cannot be identified by any experimental means, and it is the interspersing in energy of their rungs that makes the spectrum complex even after the process of level separation into polyads. Van Vleck perturbation theory is used to create polyad constants of the motion and a spectroscopic Hamiltonian from a potential fitted to experimental data. The eigen functions of this spectroscopic Hamiltonian are rewritten as semiclassical wave functions and transformed to a representation that allows us to analyze and assign the spectra with no other work other than to utilize concepts from nonlinear dynamics.     

On the relationship between the local tracking procedures and monotonic schemes in quantum optimal control

Numerical simulations of (bilinear) quantum control often rely on either monotonically convergent algorithms or tracking schemes. However, despite their mathematical simplicity, very limited intuitive understanding exists at this time to explain the former type of algorithms. Departing from the usual mathematical formalization, we present in this paper an interpretation of the monotonic algorithms as finite horizon, local in time, tracking schemes. Our purpose is not to present a new class of procedures but rather to introduce the necessary rigorous framework that supports this interpretation. As a by-product we show that at each instant, estimates of the future quality of the current control field are available and used in the optimization. When the target is expressed as reaching a prescribed final state, we also present an intuitive geometrical interpretation as the minimization of the distance between two correlated trajectories: one starting from the given initial state and the other backward in time from the target state. As an illustration, a stochastic monotonic algorithm is introduced. Numerical discretizations of the two procedures are also presented.     

Globally convergent trust-region methods for self-consistent field electronic structure calculations

As far as more complex systems are being accessible for quantum chemical calculations, the reliability of the algorithms used becomes increasingly important. Trust-region strategies comprise a large family of optimization algorithms that incorporates both robustness and applicability for a great variety of problems. The objective of this work is to provide a basic algorithm and an adequate theoretical framework for the application of globally convergent trust-region methods to electronic structure calculations. Closed shell restricted Hartree-Fock calculations are addressed as finite-dimensional nonlinear programming problems with weighted orthogonality constraints. A Levenberg-Marquardt-like modification of a trust-region algorithm for constrained optimization is developed for solving this problem. It is proved that this algorithm is globally convergent. The subproblems that ensure global convergence are easy-to-compute projections and are dependent only on the structure of the constraints, thus being extendable to other problems. Numerical experiments are presented, which confirm the theoretical predictions. The structure of the algorithm is such that accelerations can be easily associated without affecting the convergence properties.     

Nonadiabatic electron wavepacket dynamics of molecules in an intense optical field: an ab initio electronic state study

A theory of quantum electron wavepacket dynamics that nonadiabatically couples with classical nuclear motions in intense optical fields is studied. The formalism is intended to track the laser-driven electron wavepackets in terms of the linear combination of configuration-state functions generated with ab initio molecular orbitals. Beginning with the total quantum Hamiltonian for electrons and nuclei in the vector potential of classical electromagnetic field, we reduce the Hamiltonian into a mixed quantum-classical representation by replacing the quantum nuclear momentum operators with the classical counterparts. This framework gives equations of motion for electron wavepackets in an intense laser field through the time dependent variational principle. On the other hand, a generalization of the Newtonian equations provides a matrix form of forces acting on the nuclei for nonadiabatic dynamics. A mean-field approximation to the force matrix reduces this higher order formalism to the semiclassical Ehrenfest theory in intense optical fields. To bring these theories into a practical quantum chemical package for general molecules, we have implemented the relevant ab initio algorithms in it. Some numerical results in the level of the semiclassical Ehrenfest-type theory with explicit use of the nuclear kinematic (derivative) coupling and the velocity form for the optical interaction are presented.     

Quantum algorithm for alchemical optimization in material design

The development of tailored materials for specific applications is an active field of research in chemistry, material science and drug discovery. The number of possible molecules obtainable from a set of atomic species grow exponentially with the size of the system, limiting the efficiency of classical sampling algorithms. On the other hand, quantum computers can provide an efficient solution to the sampling of the chemical compound space for the optimization of a given molecular property. In this work, we propose a quantum algorithm for addressing the material design problem with a favourable scaling. The core of this approach is the representation of the space of candidate structures as a linear superposition of all possible atomic compositions. The corresponding ‘alchemical’ Hamiltonian drives the optimization in both the atomic and electronic spaces leading to the selection of the best fitting molecule, which optimizes a given property of the system, e.g., the interaction with an external potential as in drug design. The quantum advantage resides in the efficient calculation of the electronic structure properties together with the sampling of the exponentially large chemical compound space. We demonstrate both in simulations and with IBM Quantum hardware the efficiency of our scheme and highlight the results in a few test cases. This preliminary study can serve as a basis for the development of further material design quantum algorithms for near-term quantum computers.

‘Alchemical’ quantum algorithm for the simultaneous optimisation of chemical composition and electronic structure for material design. By exploiting quantum mechanical principles this approach will boost drug discovery in the near future.     

Unusual ground states via monotonic convex pair potentials

We have previously shown that inverse statistical-mechanical techniques allow the determination of optimized isotropic pair interactions that self-assemble into low-coordinated crystal configurations in the d-dimensional Euclidean space R(d). In some of these studies, pair interactions with multiple extrema were optimized. In the present work, we attempt to find pair potentials that might be easier to realize experimentally by requiring them to be monotonic and convex. Encoding information in monotonic convex potentials to yield low-coordinated ground-state configurations in Euclidean spaces is highly nontrivial. We adapt a linear programming method and apply it to optimize two repulsive monotonic convex pair potentials, whose classical ground states are counterintuitively the square and honeycomb crystals in R(2). We demonstrate that our optimized pair potentials belong to two wide classes of monotonic convex potentials whose ground states are also the square and honeycomb crystal. We show that these unexpected ground states are stable over a nonzero number density range by checking their (i) phonon spectra, (ii) defect energies and (iii) self assembly by numerically annealing liquid-state configurations to their zero-temperature ground states.     

Real-time time-dependent density functional theory approach for frequency-dependent nonlinear optical response in photonic molecules

We present ab initio calculations of frequency-dependent linear and nonlinear optical responses based on real-time time-dependent density functional theory for arbitrary photonic molecules. This approach is based on an extension of an approach previously implemented for a linear response using the electronic structure program SIESTA. Instead of calculating excited quantum states, which can be a bottleneck in frequency-space calculations, the response of large molecular systems to time-varying electric fields is calculated in real time. This method is based on the finite field approach generalized to the dynamic case. To speed the nonlinear calculations, our approach uses Gaussian enveloped quasimonochromatic external fields. We thereby obtain the frequency-dependent second harmonic generation beta(-2omega;omega,omega), the dc nonlinear rectification beta(0;-omega,omega), and the electro-optic effect beta(-omega;omega,0). The method is applied to nanoscale photonic nonlinear optical molecules, including p-nitroaniline and the FTC chromophore, i.e., 2-[3-Cyano-4-(2-{5-[2-(4-diethylamino-phenyl)-vinyl]-thiophen-2-yl}-vinyl)-5,5-dimethyl-5H-furan-2-ylidene]-malononitrile, and yields results in good agreement with experiment.