A complementary thermodynamic limit for classical Coulomb matter

The canonical equilibrium measure of classical two-component Coulomb matter with regularized interactions is analyzed in a finite volume. It is shown that, in the mean-field regime, the one-particle density is inhomogeneous on a new characteristic length scale _inh. For a system ofN positive andN negative particles, _inh and the characteristic length scale of correlations _corr (=Debye screening length) are related via _inh=(2N)^1/2 _corr. The major conceptual conclusion that is drawn from this is that one needs two nontrivial complementary thermodynamic limits to define the equilibrium thermodynamics of two-component Coulomb systems. One of them is the standard thermodynamic limit (infinite volume), where one takesN, _corr fixed. Its complementary limit is characterized byN, _inh fixed, and is a finite-volume inhomogeneous mean-field limit. The most prominent new feature in the mean-field thermodynamic limit, which is absent in the standard thermodynamic limit, is an anomalous first-order phase transition where the Coulomb system explodes or implodes, respectively. The phase transition is connected with the existence of a metastable plasma phase far below the ionization temperature.     

Likelihood estimation of soft-core interaction potentials for Gibbsian point patterns

The likelihood method is developed for the analysis of socalled regular point patterns. Approximating the normalizing factor of Gibbs canonical distribution, we simultaneously estimate two parameters, one for the scale and the other which measures the softness (or hardness), of repulsive interactions between points. The approximations are useful up to a considerably high density. Some real data are analyzed to illustrate the utility of the parameters for characterizing the regular point pattern.     

Relativistic atomic natural orbital type basis sets for the alkaline and alkaline-earth atoms applied to the ground-state potentials for the corresponding dimers

New basis sets of the atomic natural orbital (ANO) type have been developed for the atoms Li–Fr and Be–Ra. The ANOs have been obtained from the average density matrix of the ground states and the lowest excited states of the atom, the positive ion, and the dimer at its equilibirium geometry. Scalar realtivisitc effects are included through the use of a Douglas–Kroll Hamiltonian. Multiconfigurational wave functions have been used with dynamic correlation included using second-order perturbation theory (CASSCF/CASPT2). The basis sets are applied in calculations of the ground-state potentials for the dimers. Computed bond energies are accurate to within 0.05 eV for the alkaline dimers and 0.02 eV for the alkaline-earth dimers (except for Be₂).Acknowledgments.ensp;B.O.R. would like to express his gratitude to Prof. Jacopo Tomasi for all the inspiration that his scientific work has given him through the years and continues to do in particular through the work on solvent effects on molecular properties. This work has been supported by a grant from the Swedish Science Research Council, VR.Contribution to the Jacopo Tomasi Honorary Issue     

Study of a high-dimensional chaotic attractor

The nature of a very high-dimensional chaotic attractor in an infinite-dimensional phase space is examined for the purpose of studying the relationships between the physical processes occurring in the real space and the characteristics of high-dimensional attractor in the phase space. We introduce two complementary bases from which the attractor is observed, one the Lyapunov basis composed of the Lyapunov vectors and the another the Fourier basis composed of the Fourier modes. We introduce the exterior subspaces on the basis of the Lyapunov vectors and observe the chaotic motion projected onto these exteriors. It is shown that a certain statistical property of the projected motion changes markedly as the exterior subspace goes out of the attractor. The origin of such a phenomenon is attributed to more fundamental features of our attractor, which become manifest when the attractor is observed from the Lyapunov basis. A counterpart of the phenomenon can be observed also on the Fourier basis because there is a statistical one-to-one correspondence between the Lyapunov vectors and the Fourier modes. In particular, a statistical property of the high-pass filtered time series reflects clearly the difference between the interior and the exterior of the attractor.     

超球坐标下原子、分子体系的变分计算：—氦原子和氢负离子基态

应用超球会标表示氦原子和氢负离子的薛定谔方程，将二电子原子在三维空间中的运动转化为单电子原子在六维空间中受广义库仑力作用的运动，我们给出了六维空间广义角动量算符的本征值与本征函数，并以此本征函数微基构造超球波函数，得到超球径微分方程，以广义Laguerre 多项式表示超球径波函数，运用密度矩阵和线性变分法得到非正交基下超球径波函数满足的久期方程，最后求得能量和波函数，计算结果与精确的计算符合良好。     

Optimization of auxiliary basis sets for the LEDO expansion and a projection technique for LEDO-DFT

We present a systematic procedure for the optimization of the expansion basis for the limited expansion of diatomic overlap density functional theory (LEDO-DFT) and report on optimized auxiliary orbitals for the Ahlrichs split valence plus polarization basis set (SVP) for the elements H, Li--F, and Na--Cl. A new method to deal with near-linear dependences in the LEDO expansion basis is introduced, which greatly reduces the computational effort of LEDO-DFT calculations. Numerical results for a test set of small molecules demonstrate the accuracy of electronic energies, structural parameters, dipole moments, and harmonic frequencies. For larger molecular systems the numerical errors introduced by the LEDO approximation can lead to an uncontrollable behavior of the self-consistent field (SCF) process. A projection technique suggested by L?wdin is presented in the framework of LEDO-DFT, which guarantees for SCF convergence. Numerical results on some critical test molecules suggest the general applicability of the auxiliary orbitals presented in combination with this projection technique. Timing results indicate that LEDO-DFT is competitive with conventional density fitting methods.     

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.     

Electron-density-dependent fused-sphere surfaces derived from pseudopotential calculations

Fused-sphere surfaces can be used to mimic a molecular boundary associated with a constant value of the electron density. The simplest of such fused-sphere models are constructed by using the atomic radii for the spherical isodensity surfaces of individual atoms. In this work, we discuss the extension of this model to molecules containing atoms beyond the second row. In these many- electron systems, the computation of electron densities is usually simplified by adopting a pseudopotential (or effective-core potential) approach. Here, we discuss the performance of large- and small-core pseudo-potential calculations as a tool to derive atomic radii. Our results provide an optimum set of variable radii that can be used to build fused-sphere surfaces. This continuum of surfaces provides a simple approximation to the low-electron-density regions around molecules with heavy atoms.     

Accurate Gaussian basis sets for the H2O molecule generated with the molecular improved generator coordinate Hartree–Fock method

The molecular improved generator coordinate Hartree–Fock (MIGCHF) method is used to generate accurate basis sets of primitive Gaussian-type functions for the H₂O molecule. Sequences of increasing size atom centered basis sets are employed to explore the accuracy that can be achieved with this method. Using the O(24s14p8d5f2g1h);H(22s9p5d2f1g) basis set, the HF and second-order electron correlation energies of the H₂O ground state at the experimental geometry are computed as −76.0674680 and −0.3491935 hartree, respectively. The HF energy is in error by 20 μhartree and the second-order correlation energy corresponds to 96.5% of an estimate of the limiting value. The relevance of the present calculations is to show the accuracy that can be achieved in studies of small polyatomic molecules with the MIGCHF method.