Finite size effects in the presence of a chemical potential: A study in the classical nonlinear O(2) sigma model

In the presence of a chemical potential, the physics of level crossings leads to singularities at zero temperature, even when the spatial volume is finite. These singularities are smoothed out at a finite temperature but leave behind nontrivial finite size effects which must be understood in order to extract thermodynamic quantities using Monte Carlo methods, particularly close to critical points. We illustrate some of these issues using the classical nonlinear O(2) sigma model with a coupling β and chemical potential μ on a 2+1-dimensional Euclidean lattice. In the conventional formulation this model suffers from a sign problem at nonzero chemical potential and hence cannot be studied with the Wolff cluster algorithm. However, when formulated in terms of the worldline of particles, the sign problem is absent, and the model can be studied efficiently with the "worm algorithm." Using this method we study the finite size effects that arise due to the chemical potential and develop an effective quantum mechanical approach to capture the effects. As a side result we obtain energy levels of up to four particles as a function of the box size and uncover a part of the phase diagram in the (β,μ) plane. © 2010 The American Physical Society.

A mathematical theory of stochastic microlensing. II. Random images, shear, and the Kac-Rice formula

Continuing our development of a mathematical theory of stochastic microlensing, we study the random shear and expected number of random lensed images of different types. In particular, we characterize the first three leading terms in the asymptotic expression of the joint probability density function (pdf) of the random shear tensor due to point masses in the limit of an infinite number of stars. Up to this order, the pdf depends on the magnitude of the shear tensor, the optical depth, and the mean number of stars through a combination of radial position and the star's mass. As a consequence, the pdf's of the shear components are seen to converge, in the limit of an infinite number of stars, to shifted Cauchy distributions, which shows that the shear components have heavy tails in that limit. The asymptotic pdf of the shear magnitude in the limit of an infinite number of stars is also presented. All the results on the random microlensing shear are given for a general point in the lens plane. Extending to the general random distributions (not necessarily uniform) of the lenses, we employ the Kac-Rice formula and Morse theory to deduce general formulas for the expected total number of images and the expected number of saddle images. We further generalize these results by considering random sources defined on a countable compact covering of the light source plane. This is done to introduce the notion of global expected number of positive parity images due to a general lensing map. Applying the result to microlensing, we calculate the asymptotic global expected number of minimum images in the limit of an infinite number of stars, where the stars are uniformly distributed. This global expectation is bounded, while the global expected number of images and the global expected number of saddle images diverge as the order of the number of stars. © 2009 American Institute of Physics.

Reformulating time-dependent density functional theory with non-orthogonal localized molecular orbitals.

Time-dependent density functional theory (TDDFT) has broad application in the study of electronic response, excitation and transport. To extend such application to large and complex systems, we develop a reformulation of TDDFT equations in terms of non-orthogonal localized molecular orbitals (NOLMOs). NOLMO is the most localized representation of electronic degrees of freedom and has been used in ground state calculations. In atomic orbital (AO) representation, the sparsity of NOLMO is transferred to the coefficient matrix of molecular orbitals (MOs). Its novel use in TDDFT here leads to a very simple form of time propagation equations which can be solved with linear-scaling effort. We have tested the method for several long-chain saturated and conjugated molecular systems within the self-consistent charge density-functional tight-binding method (SCC-DFTB) and demonstrated its accuracy. This opens up pathways for TDDFT applications to large bio- and nano-systems.

Water exchange rates of water-soluble manganese(III) porphyrins of therapeutical potential.

The activation parameters and the rate constants of the water-exchange reactions of Mn(III)TE-2-PyP(5+) (meso-tetrakis(N-ethylpyridinium-2-yl)porphyrin) as cationic, Mn(III)TnHex-2-PyP(5+) (meso-tetrakis(N-n-hexylpyridinium-2-yl)porphyrin) as sterically shielded cationic, and Mn(III)TSPP(3-) (meso-tetrakis(4-sulfonatophenyl)porphyrin) as anionic manganese(iii) porphyrins were determined from the temperature dependence of (17)O NMR relaxation rates. The rate constants at 298 K were obtained as 4.12 x 10(6) s(-1), 5.73 x 10(6) s(-1), and 2.74 x 10(7) s(-1), respectively. On the basis of the determined entropies of activation, an interchange-dissociative mechanism (I(d)) was proposed for the cationic complexes (DeltaS(double dagger) = approximately 0 J mol(-1) K(-1)) whereas a limiting dissociative mechanism (D) was proposed for Mn(III)TSPP(3-) complex (DeltaS(double dagger) = +79 J mol(-1) K(-1)). The obtained water exchange rate of Mn(III)TSPP(3-) corresponded well to the previously assumed value used by Koenig et al. (S. H. Koenig, R. D. Brown and M. Spiller, Magn. Reson. Med., 1987, 4, 52-260) to simulate the (1)H NMRD curves, therefore the measured value supports the theory developed for explaining the anomalous relaxivity of Mn(III)TSPP(3-) complex. A magnitude of the obtained water-exchange rate constants further confirms the suggested inner sphere electron transfer mechanism for the reactions of the two positively charged Mn(iii) porphyrins with the various biologically important oxygen and nitrogen reactive species. Due to the high biological and clinical relevance of the reactions that occur at the metal site of the studied Mn(iii) porphyrins, the determination of water exchange rates advanced our insight into their efficacy and mechanism of action, and in turn should impact their further development for both diagnostic (imaging) and therapeutic purposes.

Accurate localized resolution of identity approach for linear-scaling hybrid density functionals and for many-body perturbation theory

© 2015 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.A key component in calculations of exchange and correlation energies is the Coulomb operator, which requires the evaluation of two-electron integrals. For localized basis sets, these four-center integrals are most efficiently evaluated with the resolution of identity (RI) technique, which expands basis-function products in an auxiliary basis. In this work we show the practical applicability of a localized RI-variant ('RI-LVL'), which expands products of basis functions only in the subset of those auxiliary basis functions which are located at the same atoms as the basis functions. We demonstrate the accuracy of RI-LVL for Hartree-Fock calculations, for the PBE0 hybrid density functional, as well as for RPA and MP2 perturbation theory. Molecular test sets used include the S22 set of weakly interacting molecules, the G3 test set, as well as the G2-1 and BH76 test sets, and heavy elements including titanium dioxide, copper and gold clusters. Our RI-LVL implementation paves the way for linear-scaling RI-based hybrid functional calculations for large systems and for all-electron many-body perturbation theory with significantly reduced computational and memory cost.