Biomedical subjects

Viktor Szalay

Publications and source records attributed to Viktor Szalay.

7 recordsLinked to original sources

PubMed2007

Adiabatic Jacobi corrections for H2+-like systems.

The Coulomb three-body problem in Jacobi coordinates was solved by treating the distance of the particles having equal charge as a parameter. This method allows computation of electronic energies with finite nuclear masses while maintaining the notion of a potential energy curve. The rotationless ground-state electronic and the so-called adiabatic Jacobi correction (AJC) energies are presented for H2+, D2+, and HD+ at fixed internuclear separations. The AJCs are defined as the difference between the results obtained from calculations using proper finite and infinite nuclear masses. Except at the united atom limit, the AJCs are smaller than the traditional first-order diagonal Born-Oppenheimer corrections. Expectation values of proton-electron, p-e, and deuteron-electron, d-e, distances for HD+ have been computed as a function of internuclear separation. Similarly to the fully nonadiabatic approach, the present method is able to follow the symmetry breaking in HD+. Exact and approximate analytical and numerical results are given for counterfactual systems as well. In these cases changes are allowed for the values of the electron rest mass or the elementary charge, as well as for the mass or charge of the unique particle (electron).

Journal Article↗

PubMed2006

Optimal grids for generalized finite basis and discrete variable representations: definition and method of calculation.

The method of optimal generalized finite basis and discrete variable representations (FBR and DVR) generalizes the standard, Gaussian quadrature grid-classical orthonormal polynomial basis-based FBR/DVR method to general sets of grid points and to general, nondirect product, and/or nonpolynomial bases. Here, it is shown how an optimal set of grid points can be obtained for an optimal generalized FBR/DVR calculation with a given truncated basis. Basis set optimized and potential optimized grids are defined. The optimized grids are shown to minimize a function of grid points derived by relating the optimal generalized FBR of a Hamiltonian operator to a non-Hermitian effective Hamiltonian matrix. Locating the global minimum of this function can be reduced to finding the zeros of a function in the case of one dimensional problems and to solving a system of D nonlinear equations repeatedly in the case of D>1 dimensional problems when there is an equal number of grid points and basis functions. Gaussian quadrature grids are shown to be basis optimized grids. It is demonstrated by a numerical example that an optimal generalized FBR/DVR calculation of the eigenvalues of a Hamiltonian operator with potential optimized grids can have orders of magnitude higher accuracy than a variational calculation employing the same truncated basis. Nevertheless, for numerical integration with the optimal generalized FBR quadrature rule basis optimized grids are the best among grids of the same number of points. The notions of Gaussian quadrature and Gaussian quadrature accuracy are extended to general, multivariable basis functions.

Journal Article↗

PubMed2006

Adiabatic approximations to internal rotation.

A number of subtle and confusing issues are addressed concerning large amplitude motion (LAM) coordinates (chi) for internal molecular motions, using the methyl rotation in acetaldehyde (CH(3)CHO) as a model problem. If the LAM coordinate is chosen to be one of the H-C-C-O dihedral angles rho(1), rho(2), or rho(3), it lacks the required 2pi3 periodicity, and its use is thus undesirable. An excellent local internal coordinate for this model problem is tau(3)=13(rho(1)+rho(2)+rho(3)-2pi). A similarly good but nonlocal coordinate for the adiabatic approximation of internal rotation is provided by the intrinsic reaction coordinate s. Comparison of the mass-independent V(0)(tau(3)) and the mass-dependent V(0)(s) internal rotation curves shows that the two are virtually identical for the parent isotopolog of acetaldehyde. A unified internal coordinate projection scheme for determining complementary vibrational frequencies and subsequently V(ZPVE)(chi) along a path for LAM has been formulated, where V(ZPVE)(chi) is the zero-point vibrational energy correction to the internal rotation curve. In addition to its simplicity, the projection scheme developed for a distinguished reaction path generated by constrained optimizations is appealing because the vibrational frequencies along the LAM path are invariant to chemically meaningful choices of the internal coordinates for the complementary modes.

Journal Article↗

PubMed2006

Finite basis representations with nondirect product basis functions having structure similar to that of spherical harmonics.

The currently most efficient finite basis representation (FBR) method [Corey et al., in Numerical Grid Methods and Their Applications to Schrodinger Equation, NATO ASI Series C, edited by C. Cerjan (Kluwer Academic, New York, 1993), Vol. 412, p. 1; Bramley et al., J. Chem. Phys. 100, 6175 (1994)] designed specifically to deal with nondirect product bases of structures phi(n) (l)(s)f(l)(u), chi(m) (l)(t)phi(n) (l)(s)f(l)(u), etc., employs very special l-independent grids and results in a symmetric FBR. While highly efficient, this method is not general enough. For instance, it cannot deal with nondirect product bases of the above structure efficiently if the functions phi(n) (l)(s) [and/or chi(m) (l)(t)] are discrete variable representation (DVR) functions of the infinite type. The optimal-generalized FBR(DVR) method [V. Szalay, J. Chem. Phys. 105, 6940 (1996)] is designed to deal with general, i.e., direct and/or nondirect product, bases and grids. This robust method, however, is too general, and its direct application can result in inefficient computer codes [Czako et al., J. Chem. Phys. 122, 024101 (2005)]. It is shown here how the optimal-generalized FBR method can be simplified in the case of nondirect product bases of structures phi(n) (l)(s)f(l)(u), chi(m) (l)(t)phi(n) (l)(s)f(l)(u), etc. As a result the commonly used symmetric FBR is recovered and simplified nonsymmetric FBRs utilizing very special l-dependent grids are obtained. The nonsymmetric FBRs are more general than the symmetric FBR in that they can be employed efficiently even when the functions phi(n) (l)(s) [and/or chi(m) (l)(t)] are DVR functions of the infinite type. Arithmetic operation counts and a simple numerical example presented show unambiguously that setting up the Hamiltonian matrix requires significantly less computer time when using one of the proposed nonsymmetric FBRs than that in the symmetric FBR. Therefore, application of this nonsymmetric FBR is more efficient than that of the symmetric FBR when one wants to diagonalize the Hamiltonian matrix either by a direct or via a basis-set contraction method. Enormous decrease of computer time can be achieved, with respect to a direct application of the optimal-generalized FBR, by employing one of the simplified nonsymmetric FBRs as is demonstrated in noniterative calculations of the low-lying vibrational energy levels of the H3+ molecular ion. The arithmetic operation counts of the Hamiltonian matrix vector products and the properties of a recently developed diagonalization method [Andreozzi et al., J. Phys. A Math. Gen. 35, L61 (2002)] suggest that the nonsymmetric FBR applied along with this particular diagonalization method is suitable to large scale iterative calculations. Whether or not the nonsymmetric FBR is competitive with the symmetric FBR in large-scale iterative calculations still has to be investigated numerically.

Journal Article↗

PubMed2005

On equilibrium structures of the water molecule.

Equilibrium structures are fundamental entities in molecular sciences. They can be inferred from experimental data by complicated inverse procedures which often rely on several assumptions, including the Born-Oppenheimer approximation. Theory provides a direct route to equilibrium geometries. A recent high-quality ab initio semiglobal adiabatic potential-energy surface (PES) of the electronic ground state of water, reported by Polyansky et al. [ ibid. 299, 539 (2003)] and called CVRQD here, is analyzed in this respect. The equilibrium geometries resulting from this direct route are deemed to be of higher accuracy than those that can be determined by analyzing experimental data. Detailed investigation of the effect of the breakdown of the Born-Oppenheimer approximation suggests that the concept of an isotope-independent equilibrium structure holds to about 3 x 10(-5) A and 0.02 degrees for water. The mass-independent [Born-Oppenheimer (BO)] equilibrium bond length and bond angle on the ground electronic state PES of water is r(e) (BO)=0.957 82 A and theta e (BO)=104.48(5) degrees , respectively. The related mass-dependent (adiabatic) equilibrium bond length and bond angle of H2 (16)O is r(e) (ad)=0.957 85 A and theta e (ad)=104.50(0) degrees , respectively, while those of D2 (16)O are r(e) (ad)=0.957 83 A and theta e (ad)=104.49(0) degrees . Pure ab initio prediction of J=1 and 2 rotational levels on the vibrational ground state by the CVRQD PESs is accurate to better than 0.002 cm(-1) for all isotopologs of water considered. Elaborate adjustment of the CVRQD PESs to reproduce all observed rovibrational transitions to better than 0.05 cm(-1) (or the lower ones to better than 0.0035 cm(-1)) does not result in noticeable changes in the adiabatic equilibrium structure parameters. The expectation values of the ground vibrational state rotational constants of the water isotopologs, computed in the Eckart frame using the CVRQD PESs and atomic masses, deviate from the experimentally measured ones only marginally, especially for A0 and B0. The small residual deviations in the effective rotational constants are due to centrifugal distortion, electronic, and non-Born-Oppenheimer effects. The spectroscopic (nonadiabatic) equilibrium structural parameters of H2 16O, obtained from experimentally determined A'0 and B'0 rotational constants corrected empirically to obtain equilibrium rotational constants, are r(e) (sp)=0.957 77 A and theta e (sp)=104.48 degrees .

Journal Article↗

PubMed2005

Treating singularities present in the Sutcliffe-Tennyson vibrational Hamiltonian in orthogonal internal coordinates.

Two methods are developed, when solving the related time-independent Schrodinger equation (TISE), to cope with the singular terms of the vibrational kinetic energy operator of a triatomic molecule given in orthogonal internal coordinates. The first method provides a mathematically correct treatment of all singular terms. The vibrational eigenfunctions are approximated by linear combinations of functions of a three-dimensional nondirect-product basis, where basis functions are formed by coupling Bessel-DVR functions, where DVR stands for discrete variable representation, depending on distance-type coordinates and Legendre polynomials depending on angle bending. In the second method one of the singular terms related to a distance-type coordinate, deemed to be unimportant for spectroscopic applications, is given no special treatment. Here the basis set is obtained by taking the direct product of a one-dimensional DVR basis with a two-dimensional nondirect-product basis, the latter formed by coupling Bessel-DVR functions and Legendre polynomials. With the basis functions defined, matrix representations of the TISE are set up and solved numerically to obtain the vibrational energy levels of H3+. The numerical calculations show that the first method treating all singularities is computationally inefficient, while the second method treating properly only the singularities having physical importance is quite efficient.

Journal Article↗

PubMed2004

Ab initio torsional potential and transition frequencies of acetaldehyde.

High-level ab initio electronic structure calculations, including extrapolations to the complete basis set limit as well as relativistic and diagonal Born-Oppenheimer corrections, resulted in a torsional potential of acetaldehyde in its electronic ground state. This benchmark-quality potential fully reflects the symmetry and internal rotation dynamics of this molecule in the energy range probed by spectroscopic experiments in the infrared and microwave regions. The torsional transition frequencies calculated from this potential and the ab initio torsional inverse effective mass function are within 2 cm(-1) of the available experimental values. Furthermore, the computed contortional parameter rho of the rho-axis system Hamiltonian is also in excellent agreement with that obtained from spectral analyses of acetaldehyde.

Journal Article↗