Normalization and Fermi–Coulomb and Coulomb hole sum rules for approximate wave functions

For approximate wave functions, we prove the theorem that there is a one‐to‐one correspondence between the constraints of normalization and of the Fermi–Coulomb and Coulomb hole charge sum rules at each electron position. This correspondence is surprising in light of the fact that normalization depends on the probability of finding an electron at some position. In contrast, the Fermi–Coulomb hole sum rule depends on the probability of two electrons staying apart because of correlations due to the Pauli exclusion principle and Coulomb repulsion, while the Coulomb hole sum rule depends on Coulomb repulsion. We demonstrate the theorem for the ground state of the He atom by the use of two different approximate wave functions that are functionals rather than functions. The first of these wave function functionals is constructed to satisfy the constraint of normalization, and the second that of the Coulomb hole sum rule for each electron position. Each is then shown to satisfy the other corresponding sum rule. The significance of the theorem for the construction of approximate “exchange‐correlation” and “correlation” energy functionals of density functional theory is also discussed. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Iterative hypervirial-scaling method for obtaining wave functions and eigenvalues

A procedure to improve trial wave functions is given in terms of the off-diagonal hypervirial theorem. This procedure is closely related to the optimum scaling method which is valid for the diagonal hypervirial theorem. The second excited state of the one-dimensional oscillator model is employed as an example, and improved wave functions and expectation values are calculated. The results are compared with previous values based on the diagonal hypervirial theorem.     

Polymer density functional theory approach based on scaling second-order direct correlation function

A second-order direct correlation function (DCF) from solving the polymer-RISM integral equation is scaled up or down by an equation of state for bulk polymer, the resultant scaling second-order DCF is in better agreement with corresponding simulation results than the un-scaling second-order DCF. When the scaling second-order DCF is imported into a recently proposed LTDFA-based polymer DFT approach, an originally associated adjustable but mathematically meaningless parameter now becomes mathematically meaningful, i.e., the numerical value lies now between 0 and 1. When the adjustable parameter-free version of the LTDFA is used instead of the LTDFA, i.e., the adjustable parameter is fixed at 0.5, the resultant parameter-free version of the scaling LTDFA-based polymer DFT is also in good agreement with the corresponding simulation data for density profiles. The parameter-free version of the scaling LTDFA-based polymer DFT is employed to investigate the density profiles of a freely jointed tangent hard sphere chain near a variable sized central hard sphere, again the predictions reproduce accurately the simulational results. Importance of the present adjustable parameter-free version lies in its combination with a recently proposed universal theoretical way, in the resultant formalism, the contact theorem is still met by the adjustable parameter associated with the theoretical way.     

Photon Equivalents as a Parameter for Scaling Photoredox Reactions in Flow: Translation of Photocatalytic C−N Cross‐Coupling from Lab Scale to Multikilogram Scale

With the development of new photocatalytic methods over recent decades, the translation of these chemical reactions to industrial‐production scales using continuous‐flow reactors has become a topic of increasing interest. In this context, we describe our studies toward elucidating an empirically derived parameter for scaling photocatalytic reactions in flow. By evaluating the performance of a photocatalytic C?N cross‐coupling reaction across multiple reactor sizes and geometries, it was demonstrated that expressing product yield as a function of the absorbed photon equivalents provides a predictive, empirical scaling parameter. Through the use of this scaling factor and characterization of the photonic flux within each reactor, the cross‐coupling was scaled successfully from the milligram scale in batch to a multi‐kilogram reaction in flow.     

Phase separation kinetics in mixtures of flexible polymers and low molecular weight liquid crystals

A dynamic Monte Carlo (MC) simulation is performed to investigate the phase behavior of mixtures of flexible polymers and low molecular weight thermotropic liquid crystals (LCs). The polymer is represented by three‐dimensional self‐avoiding lattice chains, while the LC is described by the Lebwohl‐Lasher nematogen model. The initially homogeneous rod‐coil mixture is, following a deep quench, separated into an isotropic phase rich in coils and a nematic phase rich in rods. The underlying spinodal decomposition (SD) process is then simulated and studied extensively. This is the first simulation of SD in a rod‐coil mixture where the nematic ordering is included. Concentration fluctuations with a conserved order parameter are thus coupled with orientation fluctuations with a nonconserved order parameter. It is found that the early stage SD in the rod‐coil mixture still exhibits the dominant spatial wavelength and that the scalar scattering functions in the late stage of SD obey the Furukawa scaling law. The kinetic difference between the so‐called isotropic and anisotropic SD regions is, however, much less pronounced than predicted recently by the mean‐field theory.     

Numerical evaluation of three- and four-center bielectronic integrals using exponential-type atomic orbitals

Three- and four-center Slater orbital bielectronic integrals are evaluated by means of a complete function set. The method provides a series to approximate the bielectronic integrals. Their corresponding partial sums are analyzed in detail for 1s orbitals. The comparison with the Fourier transform–based method brings forth encouraging perspectives for the present approach. © 1995 John Wiley & Sons, Inc.     

On the hypervirial theorem and the scaling problem

It is shown that wave-functions obtained by a limited number of off-diagonal hypervirial relations are often out of scale. Optimum scaling of these functions so as to satisfy the virial theorem gives highly improved wave-functions and expectation values. The harmonic oscillator problem is treated as an example.     

A self‐guided search for good local minima of the sum‐of‐squared‐error in nonlinear least squares regression

Hard modeling of nonlinear chemical or biological systems is highly relevant as a model function together with values for model parameters provides insights in the systems' functionalities. Deriving values for said model parameters via nonlinear regression, however, is challenging as usually one of the numerous local minima of the sum‐of‐squared errors (SSEs) is determined; furthermore, for different starting points, different minima may be found. Thus, nonlinear regression is prone to low accuracy and low reproducibility. Therefore, there is a need for a generally applicable, automated initialization of nonlinear least squares algorithms, which reaches a good, reproducible solution after spending a reasonable computation time probing the SSE‐hypersurface. For this purpose, a three‐step methodology is presented in this study. First, the SSE‐hypersurface is randomly probed in order to estimate probability density functions of initial model parameter that generally lead to an accurate fit solution. Second, these probability density functions then guide a high‐resolution sampling of the SSE‐hypersurface. This second probing focuses on those model parameter ranges that are likely to produce a low SSE. As the probing continues, the most appropriate initial guess is retained and eventually utilized in a subsequent nonlinear regression. It is shown that this “guided random search” derives considerably better regression solutions than linearization of model functions, which has so far been considered the best‐case scenario. Examples from infrared spectroscopy, cell culture monitoring, reaction kinetics, and image analyses demonstrate the broad and successful applicability of this novel method. Copyright © 2014 John Wiley & Sons, Ltd.     

Hellmann–Feynman theorem near the threshold

The features of the spectrum structure are considered for situations where some parameter λ of the N‐electron Hamiltonian reaches the threshold value η under which the discrete energy level falls into the continuous spectrum. The electron density properties are also studied. It is proved that for a sequence of the wave functions converging in energy to the lower bound of the continuous spectrum as λ approaches η the corresponding sequence of the electron densities converges to the density of the (N ? 1)‐electron ground state. The results generalize the Hellmann–Feynman theorem for the cases where only the one‐side energy derivatives exist or there is no limiting wave function. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Scaled MP3 Non‐Covalent Interaction Energies Agree Closely with Accurate CCSD(T) Benchmark Data

Scaled MP3 interaction energies calculated as a sum of MP2/CBS (complete basis set limit) interaction energies and scaled third‐order energy contributions obtained in small or medium size basis sets agree very closely with the estimated CCSD(T)/CBS interaction energies for the 22 H‐bonded, dispersion‐controlled and mixed non‐covalent complexes from the S22 data set. Performance of this so‐called MP2.5 (third‐order scaling factor of 0.5) method has also been tested for 33 nucleic acid base pairs and two stacked conformers of porphine dimer. In all the test cases, performance of the MP2.5 method was shown to be superior to the scaled spin‐component MP2 based methods, e.g. SCS–MP2, SCSN–MP2 and SCS(MI)–MP2. In particular, a very balanced treatment of hydrogen‐bonded compared to stacked complexes is achieved with MP2.5. The main advantage of the approach is that it employs only a single empirical parameter and is thus biased by two rigorously defined, asymptotically correct ab‐initio methods, MP2 and MP3. The method is proposed as an accurate but computationally feasible alternative to CCSD(T) for the computation of the properties of various kinds of non‐covalently bound systems.     

Réductions chimique et électrochimique de 5H-benzopyranno[4,3-d]pyrimidines

The chemical reduction of 5H-[1]benzopyranno[4,3-d]pyrimidines with lithium aluminum hydrrde leads to 3,4-dihydro derivations. The electro-chemical reducation in acidic medium shows two monoelectronic cathodic waves. In netural or basic medium, substituted compounds in 2 position show a single bielectronic wave while two bielectronic waves are observed for unsubstituted compounds In all cases, preparative electrolysis lead to a hydrodimer in the 4,4′-positiom.     

Application of the Stabilized Koopmans’ Theorem to the Temporary Anion States of Chlorosilanes in Long‐Range Corrected Density Functional Theory

The stabilized Koopmans’ theorem (SKT) is very successful in predicting relative vertical electron attachment energies in the Hartree‐Fock theory. It is mainly accomplished by systematically scaling the most diffuse functions in the basis set. Recently, the SKT has been extended to the temporary anion states (TASs) of polyatomic molecules in the long‐range corrected density functional theory. In this paper, this method will be further applied to chlorosilanes for their importance in the chemical processes of the semiconductor industry. Their resonance energies and lifetimes are determined by computing the density of resonance states via SKT. The detailed characteristics of resonance orbitals are then analyzed. It turns out that the lowest unfilled orbitals of chlorosilanes are essentially σ*_Si‐Cl in character. Moreover, several TASs with strong Si/Cl “d” character have been identified. These results, definitely, will help us in understanding the peculiar bonding and chemical properties of chlorosilanes.     

Efficient global optimization of reactive force‐field parameters

Reactive force fields make low‐cost simulations of chemical reactions possible. However, optimizing them for a given chemical system is difficult and time‐consuming. We present a high‐performance implementation of global force‐field parameter optimization, which delivers parameter sets of the same quality with much less effort and in far less time than before, and also offers excellent parallel scaling. We demonstrate these features with example applications targeting the ReaxFF force field. © 2015 Wiley Periodicals, Inc.     

Tandem CH Activation/Arylation Catalyzed by Low‐Valent Iron Complexes with Bisiminopyridine Ligands

Tandem C?H activation/arylation between unactivated arenes and aryl halides catalyzed by iron complexes that bear redox‐active non‐innocent bisiminopyridine ligands is reported. Similar reactions catalyzed by first‐row transition metals have been shown to involve substrate‐based aryl radicals, whereas our catalytic system likely involves ligand‐centered radicals. Preliminary mechanistic investigations based on spectroscopic and reactivity studies, in conjunction with DFT calculations, led us to propose that the reaction could proceed through an inner‐sphere C?H activation pathway, which is rarely observed in the case of iron complexes. This bielectronic noble‐metal‐like behavior could be sustained by the redox‐active non‐innocent bisiminopyridine ligands.     

Dimensionless scaling methods for capillary rise

In this article the different dimensionless scaling methods for capillary rise of liquids in a tube or a porous medium are discussed. A systematic approach is taken, and the possible options are derived by means of the Buckingham π theorem. It is found that three forces (inertial, viscous and hydrostatic forces) can be used to obtain three different scaling sets, each consisting of two dimensionless variables and one dimensionless basic parameter. From a general point of view the three scaling options are all equivalent and valid for describing the problem of capillary rise. Contrary to this we find that for certain cases (depending on the time scale and the dominant forces) one of the options can be favorable. Individually the different scalings have been discussed and used in literature previously, however, we intend to discuss the three different sets systematically in a single paper and try to evaluate when which scaling is most useful. Furthermore we investigate previous analytic solutions and determine their ranges of applicability when compared to numerical solutions of the differential equation of motion (momentum balance).     

Scaling relations,virial theorem,and energy densities for long‐range and short‐range density functionals

Decomposition of the Coulomb electron–electron interaction into a long‐range and a short‐range part is described within the framework of density functional theory, deriving some scaling relations and the corresponding virial theorem. We study the behavior of the local density approximation in the high‐density limit for the long‐range and the short‐range functionals by carrying out a detailed analysis of the correlation energy of a uniform electron gas interacting via a long‐range‐only electron–electron repulsion. Possible definitions of exchange and correlation energy densities are discussed and clarified with some examples. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Benchmark values of Shannon entropy for spherically confined hydrogen atom

The information‐theoretic measure of confined hydrogen atom has been investigated extensively in the literature. However, most of them were focused on the ground state and accurate values of information entropies, such as Shannon entropy, for confined hydrogen are still not determined. In this work, we establish the benchmark results of the Shannon entropy for confined hydrogen atom in a spherical impenetrable sphere, in both position and momentum spaces. This is done by examining the bound state energies, the normalization of wave functions, and the scaling property with respect to isoelectronic hydrogenic ions. The angular and radial parts of Shannon entropy in two conjugate spaces are provided in detail for both free and confined hydrogen atom in ground and several excited states. The entropies in position space decrease logarithmically with decreasing the size of confinement, while those in momentum space increase logarithmically. The Shannon entropy sum, however, approaches to finite values when the confinement radius closes to zero. It is also found that the Shannon entropy sum shares same trend for states with similar density distributions. Variations of entropy for nodeless bound states are significantly distinct form those owning nodes when changing the confinement radius.     

Scaling in second-order electron correlation calculations using systematic sequences of even-tempered basis sets

Improved results can often be obtained from second-order Rayleigh-Schrödinger perturbation calculations of electron correlation energies using large basis sets by introducing a scaling factor in the zero-order Hamiltonian. The scaling parameter may be determined from full third-order calculations using a smaller basis set. This scaling procedure can be applied in a systematic fashion by employing a sequence of even-tempered basis sets. Calculations illustrating this approach for the beryllium atom and the neon atom are presented. The scaling procedure is also employed in conjunction with a universal systematic sequence of basis functions. Calculations illustrating this Correlation energy — Mang-body perturbation theory.Work supported in part by S.R.C. Research Grant GR/B/4738.6.S.R.C. Advanced Fellow.