Effective Floquet Hamiltonians for dipolar and quadrupolar coupled N-spin systems in solid state nuclear magnetic resonance under magic angle spinning

Spin dynamics under magic angle spinning has been studied using different theoretical approaches and also by extensive numerical simulation programs. In this article we present a general theoretical approach that leads to analytic forms for effective Hamiltonians for an N-spin dipolar and quadrupolar coupled system under magic angle spinning (MAS) conditions, using a combination of Floquet theory and van Vleck (contact) transformation. The analytic forms presented are shown to be useful for the study of MAS spin dynamics in solids with the help of a number of simulations in two, three, and four coupled, spin-1/2 systems as well as spins in which quadrupolar interactions are also present.     

A second-quantization framework for the unified treatment of relativistic and nonrelativistic molecular perturbations by response theory

A formalism is presented for the calculation of relativistic corrections to molecular electronic energies and properties. After a discussion of the Dirac and Breit equations and their first-order Foldy-Wouthuysen [Phys. Rev. 78, 29 (1950)] transformation, we construct a second-quantization electronic Hamiltonian, valid for all values of the fine-structure constant alpha. The resulting alpha-dependent Hamiltonian is then used to set up a perturbation theory in orders of alpha(2), using the general framework of time-independent response theory, in the same manner as for geometrical and magnetic perturbations. Explicit expressions are given to second order in alpha(2) for the Hartree-Fock model. However, since all relativistic considerations are contained in the alpha-dependent Hamiltonian operator rather than in the wave function, the same approach may be used for other wave-function models, following the general procedure of response theory. In particular, by constructing a variational Lagrangian using the alpha-dependent electronic Hamiltonian, relativistic corrections can be calculated for nonvariational methods as well.     

Constrained adiabatic trajectory method: a global integrator for explicitly time-dependent Hamiltonians

The constrained adiabatic trajectory method (CATM) is reexamined as an integrator for the Schro?dinger equation. An initial discussion places the CATM in the context of the different integrators used in the literature for time-independent or explicitly time-dependent Hamiltonians. The emphasis is put on adiabatic processes and within this adiabatic framework the interdependence between the CATM, the wave operator, the Floquet, and the (t, t') theories is presented in detail. Two points are then more particularly analyzed and illustrated by a numerical calculation describing the H(2)(+) ion submitted to a laser pulse. The first point is the ability of the CATM to dilate the Hamiltonian spectrum and thus to make the perturbative treatment of the equations defining the wave function possible, possibly by using a Krylov subspace approach as a complement. The second point is the ability of the CATM to handle extremely complex time-dependencies, such as those which appear when interaction representations are used to integrate the system.     

Nonsingular two/one-component relativistic Hamiltonians accurate through arbitrary high order in α2

A series of nonsingular two-component relativistic Hamiltonians is derived from the Dirac Hamiltonian by first performing the free-particle Foldy–Wouthuysen transformation and then a block-diagonalizing transformation. The latter is defined in terms of operators which can be determined iteratively through arbitrary order in α, leading to transformed Hamiltonians with the two-component block accurate through α^2k, k=1, 2, 3,… . These Hamiltonians give relativistic energies which differ from Dirac's energies only in terms higher than α^2k. Their relation to other nonsingular methods of relativistic quantum chemistry (the Douglas–Kroll method, the regular Hamiltonian schemes) is discussed. By removing the spin-dependent operators, the derived Hamiltonians can be written in spin-free one-component form. The computational effort involved is essentially the same as in the case of the Douglas–Kroll scheme and amounts to relatively easy modification of the core Hamiltonian. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 65 : 225–239, 1997     

Combination of perturbative and variational methods for calculating molecular spectra: calculation of the upsilon=3-5 CH stretch overtone spectrum of CHF3

Highly excited states of the CHF3 molecule belonging to the third, fourth, and fifth Fermi polyad are calculated using a combination of the Van Vleck perturbation theory and a variational treatment. The perturbation theory preconditions the Hamiltonian matrix by transforming away all couplings except those between nearly degenerate states. This transformation is implemented so that eigenvalues can be found with significantly smaller matrices than that which would be needed in the original normal mode representation. Even with preconditioning, at the energies as high as 3-5 quanta in the CH stretch, it is not possible to directly diagonalize the Hamiltonian matrix due to the large basis sets required. Iterative methods, particularly the block-Davidson method, are explored for finding the eigenvalues. The methods are compared and the advantages discussed.     

Harmonic oscillator tensors in contact transformation theory

The application of contact transformation theory to the perturbed harmonic oscillator is reexamined in the light of the harmonic oscillator tensors previously presented. It is found that the recasting of the formalism of this problem in terms of harmonic oscillator tensors results in great simplifications, most of which stem from the introduction of the additional algebraic quantum numbers (l, m). The order of magnitude of each fragment of the Hamiltonian is easily recognizable, and the diagonal and nondiagonal parts contained therein are readily identifiable. The determination of the contact transformation operator is reduced to a simple formula. First, an analysis is made for a single mode of vibration, and it is subsequently extended to a multimode case. The perturbed diatomic vibrator is presented as an example.     

A paradox of grid-based representation techniques: accurate eigenvalues from inaccurate matrix elements

Several approximately variational grid-based representation techniques devised to solve the time-independent nuclear-motion Schrödinger equation share a similar behavior: while the computed eigenpairs, the only results which are of genuine interest, are accurate, many of the underlying Hamiltonian matrix elements are inaccurate, deviating substantially from their values in a variational basis representation. Examples are presented for the discrete variable representation and the Lagrange-mesh approaches, demonstrating that highly accurate eigenvalues and eigenfunctions can be obtained even if some or even all of the Hamiltonian matrix elements in these grid-based representations are inaccurate. It is shown how the apparent contradiction of obtaining accurate eigenpairs with far less accurate individual matrix elements can be resolved by considering the unitary transformation between the representations. Furthermore, the relations connecting orthonormal bases and the corresponding Lagrange bases are generalized to relations connecting nonorthogonal, regularized bases and the corresponding nonorthogonal, regularized Lagrange bases.     

Eigenvalue correlation diagrams for effective hamiltonians: The ESR spectrum of axial Gd(III)

Correlation diagrams depicting the behavior of effective Hamiltonian eigenvalues over the ranges of its variables may reveal important properties of the models. The case in which a Hamiltonian is a sum of two terms, one going to zero in one limit and the other effectively zero in the other limit, is considered here. The electron spin resonance spectrum calculated from the spin Hamiltonian of axial gadolinium In is used here as an example. The zerofield splitting term of the spin Hamiltonian is expanded in terms of normalized irreducible tensorial matrices in order to take advantage of their transformation properties under rotations. Its eigenvalues are plotted in a correlation diagram from the zero-field to the high-field limit. A similar correlation diagram for the principal transitions is used to predict a resonance spectrum.     

A Continuum Reaction Field Theory of Polarizable, Nondipolar, Quadrupolar Solvents: Ab Initio Study of Equilibrium Solvation in Benzene

A continuum theory to describe solvation in nondipolar quadrupolar solvents is developed by accounting for electronic polarizability. A general Hamiltonian for a solute–solvent system in an arbitrary nonequilibrium configuration is obtained in terms of two field variables—densities of the solvent quadrupole and induced dipole moments. Equilibrium solvation is studied by optimizing this Hamiltonian with account of cavity boundaries. As an application, electronic structures and free energies of small molecules in benzene are examined with ab initio methods. Solvation stabilization due to solvent quadrupole moments is found to be substantial; for the solutes considered here, it is comparable to and often in excess of that arising from solvent-induced dipole moments.     

Separation of quadrupolar and paramagnetic shift interactions with TOP-STMAS/MQMAS in solid-state lighting phosphors

A new approach for processing satellite-transition magic-angle spinning (STMAS) and multiple-quantum magic-angle spinning (MQMAS) data, based on the two-dimensional one-pulse (TOP) method, which separates the second-rank quadrupolar anisotropy and paramagnetic shift interactions via a double shearing transformation, is described. This method is particularly relevant in paramagnetic systems, where substantial inhomogeneous broadening may broaden the lineshapes. Furthermore, it possesses an advantage over the conventional processing of MQMAS and STMAS spectra because it overcomes the limitation on the spectral width in the indirect dimension imposed by rotor synchronization of the sampling interval. This method was applied experimentally to the Al solid-state nuclear magnetic resonance of a series of yttrium aluminum garnets (YAGs) doped with different lanthanide ions, from which the quadrupolar parameters of paramagnetically shifted and bulk unshifted sites were extracted. These parameters were then compared with density functional theory calculations, which permitted a better understanding of the local structure of Ln substituent ions in the YAG lattice.     

The generalized Hulthén potential

A generalization of the Hulthén potential is presented on the basis of an approach that uses the factorization of a general Hamiltonian by means of a specific model of operational equations with the structure ∼β(r)∓(d/dr). To achieve this goal, the treatment of the V _Hs(r) standard Hulthén potential for bound s states is carried out by proposing a particular β_p(r) ansatz to identify V _Hs(r) by means of a particular Riccati-type relationship. Once the identification has been achieved, the generalized Hulthén potential is obtained straightforwardly with the aid of a general Riccati formula. As expected, the Hamiltonian of the generalized Hulthén potential is isospectral when compared with the corresponding standard Hamiltonian. Moreover, according to the Darboux transform there exists a modified Hulthén potential which is also isospectral. We show that the latter is just a particular case of the generalized Hulthén interaction model. Received: 14 September 1999 / Accepted: 3 February 2000 / Published online: 19 April 2000     

The calculation of vibrational energy levels of polyatomic molecules including anharmonic effect using contact transformation perturbation method

Using contact transformation perturbation method based on the Taylor expansion of the potential energy function in terms of dimensionless normal coordinates up to sixth‐order, the vibrational energy levels in terms of force constants are derived. The contact transformation theory has been applied to simplify the calculation of perturbation effects. To calculate the second‐order vibrational energy correction, the third and fourth‐order terms of potential function have been placed in the first‐order perturbation Hamiltonian and the second‐order Hamiltonian contains hexatic ones. We present expressions which give relations between the fourth‐ and sixth‐order terms in dimensionless normal coordinates of the potential and the anharmonicity coefficients. For illustration, a set of vibrational energies levels of SO₂, and H₂O molecules including anharmonic effects has been calculated. © 2013 Wiley Periodicals, Inc.     

Application of the Sakurai-Sugiura projection method to core-excited-state calculation by time-dependent density functional theory

The Sakurai-Sugiura projection (SS) method was implemented and numerically assessed for diagonalization of the Hamiltonian in time-dependent density functional theory (TDDFT). Since the SS method can be used to specify the range in which the eigenvalues are computed, it may be an efficient tool for use with eigenvalues in a particular range. In this article, the SS method is applied to core excited calculations for which the eigenvalues are located within a particular range, since the eigenvalues are unique to atomic species in molecules. The numerical assessment of formaldehyde molecule by TDDFT with core-valence Becke's three-parameter exchange (B3) plus Lee-Yang-Parr (LYP) correlation (CV-B3LYP) functional demonstrates that the SS method can be used to selectively obtain highly accurate eigenvalues and eigenvectors. Thus, the SS method is a new and powerful alternative for calculating core-excitation energies without high computation costs.     

New isospectral generalized potentials

In quantum chemistry, supersymmetry, shape invariance and intertwining techniques are used to determine the class of potentials that are solvable as well as to find their isospectral and generalized partners. To do that, it is necessary to have the corresponding Witten superpotential defined by W(x)=/ where is a particular wavefunction of the Hamiltonian under study. In this work, we propose an alternative way to express the Witten superpotential in terms of reciprocal wavefunctions. Thus, when this new definition of W(x) is used as an ansatz in the Riccati equation involved, one is led to a potential identical to that resulting from the use of the standard Darboux transform, which means that it is possibly the generalization of it. Moreover, the generalization of the new Witten superpotential gives rise to a new generalized isospectral potential other than that obtained from the generalized Darboux transform. As an example of an application of the proposed approach, we found the new generalized isospectral potentials that correspond to the one-dimensional free particle, harmonic oscillator and Morse potential models. Also, owing to the fact that the proposed method is general our proposal can be used straightforwardly to obtain new, exactly solvable potentials as well as to find their isospectral generalized partners which can be used advantageously in the modeling of important quantum chemical applications.From the Proceedings of the 28th Congreso de Químicos Teóricos de Expresíon Latina (QUITEL 2002) Acknowledgement.This work was supported by CONACYT-Mexico, under scientific project No. 32762-E.     

Probing the Local Structure of Pure Ionic Liquid Salts with Solid‐ and Liquid‐State NMR

Room‐temperature ionic liquids (RTILs) are gaining increasing interest and are considered part of the green chemistry paradigm due to their negligible vapour pressure and ease of recycling. Evidence of liquid‐state order, observed by IR and Raman spectroscopy, diffraction studies, and simulated by ab initio methods, has been reported in the literature. Here, quadrupolar nuclei are used as NMR probes to extract information about the solid and possible residual order in the liquid state of RTILs. To this end, the anisotropic nature and field dependence of quadrupolar and chemical shift interactions are exploited. Relaxation time measurements and a search for residual second‐order quadrupolar coupling were employed to investigate the molecular motions present in the liquid state and infer what kind of order is present. The results obtained indicate that on a timescale of ～10^?8 sec or longer, RTILs behave as isotropic liquids without residual order.     

Generalization of coupled-cluster response theory to multireference expansion spaces: application of the coupled-cluster singles and doubles effective Hamiltonian

Employing separate cluster ansatz in time-independent and time-dependent wave-operators, coupled-cluster (CC) response theory is generalized to multireference (MR) expansion spaces. For state energies, this corresponds to the MR secular problem with an arbitrary similarity-transformed effective Hamiltonian, H˜=Ω⁻¹ HΩ. The effective Hamiltonian can be generated via size-extensive CC methods. Thus the states in MR linear response theory (MRLRT) maintain the usual CC core-extensive properties. We have used the Gelfand unitary group basis of the spin-adapted configurations to construct the matrix of H˜ in the MR excitation space. As a preliminary application, the CC singles and doubles effective Hamiltonian is applied to excitation and photoionization energies of the CH⁺ and N₂ molecules, and is compared with experimental results and results from other numerical procedures including conventional CC linear response theory (CC-LRT), MR and full configuration interaction (MRCI and FCI) methods. The numerical results indicate that MRLRT reproduces valence and external excited states quantitatively, combining the best features of CC-LRT and MRCI. Received: 2 July 1998 / Accepted: 28 August 1998 / Published online: 11 November 1998     

Approximate Eigenvalue and Eigenfunction Solutions for the Generalized Hulthén Potential with any Angular Momentum

An approximate solution of the Schr?dinger equation for the generalized Hulthén potential with non-zero angular quantum number is solved. The bound state energy eigenvalues and eigenfunctions are obtained in terms of Jacobi polynomials. The Nikiforov–Uvarov method is used in the computations. We have considered the time-independent Schr?dinger equation with the associated form of Hulthén potential which simulate the effect of the centrifugal barrier for any l-state. The energy levels of the used Hulthén potential gives satisfactory values for the non-zero angular momentum as the generalized Hulthén effective potential.