Basis set convergence on static electric dipole polarizability calculations of alkali-metal clusters
Energy Technology Data Exchange (ETDEWEB)
Souza, Fabio A. L. de; Jorge, Francisco E., E-mail: jorge@cce.ufes.br [Departamento de Fisica, Universidade Federal do Espirito Santo, 29060-900 Vitoria-ES (Brazil)
2013-07-15
A hierarchical sequence of all-electron segmented contracted basis sets of double, triple and quadruple zeta valence qualities plus polarization functions augmented with diffuse functions for the atoms from H to Ar was constructed. A systematic study of basis sets required to obtain reliable and accurate values of static dipole polarizabilities of lithium and sodium clusters (n = 2, 4, 6 and 8) at their optimized equilibrium geometries is reported. Three methods are examined: Hartree-Fock (HF), second-order Moeller-Plesset perturbation theory (MP2), and density functional theory (DFT). By direct calculations or by fitting the directly calculated values through one extrapolation scheme, estimates of the HF, MP2 and DFT complete basis set limits were obtained. Comparison with experimental and theoretical data reported previously in the literature is done (author)
Witte, Jonathon; Neaton, Jeffrey B; Head-Gordon, Martin
2016-05-21
With the aim of systematically characterizing the convergence of common families of basis sets such that general recommendations for basis sets can be made, we have tested a wide variety of basis sets against complete-basis binding energies across the S22 set of intermolecular interactions-noncovalent interactions of small and medium-sized molecules consisting of first- and second-row atoms-with three distinct density functional approximations: SPW92, a form of local-density approximation; B3LYP, a global hybrid generalized gradient approximation; and B97M-V, a meta-generalized gradient approximation with nonlocal correlation. We have found that it is remarkably difficult to reach the basis set limit; for the methods and systems examined, the most complete basis is Jensen's pc-4. The Dunning correlation-consistent sequence of basis sets converges slowly relative to the Jensen sequence. The Karlsruhe basis sets are quite cost effective, particularly when a correction for basis set superposition error is applied: counterpoise-corrected def2-SVPD binding energies are better than corresponding energies computed in comparably sized Dunning and Jensen bases, and on par with uncorrected results in basis sets 3-4 times larger. These trends are exhibited regardless of the level of density functional approximation employed. A sense of the magnitude of the intrinsic incompleteness error of each basis set not only provides a foundation for guiding basis set choice in future studies but also facilitates quantitative comparison of existing studies on similar types of systems. PMID:27208948
Petersson, George A; Malick, David K; Frisch, Michael J; Braunstein, Matthew
2006-07-28
Examination of the convergence of full valence complete active space self-consistent-field configuration interaction including all single and double excitation (CASSCF-CISD) energies with expansion of the one-electron basis set reveals a pattern very similar to the convergence of single determinant energies. Calculations on the lowest four singlet states and the lowest four triplet states of N(2) with the sequence of n-tuple-zeta augmented polarized (nZaP) basis sets (n=2, 3, 4, 5, and 6) are used to establish the complete basis set limits. Full configuration-interaction (CI) and core electron contributions must be included for very accurate potential energy surfaces. However, a simple extrapolation scheme that has no adjustable parameters and requires nothing more demanding than CAS(10e(-),8orb)-CISD/3ZaP calculations gives the R(e), omega(e), omega(e)X(e), T(e), and D(e) for these eight states with rms errors of 0.0006 Angstrom, 4.43 cm(-1), 0.35 cm(-1), 0.063 eV, and 0.018 eV, respectively. PMID:16942134
International Nuclear Information System (INIS)
Systematic sequences of basis sets are used to calculate the spin-orbit splittings of the halogen atoms F, Cl, and Br in the framework of first-order perturbation theory with the Breit-Pauli operator and internally contracted configuration interaction wave functions. The effects of both higher angular momentum functions and the presence of tight functions are studied. By systematically converging the one-particle basis set, an unambiguous evaluation of the effects of correlating different numbers of electrons in the Cl treatment is carried out. Correlation of the 2p-electrons in chlorine increases the spin-orbit splitting by ∼80 cm-1, while in bromine we observe incremental increases of 130, 145, and 93 cm-1, when adding the 3d, 3p, and 2p electrons to the set of explicitly correlated electrons, respectively. For fluorine and chlorine the final basis set limit, all-electrons correlated results match the experimentally observed spin-orbit splittings to within ∼5 cm-1, while for bromine the Breit-Pauli operator underestimates the splitting by about 100 cm-1. More extensive treatment of electron correlation results in only a slight lowering of the spin-orbit matrix elements. Thus, the discrepancy for bromine is proposed to arise from the nonrelativistic character of the underlying wave function. (c) 2000 American Institute of Physics
DEFF Research Database (Denmark)
Avery, John Scales; Rettrup, Sten; Avery, James Emil
In theoretical physics, theoretical chemistry and engineering, one often wishes to solve partial differential equations subject to a set of boundary conditions. This gives rise to eigenvalue problems of which some solutions may be very difficult to find. For example, the problem of finding...... such problems can be much reduced by making use of symmetry-adapted basis functions. The conventional method for generating symmetry-adapted basis sets is through the application of group theory, but this can be difficult. This book describes an easier method for generating symmetry-adapted basis sets...
NIMET PANCAROGLU; FATIH NURAY
2013-01-01
In this paper, we define invariant convergence, lacunary invariant convergence, invariant statistical convergence, lacunary invariant statistical convergence for sequences of sets. We investigate some relations between lacunary invariant statistical convergence and invariant statistical convergence for sequences of sets.
DEFF Research Database (Denmark)
Faber, Rasmus; Buczek, Aneta; Kupka, Teobald;
2016-01-01
The method and basis set dependence of zero-point vibrational corrections (ZPVC) to NMR shielding constants and anisotropies has been investigated using water as a test system. A systematic comparison has been made using the Hartree-Fock (HF), second-order Møller-Plesset perturbation theory (MP2......), coupled cluster singles and doubles (CCSD), coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)) and Kohn-Sham density functional theory (DFT) with the B3LYP exchange-correlation functional methods in combination with the second order vibrational perturbation theory (VPT2......) approach for the vibrational corrections. As basis sets the correlation consistent basis sets cc-pVXZ, aug-ccpVXZ, cc-pCVXZ and aug-cc-pCVXZ with X = D, T, Q, 5, 6 and the polarization consistent basis sets aug-pc-n and aug-pcS-n with n = 1, 2, 3, 4 were employed. Our results show, that basis set...
Energy Technology Data Exchange (ETDEWEB)
Spackman, Peter R.; Karton, Amir, E-mail: amir.karton@uwa.edu.au [School of Chemistry and Biochemistry, The University of Western Australia, Perth, WA 6009 (Australia)
2015-05-15
Coupled cluster calculations with all single and double excitations (CCSD) converge exceedingly slowly with the size of the one-particle basis set. We assess the performance of a number of approaches for obtaining CCSD correlation energies close to the complete basis-set limit in conjunction with relatively small DZ and TZ basis sets. These include global and system-dependent extrapolations based on the A + B/L{sup α} two-point extrapolation formula, and the well-known additivity approach that uses an MP2-based basis-set-correction term. We show that the basis set convergence rate can change dramatically between different systems(e.g.it is slower for molecules with polar bonds and/or second-row elements). The system-dependent basis-set extrapolation scheme, in which unique basis-set extrapolation exponents for each system are obtained from lower-cost MP2 calculations, significantly accelerates the basis-set convergence relative to the global extrapolations. Nevertheless, we find that the simple MP2-based basis-set additivity scheme outperforms the extrapolation approaches. For example, the following root-mean-squared deviations are obtained for the 140 basis-set limit CCSD atomization energies in the W4-11 database: 9.1 (global extrapolation), 3.7 (system-dependent extrapolation), and 2.4 (additivity scheme) kJ mol{sup –1}. The CCSD energy in these approximations is obtained from basis sets of up to TZ quality and the latter two approaches require additional MP2 calculations with basis sets of up to QZ quality. We also assess the performance of the basis-set extrapolations and additivity schemes for a set of 20 basis-set limit CCSD atomization energies of larger molecules including amino acids, DNA/RNA bases, aromatic compounds, and platonic hydrocarbon cages. We obtain the following RMSDs for the above methods: 10.2 (global extrapolation), 5.7 (system-dependent extrapolation), and 2.9 (additivity scheme) kJ mol{sup –1}.
International Nuclear Information System (INIS)
Coupled cluster calculations with all single and double excitations (CCSD) converge exceedingly slowly with the size of the one-particle basis set. We assess the performance of a number of approaches for obtaining CCSD correlation energies close to the complete basis-set limit in conjunction with relatively small DZ and TZ basis sets. These include global and system-dependent extrapolations based on the A + B/Lα two-point extrapolation formula, and the well-known additivity approach that uses an MP2-based basis-set-correction term. We show that the basis set convergence rate can change dramatically between different systems(e.g.it is slower for molecules with polar bonds and/or second-row elements). The system-dependent basis-set extrapolation scheme, in which unique basis-set extrapolation exponents for each system are obtained from lower-cost MP2 calculations, significantly accelerates the basis-set convergence relative to the global extrapolations. Nevertheless, we find that the simple MP2-based basis-set additivity scheme outperforms the extrapolation approaches. For example, the following root-mean-squared deviations are obtained for the 140 basis-set limit CCSD atomization energies in the W4-11 database: 9.1 (global extrapolation), 3.7 (system-dependent extrapolation), and 2.4 (additivity scheme) kJ mol–1. The CCSD energy in these approximations is obtained from basis sets of up to TZ quality and the latter two approaches require additional MP2 calculations with basis sets of up to QZ quality. We also assess the performance of the basis-set extrapolations and additivity schemes for a set of 20 basis-set limit CCSD atomization energies of larger molecules including amino acids, DNA/RNA bases, aromatic compounds, and platonic hydrocarbon cages. We obtain the following RMSDs for the above methods: 10.2 (global extrapolation), 5.7 (system-dependent extrapolation), and 2.9 (additivity scheme) kJ mol–1
Conductance calculations with a wavelet basis set
DEFF Research Database (Denmark)
Thygesen, Kristian Sommer; Bollinger, Mikkel; Jacobsen, Karsten Wedel
2003-01-01
. The linear-response conductance is calculated from the Green's function which is represented in terms of a system-independent basis set containing wavelets with compact support. This allows us to rigorously separate the central region from the contacts and to test for convergence in a systematic way...
Basis convergence of range-separated density-functional theory
International Nuclear Information System (INIS)
Range-separated density-functional theory (DFT) is an alternative approach to Kohn-Sham density-functional theory. The strategy of range-separated density-functional theory consists in separating the Coulomb electron-electron interaction into long-range and short-range components and treating the long-range part by an explicit many-body wave-function method and the short-range part by a density-functional approximation. Among the advantages of using many-body methods for the long-range part of the electron-electron interaction is that they are much less sensitive to the one-electron atomic basis compared to the case of the standard Coulomb interaction. Here, we provide a detailed study of the basis convergence of range-separated density-functional theory. We study the convergence of the partial-wave expansion of the long-range wave function near the electron-electron coalescence. We show that the rate of convergence is exponential with respect to the maximal angular momentum L for the long-range wave function, whereas it is polynomial for the case of the Coulomb interaction. We also study the convergence of the long-range second-order Møller-Plesset correlation energy of four systems (He, Ne, N2, and H2O) with cardinal number X of the Dunning basis sets cc − p(C)V XZ and find that the error in the correlation energy is best fitted by an exponential in X. This leads us to propose a three-point complete-basis-set extrapolation scheme for range-separated density-functional theory based on an exponential formula
Smooth Convergence Away from Singular Sets
Lakzian, Sajjad
2012-01-01
We consider sequences of metrics, $g_j$, on a Riemannian manifold, $M$, which converge smoothly on compact sets away from a singular set $S\\subset M$, to a metric, $g_\\infty$, on $M\\setminus S$. We prove theorems which describe when $M_j=(M, g_j)$ converge in the Gromov-Hausdorff sense to the metric completion, $(M_\\infty,d_\\infty)$, of $(M\\setminus S, g_\\infty)$. To obtain these theorems, we study the intrinsic flat limits of the sequences. A new method, we call hemispherical embedding, is applied to obtain explicit estimates on the Gromov-Hausdorff and Intrinsic Flat distances between Riemannian manifolds with diffeomorphic subdomains.
Property-optimized Gaussian basis sets for molecular response calculations
Rappoport, Dmitrij; Furche, Filipp
2010-10-01
With recent advances in electronic structure methods, first-principles calculations of electronic response properties, such as linear and nonlinear polarizabilities, have become possible for molecules with more than 100 atoms. Basis set incompleteness is typically the main source of error in such calculations since traditional diffuse augmented basis sets are too costly to use or suffer from near linear dependence. To address this problem, we construct the first comprehensive set of property-optimized augmented basis sets for elements H-Rn except lanthanides. The new basis sets build on the Karlsruhe segmented contracted basis sets of split-valence to quadruple-zeta valence quality and add a small number of moderately diffuse basis functions. The exponents are determined variationally by maximization of atomic Hartree-Fock polarizabilities using analytical derivative methods. The performance of the resulting basis sets is assessed using a set of 313 molecular static Hartree-Fock polarizabilities. The mean absolute basis set errors are 3.6%, 1.1%, and 0.3% for property-optimized basis sets of split-valence, triple-zeta, and quadruple-zeta valence quality, respectively. Density functional and second-order Møller-Plesset polarizabilities show similar basis set convergence. We demonstrate the efficiency of our basis sets by computing static polarizabilities of icosahedral fullerenes up to C720 using hybrid density functional theory.
Multiple Spawning with Optimal Basis Set Expansion
Yang, Sandy; Kaduk, Benjamin; Martinez, Todd J
2008-01-01
The Full Multiple Spawning (FMS) method is designed to simulate quantum dynamics in the multi-state electronic problem. The FMS nuclear wavefunction is represented in a basis of coupled, frozen Gaussians, and the spawning procedure prescribes a means of adaptively increasing the size of the basis in order to capture population transfer between electronic states. Parent trajectories create children when passing through regions of significant nonadiabatic coupling. In order to converge branching ratios without allowing the basis to reach an impractical size, population transfer at individual spawning events should be made as effective as possible. Herein we detail a new algorithm for specifying the initial conditions of freshly spawned basis functions, one that minimizes the number of spawns needed for convergence by maximizing the efficiency of individual spawning events. Optimization is achieved by maximizing the coupling between parent and child trajectories, as a function of child position and momentum, at ...
Accurate basis set truncation for wavefunction embedding
Barnes, Taylor A.; Goodpaster, Jason D.; Manby, Frederick R.; Miller, Thomas F.
2013-07-01
Density functional theory (DFT) provides a formally exact framework for performing embedded subsystem electronic structure calculations, including DFT-in-DFT and wavefunction theory-in-DFT descriptions. In the interest of efficiency, it is desirable to truncate the atomic orbital basis set in which the subsystem calculation is performed, thus avoiding high-order scaling with respect to the size of the MO virtual space. In this study, we extend a recently introduced projection-based embedding method [F. R. Manby, M. Stella, J. D. Goodpaster, and T. F. Miller III, J. Chem. Theory Comput. 8, 2564 (2012)], 10.1021/ct300544e to allow for the systematic and accurate truncation of the embedded subsystem basis set. The approach is applied to both covalently and non-covalently bound test cases, including water clusters and polypeptide chains, and it is demonstrated that errors associated with basis set truncation are controllable to well within chemical accuracy. Furthermore, we show that this approach allows for switching between accurate projection-based embedding and DFT embedding with approximate kinetic energy (KE) functionals; in this sense, the approach provides a means of systematically improving upon the use of approximate KE functionals in DFT embedding.
Mechanical basis of morphogenesis and convergent evolution of spiny seashells
Chirat, R.
2013-03-25
Convergent evolution is a phenomenon whereby similar traits evolved independently in not closely related species, and is often interpreted in functional terms. Spines in mollusk seashells are classically interpreted as having repeatedly evolved as a defense in response to shell-crushing predators. Here we consider the morphogenetic process that shapes these structures and underlies their repeated emergence. We develop a mathematical model for spine morphogenesis based on the mechanical interaction between the secreting mantle edge and the calcified shell edge to which the mantle adheres during shell growth. It is demonstrated that a large diversity of spine structures can be accounted for through small variations in control parameters of this natural mechanical process. This physical mechanism suggests that convergent evolution of spines can be understood through a generic morphogenetic process, and provides unique perspectives in understanding the phenotypic evolution of this second largest phylum in the animal kingdom.
Many-Body Basis Set Superposition Effect.
Ouyang, John F; Bettens, Ryan P A
2015-11-10
The basis set superposition effect (BSSE) arises in electronic structure calculations of molecular clusters when questions relating to interactions between monomers within the larger cluster are asked. The binding energy, or total energy, of the cluster may be broken down into many smaller subcluster calculations and the energies of these subsystems linearly combined to, hopefully, produce the desired quantity of interest. Unfortunately, BSSE can plague these smaller fragment calculations. In this work, we carefully examine the major sources of error associated with reproducing the binding energy and total energy of a molecular cluster. In order to do so, we decompose these energies in terms of a many-body expansion (MBE), where a "body" here refers to the monomers that make up the cluster. In our analysis, we found it necessary to introduce something we designate here as a many-ghost many-body expansion (MGMBE). The work presented here produces some surprising results, but perhaps the most significant of all is that BSSE effects up to the order of truncation in a MBE of the total energy cancel exactly. In the case of the binding energy, the only BSSE correction terms remaining arise from the removal of the one-body monomer total energies. Nevertheless, our earlier work indicated that BSSE effects continued to remain in the total energy of the cluster up to very high truncation order in the MBE. We show in this work that the vast majority of these high-order many-body effects arise from BSSE associated with the one-body monomer total energies. Also, we found that, remarkably, the complete basis set limit values for the three-body and four-body interactions differed very little from that at the MP2/aug-cc-pVDZ level for the respective subclusters embedded within a larger cluster. PMID:26574311
A *-mixing convergence theorem for convex set valued processes
Directory of Open Access Journals (Sweden)
R. Kleyle
1987-03-01
Full Text Available In this paper the concept of a *-mixing process is extended to multivalued maps from a probability space into closed, bounded convex sets of a Banach space. The main result, which requires that the Banach space be separable and reflexive, is a convergence theorem for *-mixing sequences which is analogous to the strong law of large numbers. The impetus for studying this problem is provided by a model from information science involving the utilization of feedback data by a decision maker who is uncertain of his goals. The main result is somewhat similar to a theorem for real valued processes and is of interest in its own right.
Basis set approach in the constrained interpolation profile method
International Nuclear Information System (INIS)
We propose a simple polynomial basis-set that is easily extendable to any desired higher-order accuracy. This method is based on the Constrained Interpolation Profile (CIP) method and the profile is chosen so that the subgrid scale solution approaches the real solution by the constraints from the spatial derivative of the original equation. Thus the solution even on the subgrid scale becomes consistent with the master equation. By increasing the order of the polynomial, this solution quickly converges. 3rd and 5th order polynomials are tested on the one-dimensional Schroedinger equation and are proved to give solutions a few orders of magnitude higher in accuracy than conventional methods for lower-lying eigenstates. (author)
Darcie, Thomas E.; Doverspike, Robert; Zirngibl, Martin; Korotky, Steven K.
2005-09-01
Call for Papers: Convergence The Journal of Optical Networking (JON) invites submissions to a special issue on Convergence. Convergence has become a popular theme in telecommunications, one that has broad implications across all segments of the industry. Continual evolution of technology and applications continues to erase lines between traditionally separate lines of business, with dramatic consequences for vendors, service providers, and consumers. Spectacular advances in all layers of optical networking-leading to abundant, dynamic, cost-effective, and reliable wide-area and local-area connections-have been essential drivers of this evolution. As services and networks continue to evolve towards some notion of convergence, the continued role of optical networks must be explored. One vision of convergence renders all information in a common packet (especially IP) format. This vision is driven by the proliferation of data services. For example, time-division multiplexed (TDM) voice becomes VoIP. Analog cable-television signals become MPEG bits streamed to digital set-top boxes. T1 or OC-N private lines migrate to Ethernet virtual private networks (VPNs). All these packets coexist peacefully within a single packet-routing methodology built on an optical transport layer that combines the flexibility and cost of data networks with telecom-grade reliability. While this vision is appealing in its simplicity and shared widely, specifics of implementation raise many challenges and differences of opinion. For example, many seek to expand the role of Ethernet in these transport networks, while massive efforts are underway to make traditional TDM networks more data friendly within an evolved but backward-compatible SDH/SONET (synchronous digital hierarchy and synchronous optical network) multiplexing hierarchy. From this common underlying theme follow many specific instantiations. Examples include the convergence at the physical, logical, and operational levels of voice and
Entrepreneurial Creativity as a Convergent Basis for Teaching Business Communication
Grenci, Richard T.
2012-01-01
Of the "21st Century" business skills of communication, collaboration, creativity, and critical thinking, creativity arguably receives among the least explicit attention in traditional business core curricula. With that in mind, the context of entrepreneurial creativity is put forth as a basis for teaching business communication. By…
Entanglement patterns in mutually unbiased basis sets
International Nuclear Information System (INIS)
A few simply stated rules govern the entanglement patterns that can occur in mutually unbiased basis sets (MUBs) and constrain the combinations of such patterns that can coexist in full complements of MUBs. We consider Hilbert spaces of prime power dimensions (D=pN), as realized by systems of N prime-state particles, where full complements of D+1 MUBs are known to exist, and we assume only that MUBs are eigenbases of generalized Pauli operators, without using any particular construction. The general rules include the following: (1) In any MUB, a given particle appears either in a pure state or totally entangled and (2) in any full MUB complement, each particle is pure in (p+1) bases (not necessarily the same ones) and totally entangled in the remaining (pN-p). It follows that the maximum number of product bases is p+1 and, when this number is realized, all remaining (pN-p) bases in the complement are characterized by the total entanglement of every particle. This ''standard distribution'' is inescapable for two-particle systems (of any p), where only product and generalized Bell bases are admissible MUB types. This and the following results generalize previous results for qubits [Phys. Rev. A 65. 032320 (2002); Phys. Rev. A 72, 062310 (2005)] and qutrits [Phys. Rev. A 70, 012302 (2004)], drawing particularly upon [Phys. Rev. A 72, 062310 (2005)]. With three particles there are three MUB types, and these may be combined in (p+2) different ways to form full complements. With N=4, there are 6 MUB types for p=2, but new MUB types become possible with larger p, and these are essential to realizing full complements. With this example, we argue that new MUB types that show new entanglement patterns should enter with every step in N and, also, when N is a prime plus 1, at a critical p value, p=N-1. Such MUBs should play critical roles in filling complements.
Completeness-optimized basis sets: application to ground-state electron momentum densities.
Lehtola, Jussi; Manninen, Pekka; Hakala, Mikko; Hämäläinen, Keijo
2012-09-14
In the current work we apply the completeness-optimization paradigm [P. Manninen and J. Vaara, J. Comput. Chem. 27, 434 (2006)] to investigate the basis set convergence of the moments of the ground-state electron momentum density at the self-consistent field level of theory. We present a black-box completeness-optimization algorithm that can be used to generate computationally efficient basis sets for computing any property at any level of theory. We show that the complete basis set (CBS) limit of the moments of the electron momentum density can be reached more cost effectively using completeness-optimized basis sets than using conventional, energy-optimized Gaussian basis sets. By using the established CBS limits, we generate a series of smaller basis sets which can be used to systematically approach the CBS and to perform calculations on larger, experimentally interesting systems. PMID:22979848
Anacker, Tony; Hill, J Grant; Friedrich, Joachim
2016-04-21
Minimal basis sets, denoted DSBSenv, based on the segmented basis sets of Ahlrichs and co-workers have been developed for use as environmental basis sets for the domain-specific basis set (DSBS) incremental scheme with the aim of decreasing the CPU requirements of the incremental scheme. The use of these minimal basis sets within explicitly correlated (F12) methods has been enabled by the optimization of matching auxiliary basis sets for use in density fitting of two-electron integrals and resolution of the identity. The accuracy of these auxiliary sets has been validated by calculations on a test set containing small- to medium-sized molecules. The errors due to density fitting are about 2-4 orders of magnitude smaller than the basis set incompleteness error of the DSBSenv orbital basis sets. Additional reductions in computational cost have been tested with the reduced DSBSenv basis sets, in which the highest angular momentum functions of the DSBSenv auxiliary basis sets have been removed. The optimized and reduced basis sets are used in the framework of the domain-specific basis set of the incremental scheme to decrease the computation time without significant loss of accuracy. The computation times and accuracy of the previously used environmental basis and that optimized in this work have been validated with a test set of medium- to large-sized systems. The optimized and reduced DSBSenv basis sets decrease the CPU time by about 15.4% and 19.4% compared with the old environmental basis and retain the accuracy in the absolute energy with standard deviations of 0.99 and 1.06 kJ/mol, respectively. PMID:27002338
Efficient basis for the Dicke model: II. Wave function convergence and excited states
International Nuclear Information System (INIS)
An extended bosonic coherent basis has been shown by Chen et al (2008 Phys. Rev. A 78 051801) and Liu T et al (2009 Phys. Rev. A 80 165308) to provide numerically exact solutions of the finite-size Dicke model. The advantages in employing this basis, as compared with the photon number (Fock) basis, are exhibited to be valid for a large region of the Hamiltonian parameter space and many excited states by analyzing the convergence in the wave functions. (paper)
Efficient basis for the Dicke model: I. Theory and convergence in energy
International Nuclear Information System (INIS)
An extended bosonic coherent basis has been shown by Chen et al (2008 Phys. Rev. A 78 051801) and Liu et al (2009 Phys. Rev. A 80 165308) to provide numerically exact solutions of the finite-size Dicke model. The advantages in employing this basis, as compared with the photon number (Fock) basis, are exhibited to be valid for a large region of the Hamiltonian parameter space by analyzing the converged values of the ground state energy. (paper)
Quantum criticality analysis by finite-size scaling and exponential basis sets
Alharbi, Fahhad H.; Kais, Sabre
2013-04-01
We combine the finite-size scaling method with the mesh-free spectral method to calculate quantum critical parameters for a given Hamiltonian. The basic idea is to expand the exact wave function in a finite exponential basis set and extrapolate the information about system criticality from a finite basis to the infinite basis set limit. The used exponential basis set, though chosen intuitively, allows handling a very wide range of exponential decay rates and calculating multiple eigenvalues simultaneously. As a benchmark system to illustrate the combined approach, we choose the Hulthen potential. The results show that the method is very accurate and converges faster when compared with other basis functions. The approach is general and can be extended to examine near-threshold phenomena for atomic and molecular systems based on even-tempered exponential and Gaussian basis functions.
Convergence to Compact Sets of Inexact Orbits of Nonexpansive Mappings in Banach and Metric Spaces
Directory of Open Access Journals (Sweden)
2009-02-01
Full Text Available We study the influence of computational errors on the convergence to compact sets of orbits of nonexpansive mappings in Banach and metric spaces. We first establish a convergence theorem assuming that the computational errors are summable and then provide examples which show that the summability of errors is necessary for convergence.
The cc-pV5Z-F12 basis set: reaching the basis set limit in explicitly correlated calculations
Peterson, Kirk A; Martin, Jan M L
2014-01-01
We have developed and benchmarked a new extended basis set for explicitly correlated calculations, namely cc-pV5Z-F12. It is offered in two variants, cc-pV5Z-F12 and cc- pV5Z-F12(rev2), the latter of which has additional basis functions on hydrogen not present in the cc-pVnZ-F12 (n=D,T,Q) sequence.A large uncontracted 'reference' basis set is used for benchmarking. cc-pVnZ-F12 (n=D, T, Q, 5) is shown to be a convergent hierarchy. Especially the cc- pV5Z-F12(rev2) basis set can yield the valence CCSD component of total atomization energies (TAEs), without any extrapolation, to an accuracy normally associated with aug-cc-pV{5,6}Z extrapolations. SCF components are functionally at the basis set limit, while the MP2 limit can be approached to as little as 0.01 kcal/mol without extrapolation. The determination of (T) appears to be the most difficult of the three components and cannot presently be accomplished without extrapolation or scaling. (T) extrapolation from cc-pV{T,Q}Z-F12 basis sets, combined with CCSD-F1...
Faster learning algorithm convergence utilizing a combined time-frequency representation as basis
Hendriks, A. J.; Uys, Hermann; du Plessis, Anton; Steenkamp, Christine
2013-10-01
Light is capable of directly manipulating and probing molecular dynamics at its most fundamental level. One versatile approach to influencing such dynamics exploits temporally shaped femtosecond laser pulses. Oftentimes the control mechanisms necessary to induce a desired reaction cannot be determined theoretically a priori. However under certain circumstances these mechanisms can be extracted experimentally through trial and error. This can be implemented systematically by using an evolutionary learning algorithm (LA) with closed loop feedback. Most frequently, pulse shaping algorithms operate within either the time or frequency domain, however seldom both. This may influence the physical insight gained due to dependence on the search basis, as well as influence the speed the algorithm takes to converge. As an alternative to the Fourier domain basis, we make use of a combined time-frequency representation known as the von Neumann basis where we observe temporal and spectral effects at the same time. We report on the numerical and experimental results obtained using the Fourier, as well as the von Neumann basis to maximize the second harmonic generation (SHG) output in a non-linear crystal. We show that the von Neumann representation converges faster than the Fourier domain when compared to searches in the Fourier domain. We also show a reduced parameter space is required for the Fourier domain to converge efficiently, but not for von Neumann domain. Finally we show the highest SHG signal is not only a consequence of the shortest pulse, but that the pulse central frequency also plays a key role. Taken together these results suggest that the von Neumann basis can be used as a viable alternative to the Fourier domain with improved convergence time and potentially deeper physical insight.
Non-Standard Rates of Convergence of Criterion-Function-Based Set Estimators
Jason R. Blevins
2013-01-01
This paper establishes conditions for consistency and potentially non-standard rates of convergence for set estimators based on contour sets of criterion functions. These conditions cover the standard parametric rate $n^{-1/2}$, non-standard polynomial rates such as $n^{-1/3}$, and an extreme case of arbitrarily fast convergence. We also establish the validity of a subsampling procedure for constructing confidence sets for the identified set. We then provide more convenient sufficient conditi...
Convergence of hausdorff approximation methods for the Edgeworth-Pareto hull of a compact set
Efremov, R. V.
2015-11-01
The Hausdorff methods comprise an important class of polyhedral approximation methods for convex compact bodies, since they have an optimal convergence rate and possess other useful properties. The concept of Hausdorff methods is extended to a problem arising in multicriteria optimization, namely, to the polyhedral approximation of the Edgeworth-Pareto hull (EPH) of a convex compact set. It is shown that the sequences of polyhedral sets generated by Hausdorff methods converge to the EPH to be approximated. It is shown that the Estimate Refinement method, which is most frequently used to approximate the EPH of convex compact sets, is a Hausdorff method and, hence, generates sequences of sets converging to the EPH.
Gidofalvi, Gergely
2014-01-01
Molecule-optimized basis sets, based on approximate natural orbitals, are developed for accelerating the convergence of quantum calculations with strongly correlated (multi-referenced) electrons. We use a low-cost approximate solution of the anti-Hermitian contracted Schr{\\"o}dinger equation (ACSE) for the one- and two-electron reduced density matrices (RDMs) to generate an approximate set of natural orbitals for strongly correlated quantum systems. The natural-orbital basis set is truncated to generate a molecule-optimized basis set whose rank matches that of a standard correlation-consistent basis set optimized for the atoms. We show that basis-set truncation by approximate natural orbitals can be viewed as a one-electron unitary transformation of the Hamiltonian operator and suggest an extension of approximate natural-orbital truncations through two-electron unitary transformations of the Hamiltonian operator, such as those employed in the solution of the ACSE. The molecule-optimized basis set from the ACS...
Point Set Denoising Using Bootstrap-Based Radial Basis Function
Ramli, Ahmad; Abd. Majid, Ahmad
2016-01-01
This paper examines the application of a bootstrap test error estimation of radial basis functions, specifically thin-plate spline fitting, in surface smoothing. The presence of noisy data is a common issue of the point set model that is generated from 3D scanning devices, and hence, point set denoising is one of the main concerns in point set modelling. Bootstrap test error estimation, which is applied when searching for the smoothing parameters of radial basis functions, is revisited. The main contribution of this paper is a smoothing algorithm that relies on a bootstrap-based radial basis function. The proposed method incorporates a k-nearest neighbour search and then projects the point set to the approximated thin-plate spline surface. Therefore, the denoising process is achieved, and the features are well preserved. A comparison of the proposed method with other smoothing methods is also carried out in this study. PMID:27315105
Real ${\\mathcal L}^2$ basis-set-size scaling of energies and widths of resonance states
Pont, Federico M; Osenda, Omar
2011-01-01
The resonance states of one-particle Hamiltonians are studied using variational expansions with real basis-set functions. The resonance energies, $E_r$, and widths, $\\Gamma$, are calculated using the density of states and an ${\\mathcal L}^2$ golden rule-like formula, respectively. We present a recipe to select adequately some solutions of the variational problem. The set of approximate energies obtained show a very regular behaviour with the basis-set size, $N$. Indeed, these particular variational eigenvalues show a quite simple scaling behaviour and convergence when $N\\rightarrow \\infty$. Following the same prescription to choose particular solutions of the variational problem we obtain a set of approximate widths. Using the scaling function that characterizes the behaviour of the approximate energies as a guide, it is possible to find a very good approximation to the actual value of the resonance width.
Localized atomic basis set in the projector augmented wave method
DEFF Research Database (Denmark)
Larsen, Ask Hjorth; Vanin, Marco; Mortensen, Jens Jørgen;
2009-01-01
is tested by calculating atomization energies and equilibrium bulk properties of a variety of molecules and solids, comparing to the grid results. Finally, it is demonstrated how a grid-quality structure optimization can be performed with significantly reduced computational effort by switching......We present an implementation of localized atomic-orbital basis sets in the projector augmented wave (PAW) formalism within the density-functional theory. The implementation in the real-space GPAW code provides a complementary basis set to the accurate but computationally more demanding grid...
DEFF Research Database (Denmark)
Prasad, Ramjee
2009-01-01
This paper presents the main conclusions which can be drawn from the discussions on Future Communication Systems and lessons on Unpredictable Future of Wireless Communication Systems. Future systems beyond the third generation are already under discussions in international bodies, such as ITU, WW...... and R&D programmes worldwide. The incoming era is characterized by the convergence of networks and access technology and the divergence of applications. Future mobile communication systems should bring something more than only faster data or wireless internet access....
International Nuclear Information System (INIS)
The authors consider the problem of computing tunneling matrix elements for bridge-mediated electron transfer reactions using the Loewdin [J. Math. Phys. 3, 969 (1962); J. Mol. Spectrosc. 13, 326 (1964)] projection-iteration technique with a nonorthogonal basis set. They compare the convergence properties of two different Loewdin projections, one containing the overlap matrix S and the other containing the inverse S-1 in the projected Hamiltonian. It was suggested in the literature that the projected Hamiltonian with S-1 has better convergence properties compared to the projected Hamiltonian with S. The authors test this proposal using a simple analytical model, and ab initio Hartree-Fock calculations on different molecules with several types of basis sets. Their calculations show that, for Gaussian-type basis sets, the projected Hamiltonian containing S has the best convergence properties, especially for diffuse basis sets and in the strong coupling limit. The limit of diffuse basis sets is relevant to tunneling matrix element calculations involving excited states and anionic electron transfer
Convergence of many-body wavefunction expansions using a plane wave basis in the thermodynamic limit
Shepherd, James J
2016-01-01
Basis set incompleteness error and finite size error can manifest concurrently in systems for which the two effects are phenomenologically well-separated in length scale. When this is true, we need not necessarily remove the two sources of error simultaneously. Instead, the errors can be found and remedied in different parts of the basis set. This would be of great benefit to a method such as coupled cluster theory since the combined cost of $n_{occ}^6 n_{virt}^4$ could be separated into $n_{occ}^6$ and $n_{virt}^4$ costs with smaller prefactors. In this Communication, we present analysis on a data set due to Baardsen and coworkers, containing coupled cluster doubles energies for the 2DEG for $r_s=$ 0.5, 1.0 and 2.0 a.u.~at a wide range of basis set sizes and particle numbers. In obtaining complete basis set limit thermodynamic limit results, we find that within a small and removable error the above assertion is correct for this simple system. This approach allows for the combination of methods which separate...
Application of discrete basis set methods to the Dirac equation
International Nuclear Information System (INIS)
The basis sets described provide a convenient method for the direct relativistic calculation of a wide variety of atomic properties. Results have been obtained for the ground state polarizability, and the 2s/sub 1/2/-1s/sub 1/2/ two-photon transition rate of hydrogenic ions. The latter calculation includes all retardation effects and contributions from higher photon multipoles
Relativistic LCAO with Minimax Principle and New Balanced Basis Sets
Zhang, Hui
2008-01-01
Relativistic density functional theory is widely applied in molecular calculations with heavy atoms, where relativistic and correlation effects are on the same footing. Variational stability of the Dirac Hamiltonian is a very important field of research from the beginning of relativistic molecular calculations on, among efforts for accuracy, efficiency, and density functional formulation, etc. Approximations of one- or two-component methods and searching for suitable basis sets are two major ...
The convergence of double Fourier-Haar series over homothetic copies of sets
International Nuclear Information System (INIS)
The paper is concerned with the convergence of double Fourier- Haar series with partial sums taken over homothetic copies of a given bounded set W⊂R+2 containing the intersection of some neighbourhood of the origin with R+2. It is proved that for a set W from a fairly broad class (in particular, for convex W) there are two alternatives: either the Fourier-Haar series of an arbitrary function f∈L([0,1]2) converges almost everywhere or Lln+ L([0,1]2) is the best integral class in which the double Fourier-Haar series converges almost everywhere. Furthermore, a characteristic property is obtained, which distinguishes which of the two alternatives is realized for a given W. Bibliography: 12 titles. (paper)
Coupled-cluster based basis sets for valence correlation calculations
Claudino, Daniel; Gargano, Ricardo; Bartlett, Rodney J.
2016-03-01
Novel basis sets are generated that target the description of valence correlation in atoms H through Ar. The new contraction coefficients are obtained according to the Atomic Natural Orbital (ANO) procedure from CCSD(T) (coupled-cluster singles and doubles with perturbative triples correction) density matrices starting from the primitive functions of Dunning et al. [J. Chem. Phys. 90, 1007 (1989); ibid. 98, 1358 (1993); ibid. 100, 2975 (1993)] (correlation consistent polarized valence X-tuple zeta, cc-pVXZ). The exponents of the primitive Gaussian functions are subject to uniform scaling in order to ensure satisfaction of the virial theorem for the corresponding atoms. These new sets, named ANO-VT-XZ (Atomic Natural Orbital Virial Theorem X-tuple Zeta), have the same number of contracted functions as their cc-pVXZ counterparts in each subshell. The performance of these basis sets is assessed by the evaluation of the contraction errors in four distinct computations: correlation energies in atoms, probing the density in different regions of space via (-3 ≤ n ≤ 3) in atoms, correlation energies in diatomic molecules, and the quality of fitting potential energy curves as measured by spectroscopic constants. All energy calculations with ANO-VT-QZ have contraction errors within "chemical accuracy" of 1 kcal/mol, which is not true for cc-pVQZ, suggesting some improvement compared to the correlation consistent series of Dunning and co-workers.
Correlation consistent basis sets for actinides. I. The Th and U atoms
Energy Technology Data Exchange (ETDEWEB)
Peterson, Kirk A., E-mail: kipeters@wsu.edu [Department of Chemistry, Washington State University, Pullman, Washington 99164-4630 (United States)
2015-02-21
New correlation consistent basis sets based on both pseudopotential (PP) and all-electron Douglas-Kroll-Hess (DKH) Hamiltonians have been developed from double- to quadruple-zeta quality for the actinide atoms thorium and uranium. Sets for valence electron correlation (5f6s6p6d), cc − pV nZ − PP and cc − pV nZ − DK3, as well as outer-core correlation (valence + 5s5p5d), cc − pwCV nZ − PP and cc − pwCV nZ − DK3, are reported (n = D, T, Q). The -PP sets are constructed in conjunction with small-core, 60-electron PPs, while the -DK3 sets utilized the 3rd-order Douglas-Kroll-Hess scalar relativistic Hamiltonian. Both series of basis sets show systematic convergence towards the complete basis set limit, both at the Hartree-Fock and correlated levels of theory, making them amenable to standard basis set extrapolation techniques. To assess the utility of the new basis sets, extensive coupled cluster composite thermochemistry calculations of ThF{sub n} (n = 2 − 4), ThO{sub 2}, and UF{sub n} (n = 4 − 6) have been carried out. After accurately accounting for valence and outer-core correlation, spin-orbit coupling, and even Lamb shift effects, the final 298 K atomization enthalpies of ThF{sub 4}, ThF{sub 3}, ThF{sub 2}, and ThO{sub 2} are all within their experimental uncertainties. Bond dissociation energies of ThF{sub 4} and ThF{sub 3}, as well as UF{sub 6} and UF{sub 5}, were similarly accurate. The derived enthalpies of formation for these species also showed a very satisfactory agreement with experiment, demonstrating that the new basis sets allow for the use of accurate composite schemes just as in molecular systems composed only of lighter atoms. The differences between the PP and DK3 approaches were found to increase with the change in formal oxidation state on the actinide atom, approaching 5-6 kcal/mol for the atomization enthalpies of ThF{sub 4} and ThO{sub 2}. The DKH3 atomization energy of ThO{sub 2} was calculated to be smaller than the DKH2
Kjaer, Hanna; Nielsen, Monia R; Pagola, Gabriel I; Ferraro, Marta B; Lazzeretti, Paolo; Sauer, Stephan P A
2012-09-01
In this article, we present the so far most extended investigation of the calculation of the coupling constant polarizability of a molecule. The components of the coupling constant polarizability are derivatives of the nuclear magnetic resonance (NMR) indirect nuclear spin-spin coupling constant with respect to an external electric field and play an important role for both chiral discrimination and solvation effects on NMR coupling constants. In this study, we illustrate the effects of one-electron basis sets and electron correlation both at the level of density functional theory as well as second-order polarization propagator approximation for the small molecule hydrogen peroxide, which allowed us to perform calculations with the largest available basis sets optimized for the calculation of NMR coupling constants. We find a systematic but rather slow convergence with the one-electron basis set and that augmentation functions are required. We observe also large and nonsystematic correlation effects with significant differences between the density functional and wave function theory methods. PMID:22618604
Petruzielo, Frank R.; Toulouse, Julien; Umrigar, C. J.
2010-01-01
A simple yet general method for constructing basis sets for molecular electronic structure calculations is presented. These basis sets consist of atomic natural orbitals from a multi-configurational self-consistent field calculation supplemented with primitive functions, chosen such that the asymptotics are appropriate for the potential of the system. Primitives are optimized for the homonuclear diatomic molecule to produce a balanced basis set. Two general features that facilitate this basis...
An efficient basis set representation for calculating electrons in molecules
Jones, Jeremiah R; Lawler, Keith V; Vecharynski, Eugene; Ibrahim, Khaled Z; Williams, Samuel; Abeln, Brant; Yang, Chao; Haxton, Daniel J; McCurdy, C William; Li, Xiaoye S; Rescigno, Thomas N
2015-01-01
The method of McCurdy, Baertschy, and Rescigno, J. Phys. B, 37, R137 (2004) is generalized to obtain a straightforward, surprisingly accurate, and scalable numerical representation for calculating the electronic wave functions of molecules. It uses a basis set of product sinc functions arrayed on a Cartesian grid, and yields 1 kcal/mol precision for valence transition energies with a grid resolution of approximately 0.1 bohr. The Coulomb matrix elements are replaced with matrix elements obtained from the kinetic energy operator. A resolution-of-the-identity approximation renders the primitive one- and two-electron matrix elements diagonal; in other words, the Coulomb operator is local with respect to the grid indices. The calculation of contracted two-electron matrix elements among orbitals requires only O(N log(N)) multiplication operations, not O(N^4), where N is the number of basis functions; N = n^3 on cubic grids. The representation not only is numerically expedient, but also produces energies and proper...
Convergent Dual Bounds Using an Aggregation of Set-Covering Constraints for Capacitated Problems
Porumbel, Daniel Cosmin; Clautiaux, François
2013-01-01
Extended formulations are now widely used to solve hard combinatorial optimization problems. Such formulations have prohibitively-many variables and are generally solved via Column Generation (CG). CG algorithms are known to have frequent convergence issues, and, up to a sometimes large number of iterations, classical Lagrangian dual bounds may be weak. This paper is devoted to set-covering problems in which all elements to cover require a given \\emph{resource consumption} and all feasible co...
Green's function multiple-scattering theory with a truncated basis set: An augmented-KKR formalism
Alam, Aftab; Khan, Suffian N.; Smirnov, A. V.; Nicholson, D. M.; Johnson, Duane D.
2014-11-01
The Korringa-Kohn-Rostoker (KKR) Green's function, multiple-scattering theory is an efficient site-centered, electronic-structure technique for addressing an assembly of N scatterers. Wave functions are expanded in a spherical-wave basis on each scattering center and indexed up to a maximum orbital and azimuthal number Lmax=(l,mmax), while scattering matrices, which determine spectral properties, are truncated at Lt r=(l,mt r) where phase shifts δl >ltr are negligible. Historically, Lmax is set equal to Lt r, which is correct for large enough Lmax but not computationally expedient; a better procedure retains higher-order (free-electron and single-site) contributions for Lmax>Lt r with δl >ltr set to zero [X.-G. Zhang and W. H. Butler, Phys. Rev. B 46, 7433 (1992), 10.1103/PhysRevB.46.7433]. We present a numerically efficient and accurate augmented-KKR Green's function formalism that solves the KKR equations by exact matrix inversion [R3 process with rank N (ltr+1 ) 2 ] and includes higher-L contributions via linear algebra [R2 process with rank N (lmax+1) 2 ]. The augmented-KKR approach yields properly normalized wave functions, numerically cheaper basis-set convergence, and a total charge density and electron count that agrees with Lloyd's formula. We apply our formalism to fcc Cu, bcc Fe, and L 1 0 CoPt and present the numerical results for accuracy and for the convergence of the total energies, Fermi energies, and magnetic moments versus Lmax for a given Lt r.
Friese, Daniel H; Ringholm, Magnus; Gao, Bin; Ruud, Kenneth
2015-10-13
We present theory, implementation, and applications of a recursive scheme for the calculation of single residues of response functions that can treat perturbations that affect the basis set. This scheme enables the calculation of nonlinear light absorption properties to arbitrary order for other perturbations than an electric field. We apply this scheme for the first treatment of two-photon circular dichroism (TPCD) using London orbitals at the Hartree-Fock level of theory. In general, TPCD calculations suffer from the problem of origin dependence, which has so far been solved by using the velocity gauge for the electric dipole operator. This work now enables comparison of results from London orbital and velocity gauge based TPCD calculations. We find that the results from the two approaches both exhibit strong basis set dependence but that they are very similar with respect to their basis set convergence. PMID:26574270
Radiobiological basis for setting neutron radiation safety standards
International Nuclear Information System (INIS)
Present neutron standards, adopted more than 20 yr ago from a weak radiobiological data base, have been in doubt for a number of years and are currently under challenge. Moreover, recent dosimetric re-evaluations indicate that Hiroshima neutron doses may have been much lower than previously thought, suggesting that direct data for neutron-induced cancer in humans may in fact not be available. These recent developments make it urgent to determine the extent to which neutron cancer risk in man can be estimated from data that are available. Two approaches are proposed here that are anchored in particularly robust epidemiological and experimental data and appear most likely to provide reliable estimates of neutron cancer risk in man. The first approach uses gamma-ray dose-response relationships for human carcinogenesis, available from Nagasaki (Hiroshima data are also considered), together with highly characterized neutron and gamma-ray data for human cytogenetics. When tested against relevant experimental data, this approach either adequately predicts or somewhat overestimates neutron tumorigenesis (and mutagenesis) in animals. The second approach also uses the Nagasaki gamma-ray cancer data, but together with neutron RBEs from animal tumorigenesis studies. Both approaches give similar results and provide a basis for setting neutron radiation safety standards. They appear to be an improvement over previous approaches, including those that rely on highly uncertain maximum neutron RBEs and unnecessary extrapolations of gamma-ray data to very low doses. Results suggest that, at the presently accepted neutron dose limit of 0.5 rad/yr, the cancer mortality risk to radiation workers is not very different from accidental mortality risks to workers in various nonradiation occupations
International Nuclear Information System (INIS)
The multireference configuration-interaction (MR-CI) method is used to calculate the binding energy of the He dimer. The convergence of the binding energy to the configuration-set limit (full-CI) is followed by progressively extending the multireference configuration set. Two variants of the Pople size-extensivity correction are applied. The distance dependence of the corrections and hence the effect upon the binding energy turns out to be very small. The effect of orbital optimization is studied and it is shown that it is sufficient to optimize the orbitals used for the multireference space in an atomic multiconfiguration self-consistent field (MCSCF) calculation. In a basis of 50 atomic orbitals, the full-CI binding energy of -9.08 K can be reproduced to 0.00 K (0.02 K) in calculations using only 37 (27) reference configurations, built from the atomic 1s, 2s, 2p, and 3s natural orbitals. Using a very large basis, the 37-reference set gives a best binding energy of -10.87 K, in satisfactory agreement with Aziz's recent semiempirical result of -10.95 K. These findings suggest that the MR-CI method can be developed into an efficient tool for calculating accurate van der Waals interaction energies for larger systems
Zhang, Gaigong; Hu, Wei; Yang, Chao; Pask, John E
2015-01-01
Recently, we have proposed the adaptive local basis set for electronic structure calculations based on Kohn-Sham density functional theory in a pseudopotential framework. The adaptive local basis set is efficient and systematically improvable for total energy calculations. In this paper, we present the calculation of atomic forces, which can be used for a range of applications such as geometry optimization and molecular dynamics simulation. We demonstrate that, under mild assumptions, the computation of atomic forces can scale nearly linearly with the number of atoms in the system using the adaptive local basis set. We quantify the accuracy of the Hellmann-Feynman forces for a range of physical systems, benchmarked against converged planewave calculations, and find that the adaptive local basis set is efficient for both force and energy calculations, requiring at most a few tens of basis functions per atom to attain accuracy required in practice. Since the adaptive local basis set has implicit dependence on a...
Nikolaev, A. V.; Lamoen, D.; Partoens, B.
2016-07-01
In order to increase the accuracy of the linearized augmented plane wave (LAPW) method, we present a new approach where the plane wave basis function is augmented by two different atomic radial components constructed at two different linearization energies corresponding to two different electron bands (or energy windows). We demonstrate that this case can be reduced to the standard treatment within the LAPW paradigm where the usual basis set is enriched by the basis functions of the tight binding type, which go to zero with zero derivative at the sphere boundary. We show that the task is closely related with the problem of extended core states which is currently solved by applying the LAPW method with local orbitals (LAPW+LO). In comparison with LAPW+LO, the number of supplemented basis functions in our approach is doubled, which opens up a new channel for the extension of the LAPW and LAPW+LO basis sets. The appearance of new supplemented basis functions absent in the LAPW+LO treatment is closely related with the existence of the u ˙ l -component in the canonical LAPW method. We discuss properties of additional tight binding basis functions and apply the extended basis set for computation of electron energy bands of lanthanum (face and body centered structures) and hexagonal close packed lattice of cadmium. We demonstrate that the new treatment gives lower total energies in comparison with both canonical LAPW and LAPW+LO, with the energy difference more pronounced for intermediate and poor LAPW basis sets.
Network Decomposition and Maximum Independent Set Part Ⅰ: Theoretic Basis
Institute of Scientific and Technical Information of China (English)
朱松年; 朱嫱
2003-01-01
The structure and characteristics of a connected network are analyzed, and a special kind of sub-network, which can optimize the iteration processes, is discovered. Then, the sufficient and necessary conditions for obtaining the maximum independent set are deduced. It is found that the neighborhood of this sub-network possesses the similar characters, but both can never be allowed incorporated together. Particularly, it is identified that the network can be divided into two parts by a certain style, and then both of them can be transformed into a pair sets network, where the special sub-networks and their neighborhoods appear alternately distributed throughout the entire pair sets network. By use of this characteristic, the network decomposed enough without losing any solutions is obtained. All of these above will be able to make well ready for developing a much better algorithm with polynomial time bound for an odd network in the the application research part of this subject.
Tweedie convergence: a mathematical basis for Taylor's power law, 1/f noise, and multifractality.
Kendal, Wayne S; Jørgensen, Bent
2011-12-01
Plants and animals of a given species tend to cluster within their habitats in accordance with a power function between their mean density and the variance. This relationship, Taylor's power law, has been variously explained by ecologists in terms of animal behavior, interspecies interactions, demographic effects, etc., all without consensus. Taylor's law also manifests within a wide range of other biological and physical processes, sometimes being referred to as fluctuation scaling and attributed to effects of the second law of thermodynamics. 1/f noise refers to power spectra that have an approximately inverse dependence on frequency. Like Taylor's law these spectra manifest from a wide range of biological and physical processes, without general agreement as to cause. One contemporary paradigm for 1/f noise has been based on the physics of self-organized criticality. We show here that Taylor's law (when derived from sequential data using the method of expanding bins) implies 1/f noise, and that both phenomena can be explained by a central limit-like effect that establishes the class of Tweedie exponential dispersion models as foci for this convergence. These Tweedie models are probabilistic models characterized by closure under additive and reproductive convolution as well as under scale transformation, and consequently manifest a variance to mean power function. We provide examples of Taylor's law, 1/f noise, and multifractality within the eigenvalue deviations of the Gaussian unitary and orthogonal ensembles, and show that these deviations conform to the Tweedie compound Poisson distribution. The Tweedie convergence theorem provides a unified mathematical explanation for the origin of Taylor's law and 1/f noise applicable to a wide range of biological, physical, and mathematical processes, as well as to multifractality. PMID:22304168
International Nuclear Information System (INIS)
New correlation consistent basis sets, cc-pVnZ-PP-F12 (n = D, T, Q), for all the post-d main group elements Ga–Rn have been optimized for use in explicitly correlated F12 calculations. The new sets, which include not only orbital basis sets but also the matching auxiliary sets required for density fitting both conventional and F12 integrals, are designed for correlation of valence sp, as well as the outer-core d electrons. The basis sets are constructed for use with the previously published small-core relativistic pseudopotentials of the Stuttgart-Cologne variety. Benchmark explicitly correlated coupled-cluster singles and doubles with perturbative triples [CCSD(T)-F12b] calculations of the spectroscopic properties of numerous diatomic molecules involving 4p, 5p, and 6p elements have been carried out and compared to the analogous conventional CCSD(T) results. In general the F12 results obtained with a n-zeta F12 basis set were comparable to conventional aug-cc-pVxZ-PP or aug-cc-pwCVxZ-PP basis set calculations obtained with x = n + 1 or even x = n + 2. The new sets used in CCSD(T)-F12b calculations are particularly efficient at accurately recovering the large correlation effects of the outer-core d electrons
Shepherd, James J
2016-07-21
Basis set incompleteness error and finite size error can manifest concurrently in systems for which the two effects are phenomenologically well-separated in length scale. When this is true, we need not necessarily remove the two sources of error simultaneously. Instead, the errors can be found and remedied in different parts of the basis set. This would be of great benefit to a method such as coupled cluster theory since the combined cost of nocc (6)nvirt (4) could be separated into nocc (6) and nvirt (4) costs with smaller prefactors. In this Communication, we present analysis on a data set due to Baardsen and co-workers, containing 2D uniform electron gas coupled cluster doubles energies for rs = 0.5, 1.0, and 2.0 a.u. at a wide range of basis set sizes and particle numbers. In obtaining complete basis set limit thermodynamic limit results, we find that within a small and removable error the above assertion is correct for this simple system. We then use this method to obtain similar results for the 3D electron gas at rs = 1.0, 2.0, and 5.0 a.u. and make comparison to the Ceperley-Alder quantum Monte Carlo results. This approach allows for the combination of methods which separately address finite size effects and basis set incompleteness error. PMID:27448865
Shepherd, James J.
2016-07-01
Basis set incompleteness error and finite size error can manifest concurrently in systems for which the two effects are phenomenologically well-separated in length scale. When this is true, we need not necessarily remove the two sources of error simultaneously. Instead, the errors can be found and remedied in different parts of the basis set. This would be of great benefit to a method such as coupled cluster theory since the combined cost of nocc 6 nvirt 4 could be separated into nocc 6 and nvirt 4 costs with smaller prefactors. In this Communication, we present analysis on a data set due to Baardsen and co-workers, containing 2D uniform electron gas coupled cluster doubles energies for rs = 0.5, 1.0, and 2.0 a.u. at a wide range of basis set sizes and particle numbers. In obtaining complete basis set limit thermodynamic limit results, we find that within a small and removable error the above assertion is correct for this simple system. We then use this method to obtain similar results for the 3D electron gas at rs = 1.0, 2.0, and 5.0 a.u. and make comparison to the Ceperley-Alder quantum Monte Carlo results. This approach allows for the combination of methods which separately address finite size effects and basis set incompleteness error.
Petruzielo, F. R.; Toulouse, Julien; Umrigar, C. J.
2011-02-01
A simple yet general method for constructing basis sets for molecular electronic structure calculations is presented. These basis sets consist of atomic natural orbitals from a multiconfigurational self-consistent field calculation supplemented with primitive functions, chosen such that the asymptotics are appropriate for the potential of the system. Primitives are optimized for the homonuclear diatomic molecule to produce a balanced basis set. Two general features that facilitate this basis construction are demonstrated. First, weak coupling exists between the optimal exponents of primitives with different angular momenta. Second, the optimal primitive exponents for a chosen system depend weakly on the particular level of theory employed for optimization. The explicit case considered here is a basis set appropriate for the Burkatzki-Filippi-Dolg pseudopotentials. Since these pseudopotentials are finite at nuclei and have a Coulomb tail, the recently proposed Gauss-Slater functions are the appropriate primitives. Double- and triple-zeta bases are developed for elements hydrogen through argon. These new bases offer significant gains over the corresponding Burkatzki-Filippi-Dolg bases at various levels of theory. Using a Gaussian expansion of the basis functions, these bases can be employed in any electronic structure method. Quantum Monte Carlo provides an added benefit: expansions are unnecessary since the integrals are evaluated numerically.
Diversity of burial rates in convergent settings decreased as Earth aged
Nicoli, Gautier; Moyen, Jean-François; Stevens, Gary
2016-05-01
The evolution and the growth of the continental crust is inextricably linked to the evolution of Earth’s geodynamic processes. The detrital zircon record within the continental crust, as well as the isotopic composition of this crust, indicates that the amount of juvenile felsic material decreased with time and that in geologically recent times, the generation of new crust is balanced by recycling of the crust back into the mantle within subduction zones. However it cannot always have been so; yet the nature of the crust and the processes of crustal reworking in the Precambrian Earth are not well constrained. Here we use both detrital zircon ages and metamorphic pressure-temperature-time (P-T-t) information from metasedimentary units deposited in proposed convergent settings from Archaean, Proterozoic and Phanerozoic terrains to characterize the evolution of minimum estimates of burial rate (km.Ma‑1) as a function of the age of the rocks. The demonstrated decrease in burial rate correlates positively with a progressive decrease in the production of juvenile felsic crust in the Archaean and Proterozoic. Burial rates are also more diverse in the Archaean than in modern times. We interpret these features to reflect a progressive decrease in the diversity of tectonic processes from Archaean to present, coupled with the emergence of the uniquely Phanerozoic modern-style collision.
Diversity of burial rates in convergent settings decreased as Earth aged.
Nicoli, Gautier; Moyen, Jean-François; Stevens, Gary
2016-01-01
The evolution and the growth of the continental crust is inextricably linked to the evolution of Earth's geodynamic processes. The detrital zircon record within the continental crust, as well as the isotopic composition of this crust, indicates that the amount of juvenile felsic material decreased with time and that in geologically recent times, the generation of new crust is balanced by recycling of the crust back into the mantle within subduction zones. However it cannot always have been so; yet the nature of the crust and the processes of crustal reworking in the Precambrian Earth are not well constrained. Here we use both detrital zircon ages and metamorphic pressure-temperature-time (P-T-t) information from metasedimentary units deposited in proposed convergent settings from Archaean, Proterozoic and Phanerozoic terrains to characterize the evolution of minimum estimates of burial rate (km.Ma(-1)) as a function of the age of the rocks. The demonstrated decrease in burial rate correlates positively with a progressive decrease in the production of juvenile felsic crust in the Archaean and Proterozoic. Burial rates are also more diverse in the Archaean than in modern times. We interpret these features to reflect a progressive decrease in the diversity of tectonic processes from Archaean to present, coupled with the emergence of the uniquely Phanerozoic modern-style collision. PMID:27216133
Bowen, J. Philip; Sorensen, Jennifer B.; Kirschner, Karl N.
2007-01-01
The analysis explains the basis set superposition error (BSSE) and fragment relaxation involved in calculating the interaction energies using various first principle theories. Interacting the correlated fragment and increasing the size of the basis set can help in decreasing the BSSE to a great extent.
Basis Set Dependence of Vibrational Raman and Raman Optical Activity Intensities.
Cheeseman, James R; Frisch, Michael J
2011-10-11
We present a systematic study of the basis set dependence of the backscattering vibrational Raman intensities and Raman Optical Activity (ROA) intensity differences. The accuracies of computed Raman intensities and ROA intensity differences for a series of commonly used basis sets are reported, relative to large reference basis sets, using the B3LYP density functional. This study attempts to separately quantify the relative accuracies obtained from particular basis set combinations: one for the geometry optimization and force field computation and the other for the computation of Raman and ROA tensors. We demonstrate here that the basis set requirements for the geometry and force fields are not similar to those of the Raman and ROA tensors. The Raman and ROA tensors require basis sets with diffuse functions, while geometry optimizations and force field computations typically do not. Eleven molecules were examined: (S)-methyloxirane, (S)-methylthirane, (R)-epichlorhydrin, (S)-CHFClBr, (1S,5S)-α-pinene, (1S,5S)-β-pinene, (1S,4S)-norborneneone, (M)-σ-[4]-helicene, an enone precursor to a cytotoxic sesquiterpene, the gauche-gauche conformer of the monosaccharide methyl-β-d-glucopyranose, and the dipeptide Ac-(alanine)2-NH2. For the molecules examined here, intensities and intensity differences obtained from Raman and ROA tensors computed using the aug-cc-pVDZ basis set are nearly equivalent to those computed with the larger aug-cc-pVTZ basis set. We find that modifying the aug-cc-pVDZ basis set by removing the set of diffuse d functions on all atoms (while keeping the diffuse s and p sets), denoted as aug(sp)-cc-pVDZ, results in a basis set which is significantly faster without much reduction in the overall accuracy. In addition, the popular rDPS basis set introduced by Zuber and Hug offers a good compromise between accuracy and efficiency. The combination of either the aug(sp)-pVDZ or rDPS basis for the computation of the Raman and ROA tensors with the 6-31G
Lehtola, Susi; Manninen, Pekka; Hakala, Mikko; Hämäläinen, Keijo
2013-01-28
Completeness-optimization is a novel method for the formation of one-electron basis sets. Contrary to conventional methods of basis set generation that optimize the basis set with respect to ground-state energy, completeness-optimization is a completely general, black-box method that can be used to form cost-effective basis sets for any wanted property at any level of theory. In our recent work [J. Lehtola, P. Manninen, M. Hakala, and K. Hämäläinen, J. Chem. Phys. 137, 104105 (2012)] we applied the completeness-optimization approach to forming primitive basis sets tuned for calculations of the electron momentum density at the Hartree-Fock (HF) level of theory. The current work extends the discussion to contracted basis sets and to the post-HF level of theory. Contractions are found to yield significant reductions in the amount of functions without compromising the accuracy. We suggest polarization-consistent and correlation-consistent basis sets for the first three rows of the periodic table, which are completeness-optimized for electron momentum density calculations. PMID:23387570
Dynamical basis sets for algebraic variational calculations in quantum-mechanical scattering theory
Sun, Yan; Kouri, Donald J.; Truhlar, Donald G.; Schwenke, David W.
1990-01-01
New basis sets are proposed for linear algebraic variational calculations of transition amplitudes in quantum-mechanical scattering problems. These basis sets are hybrids of those that yield the Kohn variational principle (KVP) and those that yield the generalized Newton variational principle (GNVP) when substituted in Schlessinger's stationary expression for the T operator. Trial calculations show that efficiencies almost as great as that of the GNVP and much greater than the KVP can be obtained, even for basis sets with the majority of the members independent of energy.
Zaleśny, Robert; Baranowska-Łączkowska, Angelika; Medveď, Miroslav; Luis, Josep M
2015-09-01
In the present work, we perform an assessment of several property-oriented atomic basis sets in computing (hyper)polarizabilities with a focus on the vibrational contributions. Our analysis encompasses the Pol and LPol-ds basis sets of Sadlej and co-workers, the def2-SVPD and def2-TZVPD basis sets of Rappoport and Furche, and the ORP basis set of Baranowska-Łączkowska and Łączkowski. Additionally, we use the d-aug-cc-pVQZ and aug-cc-pVTZ basis sets of Dunning and co-workers to determine the reference estimates of the investigated electric properties for small- and medium-sized molecules, respectively. We combine these basis sets with ab initio post-Hartree-Fock quantum-chemistry approaches (including the coupled cluster method) to calculate electronic and nuclear relaxation (hyper)polarizabilities of carbon dioxide, formaldehyde, cis-diazene, and a medium-sized Schiff base. The primary finding of our study is that, among all studied property-oriented basis sets, only the def2-TZVPD and ORP basis sets yield nuclear relaxation (hyper)polarizabilities of small molecules with average absolute errors less than 5.5%. A similar accuracy for the nuclear relaxation (hyper)polarizabilites of the studied systems can also be reached using the aug-cc-pVDZ basis set (5.3%), although for more accurate calculations of vibrational contributions, i.e., average absolute errors less than 1%, the aug-cc-pVTZ basis set is recommended. It was also demonstrated that anharmonic contributions to first and second hyperpolarizabilities of a medium-sized Schiff base are particularly difficult to accurately predict at the correlated level using property-oriented basis sets. For instance, the value of the nuclear relaxation first hyperpolarizability computed at the MP2/def2-TZVPD level of theory is roughly 3 times larger than that determined using the aug-cc-pVTZ basis set. We link the failure of the def2-TZVPD basis set with the difficulties in predicting the first-order field
DEFF Research Database (Denmark)
Hedegård, Erik D.; Jensen, Frank; Kongsted, Jacob
2012-01-01
Even for pure substances, the deduction of the absolute configuration is not always straightforward since there is no direct link between the magnitude and sign of the optical rotation and the absolute configuration. It would be very useful to use computations of the optical rotation to link...... experimentally measured optical rotations to an absolute configuration. Such electronic structure calculations of the optical rotation typically employ regular energy optimized basis sets from wave function theory, and especially the aug-cc-pVDZ basis set has been popular. Here, we have carried out extrapolation...... of the optical rotation to the basis set limits for nine small or medium sized molecules, using basis sets developed specifically for DFT and magnetic properties (aug-pcS-n series). We suggest that assignment of absolute configuration by comparisons between theoretical and experimental optical...
Gaussian Basis Set Optimization for Excited Helium and Helium-like Ions
Czech Academy of Sciences Publication Activity Database
Šmydke, Jan; Kaprálová-Žďánská, Petra Ruth
Mariapfarr : Univerzity of Graz, 2012. [Central European Symposium on Theoretical Chemistry /11./. 25.09.2011-28.09.2011, Mariapfarr] Institutional support: RVO:61388955 Keywords : helium * Gaussian basis set oprimization Subject RIV: CF - Physical ; Theoretical Chemistry
Boese, A D; Martin, J M L; Klopper, Wim; Martin, Jan M. L.
2005-01-01
In a previous contribution (Mol. Phys. {\\bf 103}, xxxx, 2005), we established the suitability of density functional theory (DFT) for the calculation of molecular anharmonic force fields. In the present work, we have assessed a wide variety of basis sets and exchange-correlation functionals for harmonic and fundamental frequencies, equilibrium and ground-state rotational constants, and thermodynamic functions beyond the RRHO (rigid rotor-harmonic oscillator) approximation. The fairly good performance of double-zeta plus polarization basis sets for frequencies results from an error compensation between basis set incompleteness and the intrinsic error of exchange-correlation functionals. Triple-zeta plus polarization basis sets are recommended, with an additional high-exponent $d$ function on second-row atoms. All conventional hybrid GGA functionals perform about equally well: high-exchange hybrid GGA and meta-GGA functionals designed for kinetics yield poor results, with the exception of of the very recently de...
Converging Higher Education Systems in a Global Setting: The Example of France and India
Pilkington, Marc
2014-01-01
We present a comparative survey between the French and Indian higher education systems. In spite of their respective idiosyncratic features, we show that the two countries have both evolved comprehensively toward a knowledge-based society, in order to ensure the prosperity of their citizens. Secondly, we single out a threefold convergence between…
Convergence and the Cauchy Property of Sequences in the Setting of Actual Infinity
Shipman, Barbara A.
2013-01-01
Traditional definitions, language, and visualizations of convergence and the Cauchy property of sequences convey a sense of the sequence as a potentially infinite process rather than an actually infinite object. This has a deep-rooted influence on how we think about and teach concepts on sequences, particularly in undergraduate calculus and…
The Bethe Sum Rule and Basis Set Selection in the Calculation of Generalized Oscillator Strengths
DEFF Research Database (Denmark)
Cabrera-Trujillo, Remigio; Sabin, John R.; Oddershede, Jens; Sauer, Stephan P. A.
1999-01-01
Fulfillment of the Bethe sum rule may be construed as a measure of basis set quality for atomic and molecular properties involving the generalized oscillator strength distribution. It is first shown that, in the case of a complete basis, the Bethe sum rule is fulfilled exactly in the random phase...
Informatics-Based Energy Fitting Scheme for Correlation Energy at Complete Basis Set Limit.
Seino, Junji; Nakai, Hiromi
2016-09-30
Energy fitting schemes based on informatics techniques using hierarchical basis sets with small cardinal numbers were numerically investigated to estimate correlation energies at the complete basis set limits. Numerical validations confirmed that the conventional two-point extrapolation models can be unified into a simple formula with optimal parameters obtained by the same test sets. The extrapolation model was extended to two-point fitting models by a relaxation of the relationship between the extrapolation coefficients or a change of the fitting formula. Furthermore, n-scheme fitting models were developed by the combinations of results calculated at several theory levels and basis sets to compensate for the deficiencies in the fitting model at one level of theory. Systematic assessments on the Gaussian-3X and Gaussian-2 sets revealed that the fitting models drastically reduced errors with equal or smaller computational effort. © 2016 Wiley Periodicals, Inc. PMID:27454327
A novel Gaussian-Sinc mixed basis set for electronic structure calculations
International Nuclear Information System (INIS)
A Gaussian-Sinc basis set methodology is presented for the calculation of the electronic structure of atoms and molecules at the Hartree–Fock level of theory. This methodology has several advantages over previous methods. The all-electron electronic structure in a Gaussian-Sinc mixed basis spans both the “localized” and “delocalized” regions. A basis set for each region is combined to make a new basis methodology—a lattice of orthonormal sinc functions is used to represent the “delocalized” regions and the atom-centered Gaussian functions are used to represent the “localized” regions to any desired accuracy. For this mixed basis, all the Coulomb integrals are definable and can be computed in a dimensional separated methodology. Additionally, the Sinc basis is translationally invariant, which allows for the Coulomb singularity to be placed anywhere including on lattice sites. Finally, boundary conditions are always satisfied with this basis. To demonstrate the utility of this method, we calculated the ground state Hartree–Fock energies for atoms up to neon, the diatomic systems H2, O2, and N2, and the multi-atom system benzene. Together, it is shown that the Gaussian-Sinc mixed basis set is a flexible and accurate method for solving the electronic structure of atomic and molecular species
An effective scheme for selecting basis sets for ab initio calculations
Institute of Scientific and Technical Information of China (English)
张瑞勤; 黄建华; 步宇翔; 韩克利; 李述汤; 何国钟
2000-01-01
An effective scheme for selecting economical basis sets for ab initio calculations has been proposed for wide-range systems based on the analysis of different functions in the currently used basis sets. Accordingly, the selection of the basis sets should be made according to the different properties and real chemical surrounding of the atoms in the systems. For normal systems, the size and level of the basis sets used for the descriptions of the constituent atoms should be increased from left to right according to the position of the atom in the periodic table. Moreover, the more the atom is negatively charged, the more the basis functions and suitable polarization functions and diffuse functions should be utilized. whereas, for the positively charged atoms, the size of basis set may be reduced. it is not necessary to use the polarization and diffuse functions for the covalently saturated atoms with normal valence states. However, for the system involving hy-drogen-bonding, weak interactions, functional
Representability of Bloch states on Projector-augmented-wave (PAW) basis sets
Agapito, Luis; Ferretti, Andrea; Curtarolo, Stefano; Buongiorno Nardelli, Marco
2015-03-01
Design of small, yet `complete', localized basis sets is necessary for an efficient dual representation of Bloch states on both plane-wave and localized basis. Such simultaneous dual representation permits the development of faster more accurate (beyond DFT) electronic-structure methods for atomistic materials (e.g. the ACBN0 method.) by benefiting from algorithms (real and reciprocal space) and hardware acceleration (e.g. GPUs) used in the quantum-chemistry and solid-state communities. Finding a `complete' atomic-orbital basis (partial waves) is also a requirement in the generation of robust and transferable PAW pseudopotentials. We have employed the atomic-orbital basis from available PAW data sets, which extends through most of the periodic table, and tested the representability of Bloch states on such basis. Our results show that PAW data sets allow systematic and accurate representability of the PAW Bloch states, better than with traditional quantum-chemistry double-zeta- and double-zeta-polarized-quality basis sets.
An effective scheme for selecting basis sets for ab initio calculations
Institute of Scientific and Technical Information of China (English)
无
2000-01-01
An effective scheme for selecting economical basis sets for ab initio calculations has been proposed for wide-range systems based on the analysis of different functions in the currently used basis sets. Accordingly, the selection of the basis sets should be made according to the different properties and real chemical surrounding of the atoms in the systems. For normal systems, the size and level of the basis sets used for the descriptions of the constituent atoms should be increased from left to right according to the position of the atom in the periodic table. Moreover, the more the atom is negatively charged, the more the basis functions and suitable polarization functions and diffuse functions should be utilized. Whereas, for the positively charged atoms, the size of basis set may be reduced. It is not necessary to use the polarization and diffuse functions for the covalently saturated atoms with normal valence states. However, for the system involving hydrogen-bonding, weak interactions, functional groups, metallic bonding with zero valence or low positive valence, and other sensitive interactions, the polarization and diffuse functions must be used. With this scheme, reliable calculation results may be obtained with suitable basis sets and smaller computational capability. By detailed analysis of a series of systems, it has been demonstrated that this scheme is very practical and effective. This scheme may be used in Hartree-Fock, M?ller-Plesset and density functional theoretical calculations. It is expected that the scheme would find important applications in the extensive calculations of large systems in chemistry, materials science, and life and biological sciences.
Basis set effects on the Hartree-Fock description of confined many-electron atoms
International Nuclear Information System (INIS)
In this work, the basis sets designed by Clementi, Bunge and Thakkar, for atomic systems, have been used to obtain the electronic structure of confined many-electron atoms by using Roothaan's approach in the Hartree-Fock context with a new code written in C, which uses the message-passing interface library. The confinement was imposed as Ludeña suggested to simulate walls with infinity potential. For closed-shell atoms, the Thakkar basis set functions give the best total energies (TE) as a function of the confinement radius, obtaining the following ordering: TE(Thakkar) < TE(Bunge) < TE(Clementi). However, for few open-shell atoms this ordering is not preserved and a trend, for the basis sets, is not observed. Although there are differences between the TE predicted by these basis set functions, the corresponding pressures are similar to each other; it means that changes in the total energy are described almost in the same way by using any of these basis sets. By analysing the total energy as a function of the inverse of the volume we propose an equation of state; for regions of small volumes, this equation predicts that the pressure is inversely proportional to the square of the volume.
Strategies for reducing basis set superposition error (BSSE) in O/AU and O/Ni
Shuttleworth, I. G.
2015-11-01
The effect of basis set superposition error (BSSE) and effective strategies for the minimisation have been investigated using the SIESTA-LCAO DFT package. Variation of the energy shift parameter ΔEPAO has been shown to reduce BSSE for bulk Au and Ni and across their oxygenated surfaces. Alternative strategies based on either the expansion or contraction of the basis set have been shown to be ineffective in reducing BSSE. Comparison of the binding energies for the surface systems obtained using LCAO were compared with BSSE-free plane wave energies.
Strategies for reducing basis set superposition error (BSSE) in O/AU and O/Ni
Shuttleworth, I.G.
2015-11-01
© 2015 Elsevier Ltd. All rights reserved. The effect of basis set superposition error (BSSE) and effective strategies for the minimisation have been investigated using the SIESTA-LCAO DFT package. Variation of the energy shift parameter ΔEPAO has been shown to reduce BSSE for bulk Au and Ni and across their oxygenated surfaces. Alternative strategies based on either the expansion or contraction of the basis set have been shown to be ineffective in reducing BSSE. Comparison of the binding energies for the surface systems obtained using LCAO were compared with BSSE-free plane wave energies.
Thin-skinned salt tectonics as a response to crustal movements in a recent convergent setting
Vendeville, Bruno; Gaullier, Virginie; Deverchere, Jacques; Sage, Françoise
2014-05-01
Published data from the seismic survey "MARADJA 1" conducted in August and September 2003 offshore the Algerian coast have imaged the Messinian salt response to tectonic activity within the basement. This helps to understand how the vertical movements of basement blocks have been recorded by thin-skinned salt tectonics. The area is undergoing crustal convergence, as attested by the Boumerdes earthquake (magnitude 6.8), which happened in 2003. The seismic data have revealed the presence of an elevated plateau, forming a 3D promontory restricted to the area offshore Algiers and is absent west and east of that area. The promontory is likely related to the contractional reactivation of the margin, as was recorded by subsalt thrusts mapped by Domzig et al. (2006). The data provided additional information on the deformation of the Messinian mobile evaporite unit and its Plio-Quaternary overburden. Margin-perpendicular profiles show mostly compressional features (anticlines and synclines) that had little or no activity during Messinian times. By contrast, margin-parallel profiles clearly show that extensional, reactive salt diapiric ridges formed early, as early as the time of deposition of the Messinian Upper Unit, as recorded by fan-shaped strata. These ridges have recorded E-W, thin-skinned gravity gliding above the Messinian salt, as a response to the rise of the basement plateau. We tested this hypothesis using analogue models. Indeed, the rise of the plateau generated preferential E-W extension above the salt, rather than N-S, which was prevented by the compressional regional tectonic stresses related to the convergence of the African and European plates.
Investigation of confined hydrogen atom in spherical cavity, using B-splines basis set
Directory of Open Access Journals (Sweden)
M Barezi
2011-03-01
Full Text Available Studying confined quantum systems (CQS is very important in nano technology. One of the basic CQS is a hydrogen atom confined in spherical cavity. In this article, eigenenergies and eigenfunctions of hydrogen atom in spherical cavity are calculated, using linear variational method. B-splines are used as basis functions, which can easily construct the trial wave functions with appropriate boundary conditions. The main characteristics of B-spline are its high localization and its flexibility. Besides, these functions have numerical stability and are able to spend high volume of calculation with good accuracy. The energy levels as function of cavity radius are analyzed. To check the validity and efficiency of the proposed method, extensive convergence test of eigenenergies in different cavity sizes has been carried out.
Borges, A; Solomon, G C
2016-05-21
Single molecule conductance measurements are often interpreted through computational modeling, but the complexity of these calculations makes it difficult to directly link them to simpler concepts and models. Previous work has attempted to make this connection using maximally localized Wannier functions and symmetry adapted basis sets, but their use can be ambiguous and non-trivial. Starting from a Hamiltonian and overlap matrix written in a hydrogen-like basis set, we demonstrate a simple approach to obtain a new basis set that is chemically more intuitive and allows interpretation in terms of simple concepts and models. By diagonalizing the Hamiltonians corresponding to each atom in the molecule, we obtain a basis set that can be partitioned into pseudo-σ and -π and allows partitioning of the Landuaer-Büttiker transmission as well as create simple Hückel models that reproduce the key features of the full calculation. This method provides a link between complex calculations and simple concepts and models to provide intuition or extract parameters for more complex model systems. PMID:27208940
On the use of locally dense basis sets in the calculation of EPR hyperfine couplings
DEFF Research Database (Denmark)
Milhøj, Birgitte Olai; Hedegård, Erik D.; Sauer, Stephan P. A.
2013-01-01
The usage of locally dense basis sets in the calculation of Electron Paramagnetic Resonance (EPR) hyperne coupling constants is investigated at the level of Density Functional Theory (DFT) for two model systems of biologically important transition metal complexes: One for the active site in the c...
On the use of Locally Dense Basis Sets in the Calculation of EPR Hyperfine Couplings
DEFF Research Database (Denmark)
Hedegård, Erik D.; Sauer, Stephan P. A.; Milhøj, Birgitte O.
2013-01-01
The usage of locally dense basis sets in the calculation of Electron Paramagnetic Resonance (EPR) hyperne coupling constants is investigated at the level of Density Functional Theory (DFT) for two model systems of biologically important transition metal complexes: One for the active site in the c...
Ca+ centers for hydrogen storage: An accurate many-body study with large basis sets
Purwanto, Wirawan; Virgus, Yudistira; Zhang, Shiwei
2011-01-01
Weak H2 physisorption energies present a significant challenge to even the best correlated theoretical many-body methods. We use the phaseless auxiliary-field quantum Monte Carlo (AFQMC) method to accurately predict the binding energy of Ca+ - 4H2. Attention has recently focused on this model chemistry to test the reliability of electronic structure methods for H2 binding on dispersed alkaline earth metal centers. A modified Cholesky decomposition is implemented to realize the Hubbard-Stratonovich transformation efficiently with large Gaussian basis sets. We employ the largest correlation-consistent Gaussian type basis sets available, up to cc-pCV5Z for Ca, to accurately extrapolate to the complete basis limit. The calculated potential energy curve exhibits binding with a double-well structure.
An editor for the maintenance and use of a bank of contracted Gaussian basis set functions
International Nuclear Information System (INIS)
A bank of basis sets to be used in ab-initio calculations has been created. The bases are sets of contracted Gaussian type orbitals to be used as input to any molecular integral package. In this communication we shall describe the organization of the bank and a portable editor program which was designed for its maintenance and use. This program is operated by commands and it may be used to obtain any kind of information about the bases in the bank as well as to produce output to be directly used as input for different integral programs. The editor may also be used to format basis sets in the conventional way utilized in publications, as well as to generate a complete, or partial, manual of the contents of the bank if so desired. (orig.)
Mougeot, Mathilde; Tribouley, Karine
2011-01-01
This paper investigates the problem of selecting variables in regression-type models for an "instrumental" setting. Our study is motivated by empirically verifying the conditional convergence hypothesis used in the economical literature concerning the growth rate. To avoid unnecessary discussion about the choice and the pertinence of instrumental variables, we embed the model in a very high dimensional setting. We propose a selection procedure with no optimization step called LOLA, for Learning Out of Leaders with Adaptation. LOLA is an auto-driven algorithm with two thresholding steps. The consistency of the procedure is proved under sparsity conditions and simulations are conducted to illustrate the practical good performances of LOLA. The behavior of the algorithm is studied when instrumental variables are artificially added without a priori significant connection to the model. Using our algorithm, we provide a solution for modeling the link between the growth rate and the initial level of the gross domest...
Truncated spherical-wave basis set for first-principles pseudopotential calculations
Monserrat, B.; Haynes, P. D.
2010-11-01
Analytic results for two- and three-centre integrals are derived for the truncated spherical-wave basis set designed for first-principles pseudopotential calculations within density-functional theory. These allow the overlap, kinetic energy and non-local pseudopotential matrix elements to be calculated efficiently and accurately. In particular, the scaling of the computational effort with maximum angular momentum component is dramatically improved and the projection method takes full account of the discontinuities in the basis functions arising from their localization within spherical regions.
Nuray, Fatih; Gök, Hafize; ULUSU, Uğur
2011-01-01
In this paper, the concepts of $sigma$-uniform density of subsets A of the set $mathbb{N}$ of positive integers and corresponding $I_{sigma}$-convergence were introduced. Furthermore, inclusion relations between $I_{sigma}$-convergence and invariant convergence also $I_{sigma}$-convergence and $[V_{sigma}]_{p}$-convergence were given.
Evarestov, R A; Losev, M V
2009-12-01
For the first time the convergence of the phonon frequencies and dispersion curves in terms of the supercell size is studied in ab initio frozen phonon calculations on LiF crystal. Helmann-Feynman forces over atomic displacements are found in all-electron calculations with the localized atomic functions (LCAO) basis using CRYSTAL06 program. The Parlinski-Li-Kawazoe method and FROPHO program are used to calculate the dynamical matrix and phonon frequencies of the supercells. For fcc lattice, it is demonstrated that use of the full supercell space group (including the supercell inner translations) enables to reduce essentially the number of the displacements under consideration. For Hartree-Fock (HF), PBE and hybrid PBE0, B3LYP, and B3PW exchange-correlation functionals the atomic basis set optimization is performed. The supercells up to 216 atoms (3 x 3 x 3 conventional unit cells) are considered. The phonon frequencies using the supercells of different size and shape are compared. For the commensurate with supercell k-points the best agreement of the theoretical results with the experimental data is found for B3PW exchange-correlation functional calculations with the optimized basis set. The phonon frequencies at the most non-commensurate k-points converged for the supercell consisting of 4 x 4 x 4 primitive cells and ensures the accuracy 1-2% in the thermodynamic properties calculated (the Helmholtz free energy, entropy, and heat capacity at the room temperature). PMID:19382176
Structural Identification of Nonlinear Static System on Basis of Analysis Sector Sets
Directory of Open Access Journals (Sweden)
Nikolay Karabutov
2013-12-01
Full Text Available Methods of structural identification of static systems with a vector input and several nonlinearities in the conditions of uncertainty are considered. We consider inputs irregular. The concept of structural space is introduced. In this space special structures (virtual portraits are analyzed. The Holder condition is applied to construction of sector set, to which belongs a virtual portrait of system of identification. Criteria of decision-making on a class of nonlinear functions on the basis of the analysis of proximity of sector sets are described. Procedures of an estimation of structural parameters of two classes of nonlinearities are stated: power and a hysteresis.
Boese, A. Daniel; Klopper, Wim; Martin, Jan M. L.
2005-01-01
In a previous contribution (Mol. Phys. {\\bf 103}, xxxx, 2005), we established the suitability of density functional theory (DFT) for the calculation of molecular anharmonic force fields. In the present work, we have assessed a wide variety of basis sets and exchange-correlation functionals for harmonic and fundamental frequencies, equilibrium and ground-state rotational constants, and thermodynamic functions beyond the RRHO (rigid rotor-harmonic oscillator) approximation. The fairly good perf...
Geminal embedding scheme for optimal atomic basis set construction in correlated calculations
International Nuclear Information System (INIS)
We introduce an efficient method to construct optimal and system adaptive basis sets for use in electronic structure and quantum Monte Carlo calculations. The method is based on an embedding scheme in which a reference atom is singled out from its environment, while the entire system (atom and environment) is described by a Slater determinant or its antisymmetrized geminal power (AGP) extension. The embedding procedure described here allows for the systematic and consistent contraction of the primitive basis set into geminal embedded orbitals (GEOs), with a dramatic reduction of the number of variational parameters necessary to represent the many-body wave function, for a chosen target accuracy. Within the variational Monte Carlo method, the Slater or AGP part is determined by a variational minimization of the energy of the whole system in presence of a flexible and accurate Jastrow factor, representing most of the dynamical electronic correlation. The resulting GEO basis set opens the way for a fully controlled optimization of many-body wave functions in electronic structure calculation of bulk materials, namely, containing a large number of electrons and atoms. We present applications on the water molecule, the volume collapse transition in cerium, and the high-pressure liquid hydrogen
Geminal embedding scheme for optimal atomic basis set construction in correlated calculations.
Sorella, S; Devaux, N; Dagrada, M; Mazzola, G; Casula, M
2015-12-28
We introduce an efficient method to construct optimal and system adaptive basis sets for use in electronic structure and quantum Monte Carlo calculations. The method is based on an embedding scheme in which a reference atom is singled out from its environment, while the entire system (atom and environment) is described by a Slater determinant or its antisymmetrized geminal power (AGP) extension. The embedding procedure described here allows for the systematic and consistent contraction of the primitive basis set into geminal embedded orbitals (GEOs), with a dramatic reduction of the number of variational parameters necessary to represent the many-body wave function, for a chosen target accuracy. Within the variational Monte Carlo method, the Slater or AGP part is determined by a variational minimization of the energy of the whole system in presence of a flexible and accurate Jastrow factor, representing most of the dynamical electronic correlation. The resulting GEO basis set opens the way for a fully controlled optimization of many-body wave functions in electronic structure calculation of bulk materials, namely, containing a large number of electrons and atoms. We present applications on the water molecule, the volume collapse transition in cerium, and the high-pressure liquid hydrogen. PMID:26723656
Geminal embedding scheme for optimal atomic basis set construction in correlated calculations
Energy Technology Data Exchange (ETDEWEB)
Sorella, S., E-mail: sorella@sissa.it [International School for Advanced Studies (SISSA), Via Beirut 2-4, 34014 Trieste, Italy and INFM Democritos National Simulation Center, Trieste (Italy); Devaux, N.; Dagrada, M., E-mail: mario.dagrada@impmc.upmc.fr [Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, Case 115, 4 Place Jussieu, 75252 Paris Cedex 05 (France); Mazzola, G., E-mail: gmazzola@phys.ethz.ch [Theoretische Physik, ETH Zurich, 8093 Zurich (Switzerland); Casula, M., E-mail: michele.casula@impmc.upmc.fr [CNRS and Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, Case 115, 4 Place Jussieu, 75252 Paris Cedex 05 (France)
2015-12-28
We introduce an efficient method to construct optimal and system adaptive basis sets for use in electronic structure and quantum Monte Carlo calculations. The method is based on an embedding scheme in which a reference atom is singled out from its environment, while the entire system (atom and environment) is described by a Slater determinant or its antisymmetrized geminal power (AGP) extension. The embedding procedure described here allows for the systematic and consistent contraction of the primitive basis set into geminal embedded orbitals (GEOs), with a dramatic reduction of the number of variational parameters necessary to represent the many-body wave function, for a chosen target accuracy. Within the variational Monte Carlo method, the Slater or AGP part is determined by a variational minimization of the energy of the whole system in presence of a flexible and accurate Jastrow factor, representing most of the dynamical electronic correlation. The resulting GEO basis set opens the way for a fully controlled optimization of many-body wave functions in electronic structure calculation of bulk materials, namely, containing a large number of electrons and atoms. We present applications on the water molecule, the volume collapse transition in cerium, and the high-pressure liquid hydrogen.
Vogel, Thomas
2015-01-01
We recently introduced a novel replica-exchange scheme in which an individual replica can sample from states encountered by other replicas at any previous time by way of a global configuration database, enabling the fast propagation of relevant states through the whole ensemble of replicas. This mechanism depends on the knowledge of global thermodynamic functions which are measured during the simulation and not coupled to the heat bath temperatures driving the individual simulations. Therefore, this setup also allows for a continuous adaptation of the temperature set. In this paper, we will review the new scheme and demonstrate its capability. The method is particularly useful for the fast and reliable estimation of the microcanonical temperature T(U) or, equivalently, of the density of states g(U) over a wide range of energies.
On fast computation of finite-time coherent sets using radial basis functions.
Froyland, Gary; Junge, Oliver
2015-08-01
Finite-time coherent sets inhibit mixing over finite times. The most expensive part of the transfer operator approach to detecting coherent sets is the construction of the operator itself. We present a numerical method based on radial basis function collocation and apply it to a recent transfer operator construction [G. Froyland, "Dynamic isoperimetry and the geometry of Lagrangian coherent structures," Nonlinearity (unpublished); preprint arXiv:1411.7186] that has been designed specifically for purely advective dynamics. The construction [G. Froyland, "Dynamic isoperimetry and the geometry of Lagrangian coherent structures," Nonlinearity (unpublished); preprint arXiv:1411.7186] is based on a "dynamic" Laplace operator and minimises the boundary size of the coherent sets relative to their volume. The main advantage of our new approach is a substantial reduction in the number of Lagrangian trajectories that need to be computed, leading to large speedups in the transfer operator analysis when this computation is costly. PMID:26328580
International Nuclear Information System (INIS)
A density functional theory (DFT) based composite electronic structure approach is proposed to efficiently compute structures and interaction energies in large chemical systems. It is based on the well-known and numerically robust Perdew-Burke-Ernzerhoff (PBE) generalized-gradient-approximation in a modified global hybrid functional with a relatively large amount of non-local Fock-exchange. The orbitals are expanded in Ahlrichs-type valence-double zeta atomic orbital (AO) Gaussian basis sets, which are available for many elements. In order to correct for the basis set superposition error (BSSE) and to account for the important long-range London dispersion effects, our well-established atom-pairwise potentials are used. In the design of the new method, particular attention has been paid to an accurate description of structural parameters in various covalent and non-covalent bonding situations as well as in periodic systems. Together with the recently proposed three-fold corrected (3c) Hartree-Fock method, the new composite scheme (termed PBEh-3c) represents the next member in a hierarchy of “low-cost” electronic structure approaches. They are mainly free of BSSE and account for most interactions in a physically sound and asymptotically correct manner. PBEh-3c yields good results for thermochemical properties in the huge GMTKN30 energy database. Furthermore, the method shows excellent performance for non-covalent interaction energies in small and large complexes. For evaluating its performance on equilibrium structures, a new compilation of standard test sets is suggested. These consist of small (light) molecules, partially flexible, medium-sized organic molecules, molecules comprising heavy main group elements, larger systems with long bonds, 3d-transition metal systems, non-covalently bound complexes (S22 and S66×8 sets), and peptide conformations. For these sets, overall deviations from accurate reference data are smaller than for various other tested DFT
Energy Technology Data Exchange (ETDEWEB)
Grimme, Stefan, E-mail: grimme@thch.uni-bonn.de; Brandenburg, Jan Gerit; Bannwarth, Christoph; Hansen, Andreas [Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn (Germany)
2015-08-07
A density functional theory (DFT) based composite electronic structure approach is proposed to efficiently compute structures and interaction energies in large chemical systems. It is based on the well-known and numerically robust Perdew-Burke-Ernzerhoff (PBE) generalized-gradient-approximation in a modified global hybrid functional with a relatively large amount of non-local Fock-exchange. The orbitals are expanded in Ahlrichs-type valence-double zeta atomic orbital (AO) Gaussian basis sets, which are available for many elements. In order to correct for the basis set superposition error (BSSE) and to account for the important long-range London dispersion effects, our well-established atom-pairwise potentials are used. In the design of the new method, particular attention has been paid to an accurate description of structural parameters in various covalent and non-covalent bonding situations as well as in periodic systems. Together with the recently proposed three-fold corrected (3c) Hartree-Fock method, the new composite scheme (termed PBEh-3c) represents the next member in a hierarchy of “low-cost” electronic structure approaches. They are mainly free of BSSE and account for most interactions in a physically sound and asymptotically correct manner. PBEh-3c yields good results for thermochemical properties in the huge GMTKN30 energy database. Furthermore, the method shows excellent performance for non-covalent interaction energies in small and large complexes. For evaluating its performance on equilibrium structures, a new compilation of standard test sets is suggested. These consist of small (light) molecules, partially flexible, medium-sized organic molecules, molecules comprising heavy main group elements, larger systems with long bonds, 3d-transition metal systems, non-covalently bound complexes (S22 and S66×8 sets), and peptide conformations. For these sets, overall deviations from accurate reference data are smaller than for various other tested DFT
Grimme, Stefan; Brandenburg, Jan Gerit; Bannwarth, Christoph; Hansen, Andreas
2015-08-01
A density functional theory (DFT) based composite electronic structure approach is proposed to efficiently compute structures and interaction energies in large chemical systems. It is based on the well-known and numerically robust Perdew-Burke-Ernzerhoff (PBE) generalized-gradient-approximation in a modified global hybrid functional with a relatively large amount of non-local Fock-exchange. The orbitals are expanded in Ahlrichs-type valence-double zeta atomic orbital (AO) Gaussian basis sets, which are available for many elements. In order to correct for the basis set superposition error (BSSE) and to account for the important long-range London dispersion effects, our well-established atom-pairwise potentials are used. In the design of the new method, particular attention has been paid to an accurate description of structural parameters in various covalent and non-covalent bonding situations as well as in periodic systems. Together with the recently proposed three-fold corrected (3c) Hartree-Fock method, the new composite scheme (termed PBEh-3c) represents the next member in a hierarchy of "low-cost" electronic structure approaches. They are mainly free of BSSE and account for most interactions in a physically sound and asymptotically correct manner. PBEh-3c yields good results for thermochemical properties in the huge GMTKN30 energy database. Furthermore, the method shows excellent performance for non-covalent interaction energies in small and large complexes. For evaluating its performance on equilibrium structures, a new compilation of standard test sets is suggested. These consist of small (light) molecules, partially flexible, medium-sized organic molecules, molecules comprising heavy main group elements, larger systems with long bonds, 3d-transition metal systems, non-covalently bound complexes (S22 and S66×8 sets), and peptide conformations. For these sets, overall deviations from accurate reference data are smaller than for various other tested DFT methods
Symmetry-adapted basis sets automatic generation for problems in chemistry and physics
Avery, John Scales; Avery, James Emil
2012-01-01
In theoretical physics, theoretical chemistry and engineering, one often wishes to solve partial differential equations subject to a set of boundary conditions. This gives rise to eigenvalue problems of which some solutions may be very difficult to find. For example, the problem of finding eigenfunctions and eigenvalues for the Hamiltonian of a many-particle system is usually so difficult that it requires approximate methods, the most common of which is expansion of the eigenfunctions in terms of basis functions that obey the boundary conditions of the problem. The computational effort needed
Corsetti, Fabiano
2014-01-01
The implementation of the orbital minimization method (OMM) for solving the self-consistent Kohn-Sham (KS) problem for electronic structure calculations in a basis of non-orthogonal numerical atomic orbitals of finite-range is reported. We explore the possibilities for using the OMM as an exact cubic-scaling solver for the KS problem, and compare its performance with that of explicit diagonalization in realistic systems. We analyze the efficiency of the method depending on the choice of line search algorithm and on two free parameters, the scale of the kinetic energy preconditioning and the eigenspectrum shift. The results of several timing tests are then discussed, showing that the OMM can achieve a noticeable speedup with respect to diagonalization even for minimal basis sets for which the number of occupied eigenstates represents a significant fraction of the total basis size (>15%). We investigate the hard and soft parallel scaling of the method on multiple cores, finding a performance equal to or better ...
Atomic self-consistent-field program by the basis set expansion method: Columbus version
Pitzer, Russell M.
2005-08-01
A revised and extended (Columbus) version of the Chicago atomic self-consistent-field (Hartree-Fock) program of 1963 is described. Its principal present use is in developing Gaussian basis sets for molecular calculations. Complete memory allocation (using Fortran 90) has been added as well as improved integral formulas and efficient and simple programming features. Energy-expression coefficients have been added sufficient to treat the ground states of all atoms to the extent that Russell-Saunders (LS) coupling applies. Excited states with large angular-momentum orbitals can be treated. Relativistic effects can be included to the extent possible with relativistic effective core potentials. A review of earlier work is included. Program summaryProgram title: atmscf Catalogue identifier: ADVR Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADVR Program available from: CPC Program Library, Queen's University of Belfast, N. Ireland Programming language: Fortran 90 Computer: Sun, SGI, PC Operating system: Solaris, Irix, Linux RAM: 10 Mbytes No. of lines in distributed program, including test data, etc.: 2113 No. of bytes in distributed program, including test data, etc.: 15 379 Distribution format: tar.gz Nature of problem: Energies and wave functions, at the Hartree-Fock level Solution method: Expansions in Gaussian or Slater functions. Iterative minimization of the total energy. Optimization of exponential parameters. Used frequently for developing Gaussian basis sets for molecular use Running time: Typical 30 s per calculation
International Nuclear Information System (INIS)
The universal basis set of de Castro and Jorge has been investigated in matrix Hartree - Fock calculations, without truncation for specific atomic species, in electron correlation studies, and in systematic schemes for basis set extension. copyright 1999 American Institute of Physics
DEFF Research Database (Denmark)
Silva-Junior, Mario R.; Sauer, Stephan P. A.; Schreiber, Marko; Thiel, Walter
Vertical electronic excitation energies and one-electron properties of 28 medium-sized molecules from a previously proposed benchmark set are revisited using the augmented correlation-consistent triple-zeta aug-cc-pVTZ basis set in CC2, CCSDR(3), and CC3 calculations. The results are compared to...... those obtained previously with the smaller TZVP basis set. For each of the three coupled cluster methods, a correlation coefficient greater than 0.994 is found between the vertical excitation energies computed with the two basis sets. The deviations of the CC2 and CCSDR(3) results from the CC3 reference...
International Nuclear Information System (INIS)
Full text: Faucher and co-workers have carried out several studies of the effect of electron correlation on the crystal-field splitting of 4fN configurations of lanthanide ions in solids. In their approach the 4fN configuration is supplemented by the 4fN-1 6p configuration. On the other hand, Reid, Burdick, and coworkers have carried out many parameter fits using pure 4fN model spaces, with the one-electron crystal field operators supplemented with two-electron correlation crystal field operators. Both approaches lead to better fits to the experimental data, but since the basis sets are different, it is difficult to compare the resulting parameters. In this work we use the matrices generated by Faucher's calculations to calculate the equivalent two-electron correlation crystal field parametrization in the 4fN model space
A two-dimensional volatility basis set – Part 3: Prognostic modeling and NOx dependence
Directory of Open Access Journals (Sweden)
W. K. Chuang
2015-06-01
Full Text Available When NOx is introduced to organic emissions, aerosol production is sometimes, but not always, reduced. Under certain conditions, these interactions will instead increase aerosol concentrations. We expanded the two-dimensional volatility basis set (2-D-VBS to include the effects of NOx on aerosol formation. This includes the formation of organonitrates, where the addition of a nitrate group contributes to a decrease of 2.5 orders of magnitude in volatility. With this refinement, we model outputs from experimental results, such as the atomic N : C ratio, organonitrate mass, and nitrate fragments in AMS measurements. We also discuss the mathematical methods underlying the implementation of the 2-D-VBS and provide the complete code in the Supplemental material. A developer version is available on Bitbucket, an online community repository.
Model for the fast estimation of basis set superposition error in biomolecular systems
Faver, John C.; Zheng, Zheng; Merz, Kenneth M.
2011-01-01
Basis set superposition error (BSSE) is a significant contributor to errors in quantum-based energy functions, especially for large chemical systems with many molecular contacts such as folded proteins and protein-ligand complexes. While the counterpoise method has become a standard procedure for correcting intermolecular BSSE, most current approaches to correcting intramolecular BSSE are simply fragment-based analogues of the counterpoise method which require many (two times the number of fragments) additional quantum calculations in their application. We propose that magnitudes of both forms of BSSE can be quickly estimated by dividing a system into interacting fragments, estimating each fragment's contribution to the overall BSSE with a simple statistical model, and then propagating these errors throughout the entire system. Such a method requires no additional quantum calculations, but rather only an analysis of the system's interacting fragments. The method is described herein and is applied to a protein-ligand system, a small helical protein, and a set of native and decoy protein folds. PMID:22010701
Directory of Open Access Journals (Sweden)
Kazimierz Włodarczyk
2010-01-01
Full Text Available In cone uniform and uniform spaces, we introduce the three kinds of dissipative set-valued dynamic systems with generalized pseudodistances and not necessarily lower semicontinuous entropies, we study the convergence of dynamic processes and generalized sequences of iterations of these dissipative dynamic systems, and we establish conditions guaranteeing the existence of periodic points and endpoints of these dissipative dynamic systems and the convergence to these periodic points and endpoints of dynamic processes and generalized sequences of iterations of these dissipative dynamic systems. The paper includes examples.
European air quality modelled by CAMx including the volatility basis set scheme
Directory of Open Access Journals (Sweden)
G. Ciarelli
2015-12-01
Full Text Available Four periods of EMEP (European Monitoring and Evaluation Programme intensive measurement campaigns (June 2006, January 2007, September–October 2008 and February–March 2009 were modelled using the regional air quality model CAMx with VBS (Volatility Basis Set approach for the first time in Europe within the framework of the EURODELTA-III model intercomparison exercise. More detailed analysis and sensitivity tests were performed for the period of February–March 2009 and June 2006 to investigate the uncertainties in emissions as well as to improve the modelling of organic aerosols (OA. Model performance for selected gas phase species and PM2.5 was evaluated using the European air quality database Airbase. Sulfur dioxide (SO2 and ozone (O3 were found to be overestimated for all the four periods with O3 having the largest mean bias during June 2006 and January–February 2007 periods (8.93 and 12.30 ppb mean biases, respectively. In contrast, nitrogen dioxide (NO2 and carbon monoxide (CO were found to be underestimated for all the four periods. CAMx reproduced both total concentrations and monthly variations of PM2.5 very well for all the four periods with average biases ranging from −2.13 to 1.04 μg m-3. Comparisons with AMS (Aerosol Mass Spectrometer measurements at different sites in Europe during February–March 2009, showed that in general the model over-predicts the inorganic aerosol fraction and under-predicts the organic one, such that the good agreement for PM2.5 is partly due to compensation of errors. The effect of the choice of volatility basis set scheme (VBS on OA was investigated as well. Two sensitivity tests with volatility distributions based on previous chamber and ambient measurements data were performed. For February–March 2009 the chamber-case reduced the total OA concentrations by about 43 % on average. On the other hand, a test based on ambient measurement data increased OA concentrations by about 47 % for the same
European air quality modelled by CAMx including the volatility basis set scheme
Ciarelli, G.; Aksoyoglu, S.; Crippa, M.; Jimenez, J. L.; Nemitz, E.; Sellegri, K.; Äijälä, M.; Carbone, S.; Mohr, C.; O'Dowd, C.; Poulain, L.; Baltensperger, U.; Prévôt, A. S. H.
2015-12-01
Four periods of EMEP (European Monitoring and Evaluation Programme) intensive measurement campaigns (June 2006, January 2007, September-October 2008 and February-March 2009) were modelled using the regional air quality model CAMx with VBS (Volatility Basis Set) approach for the first time in Europe within the framework of the EURODELTA-III model intercomparison exercise. More detailed analysis and sensitivity tests were performed for the period of February-March 2009 and June 2006 to investigate the uncertainties in emissions as well as to improve the modelling of organic aerosols (OA). Model performance for selected gas phase species and PM2.5 was evaluated using the European air quality database Airbase. Sulfur dioxide (SO2) and ozone (O3) were found to be overestimated for all the four periods with O3 having the largest mean bias during June 2006 and January-February 2007 periods (8.93 and 12.30 ppb mean biases, respectively). In contrast, nitrogen dioxide (NO2) and carbon monoxide (CO) were found to be underestimated for all the four periods. CAMx reproduced both total concentrations and monthly variations of PM2.5 very well for all the four periods with average biases ranging from -2.13 to 1.04 μg m-3. Comparisons with AMS (Aerosol Mass Spectrometer) measurements at different sites in Europe during February-March 2009, showed that in general the model over-predicts the inorganic aerosol fraction and under-predicts the organic one, such that the good agreement for PM2.5 is partly due to compensation of errors. The effect of the choice of volatility basis set scheme (VBS) on OA was investigated as well. Two sensitivity tests with volatility distributions based on previous chamber and ambient measurements data were performed. For February-March 2009 the chamber-case reduced the total OA concentrations by about 43 % on average. On the other hand, a test based on ambient measurement data increased OA concentrations by about 47 % for the same period bringing model
Complete Basis Set Extrapolation and Hybrid Schemes for Geometry Gradients of Noncovalent Complexes
Czech Academy of Sciences Publication Activity Database
Černý, Jiří; Pitoňák, M.; Riley, Kevin Eugene; Hobza, Pavel
2011-01-01
Roč. 7, č. 12 (2011), s. 3924-3934. ISSN 1549-9618 R&D Projects: GA MŠk LC512 Institutional research plan: CEZ:AV0Z50520701; CEZ:AV0Z40550506 Keywords : potential-energy surface * benzene dimer * stacking * convergence Subject RIV: CF - Physical ; Theoretical Chemistry Impact factor: 5.215, year: 2011
Assessment of quantum chemical methods and basis sets for excitation energy transfer
International Nuclear Information System (INIS)
The validity of several standard quantum chemical approaches and other models for the prediction of exciton energy transfer is investigated using the HOMO-LUMO excited states of benzene dimer as an example. The configuration interaction singles (CIS), time-dependent Hartree-Fock (TD-HF), time dependent density functional theroy (TD-DFT), and complete-active-space self-consistent-field (CASSCF) methods are applied with a supermolecule approach and compared to the previously established monomer transition density method and the ideal dipole approximation. Strong and physically incorrect admixture of charge-transfer states makes TD-DFT inappropriate for investigations of potential energy surfaces in such dimer systems. CIS, TD-HF and CASSCF perform qualitatively correct. TD-HF seems to be a particularly appropriate method due to its general applicability and overall good performance for the excited state and for transition properties. Double-zeta basis sets with polarisation functions are found to be sufficient to predict transfer rates of dipole allowed excitations. Efficient excitation energy transfer is predicted between degenerate excited states while avoided curve crossings of nearly spaced π-aggregates are identified as a possible trapping mechanism
Energy Technology Data Exchange (ETDEWEB)
Holden, Zachary C.; Richard, Ryan M.; Herbert, John M., E-mail: herbert@chemistry.ohio-state.edu [Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210 (United States)
2013-12-28
An implementation of Ewald summation for use in mixed quantum mechanics/molecular mechanics (QM/MM) calculations is presented, which builds upon previous work by others that was limited to semi-empirical electronic structure for the QM region. Unlike previous work, our implementation describes the wave function's periodic images using “ChElPG” atomic charges, which are determined by fitting to the QM electrostatic potential evaluated on a real-space grid. This implementation is stable even for large Gaussian basis sets with diffuse exponents, and is thus appropriate when the QM region is described by a correlated wave function. Derivatives of the ChElPG charges with respect to the QM density matrix are a potentially serious bottleneck in this approach, so we introduce a ChElPG algorithm based on atom-centered Lebedev grids. The ChElPG charges thus obtained exhibit good rotational invariance even for sparse grids, enabling significant cost savings. Detailed analysis of the optimal choice of user-selected Ewald parameters, as well as timing breakdowns, is presented.
Holden, Zachary C.; Richard, Ryan M.; Herbert, John M.
2013-12-01
An implementation of Ewald summation for use in mixed quantum mechanics/molecular mechanics (QM/MM) calculations is presented, which builds upon previous work by others that was limited to semi-empirical electronic structure for the QM region. Unlike previous work, our implementation describes the wave function's periodic images using "ChElPG" atomic charges, which are determined by fitting to the QM electrostatic potential evaluated on a real-space grid. This implementation is stable even for large Gaussian basis sets with diffuse exponents, and is thus appropriate when the QM region is described by a correlated wave function. Derivatives of the ChElPG charges with respect to the QM density matrix are a potentially serious bottleneck in this approach, so we introduce a ChElPG algorithm based on atom-centered Lebedev grids. The ChElPG charges thus obtained exhibit good rotational invariance even for sparse grids, enabling significant cost savings. Detailed analysis of the optimal choice of user-selected Ewald parameters, as well as timing breakdowns, is presented.
International Nuclear Information System (INIS)
Fully nonadiabatic variational computations of energies of the H2+, HD+, and D2+ ions are carried out, making detailed allowance for the motions of all particles. Computations are performed using the basis set of exponential-trigonometric functions dependent on all the interparticle distances. The method of determining the nonlinear adjustable parameters that are responsible for the vibration-modulated exponential decrease of the basis functions with increasing the interparticle distances is developed in detail. Also examined is how an accuracy of the variational energy computations rises with an increase in the number N of the basis functions. For the maximal basis-set size (N = 300), the calculated ground-state energies of the HD+ and D2+ ions deviate from their exact values by 2 x 10-10 and 1 x 10-10 a.u., respectively; whereas the energy of the 2+ ion coincides with 2 its exact value to ten significant digits. The results obtained with the exponential-trigonometric basis set are eight orders of magnitude more precise than the computations with an even larger-sized (N = 350) basis set of monotonically decreasing purely exponential functions of the interatomic distances. This indicates that the exponential-trigonometric basis functions have considerable promise for accurate fully nonadiabatic variational computations of various three-body Coulomb systems with arbitrary masses of their particles. 18 refs., 4 tabs
Indian Academy of Sciences (India)
Sarvesh Kumar Pandey; Prasanta Das; Puspendu K Das; Elangannan Arunan; Sadasivam Manogaran
2015-06-01
It has been shown earlier1 that the relaxed force constants (RFCs) could be used as a measure of bond strength only when the bonds form a part of the complete valence internal coordinates (VIC) basis. However, if the bond is not a part of the complete VIC basis, its RFC is not necessarily a measure of bond strength. Sometimes, it is possible to have a complete VIC basis that does not contain the intramolecular hydrogen bond (IMHB) as part of the basis. This means the RFC of IMHB is not necessarily a measure of bond strength. However, we know that IMHB is a weak bond and hence its RFC has to be a measure of bond strength. We resolve this problem of IMHB not being part of the complete basis by postulating `equivalent’ basis sets where IMHB is part of the basis at least in one of the equivalent sets of VIC. As long as a given IMHB appears in one of the equivalent complete VIC basis sets, its RFC could be used as a measure of bond strength parameter.
The study of basis sets for the calculation of the structure and dynamics of the benzene-Kr complex
International Nuclear Information System (INIS)
An ab initio intermolecular potential energy surface (PES) has been constructed for the benzene-krypton (BKr) van der Waals (vdW) complex. The interaction energy has been calculated at the coupled cluster level of theory with single, double, and perturbatively included triple excitations using different basis sets. As a result, a few analytical PESs of the complex have been determined. They allowed a prediction of the complex structure and its vibrational vdW states. The vibrational energy level pattern exhibits a distinct polyad structure. Comparison of the equilibrium structure, the dipole moment, and vibrational levels of BKr with their experimental counterparts has allowed us to design an optimal basis set composed of a small Dunning’s basis set for the benzene monomer, a larger effective core potential adapted basis set for Kr and additional midbond functions. Such a basis set yields vibrational energy levels that agree very well with the experimental ones as well as with those calculated from the available empirical PES derived from the microwave spectra of the BKr complex. The basis proposed can be applied to larger complexes including Kr because of a reasonable computational cost and accurate results
The study of basis sets for the calculation of the structure and dynamics of the benzene-Kr complex
Energy Technology Data Exchange (ETDEWEB)
Shirkov, Leonid; Makarewicz, Jan, E-mail: jama@amu.edu.pl [Faculty of Chemistry, Adam Mickiewicz University, Umultowska 89b, 61-614 Poznań (Poland)
2015-05-28
An ab initio intermolecular potential energy surface (PES) has been constructed for the benzene-krypton (BKr) van der Waals (vdW) complex. The interaction energy has been calculated at the coupled cluster level of theory with single, double, and perturbatively included triple excitations using different basis sets. As a result, a few analytical PESs of the complex have been determined. They allowed a prediction of the complex structure and its vibrational vdW states. The vibrational energy level pattern exhibits a distinct polyad structure. Comparison of the equilibrium structure, the dipole moment, and vibrational levels of BKr with their experimental counterparts has allowed us to design an optimal basis set composed of a small Dunning’s basis set for the benzene monomer, a larger effective core potential adapted basis set for Kr and additional midbond functions. Such a basis set yields vibrational energy levels that agree very well with the experimental ones as well as with those calculated from the available empirical PES derived from the microwave spectra of the BKr complex. The basis proposed can be applied to larger complexes including Kr because of a reasonable computational cost and accurate results.
Basis set dependence using DFT/B3LYP calculations to model the Raman spectrum of thymine.
Bielecki, Jakub; Lipiec, Ewelina
2016-02-01
Raman spectroscopy (including surface enhanced Raman spectroscopy (SERS) and tip enhanced Raman spectroscopy (TERS)) is a highly promising experimental method for investigations of biomolecule damage induced by ionizing radiation. However, proper interpretation of changes in experimental spectra for complex systems is often difficult or impossible, thus Raman spectra calculations based on density functional theory (DFT) provide an invaluable tool as an additional layer of understanding of underlying processes. There are many works that address the problem of basis set dependence for energy and bond length consideration, nevertheless there is still lack of consistent research on basis set influence on Raman spectra intensities for biomolecules. This study fills this gap by investigating of the influence of basis set choice for the interpretation of Raman spectra of the thymine molecule calculated using the DFT/B3LYP framework and comparing these results with experimental spectra. Among 19 selected Pople's basis sets, the best agreement was achieved using 6-31[Formula: see text](d,p), 6-31[Formula: see text](d,p) and 6-11[Formula: see text]G(d,p) sets. Adding diffuse functions or polarized functions for small basis set or use of a medium or large basis set without diffuse or polarized functions is not sufficient to reproduce Raman intensities correctly. The introduction of the diffuse functions ([Formula: see text]) on hydrogen atoms is not necessary for gas phase calculations. This work serves as a benchmark for further research on the interaction of ionizing radiation with DNA molecules by means of ab initio calculations and Raman spectroscopy. Moreover, this work provides a set of new scaling factors for Raman spectra calculation in the framework of DFT/B3LYP method. PMID:26508426
Turovtsev, V. V.; Orlov, Yu. D.; Tsirulev, A. N.
2015-08-01
The advantages of the orthonormal basis set of 2π-periodic Mathieu functions compared to the trigonometric basis set in calculations of torsional states of molecules are substantiated. Explicit expressions are derived for calculating the Hamiltonian matrix elements of a one-dimensional torsional Schrödinger equation with a periodic potential of the general form in the basis set of Mathieu functions. It is shown that variation of a parameter of Mathieu functions allows the rotation potential and the structural function to be approximated with a good accuracy by a small number of series terms. The conditions for the best choice of this parameter are specified, and approximations are obtained for torsional potentials of n-butane upon rotation about the central C-C bond and of its univalent radical n-butyl C2H5C·H2 upon rotation of the C·H2 group. All algorithms are implemented in the Maple package.
The convergence space of minimal usco mappings
Anguelov, R; Kalenda, O. F. K.
2006-01-01
A convergence structure generalizing the order convergence structure on the set of Hausdorff continuous interval functions is defined on the set of minimal usco maps. The properties of the obtained convergence space are investigated and essential links with the pointwise convergence and the order convergence are revealed. The convergence structure can be extended to a uniform convergence structure so that the convergence space is complete. The important issue of the denseness of the subset of...
Heaps, Charles W
2016-01-01
Quantum molecular dynamics requires an accurate representation of the molecular potential energy surface from a minimal number of electronic structure calculations, particularly for nonadiabatic dynamics where excited states are required. In this paper, we employ pseudospectral sampling of time-dependent Gaussian basis functions for the simulation of non-adiabatic dynamics. Unlike other methods, the pseudospectral Gaussian molecular dynamics tests the Schr\\"{o}dinger equation with $N$ Dirac delta functions located at the centers of the Gaussian functions reducing the scaling of potential energy evaluations from $\\mathcal{O}(N^2)$ to $\\mathcal{O}(N)$. By projecting the Gaussian basis onto discrete points in space, the method is capable of efficiently and quantitatively describing nonadiabatic population transfer and intra-surface quantum coherence. We investigate three model systems; the photodissociation of three coupled Morse oscillators, the bound state dynamics of two coupled Morse oscillators, and a two-d...
Brandenburg, Jan Gerit; Alessio, Maristella; Civalleri, Bartolomeo; Peintinger, Michael F; Bredow, Thomas; Grimme, Stefan
2013-09-26
We extend the previously developed geometrical correction for the inter- and intramolecular basis set superposition error (gCP) to periodic density functional theory (DFT) calculations. We report gCP results compared to those from the standard Boys-Bernardi counterpoise correction scheme and large basis set calculations. The applicability of the method to molecular crystals as the main target is tested for the benchmark set X23. It consists of 23 noncovalently bound crystals as introduced by Johnson et al. (J. Chem. Phys. 2012, 137, 054103) and refined by Tkatchenko et al. (J. Chem. Phys. 2013, 139, 024705). In order to accurately describe long-range electron correlation effects, we use the standard atom-pairwise dispersion correction scheme DFT-D3. We show that a combination of DFT energies with small atom-centered basis sets, the D3 dispersion correction, and the gCP correction can accurately describe van der Waals and hydrogen-bonded crystals. Mean absolute deviations of the X23 sublimation energies can be reduced by more than 70% and 80% for the standard functionals PBE and B3LYP, respectively, to small residual mean absolute deviations of about 2 kcal/mol (corresponding to 13% of the average sublimation energy). As a further test, we compute the interlayer interaction of graphite for varying distances and obtain a good equilibrium distance and interaction energy of 6.75 Å and -43.0 meV/atom at the PBE-D3-gCP/SVP level. We fit the gCP scheme for a recently developed pob-TZVP solid-state basis set and obtain reasonable results for the X23 benchmark set and the potential energy curve for water adsorption on a nickel (110) surface. PMID:23947824
... Conditions Frequently Asked Questions Español Condiciones Chinese Conditions Convergence Insufficiency En Español Read in Chinese What is convergence insufficiency? Convergence insufficiency is the inability to maintain ...
International Nuclear Information System (INIS)
A large-aperture set-up designed for nuclear physics experiments on beams of radioactive nuclei is described. The set-up includes Multiwire Proportional Chamber (MWPC) for measuring the beam profile, MWPC for measuring reaction product angular distributions, a CsI (Tl)-crystal detector and a longitudinal drift ionization chamber for identifying scattered particles and measuring their energy. The results of tests of coordinate MWPCs, particle identification on photon and ion beams, and preliminary measurements of the elastic scattering and the charge exchange reaction of 170 MeV 6He on a CH2 target are presented
Rainwater-Simons-type convergence theorems for generalized convergence methods
Hardtke, Jan-David
2010-01-01
We extend the well-known Rainwater-Simons convergence theorem to various generalized convergence methods such as strong matrix summability, statistical convergence and almost convergence. In fact we prove these theorems not only for boundaries but for the more general notion of (I)-generating sets introduced by Fonf and Lindenstrauss.
Heaps, Charles W.; Mazziotti, David A.
2016-08-01
Quantum molecular dynamics requires an accurate representation of the molecular potential energy surface from a minimal number of electronic structure calculations, particularly for nonadiabatic dynamics where excited states are required. In this paper, we employ pseudospectral sampling of time-dependent Gaussian basis functions for the simulation of non-adiabatic dynamics. Unlike other methods, the pseudospectral Gaussian molecular dynamics tests the Schrödinger equation with N Dirac delta functions located at the centers of the Gaussian functions reducing the scaling of potential energy evaluations from O ( N 2 ) to O ( N ) . By projecting the Gaussian basis onto discrete points in space, the method is capable of efficiently and quantitatively describing the nonadiabatic population transfer and intra-surface quantum coherence. We investigate three model systems: the photodissociation of three coupled Morse oscillators, the bound state dynamics of two coupled Morse oscillators, and a two-dimensional model for collinear triatomic vibrational dynamics. In all cases, the pseudospectral Gaussian method is in quantitative agreement with numerically exact calculations. The results are promising for nonadiabatic molecular dynamics in molecular systems where strongly correlated ground or excited states require expensive electronic structure calculations.
International Nuclear Information System (INIS)
We present the self-consistent numerical-basis-set linear combination of atomic orbitals (LCAO) discrete variational method for treating the electronic structure of thin films. As in the case of bulk solids, this method provides for thin films accurate solutions of the one-particle local density equations with a non-muffin-tin potential. Hamiltonian and overlap matrix elements are evaluated accurately by means of a three-dimensional numerical Diophantine integration scheme. Application of this method is made to the self-consistent solution of one-, three-, and five-layer Ni(001) unsupported films. The LCAO Bloch basis set consists of valence orbitals (3d, 4s, and 4p states for transition metals) orthogonalized to the frozen-core wave functions. The self-consistent potential is obtained iteratively within the superposition of overlapping spherical atomic charge density model with the atomic configurations treated as adjustable parameters. Thus the crystal Coulomb potential is constructed as a superposition of overlapping spherically symmetric atomic potentials and, correspondingly, the local density Kohn-Sham (α = 2/3) potential is determined from a superposition of atomic charge densities. At each iteration in the self-consistency procedure, the crystal charge density is evaluated using a sampling of 15 independent k points in (1/8)th of the irreducible two-dimensional Brillouin zone. The total density of states (DOS) and projected local DOS (by layer plane) are calculated using an analytic linear energy triangle method (presented as an Appendix) generalized from the tetrahedron scheme for bulk systems. Distinct differences are obtained between the surface and central plane local DOS. The central plane DOS is found to converge rapidly to the DOS of bulk paramagnetic Ni obtained by Wang and Callaway. Only a very small surplus charge (0.03 electron/atom) is found on the surface planes, in agreement with jellium model calculations
A stochastic convergence analysis for Tikhonov regularization with sparsity constraints
International Nuclear Information System (INIS)
In this paper we investigate convergence properties of Tikhonov regularization for linear ill-posed problems under a stochastic error model. Namely, we assume that we are given a finite amount of measurements, each contaminated by Gaussian noise with zero mean and known finite variance. Using Besov-space penalty terms to promote sparse solutions with respect to a preassigned wavelet basis, the Ky Fan metric allows us to lift deterministic convergence results into the stochastic setting. In particular, we formulate a general convergence theorem and propose a formula to directly calculate a suitable regularization parameter. This immediately leads to convergence rates. Numerical examples are presented to verify the theoretical results. (paper)
Topological Dichotomy and Unconditional Convergence
Lefevre, Pascal
1999-01-01
In this paper, we give a criterion for unconditional convergence with respect to some summability methods, dealing with the topological size of the set of choices of sign providing convergence. We obtain similar results for boundedness. In particular, quasi-sure unconditional convergence implies unconditional convergence.
International Nuclear Information System (INIS)
Prolapse-free basis sets suitable for four-component relativistic quantum chemical calculations are presented for the superheavy elements up to 118Uuo (104Rf, 105Db, 106Sg, 107Bh, 108Hs, 109Mt, 110Ds, 111Rg, 112Uub, 113Uut, 114Uuq, 115Uup, 116Uuh, 117Uus, 118Uuo) and 103Lr. These basis sets were optimized by minimizing the absolute values of the energy difference between the Dirac-Fock-Roothaan total energy and the corresponding numerical value at a milli-Hartree order of magnitude, resulting in a good balance between cost and accuracy. Parameters for generating exponents and new numerical data for some superheavy elements are also presented
Buimaga-Iarinca, Luiza
2014-01-01
We compare the density functional theory (DFT) results on the adsorption of small aromatic molecules (benzene, pyridine and thiophene) on gold surfaces obtained by using three types of van der Waals exchange-correlation functionals and localized basis set calculations. We show that the value of the molecule surface binding energy depends on the interplay between the BSSE effect and the tendency of the exchange-correlation functionals to overestimate both the molecule-surface as well as the gold-gold distances within the relaxed systems. Consequently, we find that by using different types of LCAO basis sets or geometric models for the adsorption of the molecules on the surface, the binding energy can vary up to 100 %. A critical analysis of the physical parameters resulting from the calculations is presented for each exchange-correlation functional.
Salvador Sedano, Pedro; Mayer, István
2004-01-01
The basis set superposition error-free second-order MØller-Plesset perturbation theory of intermolecular interactions was studied. The difficulties of the counterpoise (CP) correction in open-shell systems were also discussed. The calculations were performed by a program which was used for testing the new variants of the theory. It was shown that the CP correction for the diabatic surfaces should be preferred to the adiabatic ones
Czech Academy of Sciences Publication Activity Database
Ran, Jiong; Hobza, Pavel
2009-01-01
Roč. 5, č. 4 (2009), s. 1180-1185. ISSN 1549-9618 R&D Projects: GA MŠk LC512; GA AV ČR IAA400550510 Institutional research plan: CEZ:AV0Z40550506 Keywords : lone pair...pi-electron complexes * nature of bonding * CCSD(T)/ complete basis set limit calculations Subject RIV: CF - Physical ; Theoretical Chemistry Impact factor: 4.804, year: 2009
de Silva, Piotr; Wesolowski, Tomasz A.
2012-03-01
Within the linear combination of atomic orbitals (LCAO) approximation, one can distinguish two different Kohn-Sham potentials. One is the potential available numerically in calculations, and the other is the exact potential corresponding to the LCAO density. The latter is usually not available, but can be obtained from the total density by a numerical inversion procedure or, as is done here, analytically using only one LCAO Kohn-Sham orbital. In the complete basis-set limit, the lowest-lying Kohn-Sham orbital suffices to perform the analytical inversion, and the two potentials differ by no more than a constant. The relation between these two potentials is investigated here for diatomic molecules and several atomic basis sets of increasing size and quality. The differences between the two potentials are usually qualitative (wrong behavior at nuclear cusps and far from the molecule even if Slater-type orbitals are used) and δ-like features at nodal planes of the lowest-lying LCAO Kohn-Sham orbital. Such nodes occur frequently in LCAO calculations and are not physical. Whereas the behavior of the potential can be systematically improved locally by the increase of the basis sets, the occurrence of nodes is not correlated with the size of the basis set. The presence of nodes in the lowest-lying LCAO orbital can be used to monitor whether the effective potential in LCAO Kohn-Sham equations can be interpreted as the potential needed for pure-state noninteracting v-representability of the LCAO density. Squares of such node-containing lowest-lying LCAO Kohn-Sham orbitals are nontrivial examples of two-electron densities which are not pure-state noninteracting v-representable.
Buimaga-Iarinca, Luiza; Morari, Cristian
2014-01-01
We compare the density functional theory (DFT) results on the adsorption of small aromatic molecules (benzene, pyridine and thiophene) on gold surfaces obtained by using three types of van der Waals exchange-correlation functionals and localized basis set calculations. We show that the value of the molecule surface binding energy depends on the interplay between the BSSE effect and the tendency of the exchange-correlation functionals to overestimate both the molecule-surface as well as the go...
De Silva, Piotr; Wesolowski, Tomasz Adam
2012-01-01
Within the linear combination of atomic orbitals (LCAO) approximation, one can distinguish two different Kohn-Sham potentials. One is the potential available numerically in calculations, and the other is the exact potential corresponding to the LCAO density. The latter is usually not available, but can be obtained from the total density by a numerical inversion procedure or, as is done here, analytically using only one LCAO Kohn-Sham orbital. In the complete basis-set limit, the lowest-lying ...
Kolmann, Stephen J.; Jordan, Meredith J. T.
2010-02-01
One of the largest remaining errors in thermochemical calculations is the determination of the zero-point energy (ZPE). The fully coupled, anharmonic ZPE and ground state nuclear wave function of the SSSH radical are calculated using quantum diffusion Monte Carlo on interpolated potential energy surfaces (PESs) constructed using a variety of method and basis set combinations. The ZPE of SSSH, which is approximately 29 kJ mol-1 at the CCSD(T)/6-31G∗ level of theory, has a 4 kJ mol-1 dependence on the treatment of electron correlation. The anharmonic ZPEs are consistently 0.3 kJ mol-1 lower in energy than the harmonic ZPEs calculated at the Hartree-Fock and MP2 levels of theory, and 0.7 kJ mol-1 lower in energy at the CCSD(T)/6-31G∗ level of theory. Ideally, for sub-kJ mol-1 thermochemical accuracy, ZPEs should be calculated using correlated methods with as big a basis set as practicable. The ground state nuclear wave function of SSSH also has significant method and basis set dependence. The analysis of the nuclear wave function indicates that SSSH is localized to a single symmetry equivalent global minimum, despite having sufficient ZPE to be delocalized over both minima. As part of this work, modifications to the interpolated PES construction scheme of Collins and co-workers are presented.
Richard, Ryan M; Marshall, Michael S; Dolgounitcheva, O; Ortiz, J V; Brédas, Jean-Luc; Marom, Noa; Sherrill, C David
2016-02-01
In designing organic materials for electronics applications, particularly for organic photovoltaics (OPV), the ionization potential (IP) of the donor and the electron affinity (EA) of the acceptor play key roles. This makes OPV design an appealing application for computational chemistry since IPs and EAs are readily calculable from most electronic structure methods. Unfortunately reliable, high-accuracy wave function methods, such as coupled cluster theory with single, double, and perturbative triples [CCSD(T)] in the complete basis set (CBS) limit are too expensive for routine applications to this problem for any but the smallest of systems. One solution is to calibrate approximate, less computationally expensive methods against a database of high-accuracy IP/EA values; however, to our knowledge, no such database exists for systems related to OPV design. The present work is the first of a multipart study whose overarching goal is to determine which computational methods can be used to reliably compute IPs and EAs of electron acceptors. This part introduces a database of 24 known organic electron acceptors and provides high-accuracy vertical IP and EA values expected to be within ±0.03 eV of the true non-relativistic, vertical CCSD(T)/CBS limit. Convergence of IP and EA values toward the CBS limit is studied systematically for the Hartree-Fock, MP2 correlation, and beyond-MP2 coupled cluster contributions to the focal point estimates. PMID:26731487
Lee, Hee-Seung; Tuckerman, Mark E
2008-12-14
An efficient computational approach to perform Car-Parrinello ab initio molecular dynamics (CPAIMD) simulations under cluster (free) boundary conditions is presented. The general approach builds upon a recent real-space CPAIMD formalism using discrete variable representation (DVR) basis sets [Y. Liu et al., Phys. Rev. B 12, 125110 (2003); H.-S. Lee and M. E. Tuckerman, J. Phys. Chem. A 110, 5549 (2006)]. In order to satisfy cluster boundary conditions, a DVR based on sinc functions is utilized to expand the Kohn-Sham orbitals and electron density. Poisson's equation is solved in order to calculate the Hartree potential via an integral representation of the 1/r singularity. Excellent convergence properties are achieved with respect to the number of grid points (or DVR functions) and the size of the simulation cell. A straightforward implementation of the present approach leads to near linear scaling [O(N(4/3))] of the computational cost with respect to the system size (N) for the solution of Poisson's equation. The accuracy and stability of CPAIMD simulations based on sinc DVR are tested for a model problem as well as for N(2) and a water dimer. PMID:19071908
Convergence of Optimization Problems
Directory of Open Access Journals (Sweden)
K. Jeyalakshmi
2012-03-01
Full Text Available In this paper we consider a general optimization problem (OP and study the convergence and approximation of optimal values and optimal solutions to changes in the cost function and the set of feasible solutions. We consider the convergence optimization problems under the familiar notion of uniform convergence. We do not assume the convexity of the functions involved. Instead we consider a class of functions whose directional derivatives are convex. They are known as locally convex functions or following Craven and Mond nearly convex functions. We given necessary preliminaries and we prove that a sequence of locally convex optimization problems converge to a locally convex problem. We also prove that uniform convergence of locally convex optimization problems implies epi-graph convergence of the problems. Even though for simplicity we have taken locally convex functions, the results given here can be proved for locally Lipchitz functions also.
The utility of many-body decompositions for the accurate basis set extrapolation of ab initio data
International Nuclear Information System (INIS)
We present a powerful new technique for the extrapolation of ab initio data based on many-body decompositions. Using the new methodology and subtle modifications of the standard correlation consistent basis sets, the H+H2 barrier height is estimated at 9.603 kcal/mol with a precision of about 0.003 kcal/mol; this extremely accurate result is all the more striking as it can be obtained using basis sets no larger than aug-cc-pVQZ. The method is also used to yield highly accurate energies for the H+H2 system on a grid of points previously calculated by quantum Monte Carlo. The three-body energy, summed with exact one- and two-body energies, is observed to yield a useful approximate lower bound for the total energy. The highly accurate energies afforded by this method can also be used to assess the accuracy of previously calculated data that has been used to construct potential energy surfaces. As an example, we make a detailed comparison between the new results and the quantum Monte Carlo results for H+H2. copyright 1999 American Institute of Physics
Directory of Open Access Journals (Sweden)
I. Riipinen
2014-11-01
Full Text Available We present a theoretical study investigating the cloud condensation nucleus (CCN activation of multicomponent organic mixtures. We modeled these complex mixtures using the solubility basis set (SBS, analogous to the volatility basis set VBS, describing the mixture as a set of surrogate compounds with varying water-solubilities in a given range. We conducted Köhler theory calculations for 144 different mixtures with varying solubility range, number of components, assumption about the organic mixture thermodynamics and the shape of the solubility distribution, yielding approximately 6000 unique CCN-activation points. The results from these comprehensive calculations were compared to three simplifying assumptions about organic aerosol solubility: (1 complete dissolution at the point of activation, (2 combining the aerosol solubility with the molar mass and density into a single hygroscopicity parameter κ, (3 assuming a fixed water-soluble fraction ϵeff. While the complete dissolution was able to reproduce the activation points with a reasonable accuracy only when the majority (70–80% of the material was dissolved at the point of activation, the single parameter representations of complex mixture solubility were confirmed to be powerful semi-empirical tools for representing the CCN activation of organic aerosol. Depending on the condensed-phase interactions between the organic molecules, material with solubilities larger than about 1–10 g L−1 could be treated as completely soluble in the CCN activation process over particle dry diameters between 20 and 500 nm and supersaturations between 0.03 and 8%. Our results indicate that understanding the details of the solubility distribution in the range of 0.1 to 100 g L−1 is critical for capturing the CCN activation, while resolution outside this solubility range will probably not add much information except in some special cases. The connection of these results to the previous
Directory of Open Access Journals (Sweden)
Roghieh Tarlani Bashiz
2015-12-01
Full Text Available Density functional theory calculations (DFT, as well as hybrid methods (B3LYP and HF method for CNT-Calixarene complexes have been carried out to study structural stability. The geometry of the Calixarene has been optimized at DFT methods such as M062x B3LYP and HF methods with 6-31G, 6-31G*and 6-31G** basis sets. According to GIAO method, NMR parameters have been evaluated. The Gaussian quantum chemical package is used for all calculations. The gauge including atomic orbital (GIAO approach was applied for chemical shielding calculations for an isolated calix aren and a complex of Calix-SWCNTs.
Andrade, Xavier; Aspuru-Guzik, Alán
2013-10-01
We discuss the application of graphical processing units (GPUs) to accelerate real-space density functional theory (DFT) calculations. To make our implementation efficient, we have developed a scheme to expose the data parallelism available in the DFT approach; this is applied to the different procedures required for a real-space DFT calculation. We present results for current-generation GPUs from AMD and Nvidia, which show that our scheme, implemented in the free code Octopus, can reach a sustained performance of up to 90 GFlops for a single GPU, representing a significant speed-up when compared to the CPU version of the code. Moreover, for some systems, our implementation can outperform a GPU Gaussian basis set code, showing that the real-space approach is a competitive alternative for DFT simulations on GPUs. PMID:26589153
Andrade, Xavier
2013-01-01
We discuss the application of graphical processing units (GPUs) to accelerate real-space density functional theory (DFT) calculations. To make our implementation efficient, we have developed a scheme to expose the data parallelism available in the DFT approach; this is applied to the different procedures required for a real-space DFT calculation. We present results for current-generation GPUs from AMD and Nvidia, which show that our scheme, implemented in the free code OCTOPUS, can reach a sustained performance of up to 90 GFlops for a single GPU, representing an important speed-up when compared to the CPU version of the code. Moreover, for some systems our implementation can outperform a GPU Gaussian basis set code, showing that the real-space approach is a competitive alternative for DFT simulations on GPUs.
The issue of basis set size in e- + H(1s → 2s, 2p) collisions
International Nuclear Information System (INIS)
The results of the multichannel eikonal theory and three recent close-coupling calculations are compared for the σ1s→2s and σ1s→2p integral cross sections at intermediate energies. Also, the predictions of the three-state close-coupling theory and the present multichannel eikonal theory results for the moduli and phase angles of the complex scattering amplitudes for the 2p0 and 2p+1 substates of H at 54.40 eV are compared. Good agreement is obtained between previous results and the present results for the phase angles of the scattering amplitudes of both substates at small scattering angles. However, differences are noted in the modulus of the scattering amplitude of the 2p+1 substate at small scattering angles between the two theories. The implications this difference has on the composition and size of the basis sets used are addressed. (author)
Cignetti, Fabien; Salvia, Emilie; Anton, Jean-Luc; Grosbras, Marie-Hélène; Assaiante, Christine
2016-01-01
Conventional analysis of functional magnetic resonance imaging (fMRI) data using the general linear model (GLM) employs a neural model convolved with a canonical hemodynamic response function (HRF) peaking 5 s after stimulation. Incorporation of a further basis function, namely the canonical HRF temporal derivative, accounts for delays in the hemodynamic response to neural activity. A population that may benefit from this flexible approach is children whose hemodynamic response is not yet mature. Here, we examined the effects of using the set based on the canonical HRF plus its temporal derivative on both first- and second-level GLM analyses, through simulations and using developmental data (an fMRI dataset on proprioceptive mapping in children and adults). Simulations of delayed fMRI first-level data emphasized the benefit of carrying forward to the second-level a derivative boost that combines derivative and nonderivative beta estimates. In the experimental data, second-level analysis using a paired t-test showed increased mean amplitude estimate (i.e., increased group contrast mean) in several brain regions related to proprioceptive processing when using the derivative boost compared to using only the nonderivative term. This was true especially in children. However, carrying forward to the second-level the individual derivative boosts had adverse consequences on random-effects analysis that implemented one-sample t-test, yielding increased between-subject variance, thus affecting group-level statistic. Boosted data also presented a lower level of smoothness that had implication for the detection of group average activation. Imposing soft constraints on the derivative boost by limiting the time-to-peak range of the modeled response within a specified range (i.e., 4–6 s) mitigated these issues. These findings support the notion that there are pros and cons to using the informed basis set with developmental data. PMID:27471441
Weak Convergence and Weak Convergence
Narita Keiko; Shidama Yasunari; Endou Noboru
2015-01-01
In this article, we deal with weak convergence on sequences in real normed spaces, and weak* convergence on sequences in dual spaces of real normed spaces. In the first section, we proved some topological properties of dual spaces of real normed spaces. We used these theorems for proofs of Section 3. In Section 2, we defined weak convergence and weak* convergence, and proved some properties. By RNS_Real Mizar functor, real normed spaces as real number spaces already defined in the article [18...
Convergence of coupled cluster perturbation theory
Eriksen, Janus Juul; Matthews, Devin A; Jørgensen, Poul; Olsen, Jeppe
2016-01-01
The convergence of a recently proposed coupled cluster (CC) family of perturbation series [Eriksen et al., J. Chem. Phys. 140, 064108 (2014)], in which the energetic difference between a parent and a target CC model is expanded in orders of the M{\\o}ller-Plesset (MP) fluctuation potential, is investigated for four prototypical closed-shell systems (Ne, singlet methylene, distorted HF, and the fluoride anion) in standard and augmented basis sets. In these investigations, energy corrections of the various series have been calculated to high orders and their convergence radii determined by probing for possible front- and back-door intruder states. In summary, we conclude how it is primarily the choice of target state, and not the choice of parent state, which ultimately governs the convergence behavior of a given series. For example, restricting the target state to, say, triple or quadruple excitations might remove intruders present in series that target the full configuration interaction (FCI) limit, such as th...
On Convergence Properties of Shannon Entropy
Piera, Francisco J.; Parada, Patricio
2007-01-01
Convergence properties of Shannon Entropy are studied. In the differential setting, it is shown that weak convergence of probability measures, or convergence in distribution, is not enough for convergence of the associated differential entropies. A general result for the desired differential entropy convergence is provided, taking into account both compactly and uncompactly supported densities. Convergence of differential entropy is also characterized in terms of the Kullback-Liebler discrimi...
Directory of Open Access Journals (Sweden)
Claudio Amovilli
2016-02-01
Full Text Available In this work, we present a simple decomposition scheme of the Kohn-Sham optimized orbitals which is able to provide a reduced basis set, made of localized polycentric orbitals, specifically designed for Quantum Monte Carlo. The decomposition follows a standard Density functional theory (DFT calculation and is based on atomic connectivity and shell structure. The new orbitals are used to construct a compact correlated wave function of the Slater–Jastrow form which is optimized at the Variational Monte Carlo level and then used as the trial wave function for a final Diffusion Monte Carlo accurate energy calculation. We are able, in this way, to capture the basic information on the real system brought by the Kohn-Sham orbitals and use it for the calculation of the ground state energy within a strictly variational method. Here, we show test calculations performed on some small selected systems to assess the validity of the proposed approach in a molecular fragmentation, in the calculation of a barrier height of a chemical reaction and in the determination of intermolecular potentials. The final Diffusion Monte Carlo energies are in very good agreement with the best literature data within chemical accuracy.
Mohr, Stephan; Genovese, Luigi; Ratcliff, Laura; Masella, Michel
The quantum mechanics/molecular mechanis (QM/MM) method is a popular approach that allows to perform atomistic simulations using different levels of accuracy. Since only the essential part of the simulation domain is treated using a highly precise (but also expensive) QM method, whereas the remaining parts are handled using a less accurate level of theory, this approach allows to considerably extend the total system size that can be simulated without a notable loss of accuracy. In order to couple the QM and MM regions we use an approximation of the electrostatic potential based on a multipole expansion. The multipoles of the QM region are determined based on the results of a linear scaling Density Functional Theory (DFT) calculation using a set of adaptive, localized basis functions, as implemented within the BigDFT software package. As this determination comes at virtually no extra cost compared to the QM calculation, the coupling between QM and MM region can be done very efficiently. In this presentation I will demonstrate the accuracy of both the linear scaling DFT approach itself as well as of the approximation of the electrostatic potential based on the multipole expansion, and show some first QM/MM applications using the aforementioned approach.
Yang, Yanchao; Jiang, Hong; Liu, Congbin; Lan, Zhongli
2013-03-01
Cognitive radio (CR) is an intelligent wireless communication system which can dynamically adjust the parameters to improve system performance depending on the environmental change and quality of service. The core technology for CR is the design of cognitive engine, which introduces reasoning and learning methods in the field of artificial intelligence, to achieve the perception, adaptation and learning capability. Considering the dynamical wireless environment and demands, this paper proposes a design of cognitive engine based on the rough sets (RS) and radial basis function neural network (RBF_NN). The method uses experienced knowledge and environment information processed by RS module to train the RBF_NN, and then the learning model is used to reconfigure communication parameters to allocate resources rationally and improve system performance. After training learning model, the performance is evaluated according to two benchmark functions. The simulation results demonstrate the effectiveness of the model and the proposed cognitive engine can effectively achieve the goal of learning and reconfiguration in cognitive radio.
Teodoro, Tiago Quevedo; da Silva, Albérico Borges Ferreira; Haiduke, Roberto Luiz Andrade
2014-09-01
This study reports a new relativistic prolapse-free Gaussian basis set series of quadruple-ζ quality, RPF-4Z, and an augmented version that includes extra diffuse functions, aug-RPF-4Z, for all the s- and p-block elements. The relativistic adapted Gaussian basis sets (RAGBSs), which are free of variational prolapse, were used as the starting primitive sets. Exponents of correlating/polarization functions were taken from a polynomial version of the generator coordinate Dirac-Fock (p-GCDF) method, in which the previously optimized RAGBS parameters are applied. By using such procedure we aimed to reduce the computational demand of these sets in comparison with fully optimized ones. The effect of these basis set increments on the correlation energy was evaluated by atomic multireference configuration interaction calculations with single and double excitations out of the valence shell. Finally, atomic and molecular calculations of fundamental properties (bond lengths, vibrational frequencies, dipole moments and electron affinities) corroborate the quadruple-ζ quality of these new sets that are also about half-time-consuming than the correspondent Dyall's v4z sets. The read-to use format of these (aug-)RPF-4Fz sets are available as Supporting Information files and can also be found at http://basis-sets.iqsc.usp.br/ . PMID:26588525
Aurel Iancu
2008-01-01
After presenting the institutional construction during the pre-accession and post-accession to the Economic and Monetary Union (EMU), the exchange rate mechanisms (ERM) in several countries and the convergence criteria, we go on with a brief analysis of the way the CEE countries cope with the convergence criteria in accordance with the Maastricht Treaty. Then, the study deals with a topic often discussed in the scientific literature and included on the agenda of decision-makers at various lev...
DEFF Research Database (Denmark)
Kupka, Teobald; Stachów, Michal; Kaminsky, Jakub; Sauer, Stephan P. A.
2013-01-01
A linear correlation between isotropic nuclear magnetic shielding constants for seven model molecules (CH2O, H2O, HF, F2, HCN, SiH4 and H2S) calculated with 37 methods (34 density functionals, RHF, MP2 and CCSD(T) ), with affordable pcS-2 basis set and corresponding complete basis set results......, estimated from calculations with the family of polarizationconsistent pcS-n basis sets is reported. This dependence was also supported by inspection of profiles of deviation between CBS estimated nuclear shieldings and obtained with significantly smaller basis sets pcS-2 and aug-cc-pVTZ-J for the selected...... set of 37 calculation methods. It was possible to formulate a practical approach of estimating the values of isotropic nuclear magnetic shielding constants at the CCSD(T)/CBS and MP2/CBS levels from affordable CCSD(T)/pcS-2, MP2/pcS-2 and DFT/CBS calculations with pcS-n basis sets. The proposed method...
2010-10-01
... physician services in a teaching setting. 415.170 Section 415.170 Public Health CENTERS FOR MEDICARE... BY PHYSICIANS IN PROVIDERS, SUPERVISING PHYSICIANS IN TEACHING SETTINGS, AND RESIDENTS IN CERTAIN SETTINGS Physician Services in Teaching Settings § 415.170 Conditions for payment on a fee schedule...
Martin, J M L; Martin, Jan M.L.; Sundermann, Andreas
2001-01-01
We propose large-core correlation-consistent pseudopotential basis sets for the heavy p-block elements Ga-Kr and In-Xe. The basis sets are of cc-pVTZ and cc-pVQZ quality, and have been optimized for use with the large-core (valence-electrons only) Stuttgart-Dresden-Bonn relativistic pseudopotentials. Validation calculations on a variety of third-row and fourth-row diatomics suggest them to be comparable in quality to the all-electron cc-pVTZ and cc-pVQZ basis sets for lighter elements. Especially the SDB-cc-pVQZ basis set in conjunction with a core polarization potential (CPP) yields excellent agreement with experiment for compounds of the later heavy p-block elements. For accurate calculations on Ga (and, to a lesser extent, Ge) compounds, explicit treatment of 13 valence electrons appears to be desirable, while it seems inevitable for In compounds. For Ga and Ge, we propose correlation consistent basis sets extended for (3d) correlation. For accurate calculations on organometallic complexes of interest to h...
Baranowska-Łączkowska, Angelika; Fernández, Berta
2014-01-30
In order to obtain efficient basis sets for the evaluation of van der Waals complex intermolecular potentials, we carry out systematic basis set studies. For this, interaction energies at representative geometries on the potential energy surfaces are evaluated using the CCSD(T) correlation method and large polarized LPol-n and augmented polarization-consistent aug-pc-2 basis sets extended with different sets of midbond functions. On the basis of the root mean square errors calculated with respect to the values for the most accurate potentials available, basis sets are selected for fitting the corresponding interaction energies and getting analytical potentials. In this work, we study the Ne-N2 van der Waals complex and after the above procedure, the aug-pc-2-3321 and the LPol-ds-33221 basis set results are fitted. The obtained potentials are characterized by T-shaped global minima at distances between the Ne atom and the N2 center of mass of 3.39 Å, with interaction energies of -49.36 cm(-1) for the aug-pc-2-3321 surface and -50.28 cm(-1) for the LPol-ds-33221 surface. Both sets of results are in excellent agreement with the reference surface. To check the potentials further microwave transition frequencies are calculated that agree well with the experimental and the aV5Z-33221 values. The success of this study suggests that it is feasible to carry out similar accurate calculations of interaction energies and ro-vibrational spectra at reduced cost for larger complexes than has been possible hitherto. PMID:24375320
Directory of Open Access Journals (Sweden)
B. N. Murphy
2012-04-01
Full Text Available Multigenerational oxidation chemistry of atmospheric organic compounds and its effects on aerosol loadings and chemical composition is investigated by implementing the Two-Dimensional Volatility Basis Set (2-D-VBS in a Lagrangian host chemical transport model. Three model formulations were chosen to explore the complex interactions between functionalization and fragmentation processes during gas-phase oxidation of organic compounds by the hydroxyl radical. The base case model employs a conservative transformation by assuming a reduction of one order of magnitude in effective saturation concentration and an increase of oxygen content by one or two oxygen atoms per oxidation generation. A second scheme simulates functionalization in more detail using group contribution theory to estimate the effects of oxygen addition to the carbon backbone on the compound volatility. Finally, a fragmentation scheme is added to the detailed functionalization scheme to create a functionalization-fragmentation parameterization. Two condensed-phase chemistry pathways are also implemented as additional sensitivity tests to simulate (1 heterogeneous oxidation via OH uptake to the particle-phase and (2 aqueous-phase chemistry of glyoxal and methylglyoxal. The model is applied to summer and winter periods at three sites where observations of organic aerosol (OA mass and O:C were obtained during the European Integrated Project on Aerosol Cloud Climate and Air Quality Interactions (EUCAARI campaigns. The base case model reproduces observed mass concentrations and O:C well, with fractional errors (FE lower than 55% and 25%, respectively. The detailed functionalization scheme tends to overpredict OA concentrations, especially in the summertime, and also underpredicts O:C by approximately a factor of 2. The detailed functionalization model with fragmentation agrees well with the observations for OA concentration, but still underpredicts O:C. Both heterogeneous oxidation and
Directory of Open Access Journals (Sweden)
B. N. Murphy
2012-11-01
Full Text Available Multigenerational oxidation chemistry of atmospheric organic compounds and its effects on aerosol loadings and chemical composition is investigated by implementing the Two-Dimensional Volatility Basis Set (2-D-VBS in a Lagrangian host chemical transport model. Three model formulations were chosen to explore the complex interactions between functionalization and fragmentation processes during gas-phase oxidation of organic compounds by the hydroxyl radical. The base case model employs a conservative transformation by assuming a reduction of one order of magnitude in effective saturation concentration and an increase of oxygen content by one or two oxygen atoms per oxidation generation. A second scheme simulates functionalization in more detail using group contribution theory to estimate the effects of oxygen addition to the carbon backbone on the compound volatility. Finally, a fragmentation scheme is added to the detailed functionalization scheme to create a functionalization-fragmentation parameterization. Two condensed-phase chemistry pathways are also implemented as additional sensitivity tests to simulate (1 heterogeneous oxidation via OH uptake to the particle-phase and (2 aqueous-phase chemistry of glyoxal and methylglyoxal. The model is applied to summer and winter periods at three sites where observations of organic aerosol (OA mass and O:C were obtained during the European Integrated Project on Aerosol Cloud Climate and Air Quality Interactions (EUCAARI campaigns. The base case model reproduces observed mass concentrations and O:C well, with fractional errors (FE lower than 55% and 25%, respectively. The detailed functionalization scheme tends to overpredict OA concentrations, especially in the summertime, and also underpredicts O:C by approximately a factor of 2. The detailed functionalization model with fragmentation agrees well with the observations for OA concentration, but still underpredicts O:C. Both heterogeneous oxidation and
International Nuclear Information System (INIS)
Monte Carlo simulations of nuclear criticality eigenvalue problems are often performed by general purpose radiation transport codes such as MCNP. MCNP performs detailed statistical analysis of the criticality calculation and provides feedback to the user with warning messages, tables, and graphs. The purpose of the analysis is to provide the user with sufficient information to assess spatial convergence of the eigenfunction and thus the validity of the criticality calculation. As a test of this statistical analysis package in MCNP, analytic criticality verification benchmark problems have been used for the first time to assess the performance of the criticality convergence tests in MCNP. The MCNP statistical analysis capability has been recently assessed using the 75 multigroup criticality verification analytic problem test set. MCNP was verified with these problems at the 10-4 to 10-5 statistical error level using 40 000 histories per cycle and 2000 active cycles. In all cases, the final boxed combined keff answer was given with the standard deviation and three confidence intervals that contained the analytic keff. To test the effectiveness of the statistical analysis checks in identifying poor eigenfunction convergence, ten problems from the test set were deliberately run incorrectly using 1000 histories per cycle, 200 active cycles, and 10 inactive cycles. Six problems with large dominance ratios were chosen from the test set because they do not achieve the normal spatial mode in the beginning of the calculation. To further stress the convergence tests, these problems were also started with an initial fission source point 1 cm from the boundary thus increasing the likelihood of a poorly converged initial fission source distribution. The final combined keff confidence intervals for these deliberately ill-posed problems did not include the analytic keff value. In no case did a bad confidence interval go undetected. Warning messages were given signaling that the
How Riemannian Manifolds Converge: A Survey
Sormani, Christina
2010-01-01
This is an intuitive survey of extrinsic and intrinsic notions of convergence of manifolds complete with pictures of key examples and a discussion of the properties associated with each notion. We begin with a description of three extrinsic notions which have been applied to study sequences of submanifolds in Euclidean space: Hausdorff convergence of sets, flat convergence of integral currents, and weak convergence of varifolds. We next describe a variety of intrinsic notions of convergence w...
International Nuclear Information System (INIS)
The Griffin-Hill-Wheeler-Hartree-Fock equations are applied to generate universal gaussian and Slater-type basis sets for the title atoms. The results for the energies compare competitively with atom-optimized calculations at a very low computational cost. (author)
International Nuclear Information System (INIS)
Accurate Gaussian basis sets (18s for Li and Be and 20s11p for the atoms from B to Ne) for the first-row atoms, generated with an improved generator coordinate Hartree-Fock method, were contracted and enriched with polarization functions. These basis sets were tested for B2, C2, BeO, Cn-, LiF, N2, CO, BF, NO+, O2, and F2. At the Hartree-Fock (HP), second-order Moeller-Plesset (MP2), fourth-order Moeller-Plesset (MP4), and density functional theory (DFT) levels, the dipole moments, bond lengths, and harmonic vibrational frequencies were studied, and at the MP2, MP4, and DFT levels, the dissociation energies were evaluated and compared with the corresponding experimental values and with values obtained using other contracted Gaussian basis sets and numerical HF calculations. For all diatomic molecules studied, the differences between the total energies, obtained with the largest contracted basis set [6s5p3d1f], and those calculated with the numerical HF methods were always less than 3.2 mhartree
A. Brouwer; M. Hoogendoorn; E. Naarding
2015-01-01
In this paper we evaluate the International Accounting Standards Board’s (IASB) efforts, in a discussion paper (DP) of 2013, to develop a new conceptual framework (CF) in the light of its stated ambition to establish a robust and consistent basis for future standard setting, thereby guiding standard
DEFF Research Database (Denmark)
Brandbyge, Mads
2014-01-01
In a recent paper Reuter and Harrison [J. Chem. Phys.139, 114104 (2013)] question the widely used mean-field electron transport theories, which employ nonorthogonal localized basis sets. They claim these can violate an “implicit decoupling assumption,” leading to wrong results for the current...
Directory of Open Access Journals (Sweden)
Abdul Hameed Q. A. Al-Tai
2011-01-01
Full Text Available The aim of this paper is to introduce and study the fuzzy neighborhood, the limit fuzzy number, the convergent fuzzy sequence, the bounded fuzzy sequence, and the Cauchy fuzzy sequence on the base which is adopted by Abdul Hameed (every real number r is replaced by a fuzzy number r¯ (either triangular fuzzy number or singleton fuzzy set (fuzzy point. And then, we will consider that some results respect effect of the upper sequence on the convergent fuzzy sequence, the bounded fuzzy sequence, and the Cauchy fuzzy sequence.
Al-Tai, Abdul Hameed Q. A.
2011-01-01
The aim of this paper is to introduce and study the fuzzy neighborhood, the limit fuzzy number, the convergent fuzzy sequence, the bounded fuzzy sequence, and the Cauchy fuzzy sequence on the base which is adopted by Abdul Hameed (every real number r is replaced by a fuzzy number $\\overline{r}$ (either triangular fuzzy number or singleton fuzzy set (fuzzy point))). And then, we will consider that some results respect effect of the upper sequence on the convergent fuzzy sequence, the bounded f...
Directory of Open Access Journals (Sweden)
Aurel Iancu
2008-11-01
Full Text Available After presenting the institutional construction during the pre-accession and post-accession to the Economic and Monetary Union (EMU, the exchange rate mechanisms (ERM in several countries and the convergence criteria, we go on with a brief analysis of the way the CEE countries cope with the convergence criteria in accordance with the Maastricht Treaty. Then, the study deals with a topic often discussed in the scientific literature and included on the agenda of decision-makers at various levels, in order to clarify the following major issues: a shorter transition to the euro, the exchange rate equilibrium versus the inflation rate diminution and the Balassa-Samuelson effect, the exchange rates and the exchange rate deviation index, evidences concerning the real exchange rate equilibrium and the appreciation of the exchange rate in the CEE countries.
Oono, Shuhei; Hatsugai, Yasuhiro
2016-01-01
We have characterized robust propagation modes of electromagnetic waves in helical structures by the section Chern number that is defined for two-dimensional (2D) section of the three-dimensional (3D) Brillouin zone. The Weyl point in the photonic bands is associated with a discontinuous jump of the section Chern number. A spatially localized gaussian basis set is used to calculate the section Chern numbers where we have implemented the divergence-free condition on the each basis function in 3D. Validity of the bulk-edge correspondence in a 3D photonic crystal is discussed in relation to the broken inversion symmetry.
Rybczyński, Józef
2011-02-01
This paper presents the results of computer simulation of bearing misalignment defects in a power turbogenerator. This malfunction is typical for great multi-rotor and multi-bearing rotating machines and very common in power turbo-sets. Necessary calculations were carried out by the computer code system MESWIR, developed and used at the IFFM in Gdansk for calculating dynamics of rotors supported on oil bearings. The results are presented in the form of a set of journal and bush trajectories of all turbo-set bearings. Our analysis focuses on the vibrational effects of displacing the two most vulnerable machine bearings in horizontal and vertical directions by the maximum acceptable range calculated with regard to bearing vibration criterion. This assumption required preliminary assessment of the maximum values for the permissible bearing dislocations. We show the relations between the attributes of the particular bearing trajectories and the bearing displacements in relation to their base design position. The shape and dimensions of bearing trajectories are interpreted based on the theory of hydrodynamic lubrication of oil bearings. It was shown that the relative journal trajectories and absolute bush trajectories carry much important information about the dynamic state of the machine, indicating also the way in which bearings are loaded. Therefore, trajectories can be a source of information about the position and direction of bearing misalignments. This article indicates the potential of using trajectory patterns for diagnosing misalignment defects in rotating machines and suggests including sets of trajectory patterns to the knowledge base of a machine diagnostic system.
Sure, Rebecca; Brandenburg, Jan Gerit
2015-01-01
Abstract In quantum chemical computations the combination of Hartree–Fock or a density functional theory (DFT) approximation with relatively small atomic orbital basis sets of double‐zeta quality is still widely used, for example, in the popular B3LYP/6‐31G* approach. In this Review, we critically analyze the two main sources of error in such computations, that is, the basis set superposition error on the one hand and the missing London dispersion interactions on the other. We review various strategies to correct those errors and present exemplary calculations on mainly noncovalently bound systems of widely varying size. Energies and geometries of small dimers, large supramolecular complexes, and molecular crystals are covered. We conclude that it is not justified to rely on fortunate error compensation, as the main inconsistencies can be cured by modern correction schemes which clearly outperform the plain mean‐field methods. PMID:27308221
Sure, Rebecca; Brandenburg, Jan Gerit; Grimme, Stefan
2016-04-01
In quantum chemical computations the combination of Hartree-Fock or a density functional theory (DFT) approximation with relatively small atomic orbital basis sets of double-zeta quality is still widely used, for example, in the popular B3LYP/6-31G* approach. In this Review, we critically analyze the two main sources of error in such computations, that is, the basis set superposition error on the one hand and the missing London dispersion interactions on the other. We review various strategies to correct those errors and present exemplary calculations on mainly noncovalently bound systems of widely varying size. Energies and geometries of small dimers, large supramolecular complexes, and molecular crystals are covered. We conclude that it is not justified to rely on fortunate error compensation, as the main inconsistencies can be cured by modern correction schemes which clearly outperform the plain mean-field methods. PMID:27308221
2002-01-01
NASA's QuikSCAT satellite has confirmed a 30-year old largely unproven theory that there are two areas near the equator where the winds converge year after year and drive ocean circulation south of the equator. By analyzing winds, QuikSCAT has found a year-round southern and northern Intertropical Convergence Zone. This find is important to climate modelers and weather forecasters because it provides more detail on how the oceans and atmosphere interact near the equator. The Intertropical Convergence Zone (ITCZ) is the region that circles the Earth near the equator, where the trade winds of both the Northern and Southern Hemispheres come together. North of the equator, strong sun and warm water of the equator heats the air in the ITCZ, drawing air in from north and south and causing the air to rise. As the air rises it cools, releasing the accumulated moisture in an almost perpetual series of thunderstorms. Satellite data, however, has confirmed that there is an ITCZ north of the equator and a parallel ITCZ south of the equator. Variation in the location of the ITCZ is important to people around the world because it affects the north-south atmospheric circulation, which redistributes energy. It drastically affects rainfall in many equatorial nations, resulting in the wet and dry seasons of the tropics rather than the cold and warm seasons of higher latitudes. Longer term changes in the ITCZ can result in severe droughts or flooding in nearby areas. 'The double ITCZ is usually only identified in the Pacific and Atlantic Oceans on a limited and seasonal basis,' said Timothy Liu, of NASA's Jet Propulsion Laboratory and California Institute of Technology, Pasadena, Calif., and lead researcher on the project. In the eastern Pacific Ocean, the southern ITCZ is usually seen springtime. In the western Atlantic Ocean, the southern ITCZ was recently clearly identified only in the summertime. However, QuikSCAT's wind data has seen the southern ITCZ in all seasons across the
Convergence of barycentric coordinates to barycentric kernels
Kosinka, Jiří
2016-02-12
We investigate the close correspondence between barycentric coordinates and barycentric kernels from the point of view of the limit process when finer and finer polygons converge to a smooth convex domain. We show that any barycentric kernel is the limit of a set of barycentric coordinates and prove that the convergence rate is quadratic. Our convergence analysis extends naturally to barycentric interpolants and mappings induced by barycentric coordinates and kernels. We verify our theoretical convergence results numerically on several examples.
Increasing dominance of IT in ICT convergence
DEFF Research Database (Denmark)
Henten, Anders; Tadayoni, Reza
The aim of the paper is to examine the increasing dominance of IT companies in the converging ICT industry and, on the basis of this development, to contribute to extending the theoretical understanding of market and industry convergence in the ICT area.......The aim of the paper is to examine the increasing dominance of IT companies in the converging ICT industry and, on the basis of this development, to contribute to extending the theoretical understanding of market and industry convergence in the ICT area....
Energy Technology Data Exchange (ETDEWEB)
Roehle, I.
1999-11-01
A Doppler Global Velocimeter was set up in the frame of a PhD thesis. This velocimeter is optimized to carry out high accuracy, three component, time averaged planar velocity measurements. The anemometer was successfully applied to wind tunnel and test rig flows, and the measurement accuracy was investigated. A volumetric data-set of the flow field inside an industrial combustion chamber was measured. This data field contained about 400.000 vectors. DGV measurements in the intake of a jet engine model were carried out applying a fibre bundle boroskope. The flow structure of the wake of a car model in a wind tunnel was investigated. The measurement accuracy of the DGV-System is {+-}0.5 m/s when operated under ideal conditions. This study can serve as a basis to evaluate the use of DGV for aerodynamic development experiments. (orig.) [German] Im Rahmen der Dissertation wurde ein auf hohe Messgenauigkeit optimiertes DGV-Geraet fuer zeitlich gemittelte Drei-Komponenten-Geschwindigkeitsmessungen entwickelt und gebaut, an Laborstroemungen, an Teststaenden und an Windkanaelen erfolgreich eingesetzt und das Potential der Messtechnik, insbesondere im Hinblick auf Messgenauigkeit, untersucht. Im Fall einer industriellen Brennkammer konnte ein Volumen-Datensatz des Stroemungsfeldes erstellt werden, dessen Umfang bei ca. 400.000 Vektoren lag. Es wurden DGV-Messungen mittels eines flexiblen Endoskops auf Basis eines Faserbuendels durchgefuehrt und damit die Stroemung in einem Flugzeugeinlauf vermessen. Es wurden DGV-Messungen im Nachlauf eines PKW-Modells in einem Windkanal durchgefuehrt. Die Messgenauigkeit des erstellten DGV-Systems betraegt unter Idealbedingungen {+-}0,5 m/s. Durch die Arbeit wurde eine Basis zur Beurteilung des Nutzens der DGV-Technik fuer aerodynamische Entwicklungsarbeiten geschaffen. (orig.)
Hetmaniuk, U.; Lehoucq, R.
2006-10-01
The purpose of our paper is to discuss basis selection for Knyazev's locally optimal block preconditioned conjugate gradient (LOBPCG) method. An inappropriate choice of basis can lead to ill-conditioned Gram matrices in the Rayleigh-Ritz analysis that can delay convergence or produce inaccurate eigenpairs. We demonstrate that the choice of basis is not merely related to computing in finite precision arithmetic. We propose a representation that maintains orthogonality of the basis vectors and so has excellent numerical properties.
Medialogy - convergence and transdisciplinarity
DEFF Research Database (Denmark)
Nordahl, Rolf
2007-01-01
Art and design have many qualities which intuitively by society are appreciated. But on daily basis they are also assessed and evaluated. Some communities also share common understandings that norms, standards, conventions evolve, both from artistic desire as well as from more physical needs for...... changes in society, developments in taste etc. However, it certainly seems fair to say, that available technology makes a great difference to the development of any art form or practice. With the up rise of new educations such as Medialogy, new aspects of convergence and different forms of...
Rangel, Tonatiuh; Genovese, Luigi; Torrent, Marc
2016-01-01
We present a Projector Augmented-Wave~(PAW) method based on a wavelet basis set. We implemented our wavelet-PAW method as a PAW library in the ABINIT package [http://www.abinit.org] and into BigDFT [http://www.bigdft.org]. We test our implementation in prototypical systems to illustrate the potential usage of our code. By using the wavelet-PAW method, we can simulate charged and special boundary condition systems with frozen-core all-electron precision. Furthermore, our work paves the way to large-scale and potentially order-N simulations within a PAW method.
International Nuclear Information System (INIS)
A self-consistent version of the discrete variational method is described based on the use of numerical LCAO basis functions obtained as solutions of the Hartree-Fock equations for free atoms. Sets of single-zeta Slater function are applied to approximate atomic densities further employed in the calculations of the Coulomb potential. The computer programs realizing this approach have been written and utilized to calculated electronic structures of molybdenum, tungsten and uranium hexafluorides. The ionization potentials calculated are in good quantitative agreement with experimental data. The deviations of the calculated valence state IP's from those determined by photoelectron spectroscopy do not exceed 1 eV. (orig.)
Indian Academy of Sciences (India)
M Talebian; E Talebian; A Abdi
2012-05-01
We obtained an approximation of the force ﬁeld of -quartz crystal using a new idea of applying density functional theory [J Purton, R Jones, C R A Catlow and M Leslie, Phys. Chem. Minerals 19, 392 (1993)]. Our calculations were based on B3LYP Hamiltonian [A N Lazarev and A P Mirgorodsky, Phys. Chem. Minerals 18, 231 (1991)] in 6−311+G(2d) basis set for H16Si7O6 cluster and included a unit cell of the lattice. The advantage of our method is the increase in the speed of calculations and the better adaption of simulation results with the experimental data.
International Nuclear Information System (INIS)
The results of LCAO DFT calculations of lattice parameters, cohesive energy and bulk modulus of the crystalline uranium nitrides UN, U2N3 and UN2 are presented and discussed. The LCAO computer codes Gaussian03 and Crystal06 are applied. The calculations are made with the uranium atom relativistic effective small core potential by Stuttgart-Cologne group (60 electrons in the core). The calculations include the U atom basis set optimization. Powell, Hooke-Jeeves, conjugated gradient and Box methods are implemented in the author's optimization package, being external to the codes for molecular and periodic calculations. The basis set optimization in LCAO calculations improves the agreement of the lattice parameter and bulk modulus of UN crystal with the experimental data, the change of the cohesive energy due to the optimization is small. The mixed metallic-covalent chemical bonding is found both in LCAO calculations of UN and U2N3 crystals; UN2 crystal has the semiconducting nature
Balabanov, Nikolai B.; Peterson, Kirk A.
2006-08-01
Recently developed correlation consistent basis sets for the first row transition metal elements Sc-Zn have been utilized to determine complete basis set (CBS) scalar relativistic electron affinities, ionization potentials, and 4s23dn -2-4s1dn -1 electronic excitation energies with single reference coupled cluster methods [CCSD(T), CCSDT, and CCSDTQ] and multireference configuration interaction with three reference spaces: 3d4s, 3d4s4p, and 3d4s4p3d'. The theoretical values calculated with the highest order coupled cluster techniques at the CBS limit, including extrapolations to full configuration interaction, are well within 1kcal/mol of the corresponding experimental data. For the early transition metal elements (Sc-Mn) the internally contracted multireference averaged coupled pair functional method yielded excellent agreement with experiment; however, the atomic properties for the late transition metals (Mn-Zn) proved to be much more difficult to describe with this level of theory, even with the largest reference function of the present work.
Acceleration of Logarithmic Convergence
Gaskin, J. G.; Ford, W. F.
1998-01-01
In this paper, we shall give a characterization of all monotonically decreasing sequence of positive terms, whose sum converge and then introduce a Transformation which can be used to accelerate the convergence of a large class of logarithmically convergent series.
Bongardt, Annette; Torres, Francisco; WUNSCH, Pierre; Hefeker, Carsten; Hermann, Christoph
2013-01-01
The EU has long viewed economic and institutional convergence as important goals, but the results thus far have been decidedly mixed, and there remain several open questions: How exactly should convergence be defined? How much convergence is necessary? What steps can be taken to improve convergence in the EU, and how can success be defined? Finally, how much convergence can be achieved by improving the economic performance in underperforming regions, and how can convergence in the form of har...
Institute of Scientific and Technical Information of China (English)
曹静
2012-01-01
本文在对零售业发展路径的相关文献进行梳理的基础上，从产业融合的视角，对我国零售业的发展路径进行分析，并指出我国现代零售业的发展是产业关联、产业协同和产业融合不断演化的结果，零售产业融合是现代零售业发展的必然趋势。在此基础上，建立现代零售业发展的产业融合理论模型，进而提出零售企业必须抓住产业融合的契机，占领未来市场竞争的制高点。%Based on literature concerning retail development, this article analyzes the development path of retailing from the perspective of industrial convergence and points out that Chinese modern retailing development is the evolution of industrial correlation, synergy and convergence. Retailing convergence is the inevitable result of the development of Chinese modern retailing. Furthermore, a theoretical model of modern retailing convergence is built. The author holds the view that retailing enterprises must catch the opportunity of the industrial convergence and gain an edge competitive advantage.
Strictly convergent analytic structures
Cluckers, Raf; Lipshitz, Leonard
2013-01-01
We give conclusive answers to some questions about definability in analytic languages that arose shortly after the work by Denef and van den Dries, [DD], on $p$-adic subanalytic sets, and we continue the study of non-archimedean fields with analytic structure of [LR3], [CLR1] and [CL1]. We show that the language $L_K$ consisting of the language of valued fields together with all strictly convergent power series over a complete, rank one valued field $K$ can be expanded, in a definitial way, t...
Schmitz, G. J.
2016-01-01
The importance of microstructure simulation in integrated computational materials engineering settings in relation to the added value provided for macroscopic process simulation, as well as the contribution this kind of simulation can make in predicting material properties, are discussed. The roles of microstructure simulation in integrating scales ranging from component/process scales down to atomistic scales, and also in integrating experimental and virtual worlds, are highlighted. The hierarchical data format (HDF5) as a basis for enhancing the interoperability of the heterogeneous range of simulation tools and experimental datasets in the area of computational materials engineering is discussed. Several ongoing developments indicate that HDF5 might evolve into a de facto standard for digital microstructure representation of all length scales.
Alam, Aftab; Johnson, D. D.
2009-09-01
Site-centered, electronic-structure methods use an expansion inside nonoverlapping “muffin-tin” (MT) spheres plus an interstitial basis set. As the boundary separating the more spherical from nonspherical density between atoms, the “saddle-point” radii (SPR) in the density provide an optimal spherical region for expanding in spherical harmonics, as used in augmented plane wave, muffin-tin orbital, and multiple-scattering [Korringa, Kohn, and Rostoker (KKR)] methods. These MT-SPR guarantee unique, convex Voronoi polyhedra at each site, in distinction to Bader topological cells. We present a numerically fast, two-center expansion to find SPR a priori from overlapping atomic charge densities, valid also for disordered alloys. We adopt this MT-SPR basis for KKR in the atomic sphere approximation and study (dis)ordered alloys with large differences in atomic size (fcc CoPt and bcc CrW). For this simple and unique improvement, we find formation energies and structural parameters in strikingly better agreement with more exact methods or experiment, and resolve issues with former results.
Li, Y. Q.; Ma, F. C.; Sun, M. T.
2013-10-01
A full three-dimensional global potential energy surface is reported first time for the title system, which is important for the photodissociation processes. It is obtained using double many-body expansion theory and an extensive set of accurate ab initio energies extrapolated to the complete basis set limit. Such a work can be recommended for dynamics studies of the N(2D) + H2 reaction, a reliable theoretical treatment of the photodissociation dynamics and as building blocks for constructing the double many-body expansion potential energy surface of larger nitrogen/hydrogen containing systems. In turn, a preliminary theoretical study of the reaction N(^2D)+H_2(X^1Σ _g^+)(ν =0,j=0)rArr NH(a^1Δ )+H(^2S) has been carried out with the method of quasi-classical trajectory on the new potential energy surface. Integral cross sections and thermal rate constants have been calculated, providing perhaps the most reliable estimate of the integral cross sections and the rate constants known thus far for such a reaction.
International Nuclear Information System (INIS)
Augmented quasiparticle (QP) mappings, as applied to indistinguishable point sets of (Liouvillian) democratic-recoupled (DR) tensors, provide for a 1:1 invariant labelling of the underlying (disjoint) dual projective map carrier subspaces, where the Liouville pattern basis set is defined via superboson unit-tensor actions on a null space, | )). The co-operative-action Liouville algebras described here imply parallel limitations to Jucys graph recoupling and its related Racah-Wigner (R-W) algebras once DR indistinguishable point tensorial sets are involved, as in non-SR Sn,n>=4 dominant (NMR) spin symmetry. The importance of Sn G-invariants, as labels for disjoint carrier subspaces in such automorphic spin symmetries, arises from their essential role in defining the quantal-completeness of indistinguishable point sets. From the established properties of augmented-QPs as super-bosons (Temme 2002 Int. J. Quantum Chem. 89 429) (i.e., beyond the earlier Hilbert-space-based Louck and Biedenharn boson pattern views), insight into Atiyah and Sutcliffe's (A-S) assertions (Atiyah and Sutcliffe 2002 Proc. R. Soc. A 458 1089) on the limitations of graph recoupling theory to distinct point sets is obtained. This clarifies the wider analytic intractible of automorphic DR spin systems-beyond the Levi-Civita cyclic-commutation (R-W) approach (Levy-Leblond and Levy-Nahas 1965 J. Math. Phys. 6 1372) which holds for a mono-invariant problem. For (rotating-frame) density matrix approaches to [A]n, [A]n(X) and [AX]n(SU(2) x Sn) (dual) NMR systems, the focus is necessarily on the specialized nature of indistinguishable point sets within multiple invariant-theoretic-based, dynamical spin physics. Here, the GI(s) (GI-cardinality) constitute an important part of the dual irrep set, {Dk(U-tilde) x Γ-tilde[λ](v-tilde)(P)}, with combinatorics, as a central facet of invariant theory, playing a crucial role in the concept of 'quantal completeness' and the impact of A-S's assertion on the
Measure theory of statistical convergence
Institute of Scientific and Technical Information of China (English)
2008-01-01
The question of establishing measure theory for statistical convergence has been moving closer to center stage, since a kind of reasonable theory is not only fundamental for unifying various kinds of statistical convergence, but also a bridge linking the studies of statistical convergence across measure theory, integration theory, probability and statistics. For this reason, this paper, in terms of subdifferential, first shows a representation theorem for all finitely additive probability measures defined on the σ-algebra A of all subsets of N, and proves that every such measure can be uniquely decomposed into a convex combination of a countably additive probability measure and a statistical measure (i.e. a finitely additive probability measure μ with μ(k) = 0 for all singletons {k}). This paper also shows that classical statistical measures have many nice properties, such as: The set S of all such measures endowed with the topology of point-wise convergence on A forms a compact convex Hausdorff space; every classical statistical measure is of continuity type (hence, atomless), and every specific class of statistical measures fits a complementation minimax rule for every subset in N. Finally, this paper shows that every kind of statistical convergence can be unified in convergence of statistical measures.
Measure theory of statistical convergence
Institute of Scientific and Technical Information of China (English)
CHENG LiXin; LIN GuoChen; LAN YongYi; LIU Hui
2008-01-01
The question of establishing measure theory for statistical convergence has been moving closer to center stage, since a kind of reasonable theory is not only fundamental for unifying various kinds of statistical convergence, but also a bridge linking the studies of statistical convergence across measure theory, integration theory, probability and statistics. For this reason, this paper, in terms of subdifferential, first shows a representation theorem for all finitely additive probability measures defined on the σ-algebra of all subsets of N, and proves that every such measure can be uniquely decomposed into a convex combination of a countably additive probability measure and a statistical measure (i.e. a finitely additive probability measure μ with μ(k) = 0 for all singletons {k}). This paper also shows that classical statistical measures have many nice properties, such as: The set of all such measures endowed with the topology of point-wise convergence on forms a compact convex Hausdorff space; every classical statistical measure is of continuity type (hence, atomless), and every specific class of statistical measures fits a complementation minimax rule for every subset in N. Finally, this paper shows that every kind of statistical convergence can be unified in convergence of statistical measures.
Topologies of (strong) uniform convergence on bornologies
Holá, Lubica; Novotný, Branislav
2012-01-01
We continue the study of topologies of strong uniform convergence on bornologies initiated in [G. Beer and S. Levi, Strong uniform continuity, J. Math Anal. Appl., 350:568-589, 2009] and [G. Beer and S. Levi, Uniform continuity, uniform convergence and shields, Set-Valued and Variational Analysis, 18:251-275, 2010]. We study cardinal invariants of topologies of (strong) uniform convergence on bornologies on the space of continuous real-valued functions and we also generalize some known result...
Converging Information and Communication Systems
DEFF Research Database (Denmark)
Øst, Alexander
2003-01-01
the future to have - significant importance to the process and consequences of the convergence. The project focuses on the appliances, i.e. the TV sets, the computers and their peripheral equipment. It also takes into account the infrastructure and signals, which contain and deliver the information...
Strong path convergence from Loewner driving convergence
Sheffield, Scott
2010-01-01
We show that, under mild assumptions on the limiting curve, a sequence of simple chordal planar curves converges uniformly whenever certain Loewner driving functions converge. We extend this result to random curves. The random version applies in particular to random lattice paths that have chordal SLE as a scaling limit, with kappa less than 8 (non-space-filling). Existing SLE convergence proofs often begin by showing that the Loewner driving functions of these paths (viewed from infinity) converge to Brownian motion. Unfortunately, this is not sufficient, and additional arguments are required to complete the proofs. We show that driving function convergence is sufficient if it can be established for both parametrization directions and a generic target.
Rodriguez-Bautista, Mariano; Díaz-García, Cecilia; Navarrete-López, Alejandra M; Vargas, Rubicelia; Garza, Jorge
2015-07-21
In this report, we use a new basis set for Hartree-Fock calculations related to many-electron atoms confined by soft walls. One- and two-electron integrals were programmed in a code based in parallel programming techniques. The results obtained with this proposal for hydrogen and helium atoms were contrasted with other proposals to study just one and two electron confined atoms, where we have reproduced or improved the results previously reported. Usually, an atom enclosed by hard walls has been used as a model to study confinement effects on orbital energies, the main conclusion reached by this model is that orbital energies always go up when the confinement radius is reduced. However, such an observation is not necessarily valid for atoms confined by penetrable walls. The main reason behind this result is that for atoms with large polarizability, like beryllium or potassium, external orbitals are delocalized when the confinement is imposed and consequently, the internal orbitals behave as if they were in an ionized atom. Naturally, the shell structure of these atoms is modified drastically when they are confined. The delocalization was an argument proposed for atoms confined by hard walls, but it was never verified. In this work, the confinement imposed by soft walls allows to analyze the delocalization concept in many-electron atoms. PMID:26203010
International Nuclear Information System (INIS)
In this report, we use a new basis set for Hartree-Fock calculations related to many-electron atoms confined by soft walls. One- and two-electron integrals were programmed in a code based in parallel programming techniques. The results obtained with this proposal for hydrogen and helium atoms were contrasted with other proposals to study just one and two electron confined atoms, where we have reproduced or improved the results previously reported. Usually, an atom enclosed by hard walls has been used as a model to study confinement effects on orbital energies, the main conclusion reached by this model is that orbital energies always go up when the confinement radius is reduced. However, such an observation is not necessarily valid for atoms confined by penetrable walls. The main reason behind this result is that for atoms with large polarizability, like beryllium or potassium, external orbitals are delocalized when the confinement is imposed and consequently, the internal orbitals behave as if they were in an ionized atom. Naturally, the shell structure of these atoms is modified drastically when they are confined. The delocalization was an argument proposed for atoms confined by hard walls, but it was never verified. In this work, the confinement imposed by soft walls allows to analyze the delocalization concept in many-electron atoms
Rodriguez-Bautista, Mariano; Díaz-García, Cecilia; Navarrete-López, Alejandra M.; Vargas, Rubicelia; Garza, Jorge
2015-07-01
In this report, we use a new basis set for Hartree-Fock calculations related to many-electron atoms confined by soft walls. One- and two-electron integrals were programmed in a code based in parallel programming techniques. The results obtained with this proposal for hydrogen and helium atoms were contrasted with other proposals to study just one and two electron confined atoms, where we have reproduced or improved the results previously reported. Usually, an atom enclosed by hard walls has been used as a model to study confinement effects on orbital energies, the main conclusion reached by this model is that orbital energies always go up when the confinement radius is reduced. However, such an observation is not necessarily valid for atoms confined by penetrable walls. The main reason behind this result is that for atoms with large polarizability, like beryllium or potassium, external orbitals are delocalized when the confinement is imposed and consequently, the internal orbitals behave as if they were in an ionized atom. Naturally, the shell structure of these atoms is modified drastically when they are confined. The delocalization was an argument proposed for atoms confined by hard walls, but it was never verified. In this work, the confinement imposed by soft walls allows to analyze the delocalization concept in many-electron atoms.
Scalar Curvature and Intrinsic Flat Convergence
Sormani, Christina
2016-01-01
Herein we present open problems and survey examples and theorems concerning sequences of Riemannian manifolds with uniform lower bounds on scalar curvature and their limit spaces. Examples of Gromov and of Ilmanen which naturally ought to have certain limit spaces do not converge with respect to smooth or Gromov-Hausdorff convergence. Thus we focus here on the notion of Intrinsic Flat convergence, developed jointly with Wenger. This notion has been applied successfully to study sequences that arise in General Relativity. Gromov has suggested it should be applied in other settings as well. We first review intrinsic flat convergence, its properties, and its compactness theorems, before presenting the applications and the open problems.
Crespo Cuaresma, Jesus; Klasen, Stephan; Konstantin M. Wacker
2016-01-01
Martin Ravallion ("Why Don't We See Poverty Convergence?" American Economic Review, 102(1): 504-23; 2012) presents evidence against the existence of convergence in global poverty rates despite convergence in household mean income levels and the close linkage between income growth and poverty reduction. We show that this finding is driven by a specification that demands more than simple convergence in poverty headcount rates and assumes a growth elasticity of poverty reduction, which is well-k...
Worldwide Educational Convergence
Directory of Open Access Journals (Sweden)
Connie L. McNeely
1994-11-01
Full Text Available We argue for an examination of the role of the transnational organizational apparatus vis-a-vis nation-states in organizing national educational systems in accordance with world level educational ideologies, structures, and practices. We propose that more analytic attention be given international organizations as an institutionalizing force in examining educational convergence and change, and suggest four primary international organization activities as potentially fruitful avenues for research in this area: 1 the exchange of information, 2 charters and constitutions, 3 standard-setting instruments, and 4 technical and financial resources. Focusing on these activities, we present and discuss evidence of international organizations as world- level agencies influencing the incorporation and diffusion of educational ideologies and practices within and among nation-states.
Pal, Sudip Kumar; Ch, Debraj; ,; Dutta, Sudipta
2013-01-01
In this paper we extend the notion of rough convergence using theconcept of ideals which automatically extends the earlier notions ofrough convergence and rough statistical convergence. We define the setof rough ideal limit points and prove several results associated with thisset.
Mačaj, Martin; Sleziak, Martin
2011-01-01
In this paper we introduce I^K-convergence which is a common generalization of the I^K-convergence of sequences, double sequences and nets. We show that many results that were shown before for these special cases are true for the I^K-convergence, too
Enyeart, Mike; Staman, E. Michael; Valdes, Jose J., Jr.
2007-01-01
The concept of convergence has evolved significantly during recent years. Today, "convergence" refers to the integration of the communications and computing resources and services that seamlessly traverse multiple infrastructures and deliver content to multiple platforms or appliances. Convergence is real. Those in higher education, and especially…
Directory of Open Access Journals (Sweden)
A. P. Tsimpidi
2010-01-01
Full Text Available New primary and secondary organic aerosol modules have been added to PMCAMx, a three dimensional chemical transport model (CTM, for use with the SAPRC99 chemistry mechanism based on recent smog chamber studies. The new modelling framework is based on the volatility basis-set approach: both primary and secondary organic components are assumed to be semivolatile and photochemically reactive and are distributed in logarithmically spaced volatility bins. This new framework with the use of the new volatility basis parameters for low-NO_{x} and high-NO_{x} conditions tends to predict 4–6 times higher anthropogenic SOA concentrations than those predicted with the older generation of models. The resulting PMCAMx-2008 was applied in Mexico City Metropolitan Area (MCMA for approximately a week during April 2003 during a period of very low regional biomass burning impact. The emission inventory, which uses as a starting point the MCMA 2004 official inventory, is modified and the primary organic aerosol (POA emissions are distributed by volatility based on dilution experiments. The predicted organic aerosol (OA concentrations peak in the center of Mexico City, reaching values above 40 μg m^{−3}. The model predictions are compared with the results of the Positive Matrix Factorization (PMF analysis of the Aerosol Mass Spectrometry (AMS observations. The model reproduces both Hydrocarbon-like Organic Aerosol (HOA and Oxygenated Organic Aerosol (OOA concentrations and diurnal profiles. The small OA underprediction during the rush-hour periods and overprediction in the afternoon suggest potential improvements to the description of fresh primary organic emissions and the formation of the oxygenated organic aerosols, respectively, although they may also be due to errors in the simulation of dispersion and vertical mixing. However, the AMS OOA data are not specific enough to prove that the model reproduces the organic aerosol
Tsimpidi, A. P.; Karydis, V. A.; Pandis, S. N.; Zavala, M.; Lei, W.; Molina, L. T.
2007-12-01
Anthropogenic air pollution is an increasingly serious problem for public health, agriculture, and global climate. Organic material (OM) contributes ~ 20-50% to the total fine aerosol mass at continental mid-latitudes. Although OM accounts for a large fraction of PM2.5 concentration worldwide, the contributions of primary and secondary organic aerosol have been difficult to quantify. In this study, new primary and secondary organic aerosol modules were added to PMCAMx, a three dimensional chemical transport model (Gaydos et al., 2007), for use with the SAPRC99 chemistry mechanism (Carter, 2000; ENVIRON, 2006) based on recent smog chamber studies (Robinson et al., 2007). The new modeling framework is based on the volatility basis-set approach (Lane et al., 2007): both primary and secondary organic components are assumed to be semivolatile and photochemically reactive and are distributed in logarithmically spaced volatility bins. The emission inventory, which uses as starting point the MCMA 2004 official inventory (CAM, 2006), is modified and the primary organic aerosol (POA) emissions are distributed by volatility based on dilution experiments (Robinson et al., 2007). Sensitivity tests where POA is considered as nonvolatile and POA and SOA as chemically reactive are also described. In all cases PMCAMx is applied in the Mexico City Metropolitan Area during March 2006. The modeling domain covers a 180x180x6 km region in the MCMA with 3x3 km grid resolution. The model predictions are compared with Aerodyne's Aerosol Mass Spectrometry (AMS) observations from the MILAGRO Campaign. References Robinson, A. L.; Donahue, N. M.; Shrivastava, M. K.; Weitkamp, E. A.; Sage, A. M.; Grieshop, A. P.; Lane, T. E.; Pandis, S. N.; Pierce, J. R., 2007. Rethinking organic aerosols: semivolatile emissions and photochemical aging. Science 315, 1259-1262. Gaydos, T. M.; Pinder, R. W.; Koo, B.; Fahey, K. M.; Pandis, S. N., 2007. Development and application of a three- dimensional aerosol
Directory of Open Access Journals (Sweden)
Q. J. Zhang
2012-11-01
Full Text Available Results of the chemistry transport model CHIMERE are compared with the measurements performed during the MEGAPOLI summer campaign in the Greater Paris Region in July, 2009. The Volatility-Basis-Set approach (VBS is implemented into this model, taking into account the volatility of primary organic aerosol (POA and the chemical aging of semi-volatile organic species. Organic aerosol is the main focus and is simulated with three different configurations related to the volatility of POA and the scheme of secondary organic aerosol (SOA formation. In addition, two types of emission inventories are used as model input in order to test the uncertainty related to the emissions. Predictions of basic meteorological parameters and primary and secondary pollutant concentrations are evaluated and four pollution regimes according to the air mass origin are defined. Primary pollutants are generally overestimated, while ozone is consistent with observations. Sulfate is generally overestimated, while ammonium and nitrate levels are well simulated with the refined emission data set. As expected, the simulation with non-volatile POA and a single-step SOA formation mechanism largely overestimates POA and underestimates SOA. Simulation of organic aerosol with the VBS approach taking into account the aging of semi-volatile organic compounds (SVOC shows the best correlation with measurements. All observed high concentration events are reproduced by the model mostly after long range transport, indicating that long range transport of SOA to Paris is well reproduced. Depending on the emission inventory used, simulated POA levels are either reasonable or underestimated, while SOA levels tend to be overestimated. Several uncertainties related to the VBS scheme (POA volatility, SOA yields, the aging parameterization, to emission input data, and to simulated OH levels can be responsible for this behavior. Despite these uncertainties, the implementation of the VBS scheme into
Knote, C. J.; Hodzic, A.; Aumont, B.; Madronich, S.
2014-12-01
Traditional understanding views secondary organic aerosol (SOA) formation in the atmosphere as continuous gas-phase oxidation of precursors such as isoprene, aromatics or alkanes. Recent research found that these oxidation products are also highly water soluble. It is further understood that the liquid-phase of cloud droplets as well as deliquesced particles could mediate SOA formation through chemistry in the aqueous-phase. While the effect of multi-phase processing has been studied in detailed for specific compounds like glyoxal or methylglyoxal, an integrated approach that considers the large number of individual compounds has been missing due to the complexity involved. In our work we explore the effects of multi-phase processing on secondary organic aerosol from an explicit modeling perspective.Volatility and solubility determine in which phase a given molecule will be found under given atmospheric conditions. Volatility has already been used to simplify the description of SOA formation in the gas-phase in what became known as the Volatility Basis Set approach (VBS). Compounds contributing to SOA formation are grouped by volatility and then treated as a whole. A number of studies extended the VBS by adding a second dimension like oxygen to carbon ratio or the mean oxidation state. In our work we use functional groups as second dimension.Using explicit oxidation chemistry modeling (GECKO-A) we derive SOA yields as well as their composition in terms of functional groups for commonly used precursors. We then investigate the effect of simply partitioning functional-group specific organic mass into cloud droplets and deliquesced aerosol based on their estimated solubility. Finally we apply simple chemistry in the aqueous-phase and relate changes in functional groups to changes in volatility and subsequent changes in partitioning between gas- and aerosol-phase.In our presentation we will explore the sensitivites of the multi-phase system in a box model setting with
Maroulis, George
1998-04-01
The electric multipole moments, dipole and quadrupole polarizability and hyperpolarizability of hydrogen chloride have been determined from an extensive and systematic study based on finite-field fourth-order many-body perturbation theory and coupled-cluster calculations. Our best values for the dipole, quadrupole, octopole and hexadecapole moment at the experimental internuclear separation of Re=2.408645a0 are μ=0.4238ea0, Θ=2.67ea02, Ω=3.94ea03, and Φ=13.37ea04, respectively. For the mean and the anisotropy of the dipole polarizability ααβ we recommend ᾱ=17.41±0.02 and Δα=1.60±0.03e2a02Eh-1. For the mean value of the first dipole hyperpolarizability βαβγ we advance β¯=-6.8±0.3e3a03Eh-2. Extensive calculations with a [8s6p6d3f/5s4p2d1f] basis set at the CCSD(T) level of theory yield the R-dependence of the Cartesian components and the mean of the second dipole hyperpolarizability γαβγδ(R)/e4a04Eh-3 around Re as γzzzz(R)=1907+1326(R-Re)+570(R-Re)2+10(R-Re)3-40(R-Re)4, γxxxx(R)=3900+747(R-Re)-65(R-Re)2-38(R-Re)3-7(R-Re)4, γxxzz(R)=962+222(R-Re)+88(R-Re)2+49(R-Re)3+5(R-Re)4, γ¯(R)=3230+841(R-Re)+151(R-Re)2+21(R-Re)3-9(R-Re)4, with z as the molecular axis. The present investigation suggests an estimate of (26.7±0.3)×102e4a04Eh-3 for the Hartree-Fock limit of the mean value γ¯ at Re. CCSD(T) calculations with basis sets of [8s6p6d3f/5s4p2d1f] and [9s7p5d4f/6s5p4d1f] size and MP4 calculations with the even larger [15s12p7d3f/12s7p2d1f] give (7.0±0.3)×102e4a04Eh-3 for the electron correlation effects for this property, thus leading to a recommended value of γ¯=(33.7±0.6)×102e4a04Eh-3. For the quadrupole polarizability Cαβ,γδ/e2a04Eh-1 at Re our best values are Czz,zz=41.68, Cxz,xz=26.11, and Cxx,xx=35.38, calculated with the [9s7p5d4f/6s5p4d1f] basis set at the CCSD(T) level of theory. The following CCSD(T) values were obtained with [8s6p6d3f/5s4p2d1f] at Re: dipole-quadrupole polarizability Aα,βγ/e2a03Eh-1, Az,zz=14.0, and
Evolution and convergence in telecommunications
International Nuclear Information System (INIS)
These lectures throw a spotlight on different aspects of the evolution of telecommunications networks, namely on the various facets of service and network convergence. The last years progress in data and telecommunications technologies, such as P-based networks, and the enormous potential of mobile communication systems and users' demands for comprehensive and network-independent have led to a convergence of data and telecommunications infrastructures in many aspects. In order to help the reader to an easier understanding of the phenomenon convergence, in the first two parts of this volume the evolution of the basic technologies is described one by one. This is done briefly and is focused on the principle topics, just to build a basis for the third part devoted to problems of convergence in telecommunications. These notes are addressed to those readers, who in a quick overview want to be informed on the future service and network landscape. The notes are equally suited for professionals with the desire to extend their horizon as well as for students looking for an introduction into telecommunications under more general aspects. The authors clearly understand the difficulties in writing a book devoted to the evolution in telecommunications. Today, telecom landscape varies at very high speed. Every few months new network technologies, new products and new services are developed. Attempts to present them in time can be accessible only for magazine publications or contributions to conferences. Therefore, in a number of areas, such as Voice over IP and new switching technologies not much more than the starting point of new paradigms is described. However, the content is up-to-date to a degree, that the phenomenon convergence can be fully understood. Partially, these notes are based on a number of lecture courses that were delivered during recent ICTP winter schools devoted to multimedia and digital communications
Jahan, Sumbul
2013-01-01
DOES CONVERGENCE EXIST? Sumbul Jahan Institute of Business Administration (IBA), 2013 Research Project Supervisor: Dr. Farooq Pasha & Dr. Heman D. Lohano The idea of convergence in economics is the hypothesis that poorer economies income will tend to grow at faster rates than richer economies. As a result, all economies should eventually converge; Developing countries have the potential to grow at a faster rate than the developed countries because of availability of better health...
Convergent Yang-Mills matrix theories
International Nuclear Information System (INIS)
We consider the partition function and correlation functions in the bosonic and supersymmetric Yang-Mills matrix models with compact semi-simple gauge group. In the supersymmetric case, we show that the partition function converges when D=4, 6 and 10, and that correlation functions of degree kc=2(D-3) are convergent independently of the group. In the bosonic case we show that the partition function is convergent when D≥Dc, and that correlation functions of degree kc are convergent, and calculate Dc and kc for each group, thus extending our previous results for SU(N). As a special case these results establish that the partition function and a set of correlation functions in the IKKT IIB string matrix model are convergent. (author)
DEFF Research Database (Denmark)
Lema, Rasmus; Zhou, Yuan; Sagar, Ambuj
2016-01-01
characteristics and therefore differ along key dimensions, the increasing roles of cross-national firm interactions amplify tendencies towards global convergence. These patterns of divergence and convergence can potentially enhance the contribution of wind power to the low-carbon technology transition but also...
Concordant Convergence Empirics
Webber, Don J; Paul White
2004-01-01
We present a new model to test the convergence hypothesis based on the ideas of concordance and then employ the model to test empirically for GDP per capita convergence across 97 countries. Our results suggest the presence of switching, while there is more ‘strong divergence’ than ‘strong convergence’.
Basis for breakup states of three identical particles
International Nuclear Information System (INIS)
A new basis for expanding three-body momentum-space states for three identical particles is studied. The basis states are simultaneously eigenstates of the total angular momentum and the total antisymmetrization operator. The total kinetic energy and two Dalitz-Fabri variables are chosen as the remaining three continuous variables. Zernike polynomials are used as a basis set for a generalized Fourier expansion in the Dalitz-Fabri variables. Born approximations to the nucleon-deuteron breakup amplitude zero total orbital angular momentum) are calculated for Malfliet-Tjon I-III potentials and displayed in a Dalitz plot that shows the global structures of the reaction probabilities. Numerical results are presented, which indicate favorable convergence properties of the generalized Fourier expansion. These results suggest that the new basis set may be attractive in more realistic calculations. (author)
Normalized Convergence in Stochastic Optimization
Y.M. Ermoliev; V.I. Norkin
1989-01-01
A new concept of (normalized) convergence of random variables is introduced. The normalized convergence in preserved under Lipschitz transformations. This convergence follows from the convergence in mean and itself implies the convergence in probability. If a sequence of random variables satisfies a Limit theorem then it in a normalized convergent sequence. The introduced concept is applied to the convergence rate study of a statistical approach in stochastic optimization.
Bor, G.; Kettle, Sidney F. A.
1999-12-01
The unifying power of the use of trigonometric basis functions in group character tables is demonstrated. Additionally, these functions provide a simple way of generating pictures of symmetry coordinates. This is illustrated for the in-plane stretch and out-of-plane deformation motions of the C-H bonds in benzene. Their application to orbital symmetry applications is also indicated.
Convergences and divergences in mortality
Directory of Open Access Journals (Sweden)
2004-04-01
Full Text Available Abdel Omran's 1971 theory of "Epidemiologic Transition" was the first attempt to account for the extraordinary advances in health care made in industrialized countries since the 18th century. In the framework of the Demographic Transition, it implied a general convergence of life expectancies toward a limit imposed by the new epidemiological features of modern societies. However, important failures, occurred in the past decades (mainly the health crisis in Eastern Europe and AIDS in Africa, seem to have stopped that process of convergence. In fact such failures do not really contradict the theory. The latter is much more ruined by the unexpected dramatic improvement in the field of cardiovascular disease experienced since the seventies, which results in a new step of a more general process. On the basis of the broader concept of "Health Transition" initiated by Julio Frenk et al., the present paper tries to rethink the full process in term of divergence/convergence sequences inferred by successive major changes in health technologies and strategies.
Quiney, HM; Glushkov, VN; Wilson, S
2002-01-01
Using basis sets of distributed s-type Gaussian functions with positions and exponents optimized so as to support Hartree-Fock total energies with an accuracy approaching the sub-muHartree level, Dirac-Hartree-Fock-Coulomb calculations are reported for the ground states of the H-2, LiH, and BH molec
Czech Academy of Sciences Publication Activity Database
Li, F.; Wang, L.; Zhao, J.; Xie, J. R. H.; Riley, Kevin Eugene; Chen, Z.
2011-01-01
Roč. 130, 2/3 (2011), s. 341-352. ISSN 1432-881X Institutional research plan: CEZ:AV0Z40550506 Keywords : water cluster * density functional theory * MP2 . CCSD(T) * basis set * relative energies Subject RIV: CF - Physical ; Theoretical Chemistry Impact factor: 2.162, year: 2011
Economic convergence and climate policy
International Nuclear Information System (INIS)
This paper addresses the relevance of the economic convergence hypotheses between the developing and the developed world in international greenhouse gas (GHG) emissions negotiations. The results are based on a two-region (the OECD and the rest of the world, ROW) neo-classical growth model with exogenous technical progress, different technological diffusion patterns, and a set of geophysical relationships that consider an environmental externality linked to GHG emissions. A game framework is taken into account in the model to capture the strategic interactions between agents. The outcome of the negotiations seems indeed to depend on the economic convergence hypotheses. Faster economic growth of the ROW countries would encourage them to further mitigate carbon emissions. (Author)
The Convergence of European Business Cycles 1980-2004
International Nuclear Information System (INIS)
The degree of convergence of the business cycles of the economies of the European Union is a key policy issue. In particular, a substantial degree of convergence is needed if the European Central Bank is to be capable of setting a monetary policy which is appropriate to the stage of the cycle of the Euro zone economies. I consider the annual rates of real GDP growth on a quarterly basis in the main economies of the EU (France, Germany, Italy, UK, Spain, Belgium and the Netherlands) over the period 1980Q1-2004Q4. An important empirical question is the degree to which the correlations between these growth rates contain true information rather than noise. The technique of random matrix theory is able to answer this question, and has been applied successfully in the physics journals to financial markets data. I find that the correlations between the growth rates of most of the core EU economies contain substantial amounts of true information, and exhibit considerable stability over time. Even in the late 1970s and early 1980s, these economies moved together closely over the course of the business cycle. There was a slight loosening at the time of German re-unification, but the economies have moved back into close synchronisation. The same result holds when Spain is added to the group of core EU countries. However, the problems of the German economy which arose from the early 1990s onwards has led to Germany becoming increasingly less synchronised with the rest of the core EU. Further, the results obtained with a data set of the converged EU core plus the UK show no real convergence between the UK and this group of economies. (author)
The convergence of corporate social responsibility practices
Misani, Nicola
2010-01-01
Purpose – This paper tries to explain why many socially-responsible firms appear to converge on a standard set of corporate social responsibility (CSR) practices instead of striving to differentiate themselves from rivals and achieve competitive advantage. Design/methodology/approach – Three explanations of this convergence are presented: herd behaviour, institutional isomorphism, and strategic cooperation. The different empirical predictions of these theories are laid down. The resulting fra...
On Conditions for Convergence to Consensus
Lorenz, Jan; Lorenz, Dirk A.
2008-01-01
A new theorem on conditions for convergence to consensus of a multiagent time-dependent time-discrete dynamical system is presented. The theorem is build up on the notion of averaging maps. We compare this theorem to results by Moreau (IEEE Transactions on Automatic Control, vol. 50, no. 2, 2005) about set-valued Lyapunov theory and convergence under switching communication topologies. We give examples that point out differences of approaches including examples where Moreau's theorem is not a...
An $L^p$ theory of sparse graph convergence II: LD convergence, quotients, and right convergence
Borgs, Christian; Chayes, Jennifer T.; Cohn, Henry; Zhao, Yufei
2014-01-01
We extend the $L^p$ theory of sparse graph limits, which was introduced in a companion paper, by analyzing different notions of convergence. Under suitable restrictions on node weights, we prove the equivalence of metric convergence, quotient convergence, microcanonical ground state energy convergence, microcanonical free energy convergence, and large deviation convergence. Our theorems extend the broad applicability of dense graph convergence to all sparse graphs with unbounded average degre...
Toeplitz Lemma, Complete Convergence and Complete Moment Convergence
Li, Jiyanglin; Hu, Ze-Chun
2014-01-01
In this paper, we study the Toeplitz lemma, the Ces\\`{a}ro mean convergence theorem and the Kronecker lemma. At first, we study "complete convergence" versions of the Toeplitz lemma, the Ces\\`{a}ro mean convergence theorem and the Kronecker lemma. Two counterexamples show that they can fail in general and some sufficient conditions for "complete convergence" version of the Ces\\`{a}ro mean convergence theorem are given. Secondly we introduce two classes of complete moment convergence, which ar...
OpenMC In Situ Source Convergence Detection
Energy Technology Data Exchange (ETDEWEB)
Aldrich, Garrett Allen [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Univ. of California, Davis, CA (United States); Dutta, Soumya [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); The Ohio State Univ., Columbus, OH (United States); Woodring, Jonathan Lee [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
2016-05-07
We designed and implemented an in situ version of particle source convergence for the OpenMC particle transport simulator. OpenMC is a Monte Carlo based-particle simulator for neutron criticality calculations. For the transport simulation to be accurate, source particles must converge on a spatial distribution. Typically, convergence is obtained by iterating the simulation by a user-settable, fixed number of steps, and it is assumed that convergence is achieved. We instead implement a method to detect convergence, using the stochastic oscillator for identifying convergence of source particles based on their accumulated Shannon Entropy. Using our in situ convergence detection, we are able to detect and begin tallying results for the full simulation once the proper source distribution has been confirmed. Our method ensures that the simulation is not started too early, by a user setting too optimistic parameters, or too late, by setting too conservative a parameter.
Convergence in pigmentation at multiple levels: mutations, genes and function
Manceau, Marie; Domingues, Vera S.; Linnen, Catherine R.; Rosenblum, Erica Bree; Hoekstra, Hopi E.
2010-01-01
Convergence—the independent evolution of the same trait by two or more taxa—has long been of interest to evolutionary biologists, but only recently has the molecular basis of phenotypic convergence been identified. Here, we highlight studies of rapid evolution of cryptic coloration in vertebrates to demonstrate that phenotypic convergence can occur at multiple levels: mutations, genes and gene function. We first show that different genes can be responsible for convergent phenotypes even among...
Convergence in variable Lebesgue spaces
Cruz-Uribe, David; SFO; Fiorenza, Alberto
2010-01-01
We consider the relationship in the variable Lebesgue space Lp(•)(Ω) between convergence in norm, convergence in modular, and convergence in measure, for both bounded and unbounded exponent functions.
Convergence in variable Lebesque spaces
Cruz-Uribe, David; Fiorenza, Alberto
2010-01-01
We consider the relationship in the variable Lebesgue space $L^{p(\\cdot)}(\\Omega)$ between convergence in norm, convergence in modular, and convergence in measure, for both bounded and unbounded exponent functions.
Convergence properties and compactifications
Leek, Robert
2014-01-01
In this paper, we will use investigate the existence of compactifications with particular convergence properties - pseudoradial, radial, sequential and Fr\\'echet-Urysohn - through the use of spoke systems.
Convergent Aeronautics Solutions Project
National Aeronautics and Space Administration — The Convergent Aeronautics Solutions (CAS) Project uses short-duration activities to establish early-stage concept and technology feasibility for high-potential...
Fixed mobile convergence handbook
Ahson, Syed A
2010-01-01
From basic concepts to future directions, this handbook provides technical information on all aspects of fixed-mobile convergence (FMC). The book examines such topics as integrated management architecture, business trends and strategic implications for service providers, personal area networks, mobile controlled handover methods, SIP-based session mobility, and supervisory and notification aggregator service. Case studies are used to illustrate technical and systematic implementation of unified and rationalized internet access by fixed-mobile network convergence. The text examines the technolo
Hefeker, Carsten; Torres, Francisco S; Bongardt, Annette; WUNSCH, Pierre; Hermann, Christoph
2013-01-01
The EU has long viewed economic and institutional convergence as important goals, but the results thus far have been decidedly mixed, and there remain several open questions: How exactly should convergence be defined? How much convergence is necessary? What steps can be taken to improve convergence in the EU, and how can success be defined? Finally, how much convergence can be achieved by improving the economic performance in underperforming regions, and how can convergence in the form of har...
International Nuclear Information System (INIS)
This book explains IT-BT convergence technology as the future technology, which includes a prolog, easy IT-BT convergence technology that has infinite potentials for new value, policy of IT-BT convergence technology showing the potential of smart Korea, IT-BT convergence opening happy future, for the new future of IT powerful nation Korea with IT-BT convergence technology and an epilogue. This book reveals the conception, policy, performance and future of IT-BT convergence technology.
Technology assessment using NBIC convergence
International Nuclear Information System (INIS)
Full text: Notwithstanding progress in the areas of nanotechnology/nanoscience, biotechnology, information technology, and cognitive sciences (NBIC), the synergy arising from convergence of these disciplines offers great potential for transformational, revolutionary, and embryonic opportunities with many technological applications. In addition, advances in synthesis and characterization methods have provided the means to study, understand, control, or even manipulate the transitional characteristics between isolated atoms and bulk material. In recent years, newly developed architectures in nanostructures and nanosystems with improved functionality, and ensuing unique characteristics have been developed with applications in chemical and biological sensors, nanobiotechnology, nanophotonics, and analysis of cellular processes. Novel convergence methodologies will integrate and advance next generation solutions to current and future technical challenges. Convergence in research methodologies transform the way research is conducted by overcoming specific barriers or filling existing knowledge gaps. NBIC Convergence and associated research methodologies has exceptionally high potential for transforming the manner in which state-of-the-art information is gathered, analyzed, and leveraged to enable future advances and applications. Examples of synergy of these disciplines include: label free, highly multiplexed over broad dynamic range, and decentralized nanotechnology based sensor platform for detection of biological and chemical agents; nucleic acid layers in conjunction with nanomaterials-based electrochemical or optical transducers as DNA Biosensor; potential targets for the next generation of vaccines, therapeutics and diagnostics to enhance human performance; basis for 'Gene Ontology' to provides an important link between gene function and systems biology to understand a global picture of host-microbe interactions - to name a few. Since the idea of 'Converging
Revisiting world energy intensity convergence for regional differences
Energy Technology Data Exchange (ETDEWEB)
Liddle, Brantley [Centre for Strategic Economic Studies, Victoria University, Melbourne, VIC 8001 (Australia)
2010-10-15
World convergence in energy intensity is revisited using two new large data sets: a 111-country sample spanning 1971-2006, and a 134-country sample spanning 1990-2006. Both data sets confirm continued convergence. However, the larger data set, which adds the former Soviet Union republics and additional Balkan countries, indicates greater convergence over its more recent time-frame. Further investigation of geographical differences reveals that the OECD and Eurasian countries have shown considerable, continued convergence, while the Sub-Saharan African countries have converged amongst themselves, but at a slower rate than the OECD and Eurasian countries; by contrast, Latin American and Caribbean and Middle East and North African countries have exhibited no convergence to divergence in energy intensity. (author)
Ensemble based convergence assessment of biomolecular trajectories
Lyman, E; Lyman, Edward; Zuckerman, Daniel M.
2006-01-01
Assessing the convergence of a biomolecular simulation is an essential part of any computational investigation. This is because many important quantities (e.g., free energy differences) depend on the relative populations of different conformers; insufficient convergence translates into systematic errors. Here we present a simple method to self-consistently assess the convergence of a simulation. Standard clustering methods first generate a set of reference structures to any desired precision. The trajectory is then classified by proximity to the reference structures, yielding a one-dimensional histogram of structurally distinct populations. Comparing ensembles of different trajectories (or different parts of the same trajectory) built with the same reference structures provides a sensitive, quantitative measure of convergence. Please note: this is a preliminary manuscript, and should be read as such. Comments are most welcome, especially regarding pertinent prior work.
"First pain" in humans: convergent and specific forebrain responses
Directory of Open Access Journals (Sweden)
Tran Tuan D
2010-11-01
Full Text Available Abstract Background Brief heat stimuli that excite nociceptors innervated by finely myelinated (Aδ fibers evoke an initial, sharp, well-localized pain ("first pain" that is distinguishable from the delayed, less intense, more prolonged dull pain attributed to nociceptors innervated by unmyelinated (C fibers ("second pain". In the present study, we address the question of whether a brief, noxious heat stimulus that excites cutaneous Aδ fibers activates a distinct set of forebrain structures preferentially in addition to those with similar responses to converging input from C fibers. Heat stimuli at two temperatures were applied to the dorsum of the left hand of healthy volunteers in a functional brain imaging (fMRI paradigm and responses analyzed in a set of volumes of interest (VOI. Results Brief 41°C stimuli were painless and evoked only C fiber responses, but 51°C stimuli were at pain threshold and preferentially evoked Aδ fiber responses. Most VOI responded to both intensities of stimulation. However, within volumes of interest, a contrast analysis and comparison of BOLD response latencies showed that the bilateral anterior insulae, the contralateral hippocampus, and the ipsilateral posterior insula were preferentially activated by painful heat stimulation that excited Aδ fibers. Conclusions These findings show that two sets of forebrain structures mediate the initial sharp pain evoked by brief cutaneous heat stimulation: those responding preferentially to the brief stimulation of Aδ heat nociceptors and those with similar responses to converging inputs from the painless stimulation of C fibers. Our results suggest a unique and specific physiological basis, at the forebrain level, for the "first pain" sensation that has long been attributed to Aδ fiber stimulation and support the concept that both specific and convergent mechanisms act concurrently to mediate pain.
Energy Technology Data Exchange (ETDEWEB)
Kobus, J [Instytut Fizyki, Uniwersytet Mikolaja Kopernika, ul. Grudziadzka 5/7, 87-100 Torun (Poland); Moncrieff, D [School of Computational Science and Information Technology, Florida State University, Tallahassee, FL 32306 (United States); Wilson, S [Department of Chemical Physics, Faculty of Mathematics and Physics, Comenius University, 842 15 Bratislava (Slovakia)
2007-03-14
We investigate the accuracy with which the electric dipole polarizability, {alpha}{sub zz}, and the hyperpolarizability, {beta}{sub zzz}, can be calculated by using the algebraic approximation, i.e. finite basis set expansions, and by means of the finite difference method in calculations for the ground states of the 14 electron systems N{sub 2}, CO and BF within the Hartree-Fock model at their respective experimental equilibrium geometries. For a well-chosen grid, the finite difference technique can provide Hartree-Fock energy and dipole moment expectation values approaching machine precision which can be used to assess the accuracy of corresponding calculations carried out within the algebraic approximation. The finite field approximation is used to determine polarizabilities and hyperpolarizabilities from finite difference Hartree-Fock dipole moment expectation values. The results are compared with finite basis set calculations of the corresponding quantities which are carried out analytically using coupled perturbed Hartree-Fock theory. For the N{sub 2} molecule, the Hartree-Fock polarizability is found to be 14.9512 au within the finite basis set approximation and 14.945 au within the finite difference approach. For the CO molecule, the corresponding results are 14.4668 au and 14.4668 au, whilst for the BF molecule the values are 16.6450 au and 16.6450 au, respectively. The Hartree-Fock hyperpolarizability of the CO molecule is found to be 31.4081 au and 31.411 au within the finite basis set and finite difference approximations, respectively. The corresponding hyperpolarizability values for the BF molecule are 63.9687 au and 63.969 au, respectively.
Extended Global Convergence Framework for Unconstrained Optimization
Institute of Scientific and Technical Information of China (English)
(A)rpád B(U)RMEN; Franc BRATKOVI(C); Janez PUHAN; Iztok FAJFAR; Tadej TUMA
2004-01-01
An extension of the global convergence framework for unconstrained derivative-free optimization methods is presented. The extension makes it possible for the framework to include optimization methods with varying cardinality of the ordered direction set. Grid-based search methods are shown to be a special case of the more general extended global convergence framework. Furthermore,the required properties of the sequence of ordered direction sets listed in the definition of grid-based methods are re]axed and simplified by removing the requirement of structural equivalence.
Fresnel lens solar beam convergence properties
Huang, Minshuang; Tian, Feng
2014-11-01
Use the law of refraction of light for theoretical analysis of sunlight incident Fresnel lens, summary of the relationship between Fresnel lens solar light beam convergence characteristics and beam incident angle,beam incident position.Through the analysis of the calculation and experimental, to obtain a variety of incident angle lens center and best gathering focal spot the size of the data.Discuss the incident angle, incident location,the relationship between convergence center and focal spot size.For the sun acquisition system of the light from the sun and optical fiber coupling design do sufficient theory design basis.
Directory of Open Access Journals (Sweden)
Johanna Egetemeir
2011-09-01
Full Text Available Many everyday life situations require two or more individuals to execute actions together. Assessing brain activation during naturalistic tasks to uncover relevant processes underlying such real-life joint action situations has remained a methodological challenge. In the present study, we introduce a novel joint action paradigm that enables the assessment of brain activation during real-life joint action tasks using functional near-infrared spectroscopy (fNIRS. We monitored brain activation of participants who coordinated complex actions with a partner sitting opposite them. Participants performed table-setting tasks, either alone (solo action or in cooperation with a partner (joint action, or they observed the partner performing the task (action observation. Comparing joint action and solo action revealed stronger activation (higher [oxy-Hb]-concentration during joint action in a number of areas. Among these were areas in the inferior parietal lobule (IPL that additionally showed an overlap of activation during action observation and solo action. Areas with such a close link between action observation and action execution have been associated with action simulation processes. The magnitude of activation in these IPL areas also varied according to joint action type and its respective demand on action simulation. The results validate fNIRS as an imaging technique for exploring the functional correlates of interindividual action coordination in real-life settings and suggest that coordinating actions in real-life situations requires simulating the actions of the partner.