A central limit theorem for age- and density-dependent population processes

Consider a population consisting of one type of individual living in a fixed region with area A. In [8], we constructed a stochastic population model in which the death rate is affected by the age of the individual and the birth rate is affected by the population density P_A(t), i.e., the population size divided by the area A of the given region. In [8], we proposed a continuous deterministic model which in general is a nonlinear Volterra type integral equation and proved that under appropriate conditions the sequence P_A(t) would converge to the solution P(t) of our integral equation in the sense that

$\text{lim}\text{→∞}\text{P}\text{sup}\text{0?s?t}\text{|P}{}_{\mathrm{A}}\text{(s) ? P(s)|>ε}\text{=0}\text{for every}\text{ε > 0}$

.In this paper, we obtain a “central limit theorem” for the random element √A(P_A(t)?P(t)). We prove that under appropriate conditions √A(P_A(t)?P(t)) will converge to a Gaussian process. (See Theorem 3.4 for the explicit formula of this Gaussian process.)     

邻位吡啶自由基多通道分解反应的动力学和反应机理研究

Utilizing Gaussian94 program package, all species involved in decomposition reactions of o-pyridyl radical were optimized fully at B3LYP/6-311＋＋G^＊＊ level. Intrinsic reaction coordinate calculations were employed to confirm the connections of the transition states and products, and transition states were ascertained by the number of imaginary frequency （0 or 1）. The reaction mechanism was elucidated by the vibrational mode analysis and electronic population analysis, and the reaction rate constants were calculated with transition state theory.     

Improving the efficiency and reliability of free energy perturbation calculations using overlap sampling methods

A challenge in free energy calculation for complex molecular systems by computer simulation is to obtain a reliable estimate within feasible computational time. In this study, we suggest an answer to this challenge by exploring a simple method, overlap sampling (OS), for producing reliable free-energy results in an efficient way. The formalism of the OS method is based on ensuring sampling of important overlapping phase space during perturbation calculations. This technique samples both forward and reverse free energy perturbation (FEP) to improve the free-energy calculation. It considers the asymmetry of the FEP calculation and features an ability to optimize both the precision and the accuracy of the measurement without affecting the simulation process itself. The OS method is tested at two optimization levels: no optimization (simple OS), and full optimization (equivalent to Bennett's method), and compared to conventional FEP techniques, including the widely used direct FEP averaging method, on three alchemical mutation systems: (a) an anion transformation in water solution, (b) mutation between methanol and ethane, and (c) alchemical change of an adenosine molecule. It is consistently shown that the reliability of free-energy estimates can be greatly improved using the OS techniques at both optimization levels, while the performance of Bennett's method is particularly striking. In addition, the efficiency of a calculation can be significantly improved because the method is able to (a) converge to the right answer quickly, and (b) work for large perturbations. The basic two-stage OS method can be extended to admit additional stages, if needed. We suggest that the OS method can be used as a general perturbation technique for computing free energy differences in molecular simulations.     

Semiempirical local spin: Theory and implementation of the ZILSH method for predicting Heisenberg exchange constants of polynuclear transition metal complexes

The local spin formalism ( 3 ) for computing expectation values 〈S_A · S_B〉 that appear in the Heisenberg spin model has been extended to semiempirical single determinant wave functions. An alternative derivation of expectation values in restricted and unrestricted cases is given that takes advantage of the zero differential overlap (ZDO) approximation. A formal connection between single determinant wave functions (which are not in general spin eigenfunctions) and the Heisenberg spin model was established by demonstrating that energies of single determinants that are eigenfunctions of the local spin operators with eigenvalues corresponding to high‐spin radical centers are given by the same Heisenberg coupling constants {J_AB} that describe the true spin states of the system. Unrestricted single determinant wave functions for transition metal complexes are good approximations of local spin eigenfunctions when the metal d orbitals are local in character and all unpaired electrons on each metal have the same spin (although spins on different metals might be reversed). Good approximations of the coupling constants can then be extracted from local spin expectation values 〈S_A · S_B〉 energies of the single determinant wave functions. Once the coupling constants are obtained, diagonalization of the Heisenberg spin Hamiltonian provides predictions of the energies and compositions of the spin states. A computational method is presented for obtaining coupling constants and spin‐state energies in this way for polynuclear transition metal complexes using the intermediate neglect of differential overlap Hamiltonian parameterized for optical spectroscopy (INDO/S) in the ZINDO program. This method is referred to as ZILSH, derived from ZINDO, Davidson's local spin formalism, and the Heisenberg spin model. Coupling constants and spin ground states obtained for 10 iron complexes containing from 2 to 6 metals are found to agree well with experimental results in most cases. In the case of the complex [Fe₆O₃(OAc)₉(OEt)₂(bpy)₂]⁺, a priori predictions of the coupling constants yield a ground‐state spin of zero, in agreement with variable‐temperature magnetization data, and corroborate spin alignments proposed earlier on the basis of structural considerations. This demonstrates the potential of the ZILSH method to aid in understanding magnetic interactions in polynuclear transition metal complexes. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

Survey of ligand field parameters of strong field d complexes obtained from the g matrix

Using a newly proposed approach involving an internally consistent set of equations, the ligand field parameters Δ/ξ, V/ξ and k are obtained from literature values of the g matrix for strong field d⁵ systems of various conformations in which |Δ/ξ|≤10. Qualitative analysis of the observed results is done using the Angular Overlap Model, AOM.     

A theoretical investigation into the conformational changes of dibenzo-p-dioxin, thianthrene, and selenanthrene

Theoretical ab initio calculations using the HF and B3LYP methods have been carried out to investigate the conformational differences of three cyclic rings, dibenzo-p-dioxin (DPD), thianthrene (THT), and selenanthrene (SET). The physical origin for the conformational preference of each molecule has been studied using the natural bond orbital (NBO) analysis. The NBO results indicate that DPD exists in a planar form due to strong electron delocalization caused by the specific orbital interaction, around the X atom. On the other hand, THT and SET exist as puckered forms with high inversion barriers due to less effective electron delocalization. The NBO analysis also shows that the conformational stabilization in DPD is caused by a more effective overlap of the orbitals, compared with the overlap of the orbitals in THT.     

Analysis of B3LYP and MP2 conformational population distributions in trans‐nicotine,acetylcholine, and ABT‐594

This work presents an analysis of the equivalence of MP2 and DFT (B3LYP functional) conformational populations. As a test case, we select three cholinergic agents (trans‐nicotine, acetylcholine, and the nicotinic analgesic ABT‐594), where the minima on the conformational energy hypersurfaces expand a large range of energies (～0–30 kJ mol^?1). From energetic and structural data obtained in vacuo at the MP2 and B3LYP/cc‐pVDZ levels, we build conformational partition functions, including the effect of the conformational kinetic energy and the rotovibrational coupling. Our results at a physiological temperature (37°C) show qualitative agreement in all cases. Quantitative agreement, however, is only found for trans‐nicotine and ABT‐594. In the first case, energy minima differ by <0.2 kJ mol^?1. Therefore, the equivalence of structural results translates in the equivalence of the conformational distribution. For ABT‐594, the minima are separated by as much as 8.0 kJ mol^?1, and the conformational energy determines the conformational distribution. In this case, the slight relative variation of conformational energy, between B3LYP and MP2, does not affect the population, since the secondary minima are high in energy and very low in population. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Estimating reaction parameters in mechanism-enabled population balance models of nanoparticle size distributions: A Bayesian inverse problem approach

In order to quantitatively predict nano- as well as other particle-size distributions, one needs to have both a mathematical model and estimates of the parameters that appear in these models. Here, we show how one can use Bayesian inversion to obtain statistical estimates for the parameters that appear in recently derived mechanism-enabled population balance models (ME-PBM) of nanoparticle growth. The Bayesian approach addresses the question of “how well do we know our parameters, along with their uncertainties?.” The results reveal that Bayesian inversion statistical analysis on an example, prototype nanoparticle formation system allows one to estimate not just the most likely rate constants and other parameter values, but also their SDs, confidence intervals, and other statistical information. Moreover, knowing the reliability of the mechanistic model's parameters in turn helps inform one about the reliability of the proposed mechanism, as well as the reliability of its predictions. The paper can also be seen as a tutorial with the additional goal of achieving a “Gold Standard” Bayesian inversion ME-PBM benchmark that others can use as a control to check their own use of this methodology for other systems of interest throughout nature. Overall, the results provide strong support for the hypothesis that there is substantial value in using a Bayesian inversion methodology for parameter estimation in particle formation systems.     

Calculation of the vibration–rotational transition intensities of water molecules trapped in an argon matrix: stretching O–H vibrations spectral region

The frequencies and intensities of vibration–rotational transitions of water molecules in an argon matrix were calculated for temperatures of 6 and 30 K. The rigid asymmetric top approximation was used with available literature values of the effective rotational constants in the ground and excited vibrational states. The calculations were carried out by taking into account the existence of a non-equilibrium population distribution between the rotational levels of ortho- and para-water isomers. It was assumed that the temperature relaxation of the population of rotational levels is independent of the ortho- and para-isomers. Comparison of the results of the theoretical calculations with experimental literature data shows good agreement for the majority of the rotational structure lines for symmetric and antisymmetric stretching vibrations both in the frequency values and in the values of the relative intensities.     

New analytical expressions,symmetry relations and numerical solutions for the rotational overlap integrals

In this article, extremely simple analytical formulas are obtained for rotational overlap integrals which occur in integrals over two reduced rotation matrix elements. The analytical derivations are based on the properties of the Jacobi polynomials and beta functions. Numerical results and special values for rotational overlap integrals are obtained by using symmetry properties and recurrence relationships for reduced rotation matrix elements. The final results are of surprisingly simple structures and very useful for practical applications. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2012