Quantum control mechanism analysis through field based Hamiltonian encoding

Optimal control of quantum dynamics in the laboratory is proving to be increasingly successful. The control fields can be complex, and the mechanisms by which they operate have often remained obscure. Hamiltonian encoding (HE) has been proposed as a method for understanding mechanisms in quantum dynamics. In this context mechanism is defined in terms of the dominant quantum pathways leading to the final state of the controlled system. HE operates by encoding a special modulation into the Hamiltonian and decoding its signature in the dynamics to determine the dominant pathway amplitudes. Earlier work encoded the modulation directly into the Hamiltonian operators. This present work introduces the alternative scheme of field based HE, where the modulation is encoded into the control field and not directly into the Hamiltonian operators. This distinct form of modulation yields a new perspective on mechanism and is computationally faster than the earlier approach. Field based encoding is also an important step towards a laboratory based algorithm for HE as it is the only form of encoding that may be experimentally executed. HE is also extended to cover systems with noise and uncertainty and finally, a hierarchical algorithm is introduced to reveal mechanism in a stepwise fashion of ever increasing detail as desired. This new hierarchical algorithm is an improvement over earlier approaches to HE where the entire mechanism was determined in one stroke. The improvement comes from the use of less complex modulation schemes, which leads to fewer evaluations of Schrodinger's equation. A number of simulations are presented on simple systems to illustrate the new field based encoding technique for mechanism assessment.     

Vibration-rotation interactions between overtone and combination levels of asymmetric-top molecules: application to the infrared spectroscopy of formaldehyde and ketene

The conventional vibration-rotation Hamiltonian for an asymmetric-top molecule is rewritten by expanding the elements of the inverse inertial tensor about the equilibrium molecular geometry. The approach allows the identification of terms in the Hamiltonian that couple states differing by two, three, or four vibrational quanta and hence the calculation of dimensioned Coriolis xi coupling coefficients for interacting fundamental, overtone, and combination levels. The matrix elements that result from the application of the expanded Hamiltonian depend upon the harmonic vibrational wave numbers, equilibrium moments of inertia, Coriolis zeta parameters, and the derivatives of the elements of the inertial tensor matrix with respect to each of the normal coordinates. The Coriolis coupling coefficients may be calculated through evaluation of the summations that result from the appropriate terms. The validity of the approach is demonstrated through the calculation of coupling coefficients for interacting levels in formaldehyde and ketene. The uncertainty in the calculated values of the coupling coefficients is typically better than +/-6%, although the values calculated for interactions that involve low-frequency vibrational modes are less reliable. Comparisons are made between the calculated values and experimental results.     

Density matrix structure for saturated molecules

Conclusion The existence of a common Hamiltonian matrix structure for saturated systems results in common structural properties of the density matrices for the whole class of molecules, such as the zero occupation of AO in the first approximation, the density matrix perturbations due to a heteroatom, etc. This fact can be taken as a quantum-mechanical foundation for viewing saturated molecules as a separate class of compounds. The endowment of this system with the transferability of electronic structure properties, relative to atoms and bonds, to high accuracy, within the framework of the effective Hamiltonian method follows from an analysis of the general expressions for the density matrix elements. The transferability of the saturated system Hamiltonian matrix elements requisite for this is supported by a comparison among the self-consistent Fock matrix elements of various hydrocarbons in a localized orbital basis [9]. Independently of the detailed structure of the actual molecules, the influence of a heteroatom on the electron density distribution in saturated systems dies off quickly with distance from the heteroatom. From an analysis of expressions for the nondiagonal elements of the density matrix corresponding to nonneighboring AO we establish a connection between the degree of electron localization in saturated systems and the size or certain Hamiltonian matrix elements.There is a consequent analogy between saturated and alternatively conjugated hydrocarbons, which, starting from the common structure of the Hamiltonians, also leads to common properties of the density matrices [14]. However, the study of the influence of heteroatoms on the density matrices in these systems by means of perturbation theory is complicated by the dependence of the matrix N(4) on the molecular structure, which makes it necessary to introduce highly simplified approaches for the solution of Eq. (2) [8]. Therefore, for alternating hydrocarbons we have succeeded in establishing only the sign of the orbital — orbital polarizability (alternating polarity theorem [14]), while, as for saturated systems, the equality N=1 permits an analytic expression for the polarizability.V. Kapsukas Vilnius State University. Translated from Zhurnal Strukturnoi Khimii, Vol. 29, No. 5, pp. 3–8, September–October, 1988.     

Quantum control by means of hamiltonian structure manipulation

A traditional quantum optimal control experiment begins with a specific physical system and seeks an optimal time-dependent field to steer the evolution towards a target observable value. In a more general framework, the Hamiltonian structure may also be manipulated when the material or molecular 'stockroom' is accessible as a part of the controls. The current work takes a step in this direction by considering the converse of the normal perspective to now start with a specific fixed field and employ the system's time-independent Hamiltonian structure as the control to identify an optimal form. The Hamiltonian structure control variables are taken as the system energies and transition dipole matrix elements. An analysis is presented of the Hamiltonian structure control landscape, defined by the observable as a function of the Hamiltonian structure. A proof of system controllability is provided, showing the existence of a Hamiltonian structure that yields an arbitrary unitary transformation when working with virtually any field. The landscape analysis shows that there are no suboptimal traps (i.e., local extrema) for controllable quantum systems when unconstrained structural controls are utilized to optimize a state-to-state transition probability. This analysis is corroborated by numerical simulations on model multilevel systems. The search effort to reach the top of the Hamiltonian structure landscape is found to be nearly invariant to system dimension. A control mechanism analysis is performed, showing a wide variety of behavior for different systems at the top of the Hamiltonian structure landscape. It is also shown that reducing the number of available Hamiltonian structure controls, thus constraining the system, does not always prevent reaching the landscape top. The results from this work lay a foundation for considering the laboratory implementation of optimal Hamiltonian structure manipulation for seeking the best control performance, especially with limited electromagnetic resources.     

Finite group theory for large systems. 3. Symmetry-generation of reduced matrix elements for icosahedral C(20) and C(60) molecules

This paper uses symmetry-generation to simplify the determination of Hamiltonian reduced matrix elements. It is part of a series on using computers to apply finite group theory to quantum mechanical calculations on large systems. Symmetry-generation is an expression of the whole molecule as a sum of symmetry transformations on a smaller reference structure. Then on a suitably-conditioned symmetry-adapted basis, the reduced matrix elements of the Hamiltonian are averages of certain elements of the simpler reference structure matrix. The smaller the reference structure, the greater is the computational savings. Single atom reference structures are used here for the Hückel treatment of icosahedral C(20) and C(60) fullerenes. The analytical power of this approach is illustrated by determining the two bond lengths of C(60) from spectral data.     

Comparative study of perturbative methods for computing electron transfer tunneling matrix elements with a nonorthogonal basis set

The authors consider the problem of computing tunneling matrix elements for bridge-mediated electron transfer reactions using the Lowdin [J. Math. Phys. 3, 969 (1962); J. Mol. Spectrosc. 13, 326 (1964)] projection-iteration technique with a nonorthogonal basis set. They compare the convergence properties of two different Lowdin projections, one containing the overlap matrix S and the other containing the inverse S-1 in the projected Hamiltonian. It was suggested in the literature that the projected Hamiltonian with S-1 has better convergence properties compared to the projected Hamiltonian with S. The authors test this proposal using a simple analytical model, and ab initio Hartree-Fock calculations on different molecules with several types of basis sets. Their calculations show that, for Gaussian-type basis sets, the projected Hamiltonian containing S has the best convergence properties, especially for diffuse basis sets and in the strong coupling limit. The limit of diffuse basis sets is relevant to tunneling matrix element calculations involving excited states and anionic electron transfer.     

Elements of irreducible tensorial matrices generated by finite groups with applications to ligand field Hamiltonians

Using symmetry to determine Hamiltonian matrix elements for quantum systems with finite group symmetry is a special case of obtaining group-generated irreducible tensorial matrices. A group-generated irreducible tensorial matrix transforms irreducibly under the group and is a linear combination of group transformations on a reference matrix. The reference matrix elements may be appropriate integrals or parameters. The methods of normalized irreducible tensorial matrices (NITM) are employed to express elements of the generated matrix in terms of those of the reference matrix without performing the actual transformations. Only NTTM components of the reference matrix with the same transformation properties as the group-generated matrix will contribute to its elements. The elements of invariant symmetry-generated matrices are proportional to simple averages of certain elements of the reference matrix. This relation is substantially more efficient than previous techniques for evaluating matrix elements of octahedral and tetragonal d-type ligand-field Hamiltonians.     

Implementation of screened hybrid density functional for periodic systems with numerical atomic orbitals: basis function fitting and integral screening

We present an efficient O(N) implementation of screened hybrid density functional for periodic systems with numerical atomic orbitals (NAOs). NAOs of valence electrons are fitted with gaussian-type orbitals, which is convenient for the calculation of electron repulsion integrals and the construction of Hartree-Fock exchange matrix elements. All other parts of Hamiltonian matrix elements are constructed directly with NAOs. The strict locality of NAOs is adopted as an efficient two-electron integral screening technique to speed up calculations.     

Spin-independent three-body effective valence-shell operators: Application to molecular oxygen

A mathematical construction is presented that uniquely defines a set of spin-independent effective valence-shell Hamiltonian (H^v) three-body matrix elements. These spin-independent H^v matrix elements separate direct and exchange portions of the three-body H^v matrix elements and therefore provide the most natural form for comparisons with parameterization schemes of semiempirical electronic structure methods in which the three-body matrix elements are incorporated into semiempirical one- and two-body Hamiltonian matrix elements in an averaged manner. Ab initio H^v three-body matrix elements of O₂ are computed through third order of quasidegenerate perturbation theory and are analyzed as a function of internuclear distance and atomic orbital overlap to aid in understanding how these three-body matrix elements may be averaged into semiempirical one- and two-body matrix elements. © 1992 John Wiley & Sons, Inc.     

Renner-Teller vibronic analysis for a tetra-atomic molecule. I. The effective Hamiltonian and matrix elements

The effective vibronic Hamiltonian for a linear tetra-atomic molecule in a Pi state has been investigated. In addition to the usual vibrational and Renner-Teller coupling terms, the bending mode anharmonicity, spin-orbit coupling, and Fermi resonance interactions have been added to the model. Terms in the Hamiltonian up to the fourth order are given explicitly for molecules of C(infinityupsilon) symmetry and simplifications for symmetric D(infinityh) molecules are discussed. The matrix elements for the HCCS free radical have been obtained and are used to analyze the observed ground-state levels of HCCS and DCCS in a companion paper. The Sears resonance vibronic interaction that couples levels with the selection rules DeltaK=+/-1, DeltaSigma=-/+1, and DeltaP=0 has also been studied and the matrix elements derived. The determinable combinations of signs for the major parameters in the model are discussed.     

Electronic structure of saturated hydrocarbons in the semiempirical equivalent orbital method

The problem of matrix elements of the Hartree-Fock Hamiltonian for saturated hydrocarbons in the EO method is considered. It is stated, that all the proposed scales of such parameters are incomplete or lead to an incorrect band structure of diamond. It is shown that the data on the band structure of diamond may be useful to obtain the full scale of the parameters for the calculations of the electronic structure of hydrocarbons.     

Triplet-triplet energy-transfer coupling: theory and calculation

Triplet-triplet (TT) energy transfer requires two molecular fragments to exchange electrons that carry different spin and energy. In this paper, we analyze and report values of the electronic coupling strengths for TT energy transfer. Two different methods were proposed and tested: (1) Directly calculating the off-diagonal Hamiltonian matrix element. This direct coupling scheme was generalized from the one used for electron transfer coupling, where two spin-localized unrestricted Hartree-Fock wave functions are used as the zero-order reactant and product states, and the off-diagonal Hamiltonian matrix elements are calculated directly. (2) From energy gaps derived from configuration-interaction-singles (CIS) scheme. Both methods yielded very similar results for the systems tested. For TT coupling between a pair of face-to-face ethylene molecules, the exponential attenuation factor is 2.59 A(-1)(CIS6-311+G(**)), which is about twice as large as typical values for electron transfer. With a series of fully stacked polyene pairs, we found that the TT coupling magnitudes and attenuation rates are very similar irrespective of their molecular size. If the polyenes were partially stacked, TT couplings were much reduced, and they decay more rapidly with distance than those of full-stacked systems. Our results showed that the TT coupling arises mainly from the region of close contact between the donor and acceptor frontier orbitals, and the exponential decay of the coupling with separation depends on the details of the molecular contacts. With our calculated results, nanosecond or picosecond time scales for TT energy-transfer rates are possible.     

Optimization of molecular electronic structure calculations. 2. Calculation of symmetry-reduced matrix elements

A suggested formalism of the local symmetricized orbitals in conjunction with the selection technique for independent blocks of integrals in an original basis is used for a construction of multielectron Hamiltonian matrix elements in the symmetry orbital basis. The optimal molecular electronic structure calculation algorithm with the Hartree–Fock–Roothaan method in the symmetricized basis was obtained as a result. The minimal number of fundamentally distinguished (symmetry attributed) elements both in original and in symmetricized basis is used in the calculations.     

On the accuracy of the potential harmonic functions

The matrix elements of the correlation function between symmetric potential harmonics are simplified into the analytical summations of the grand angular momenta by smartly using the recurrence and coupling relations of the potential harmonics. The correlation-function potential-harmonic and generalized-Laguerre-function method (CFPHGLF) , recently developed by us, is applied to the S states of the helium-like systems for Z = 2 to 9. The results exhibit good convergence with the bases in terms of both the angular and radial directions. The final eigen-energies agree excellently with the best s-limits of the variational configuration interaction (CI) method for the involved low-lying S states. The accuracy of the potential harmonic (PH) expansion scheme is discussed relative to the exact Hylleraas CI results (HCI), and Hartree-Fock results. Moreover, suggestion is given for the future improvement of the PH scheme.     

A diagrammatic valence bond method for configuration interaction calculations in atoms and molecules

A diagrammatic valence bond method based on Rumer-Pauling rules for configuration interaction calculations is described. The advantages of this method are that it is simple and flexible and is expected to be computationally efficient as the basis functions can be coded as increasing integers. Evaluation of Hamiltonian matrix elements involves simple bit manipulations and binary searches. The basis, being represented pictorially, should also help in utilizing spatial symmetries for further block-diagonalizing the Hamiltonian matrix. The eigenfunctions of the Hamiltonian can also be used to compute matrix elements between different electronic states.     

Electronic coupling matrix elements from charge constrained density functional theory calculations using a plane wave basis set

We present a plane wave basis set implementation for the calculation of electronic coupling matrix elements of electron transfer reactions within the framework of constrained density functional theory (CDFT). Following the work of Wu and Van Voorhis [J. Chem. Phys. 125, 164105 (2006)], the diabatic wavefunctions are approximated by the Kohn-Sham determinants obtained from CDFT calculations, and the coupling matrix element calculated by an efficient integration scheme. Our results for intermolecular electron transfer in small systems agree very well with high-level ab initio calculations based on generalized Mulliken-Hush theory, and with previous local basis set CDFT calculations. The effect of thermal fluctuations on the coupling matrix element is demonstrated for intramolecular electron transfer in the tetrathiafulvalene-diquinone (Q-TTF-Q(-)) anion. Sampling the electronic coupling along density functional based molecular dynamics trajectories, we find that thermal fluctuations, in particular the slow bending motion of the molecule, can lead to changes in the instantaneous electron transfer rate by more than an order of magnitude. The thermal average, (<|H(ab)|(2)>)(1/2)=6.7 mH, is significantly higher than the value obtained for the minimum energy structure, |H(ab)|=3.8 mH. While CDFT in combination with generalized gradient approximation (GGA) functionals describes the intermolecular electron transfer in the studied systems well, exact exchange is required for Q-TTF-Q(-) in order to obtain coupling matrix elements in agreement with experiment (3.9 mH). The implementation presented opens up the possibility to compute electronic coupling matrix elements for extended systems where donor, acceptor, and the environment are treated at the quantum mechanical (QM) level.     

Calculation of nonadiabatic couplings with restricted open-shell Kohn-Sham density-functional theory

This paper generalizes the recently proposed approaches for calculating the derivative couplings between adiabatic states in density-functional theory (DFT) based on a Slater transition-state density to transitions such as singlet-singlet excitations, where a single-determinant ansatz is insufficient. The proposed approach is based on restricted open-shell Frank et al. [J. Chem. Phys. 108, 4060 (1998)] theory used to describe a spin-adapted Slater transition state. To treat the dependence of electron-electron interactions on the nuclear positions, variational linear-response density-functional perturbation theory is generalized to reference states with an orbital-dependent Kohn-Sham Hamiltonian and nontrivial occupation patterns. The methods proposed in this paper are not limited to the calculation of derivative coupling vectors, but can also be used for the calculation of other transition matrix elements. Moreover, they can be used to calculate the linear response of open-shell systems to arbitrary external perturbations in DFT.     

One-electron spin-orbit contribution by effective nuclear charges

Effective nuclear charges of the main group elements from the second up to the fifth row have been developed for the one-electron part of the spin-orbit (SO) coupling Hamiltonian. These parameters, suitable to be used for SO calculations of large molecular systems, provide a useful and remarkably good approximation to the full SO Hamiltonian. We have derived atomic effective nuclear charges by fitting procedure. Computed fine-structure splitting (FSS) of the doublet and triplet II states of AH species (A is one of the abovementioned elements) have been chosen for this purpose. We have adopted the noniterative scheme, previously reported, according to which SO contributions can be calculated through direct coupling between the II states. The latter have been optimized at B3LYP level using DZVP basis sets. As surrogates for a large number of possible applications, we have widely employed the empirical parameters to compute II-FSSs of diatomic species for which experimental data are available.     

Explicit formulas for the Hamiltonian matrix elements

The explicit formulas for the evaluation of the Hamiltonian matrix elements are presented. The calculation of the integral coefficients is independent of both the nature of the orbitals and th spin coupling schemes. It is fully automatic and only dependent on the number of doubly and singly occupied orbitals. Further-more, the symmetric group representation matrices are not needed, and the N! problem can be avoided.     

Gradients for two-component quasirelativistic methods. Application to dihalogenides of element 116

The authors report the implementation of geometry gradients for quasirelativistic two-component Hartree-Fock and density functional methods using either the zero-order regular approximation Hamiltonian or spin-dependent effective core potentials. The computational effort of the resulting program is comparable to that of corresponding nonrelativistic calculations, as it is dominated by the evaluation of derivative two-electron integrals, which is the same for both types of calculations. Besides the implementation of derivatives of matrix elements of the one-particle Hamiltonian with respect to nuclear displacements, the calculation of the derivative exchange-correlation energy for the open shell case involves complicated expressions because of the noncollinear approach chosen to define the spin density. A pilot application to dihalogenides of element 116 shows how spin-orbit coupling strongly affects the chemistry of the superheavy p-block elements. While these molecules are bent at a scalar-relativistic level, spin-orbit coupling is so strong that only the 7p3/2 atomic orbitals of element 116 are involved in bonding, which favors linear molecular geometries for dihalogenides with heavy terminal halogen atoms.