Entangled trajectory dynamics in the Husimi representation

We solve quantum dynamical equations of simple systems by propagating ensembles of interacting trajectories. A scheme is proposed which uses adaptive kernel density estimation for representing probability distribution functions and their derivatives. The formulation is carried on in the Husimi representation to ensure the positiveness of the distribution functions. By comparing to previous work, the effect of changing representations is studied as well as the advantage of using adaptive kernels for the estimation of probability distributions. We found significant improvement in the accuracy of the results.     

Constraints in the variational principle for stationary states

Constraints in the variational principle for stationary states (VPSS) are classified in accordance with Dirac’s constrained classical mechanics and the time-dependent variational principle (TDVP). All of the VPSS constraints are required to belong to the first-class TDVP as constants of motion to ensure the real-valuedness of the Lagrange multipliers. The VPSS constraints are further classified as either first-class or second-class. The first-class VPSS constraints are constants of variation with symmetry-adapted wave functions. If the representation basis for the constraint operators is incomplete, however, the first-class VPSS constraints lead to broken-symmetry solutions. The nondegenerate energies of \({}^2E'\) at the \(D_{3h}\) geometry in the Jahn–Teller distortion of H\(_3\) are presented as an example of a problem with broken-symmetry. An ad hoc adjustment is suggested by considering the second-class pseudo-VPSS constraints with new adiabatic parameters.     

Application of Wigner and Husimi intracule based electron correlation models to excited states

A new approach to the electron correlation problem based on phase space intracules derived from the Wigner distribution is applied to excited states. The computed electron correlation energy reduces the mean absolute error in the prediction of the excitation energies of 55 atomic excited states from 0.65 eV for unrestricted Hartree-Fock to 0.32 eV. This compares favorably to a mean absolute deviation of 0.52 eV for second order Moller-Plesset perturbation theory and 0.35 eV for the Lee-Yang-Parr functional. An analogous correlation model based on the Husimi distribution is developed. Predicted correlation energies and excitation energies from this model are significantly worse than for the Wigner intracule based model. Alternative correlation kernels may be more suitable for the Husimi intracule based approach.     

Construction of stationary quantum states with targeted energies

Given an exact stationary state and its energy-eigenvalue, we devise a simple way of generating new families of exact eigenstates with the same eigenenergy, but for different potentials. While this recipe of designing isoergic states is quite general in the context of arriving at new exact solutions from a known premise, relevance of this route with the scheme of construction of isospectral potentials is noted under specific situations. The idea is extended to finding new exact eigenstates with eigenvalues as the sum of energies of two or more separate but known stationary states in dissimilar potentials. Implication of the latter endeavor in the contexts of manifold energy minimization is discussed.     

Intermediate state representation approach to physical properties of dicationic states

The second-order algebraic construction (ADC(2)) approach to the two-particle (pp) propagator, devised to compute double ionization energies and associated spectroscopic amplitudes, is reformulated and extended using the concept of intermediate state representations (ISR). The ISR formulation allows one to go beyond the general limitations inherent to the propagator approach, as here (N-2)-electron wave functions and properties become directly accessible. The (N-2)-electron ISR(2) equations for a general one-particle operator have been derived and implemented in a recent version of the double ionization ADC(2) program. As a first test of the method, the dipole moments of a series of 2h states of LiH, HF, and H(2)O were computed and compared to the results of a full configuration interaction (FCI) treatment. The dipole moments obtained at the ADC(2)/ISR(2) computational level are in good agreement with the FCI results.     

Kinetic insights over a PEMFC operating on stationary and oscillatory states

Kinetic investigations in the oscillatory state have been carried out in order to shed light on the interplay between the complex kinetics exhibited by a proton exchange membrane fuel cell fed with poisoned H(2) (108 ppm of CO) and the other in serie process. The apparent activation energy (E(a)) in the stationary state was investigated in order to clarify the E(a) observed in the oscillatory state. The apparent activation energy in the stationary state, under potentiostatic control, rendered (a) E(a) ≈ 50-60 kJ mol(-1) over 0.8 V < E < 0.6 V and (b) E(a) ≈ 10 kJ mol(-1) at E = 0.3 V. The former is related to the H(2) adsorption in the vacancies of the surface poisoned by CO and the latter is correlated to the process of proton conductivity in the membrane. The dependence of the period-one oscillations on the temperature yielded a genuine Arrhenius dependence with two E(a) values: (a) E(a) around 70 kJ mol(-1), at high temperatures, and (b) E(a) around 10-15 kJ mol(-1), at lower temperatures. The latter E(a) indicates the presence of protonic mass transport coupled to the essential oscillatory mechanism. These insights point in the right direction to predict spatial couplings between anode and cathode as having the highest strength as well as to speculate the most likely candidates to promote spatial inhomogeneities.     

A self-consistent field procedure for stationary states using an algebraic approach and the maximum entropy formalism

A stationary state of maximal entropy is derived as a solution of a variational procedure. Generators of a continuous group are used as the constraints. The self-consistent hamiltonian is linear in these generators so that the solution of the self-consistency problem is replaced by a solution of an algebraic equation. The familiar Hartree-Fock procedure is a special case.     

One dimensional stationary states with constant electrical induction in ferroelectric liquid crystals

In a ferroelectric liquid crystal, it is the electrical induction, D, and not the electric field, E, which is independent of the space coordinate in a one dimensional analysis. This causes, at small applied voltages, stationary states which have a homogeneous orientation and are field free over most of the liquid crystal. At higher voltages homogeneous field states with zero, one or two surface layers may exist.     

Spectra of Husimi cacti: exact results and applications

Starting from exact relations for finite Husimi cacti we determine their complete spectra to very high accuracy. The Husimi cacti are dual structures to the dendrimers but, distinct from these, contain loops. Our solution makes use of a judicious analysis of the normal modes. Although close to those of dendrimers, the spectra of Husimi cacti differ. From the wealth of applications for measurable quantities which depend only on the spectra, we display for Husimi cacti the behavior of the fluorescence depolarization under quasiresonant Forster energy transfer.     

Ab initio many-electron study for the low-lying states of the alkali hydride cations in the adiabatic representation

An ab initio multireference single- and double-excitation configuration interaction (CI) study is carried out for the ground and excited electronic states of alkali-hydride cations (LiH(+), NaH(+), KH(+), RbH(+), and CsH(+)). For all alkali-metal atoms, the first inner-shell and valence electrons (nine active electrons, three for Li) are considered explicitly in the ab initio self-consistent-field and CI calculations. The adiabatic potential energy curves, radial and rotational couplings are calculated and presented. Short-range (～3 a.u.) potential wells produced by the excitation of the inner-shell electrons are found. The depths of the inner potential wells are much greater than those of the outer wells for the CsH(+) system. The computed spectroscopic constants for the long-range potential well of the 2 (2)Σ(+) state are very close to the available theoretical and experimental data. The electronic states of alkali-hydrogen cations are also compared with each other, it is found that the positions of the potential wells shift to larger internuclear distances gradually, and the depths of these potential wells become greater with increasing alkali-metal atomic number. The relationships between structures of the radial coupling matrix elements and the avoiding crossings of the potential curves are analyzed. From NaH(+) to CsH(+), radial coupling matrix elements display more and more complex structures due to the gradual decrease of energy separations for avoided crossings. Finally, the behavior of some rotational couplings is also shown.     

Grassmann coherent states representation of the path integral: Evaluation of the generating function for spin systems

An interacting spin system is investigated within the scenario of the Feynman path integral representation of quantum mechanics. Short‐time propagator algorithms and a discrete time formalism are used in combination with a basis set involving Grassmann variables coherent states to get a many‐body analytic propagator. The generating function thus obtained leads, after an adequate tracing over Grassmann variables in the imaginary time domain, to the partition function. A spin 1/2 Hamiltonian involving the whole set of interactions is considered. Fermion operators satisfying the standard anticommutation relations are constructed from the raising and lowering spin operators via the Jordan–Wigner transformation. The partition function obtained is more general than the partition function of the traditional Ising model involving only first‐neighbor interactions. Computations were performed assuming that the coupling as a function of the distance can be reasonably well represented by an Airy function. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Comments about the representation of Rydberg and ionic excited states and their photochemistry

This paper discusses the conditions for representing Rydberg and ionic excited states of molecules. It especially shows the intrinsic difficulties of MO methods to treat the weak resonances between strongly polarized situations in highly polarizable symmetric systems; such situations occur in the long distance region for Rydberg excited states of homonuclear molecules, and for the 90 ° twisted singlet excited states of polyenes. The valence/Rydberg mixing is discussed, and some principles for the understanding of Rydberg photochemistry are proprosed, based on a few examples. The present knowledge of the photochemistry of zwitterionic excited states of polyenes is summarized.     

Explicit representation of Gel'fand–Tsetlin states in clifford algebra unitary group approach

A explicit expression for the unitary group Clebsch–Gordan coefficients, which couple two fully antisymmetric single-column states into the two-column Gel'fand–Tsetlin states, is given in terms of isoscalar factors for the canonical subgroup chain U(n) ? U(n – 1) ? …? ? U(1). The isoscalar factors are expressed through the step numbers labeling canonical basis states and enable a straightforward construction of Gel'fand–Tsetlin states in the Clifford algebra unitary group approach, without the use of the tables for the symmetric group outer-product reduction coefficients.     

Geometric constraints in semiclassical initial value representation calculations in Cartesian coordinates: excited states

The authors show that a recently proposed approach [J. Chem. Phys. 123, 084103 (2005)] for the inclusion of geometric constraints in semiclassical initial value representation calculations can be used to obtain excited states of weakly bound complexes. Sample calculations are performed for free and constrained rare gas clusters. The results show that the proposed approach allows the evaluation of excited states with reasonable accuracy when compared to exact basis set calculations.     

Exact solution of a RNA-like polymer model on the Husimi lattice

We investigate a two-tolerant polymer model on the square Husimi lattice, which aims at describing the properties of RNA-like macromolecules. We solve the model in a numerically exact way, working out the grand-canonical phase diagram, both with and without taking into account the stacking effect. Besides a nonpolymerized phase, we observe two different polymerized phases characterized by a lower or higher density of doubly visited lattice bonds. The system exhibits three qualitatively different regimes, as a function of the monomer chemical potential. Below some T1 temperature and above some T2 temperature, the transition to the nonpolymerized phase is continuous, whereas, in the (T1,T2) temperature range, the transition is first order. In the dilute-solution limit, the high temperature regime corresponds to a swollen ("coil") state, the intermediate regime to a moderately collapsed ("molten") state, with a small fraction of paired segments, and the low temperature regime to an almost fully paired ("native") state. The molten state ends in a tricritical (Theta-like) transition at high temperature and in a critical end point at low temperature. Upon increasing the stacking energy parameter, the temperature range of the molten state turns out to be progressively reduced but never completely removed.     

Discrete variable representation method applied to the determination of rotation-vibration bound states of NO2

The discrete variable representation method is applied to the determination of the rotation-vibration energy levels of the fundamental electronic state of NO₂. The Hamiltonian is expressed in Johnson hyperspherical coordinates and developed on a DVR basis for each internal coordinate, while parity-adapted linear combinations of Wigner functions are used to describe the rotational motion. The diagonalization of the Hamiltonian matrix is performed using the Lanczos algorithm for large symmetric and Hermitian matrices. Results for rovibrational states up to J = 11 for the first five vibrational energy levels are presented. © 1997 John Wiley & Sons, Inc.