Electronic structure modulation of metal–organic frameworks (MOFs) through the connection of linker “wires” as a function of an external stimulus is reported for the first time. The established correlation between MOF electronic properties and photoisomerization kinetics as well as changes in an absorption profile is unprecedented for extended well-defined structures containing coordinatively integrated photoresponsive linkers. The presented studies were carried out on both single crystal and bulk powder with preservation of framework integrity. An LED-containing electric circuit, in which the switching behavior was driven by the changes in MOF electronic profile, was built for visualization of experimental findings. The demonstrated concept could be used as a blueprint for development of stimuli-responsive materials with dynamically controlled electronic behavior.

Solution of the Schrödinger 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.

We introduce quasirandom distributed Gaussian bases ͑QDGB͒ that are well suited for bound problems. The positions of the basis functions are chosen quasirandomly while their widths and density are functions of the potential. The basis function overlap and kinetic energy matrix elements are analytical. The potential energy matrix elements are accurately evaluated using few-point quadratures, since the Gaussian basis functions are localized. The resulting QDGB can be easily constructed and is shown to be accurate and efficient for eigenvalue calculation for several multidimensional model vibrational problems. As more demanding examples, we used a 2D QDGB-DVR basis to calculate the lowest 400 or so energy levels of the water molecule for zero total angular momentum to sub-wave-number precision. Finally, the lower levels of Ar 3 and Ne 3 were calculated using a symmetrized QDGB. The QDGB was shown to be accurate with a small basis.

Calculation of chemical reaction dynamics is central to theoretical chemistry. The majority of calculations use either classical mechanics, which is computationally inexpensive but misses quantum effects, such as tunneling and interference, or quantum mechanics, which is computationally expensive and often conceptually opaque. An appealing middle ground is the use of semiclassical mechanics. Indeed, since the early 1970s there has been great interest in using semiclassical methods to calculate reaction probabilities. However, despite the elegance of classical S-matrix theory, numerical results on even the simplest reactive systems remained out of reach. Recently, with advances both in correlation function formulations of reactive scattering as well as in semiclassical methods, it has become possible for the first time to calculate reaction probabilities semiclassically. The correlation function methods are contrasted with recent flux-based methods, which, although providing somewhat more compact expressions for the cumulative reactive probability, are less compatible with semiclassical implementation.

No abstract

The quantum trajectory approach is generalized to arbitrary coordinate systems, including curvilinear coordinates. This allows one to perform an approximate quantum trajectory propagation, which scales favorably with system size, in the same framework as standard quantum wave packet dynamics. The trajectory formulation is implemented in Jacobi coordinates for a nonrotating triatomic molecule. Wave packet reaction probabilities are computed for the O(3P) + H2 --> OH + H reaction using the approximate quantum potential. The latter is defined by the nonclassical component of the momentum operator expanded in terms of linear and exponential functions. Unlike earlier implementations with linear functions, the introduction of the exponential function provides an accurate description of asymptotic dynamics for this system and gives good agreement of approximate reaction probabilities with accurate quantum calculations.

In the de Broglie-Bohm formulation of quantum mechanics the time-dependent Schrödinger 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.

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.