Biomedical subjects

Sophya Garashchuk

Publications and source records attributed to Sophya Garashchuk.

4 recordsLinked to original sources

PubMed2004

Modified quantum trajectory dynamics using a mixed wave function representation.

Dynamics of quantum trajectories provides an efficient framework for description of various quantum effects in large systems, but it is unstable near the wave function density nodes where the quantum potential becomes singular. A mixed coordinate space/polar representation of the wave function is used to circumvent this problem. The resulting modified trajectory dynamics associated with the polar representation is nonsingular and smooth. The interference structure and the nodes of the wave function density are described, in principle, exactly in the coordinate representation. The approximate version of this approach is consistent with the semiclassical linearized quantum force method [S. Garashchuk and V. A. Rassolov, J. Chem. Phys. 120, 1181 (2004)]. This approach is exact for general wave functions with the density nodes in a locally quadratic potential.

Journal Article↗

PubMed2004

Geminal model chemistry II. Perturbative corrections.

We introduce and investigate a chemical model based on perturbative corrections to the product of singlet-type strongly orthogonal geminals wave function. Two specific points are addressed (i) Overall chemical accuracy of such a model with perturbative corrections at a leading order; (ii) Quality of strong orthogonality approximation of geminals in diverse chemical systems. We use the Epstein-Nesbet form of perturbation theory and show that its known shortcomings disappear when it is used with the reference Hamiltonian based on strongly orthogonal geminals. Application of this model to various chemical systems reveals that strongly orthogonal geminals are well suited for chemical models, with dispersion interactions between the geminals being the dominant effect missing in the reference wave functions.

Journal Article↗

PubMed2004

Bohmian dynamics on subspaces using linearized quantum force.

In the de Broglie-Bohm formulation of quantum mechanics the time-dependent Schrodinger equation is solved in terms of quantum trajectories evolving under the influence of quantum and classical potentials. For a practical implementation that scales favorably with system size and is accurate for semiclassical systems, we use approximate quantum potentials. Recently, we have shown that optimization of the nonclassical component of the momentum operator in terms of fitting functions leads to the energy-conserving approximate quantum potential. In particular, linear fitting functions give the exact time evolution of a Gaussian wave packet in a locally quadratic potential and can describe the dominant quantum-mechanical effects in the semiclassical scattering problems of nuclear dynamics. In this paper we formulate the Bohmian dynamics on subspaces and define the energy-conserving approximate quantum potential in terms of optimized nonclassical momentum, extended to include the domain boundary functions. This generalization allows a better description of the non-Gaussian wave packets and general potentials in terms of simple fitting functions. The optimization is performed independently for each domain and each dimension. For linear fitting functions optimal parameters are expressed in terms of the first and second moments of the trajectory distribution. Examples are given for one-dimensional anharmonic systems and for the collinear hydrogen exchange reaction.

Journal Article↗

PubMed2004

Energy conserving approximations to the quantum potential: dynamics with linearized quantum force.

Solution of the Schrodinger equation within the de Broglie-Bohm formulation is based on propagation of trajectories in the presence of a nonlocal quantum potential. We present a new strategy for defining approximate quantum potentials within a restricted trial function by performing the optimal fit to the log-derivatives of the wave function density. This procedure results in the energy-conserving dynamics for a closed system. For one particular form of the trial function leading to the linear quantum force, the optimization problem is solved analytically in terms of the first and second moments of the weighted trajectory distribution. This approach gives exact time-evolution of a correlated Gaussian wave function in a locally quadratic potential. The method is computationally cheap in many dimensions, conserves total energy and satisfies the criterion on the average quantum force. Expectation values are readily found by summing over trajectory weights. Efficient extraction of the phase-dependent quantities is discussed. We illustrate the efficiency and accuracy of the linear quantum force approximation by examining a one-dimensional scattering problem and by computing the wavepacket reaction probability for the hydrogen exchange reaction and the photodissociation spectrum of ICN in two dimensions.

Journal Article↗