Rossi, Tuomas P. Sakko, Arto; Puska, Martti J.; Lehtola, Susi; Nieminen, Risto M.
2015-03-07
We present an approach for generating local numerical basis sets of improving accuracy for first-principles nanoplasmonics simulations within time-dependent density functional theory. The method is demonstrated for copper, silver, and gold nanoparticles that are of experimental interest but computationally demanding due to the semi-core d-electrons that affect their plasmonic response. The basis sets are constructed by augmenting numerical atomic orbital basis sets by truncated Gaussian-type orbitals generated by the completeness-optimization scheme, which is applied to the photoabsorption spectra of homoatomic metal atom dimers. We obtain basis sets of improving accuracy up to the complete basis set limit and demonstrate that the performance of the basis sets transfers to simulations of larger nanoparticles and nanoalloys as well as to calculations with various exchange-correlation functionals. This work promotes the use of the local basis set approach of controllable accuracy in first-principles nanoplasmonics simulations and beyond.
Rossi, Tuomas P; Lehtola, Susi; Sakko, Arto; Puska, Martti J; Nieminen, Risto M
2015-03-01
We present an approach for generating local numerical basis sets of improving accuracy for first-principles nanoplasmonics simulations within time-dependent density functional theory. The method is demonstrated for copper, silver, and gold nanoparticles that are of experimental interest but computationally demanding due to the semi-core d-electrons that affect their plasmonic response. The basis sets are constructed by augmenting numerical atomic orbital basis sets by truncated Gaussian-type orbitals generated by the completeness-optimization scheme, which is applied to the photoabsorption spectra of homoatomic metal atom dimers. We obtain basis sets of improving accuracy up to the complete basis set limit and demonstrate that the performance of the basis sets transfers to simulations of larger nanoparticles and nanoalloys as well as to calculations with various exchange-correlation functionals. This work promotes the use of the local basis set approach of controllable accuracy in first-principles nanoplasmonics simulations and beyond. PMID:25747068
Property-optimized Gaussian basis sets for molecular response calculations
NASA Astrophysics Data System (ADS)
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.
Optimization of selected molecular orbitals in group basis sets.
Ferenczy, György G; Adams, William H
2009-04-01
We derive a local basis equation which may be used to determine the orbitals of a group of electrons in a system when the orbitals of that group are represented by a group basis set, i.e., not the basis set one would normally use but a subset suited to a specific electronic group. The group orbitals determined by the local basis equation minimize the energy of a system when a group basis set is used and the orbitals of other groups are frozen. In contrast, under the constraint of a group basis set, the group orbitals satisfying the Huzinaga equation do not minimize the energy. In a test of the local basis equation on HCl, the group basis set included only 12 of the 21 functions in a basis set one might ordinarily use, but the calculated active orbital energies were within 0.001 hartree of the values obtained by solving the Hartree-Fock-Roothaan (HFR) equation using all 21 basis functions. The total energy found was just 0.003 hartree higher than the HFR value. The errors with the group basis set approximation to the Huzinaga equation were larger by over two orders of magnitude. Similar results were obtained for PCl(3) with the group basis approximation. Retaining more basis functions allows an even higher accuracy as shown by the perfect reproduction of the HFR energy of HCl with 16 out of 21 basis functions in the valence basis set. When the core basis set was also truncated then no additional error was introduced in the calculations performed for HCl with various basis sets. The same calculations with fixed core orbitals taken from isolated heavy atoms added a small error of about 10(-4) hartree. This offers a practical way to calculate wave functions with predetermined fixed core and reduced base valence orbitals at reduced computational costs. The local basis equation can also be used to combine the above approximations with the assignment of local basis sets to groups of localized valence molecular orbitals and to derive a priori localized orbitals. An
Gidofalvi, Gergely; Mazziotti, David A
2014-01-16
Molecule-optimized basis sets, based on approximate natural orbitals, are developed for accelerating the convergence of quantum calculations with strongly correlated (multireferenced) electrons. We use a low-cost approximate solution of the anti-Hermitian contracted Schrö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 ACSE improves the accuracy of the equivalent standard atom-optimized basis set at little additional computational cost. We illustrate the method with the potential energy curves of hydrogen fluoride and diatomic nitrogen. Relative to the hydrogen fluoride potential energy curve from the ACSE in a polarized triple-ζ basis set, the ACSE curve in a molecule-optimized basis set, equivalent in size to a polarized double-ζ basis, has a nonparallelity error of 0.0154 au, which is significantly better than the nonparallelity error of 0.0252 au from the polarized double-ζ basis set. PMID:24387056
Pulay forces from localized orbitals optimized in situ using a psinc basis set.
Ruiz-Serrano, Álvaro; Hine, Nicholas D M; Skylaris, Chris-Kriton
2012-06-21
In situ optimization of a set of localized orbitals with respect to a systematically improvable basis set independent of the position of the atoms, such as psinc functions, would theoretically eliminate the correction due to Pulay forces from the total ionic forces. We demonstrate that for strict localization constraints, especially with small localization regions, there can be non-negligible Pulay forces that must be calculated as a correction to the Hellmann-Feynman forces in the ground state. Geometry optimization calculations, which rely heavily upon accurate evaluation of the total ionic forces, show much better convergence when Pulay forces are included. The more conventional case, where the local orbitals remain fixed to pseudo-atomic orbital multiple-ζ basis sets, also benefits from this implementation. We have validated the method on several test cases, including a DNA fragment with 1045 atoms. PMID:22779575
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.
Geminal embedding scheme for optimal atomic basis set construction in correlated calculations
NASA Astrophysics Data System (ADS)
Sorella, S.; Devaux, N.; Dagrada, M.; Mazzola, G.; Casula, M.
2015-12-01
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
NASA Astrophysics Data System (ADS)
Gidopoulos, Nikitas I.; Lathiotakis, Nektarios N.
2013-10-01
The Comment by Friedrich does not dispute the central result of our paper [Phys. Rev. APLRAAN1050-294710.1103/PhysRevA.85.052508 85, 052508 (2012)] that nonanalytic behavior is present in long-established mathematical pathologies arising in the solution of finite basis optimized effective potential (OEP) equations. In the Comment, the terms “balancing of basis sets” and “basis-set convergence” imply a particular order towards the limit of a large orbital basis sets where the large-orbital-base limit is always taken first, before the large-auxiliary-base limit, until overall convergence is achieved, at a high computational cost. The authors claim that, on physical grounds, this order of limits is not only sufficient, but also necessary in order to avoid the mathematical pathologies. In response to the Comment, we remark that it is already written in our paper that the nonanalyticity trivially disappears with large orbital basis sets. We point out that the authors of the Comment give an incorrect proof of this statement. We also show that the order of limits towards convergence of the potential is immaterial. A recent paper by the authors of the Comment proposes a partial correction for the incomplete orbital basis error in the full-potential linearized augmented-plane-wave method. Similar to the correction developed in our paper, this correction also benefits from an effectively complete orbital basis, even though only a finite orbital basis is employed in the calculation. This shows that it is unnecessary to take, in practice, the limit of an infinite orbital basis in order to avoid mathematical pathologies in the OEP. Our paper is a significant contribution in that direction with general applicability to any choice of basis sets. Finally, contrary to an allusion in the abstract and assertions in the main text of the Comment that unphysical oscillations of the OEP are supposedly attributed to the common energy denominator approximation, in fact, such
An 'optimal' spawning algorithm for adaptive basis set expansion in nonadiabatic dynamics
Yang, Sandy; Coe, Joshua D.; Kaduk, Benjamin; Martinez, Todd J.
2009-04-07
The full multiple spawning (FMS) method has been developed to simulate quantum dynamics in the multistate electronic problem. In FMS, the nuclear wave function is represented in a basis of coupled, frozen Gaussians, and a 'spawning' procedure prescribes a means of adaptively increasing the size of this basis in order to capture population transfer between electronic states. Herein we detail a new algorithm for specifying the initial conditions of newly spawned basis functions that minimizes the number of spawned basis functions needed for convergence. 'Optimally' spawned basis functions are placed to maximize the coupling between parent and child trajectories at the point of spawning. The method is tested with a two-state, one-mode avoided crossing model and a two-state, two-mode conical intersection model.
Feller, D; Schuchardt, Karen L.; Didier, Brett T.; Elsethagen, Todd; Sun, Lisong; Gurumoorthi, Vidhya; Chase, Jared; Li, Jun
The Basis Set Exchange (BSE) provides a web-based user interface for downloading and uploading Gaussian-type (GTO) basis sets, including effective core potentials (ECPs), from the EMSL Basis Set Library. It provides an improved user interface and capabilities over its predecessor, the EMSL Basis Set Order Form, for exploring the contents of the EMSL Basis Set Library. The popular Basis Set Order Form and underlying Basis Set Library were originally developed by Dr. David Feller and have been available from the EMSL webpages since 1994. BSE not only allows downloading of the more than 500 Basis sets in various formats; it allows users to annotate existing sets and to upload new sets. (Specialized Interface)
Spackman, Peter R.; Karton, Amir
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}.
NASA Astrophysics Data System (ADS)
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.
Atomic orbital basis sets for use with effective core potentials
Blaudeau, J.P.; Brozell, S.R.; Matsika, S.; Zhang, Z.; Pitzer, R.M.
2000-03-15
Basis sets developed for use with effective core potentials describe pseudo-orbitals rather than orbitals. The primitive Gaussian functions and the contraction coefficients in the basis set must therefore both describe the valence region effectively and allow the pseudo-orbital to be small in the core region. The latter is particularly difficult using 1s primitive functions, which have their maxima at the nucleus. Several methods of choosing contraction coefficients are tried, and it is found that natural orbitals give the best results. The number and optimization of primitive functions are done following Dunning's correlation-consistent procedure. Optimization of orbital exponents for larger atoms frequently results in coalescence of adjacent exponents; use of orbitals with higher principal quantum number is one alternative. Actinide atoms or ions provide the most difficult cases in that basis sets must be optimized for valence shells of different radial size simultaneously considering correlation energy and spin-orbit energy.
High quality Gaussian basis sets for fourth-row atoms
NASA Technical Reports Server (NTRS)
Partridge, Harry; Faegri, Knut, Jr.
1992-01-01
Energy optimized Gaussian basis sets of triple-zeta quality for the atoms Rb-Xe have been derived. Two series of basis sets are developed: (24s 16p 10d) and (26s 16p 10d) sets which were expanded to 13d and 19p functions as the 4d and 5p shells become occupied. For the atoms lighter than Cd, the (24s 16p 10d) sets with triple-zeta valence distributions are higher in energy than the corresponding double-zeta distribution. To ensure a triple-zeta distribution and a global energy minimum, the (26s 16p 10d) sets were derived. Total atomic energies from the largest basis sets are between 198 and 284 (mu)E(sub H) above the numerical Hartree-Fock energies.
Accurate basis set truncation for wavefunction embedding
NASA Astrophysics Data System (ADS)
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.
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
Chopped random-basis quantum optimization
Caneva, Tommaso; Calarco, Tommaso; Montangero, Simone
2011-08-15
In this work, we describe in detail the chopped random basis (CRAB) optimal control technique recently introduced to optimize time-dependent density matrix renormalization group simulations [P. Doria, T. Calarco, and S. Montangero, Phys. Rev. Lett. 106, 190501 (2011)]. Here, we study the efficiency of this control technique in optimizing different quantum processes and we show that in the considered cases we obtain results equivalent to those obtained via different optimal control methods while using less resources. We propose the CRAB optimization as a general and versatile optimal control technique.
Gravitational Lens Modeling with Basis Sets
NASA Astrophysics Data System (ADS)
Birrer, Simon; Amara, Adam; Refregier, Alexandre
2015-11-01
We present a strong lensing modeling technique based on versatile basis sets for the lens and source planes. Our method uses high performance Monte Carlo algorithms, allows for an adaptive build up of complexity, and bridges the gap between parametric and pixel based reconstruction methods. We apply our method to a Hubble Space Telescope image of the strong lens system RX J1131-1231 and show that our method finds a reliable solution and is able to detect substructure in the lens and source planes simultaneously. Using mock data, we show that our method is sensitive to sub-clumps with masses four orders of magnitude smaller than the main lens, which corresponds to about {10}8{M}⊙ , without prior knowledge of the position and mass of the sub-clump. The modeling approach is flexible and maximizes automation to facilitate the analysis of the large number of strong lensing systems expected in upcoming wide field surveys. The resulting search for dark sub-clumps in these systems, without mass-to-light priors, offers promise for probing physics beyond the standard model in the dark matter sector.
Correlation consistent basis sets for the atoms In–Xe
Mahler, Andrew; Wilson, Angela K.
2015-02-28
In this work, the correlation consistent family of Gaussian basis sets has been expanded to include all-electron basis sets for In–Xe. The methodology for developing these basis sets is described, and several examples of the performance and utility of the new sets have been provided. Dissociation energies and bond lengths for both homonuclear and heteronuclear diatomics demonstrate the systematic convergence behavior with respect to increasing basis set quality expected by the family of correlation consistent basis sets in describing molecular properties. Comparison with recently developed correlation consistent sets designed for use with the Douglas-Kroll Hamiltonian is provided.
Entanglement patterns in mutually unbiased basis sets
Lawrence, Jay
2011-08-15
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=p{sup N}), 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 (p{sup N}-p). It follows that the maximum number of product bases is p+1 and, when this number is realized, all remaining (p{sup N}-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.
QUALITY: A program to assess basis set quality
NASA Astrophysics Data System (ADS)
Sordo, J. A.
1998-09-01
A program to analyze in detail the quality of basis sets is presented. The information provided by the application of a wide variety of (atomic and/or molecular) quality criteria is processed by using a methodology that allows one to determine the most appropriate quality test to select a basis set to compute a given (atomic or molecular) property. Fuzzy set theory is used to choose the most adequate basis set to compute simultaneously a set of properties.
Tests for Wavelets as a Basis Set
NASA Astrophysics Data System (ADS)
Baker, Thomas; Evenbly, Glen; White, Steven
A wavelet transformation is a special type of filter usually reserved for image processing and other applications. We develop metrics to evaluate wavelets for general problems on test one-dimensional systems. The goal is to eventually use a wavelet basis in electronic structure calculations. We compare a variety of orthogonal wavelets such as coiflets, symlets, and daubechies wavelets. We also evaluate a new type of orthogonal wavelet with dilation factor three which is both symmetric and compact in real space. This work was supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences under Award #DE-SC008696.
Basis Set Exchange: A Community Database for Computational Sciences
Schuchardt, Karen L.; Didier, Brett T.; Elsethagen, Todd O.; Sun, Lisong; Gurumoorthi, Vidhya; Chase, Jared M.; Li, Jun; Windus, Theresa L.
2007-05-01
Basis sets are one of the most important input data for computational models in the chemistry, materials, biology and other science domains that utilize computational quantum mechanics methods. Providing a shared, web accessible environment where researchers can not only download basis sets in their required format, but browse the data, contribute new basis sets, and ultimately curate and manage the data as a community will facilitate growth of this resource and encourage sharing both data and knowledge. We describe the Basis Set Exchange (BSE), a web portal that provides advanced browsing and download capabilities, facilities for contributing basis set data, and an environment that incorporates tools to foster development and interaction of communities. The BSE leverages and enables continued development of the basis set library originally assembled at the Environmental Molecular Sciences Laboratory.
Near Hartree-Fock quality GTO basis sets for the second-row atoms
NASA Technical Reports Server (NTRS)
Partridge, Harry
1987-01-01
Energy optimized, near Hartree-Fock quality Gaussian basis sets ranging in size from (17s12p) to (20s15p) are presented for the ground states of the second-row atoms for Na(2P), Na(+), Na(-), Mg(3P), P(-), S(-), and Cl(-). In addition, optimized supplementary functions are given for the ground state basis sets to describe the negative ions, and the excited Na(2P) and Mg(3P) atomic states. The ratios of successive orbital exponents describing the inner part of the 1s and 2p orbitals are found to be nearly independent of both nuclear charge and basis set size. This provides a method of obtaining good starting estimates for other basis set optimizations.
Hill, J. Grant E-mail: kipeters@wsu.edu; Peterson, Kirk A. E-mail: kipeters@wsu.edu
2014-09-07
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.
Segmented contracted basis sets for one- and two-component Dirac-Fock effective core potentials.
Weigend, Florian; Baldes, Alexander
2010-11-01
Segmented contracted basis sets for 4d, 5d, 5s, and 6s elements of split (double zeta) valence to quadruple zeta valence quality optimized for Dirac-Fock effective core potentials (ECPs) are presented. They were obtained from previous bases optimized for Wood-Boring ECPs by comparably small modifications and reoptimizations. Additionally extensions for two-component self-consistent-field treatments accounting for spin-orbit (SO) coupling were designed and optimized. Reliability for chemical applications was assessed by comparing results to those obtained with a very large (19s16p17d7f6g) reference basis for a set of more than 80 representatively chosen 5s-5d compounds. Moreover, the effect of different types of ECPs and that of the SO-coupling at the basis set limit of density functional theory is documented for the above set of molecules extended by 40 5p-6p compounds. PMID:21054001
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
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
Localized orbitals from basis sets augmented with diffuse functions
NASA Astrophysics Data System (ADS)
Høyvik, Ida-Marie; Jørgensen, Poul
2013-05-01
Orbital localization of occupied and virtual Hartree-Fock orbitals generated from basis sets augmented with diffuse functions is performed using the Pipek-Mezey, Boys, powers of the second central moment, and powers of the fourth central moment localizations. The locality of the obtained orbital sets are presented in terms of second and fourth moment orbital spreads. The results show that both local occupied and virtual orbitals may be obtained when using powers of the second central moment and powers of the fourth central moment localizations, while the Pipek-Mezey and Boys localizations fail to produce sets of local virtual orbitals. The locality of the fourth central moment virtual orbitals exhibits a locality similar to the locality of a Boys localization for non-augmented basis sets.
NASA Astrophysics Data System (ADS)
Goldey, Matthew; Head-Gordon, Martin
2015-03-01
Second order Møller-Plesset perturbation theory (MP2) in finite basis sets describes several classes of noncovalent interactions poorly due to basis set superposition error (BSSE) and underlying inaccurate physics for dispersion interactions. Attenuation of the Coulomb operator provides a direct path toward improving MP2 for noncovalent interactions. In limited basis sets, we demonstrate improvements in accuracy for intermolecular interactions with a three to five-fold reduction in RMS errors. For a range of inter- and intermolecular test cases, attenuated MP2 even outperforms complete basis set estimates of MP2. Finite basis attenuated MP2 is useful for inter- and intramolecular interactions where higher cost approaches are intractable. Extending this approach, recent research pairs attenuated MP2 with long-range correction to describe potential energy landscapes, and further results for large systems with noncovalent interactions are shown. This work was supported by the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. We acknowledge computational resources obtained under NSF Award CHE-1048789.
Basis-set extensions for two-component spin-orbit treatments of heavy elements.
Armbruster, Markus K; Klopper, Wim; Weigend, Florian
2006-11-14
The accuracy of standard basis sets of quadruple-zeta and lower quality for the use in two-component self-consistent field procedures including spin-orbit coupling is investigated for the elements In-I and Au-At. Spin-orbit coupling leads to energetic and spatial splittings of inner shells, which are not described accurately with standard basis sets optimized for scalar relativistic calculations. This results in large errors in total atomic energies and significant errors in atomization energies of compounds containing these atoms. We show how these errors can be corrected by adding just a few steep sets of basis functions and demonstrate the quality of the resulting extended basis sets. PMID:17066175
Accurate Complete Basis Set Extrapolation of Direct Random Phase Correlation Energies.
Mezei, Pál D; Csonka, Gábor I; Ruzsinszky, Adrienn
2015-08-11
The direct random phase approximation (dRPA) is a promising way to obtain improvements upon the standard semilocal density functional results in many aspects of computational chemistry. In this paper, we address the slow convergence of the calculated dRPA correlation energy with the increase of the quality and size of the popular Gaussian-type Dunning's correlation consistent aug-cc-pVXZ split valence atomic basis set family. The cardinal number X controls the size of the basis set, and we use X = 3-6 in this study. It is known that even the very expensive X = 6 basis sets lead to large errors for the dRPA correlation energy, and thus complete basis set extrapolation is necessary. We study the basis set convergence of the dRPA correlation energies on a set of 65 hydrocarbon isomers from CH4 to C6H6. We calculate the iterative density fitted dRPA correlation energies using an efficient algorithm based on the CC-like form of the equations using the self-consistent HF orbitals. We test the popular inverse cubic, the optimized exponential, and inverse power formulas for complete basis set extrapolation. We have found that the optimized inverse power based extrapolation delivers the best energies. Further analysis showed that the optimal exponent depends on the molecular structure, and the most efficient two-point energy extrapolations that use X = 3 and 4 can be improved considerably by considering the atomic composition and hybridization states of the atoms in the molecules. Our results also show that the optimized exponents that yield accurate X = 3 and 4 extrapolated dRPA energies for atoms or small molecules might be inaccurate for larger molecules. PMID:26574475
NASA Astrophysics Data System (ADS)
Balabin, Roman M.
2011-03-01
The quantum chemistry of conformation equilibrium is a field where great accuracy (better than 100 cal mol-1) is needed because the energy difference between molecular conformers rarely exceeds 1000-3000 cal mol-1. The conformation equilibrium of straight-chain (normal) alkanes is of particular interest and importance for modern chemistry. In this paper, an extra error source for high-quality ab initio (first principles) and DFT calculations of the conformation equilibrium of normal alkanes, namely the intramolecular basis set superposition error (BSSE), is discussed. In contrast to out-of-plane vibrations in benzene molecules, diffuse functions on carbon and hydrogen atoms were found to greatly reduce the relative BSSE of n-alkanes. The corrections due to the intramolecular BSSE were found to be almost identical for the MP2, MP4, and CCSD(T) levels of theory. Their cancelation is expected when CCSD(T)/CBS (CBS, complete basis set) energies are evaluated by addition schemes. For larger normal alkanes (N > 12), the magnitude of the BSSE correction was found to be up to three times larger than the relative stability of the conformer; in this case, the basis set superposition error led to a two orders of magnitude difference in conformer abundance. No error cancelation due to the basis set superposition was found. A comparison with amino acid, peptide, and protein data was provided.
Full Waveform Inversion with Optimal Basis Functions
NASA Astrophysics Data System (ADS)
Sun, Gang; Chang, Qianshun; Sheng, Ping
2003-03-01
Based on the approach suggested by Tarantola, and Gauthier etal., we show that the alternate use of the step (linear) function basis and the block function (quasi-δ function) basis can give accurate full waveform inversion results for the layered acoustic systems, starting from a uniform background. Our method is robust against additive white noise (up to 20% of the signal) and can resolve layers that are comparable to or smaller than a wavelength in thickness. The physical reason for the success of our approach is illustrated through a simple example.
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
Universal perturbative explicitly correlated basis set incompleteness correction
NASA Astrophysics Data System (ADS)
Torheyden, Martin; Valeev, Edward F.
2009-11-01
Basis set incompleteness error for an arbitrary approximate electronic wave function is robustly reduced using a second-order perturbative correction into a basis of explicitly correlated, internally contracted geminal functions. The Hylleraas functional for the second-order energy correction is evaluated algebraically involving at most a four-electron reduced density matrix and four-electron integrals. By using the R12 technology in combination with screening approximations such a correction only requires a two-electron reduced density matrix and two-electron integrals. Preliminary investigations of potential energy surfaces of hydrogen fluoride and nitrogen molecules at the multireference configuration interaction singles and doubles indicate that with the perturbative correction only an aug-cc-pVDZ basis is necessary to compute correlation energies of an aug-cc-pVQZ quality, or better. The proposed correction, dubbed [2]R12, can in principle be combined with any single reference and multireference method in use today.
The Neural Basis of Optimism and Pessimism
2013-01-01
Our survival and wellness require a balance between optimism and pessimism. Undue pessimism makes life miserable; however, excessive optimism can lead to dangerously risky behaviors. A review and synthesis of the literature on the neurophysiology subserving these two worldviews suggests that optimism and pessimism are differentially associated with the two cerebral hemispheres. High self-esteem, a cheerful attitude that tends to look at the positive aspects of a given situation, as well as an optimistic belief in a bright future are associated with physiological activity in the left-hemisphere (LH). In contrast, a gloomy viewpoint, an inclination to focus on the negative part and exaggerate its significance, low self-esteem as well as a pessimistic view on what the future holds are interlinked with neurophysiological processes in the right-hemisphere (RH). This hemispheric asymmetry in mediating optimistic and pessimistic outlooks is rooted in several biological and functional differences between the two hemispheres. The RH mediation of a watchful and inhibitive mode weaves a sense of insecurity that generates and supports pessimistic thought patterns. Conversely, the LH mediation of an active mode and the positive feedback it receives through its motor dexterity breed a sense of confidence in one's ability to manage life's challenges, and optimism about the future. PMID:24167413
Coupled-cluster based basis sets for valence correlation calculations
NASA Astrophysics Data System (ADS)
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
Coupled-cluster based basis sets for valence correlation calculations.
Claudino, Daniel; Gargano, Ricardo; Bartlett, Rodney J
2016-03-14
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 ⟨r(n)⟩ (-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. PMID:26979680
Scar Functions, Barriers for Chemical Reactivity, and Vibrational Basis Sets.
Revuelta, F; Vergini, E; Benito, R M; Borondo, F
2016-07-14
The performance of a recently proposed method to efficiently calculate scar functions is analyzed in problems of chemical interest. An application to the computation of wave functions associated with barriers relevant for the LiNC ⇄ LiCN isomerization reaction is presented as an illustration. These scar functions also constitute excellent elements for basis sets suitable for quantum calculation of vibrational energy levels. To illustrate their efficiency, a calculation of the LiNC/LiCN eigenfunctions is also presented. PMID:26905100
Global and Local Optimization Algorithms for Optimal Signal Set Design
Kearsley, Anthony J.
2001-01-01
The problem of choosing an optimal signal set for non-Gaussian detection was reduced to a smooth inequality constrained mini-max nonlinear programming problem by Gockenbach and Kearsley. Here we consider the application of several optimization algorithms, both global and local, to this problem. The most promising results are obtained when special-purpose sequential quadratic programming (SQP) algorithms are embedded into stochastic global algorithms.
Fuzzy set applications in engineering optimization: Multilevel fuzzy optimization
NASA Technical Reports Server (NTRS)
Diaz, Alejandro R.
1989-01-01
A formulation for multilevel optimization with fuzzy objective functions is presented. With few exceptions, formulations for fuzzy optimization have dealt with a one-level problem in which the objective is the membership function of a fuzzy set formed by the fuzzy intersection of other sets. In the problem examined here, the goal set G is defined in a more general way, using an aggregation operator H that allows arbitrary combinations of set operations (union, intersection, addition) on the individual sets Gi. This is a straightforward extension of the standard form, but one that makes possible the modeling of interesting evaluation strategies. A second, more important departure from the standard form will be the construction of a multilevel problem analogous to the design decomposition problem in optimization. This arrangement facilitates the simulation of a system design process in which different components of the system are designed by different teams, and different levels of design detail become relevant at different time stages in the process: global design features early, local features later in the process.
Near Hartree-Fock quality Gaussian type orbital basis sets for the first- and third-row atoms
NASA Technical Reports Server (NTRS)
Partridge, Harry
1989-01-01
Energy-optimized, near Hartree-Fock (NHF) quality Gaussian type orbital (GTO) basis sets are reported for the second-row (Li to Ne) and fourth-row (K to Kr) atoms. The most accurate basis sets reported for the second row are (18s 13p) sets which are with 4 micro E(H) of the numerical Hartree-Fock (NHF) results. For B to Ne basis sets with more than 15s functions are quadruple zeta in the valence space. For the second-row transition metal atoms the (20s 12p 9d) basis sets are triple zeta in the valence space and are approximately equivalent to Clementi and Roetti's accurate Slater type orbital sets. Supplementing the (20s 12p 9d) basis sets optimized for the lowest state with the 4s(2)3d(n) occupation with a diffuse d function gives self-consistent-field energy separations to the 4s(1)3d(n+1) and 3d(n+2) states which are within 100 micro E(H) of the NHF results. The most accurate basis sets for the transition metal atoms are with 30 micro E(H) of the NHF results. In addition, energy optimized sets are reported for He(3P), Li(2P) and Be(3P).
Optimal Piecewise Linear Basis Functions in Two Dimensions
Brooks III, E D; Szoke, A
2009-01-26
We use a variational approach to optimize the center point coefficients associated with the piecewise linear basis functions introduced by Stone and Adams [1], for polygonal zones in two Cartesian dimensions. Our strategy provides optimal center point coefficients, as a function of the location of the center point, by minimizing the error induced when the basis function interpolation is used for the solution of the time independent diffusion equation within the polygonal zone. By using optimal center point coefficients, one expects to minimize the errors that occur when these basis functions are used to discretize diffusion equations, or transport equations in optically thick zones (where they approach the solution of the diffusion equation). Our optimal center point coefficients satisfy the requirements placed upon the basis functions for any location of the center point. We also find that the location of the center point can be optimized, but this requires numerical calculations. Curiously, the optimum center point location is independent of the values of the dependent variable on the corners only for quadrilaterals.
Basis set limit geometries for ammonia at the SCF and MP2 levels of theory
NASA Technical Reports Server (NTRS)
Defrees, D. J.; Mclean, A. D.
1984-01-01
The controversy over the Hartree-Fock bond angle of NH3 is resolved and the convergence of the geometry for the molecule as the basis set is systematically improved with both SCF and correlated MP2 wave functions. The results of the geometrical optimizations, carried out in four stages with a series of uncontracted bases sets, are shown. The obtained structure for NH3 supports the results of Radom and Rodwell (1980) that the Hartree-Fock limit angle is significantly greater than was previously believed.
Mao, Yuezhi; Horn, Paul R; Mardirossian, Narbe; Head-Gordon, Teresa; Skylaris, Chris-Kriton; Head-Gordon, Martin
2016-07-28
Recently developed density functionals have good accuracy for both thermochemistry (TC) and non-covalent interactions (NC) if very large atomic orbital basis sets are used. To approach the basis set limit with potentially lower computational cost, a new self-consistent field (SCF) scheme is presented that employs minimal adaptive basis (MAB) functions. The MAB functions are optimized on each atomic site by minimizing a surrogate function. High accuracy is obtained by applying a perturbative correction (PC) to the MAB calculation, similar to dual basis approaches. Compared to exact SCF results, using this MAB-SCF (PC) approach with the same large target basis set produces <0.15 kcal/mol root-mean-square deviations for most of the tested TC datasets, and <0.1 kcal/mol for most of the NC datasets. The performance of density functionals near the basis set limit can be even better reproduced. With further improvement to its implementation, MAB-SCF (PC) is a promising lower-cost substitute for conventional large-basis calculations as a method to approach the basis set limit of modern density functionals. PMID:27475350
NASA Astrophysics Data System (ADS)
Mao, Yuezhi; Horn, Paul R.; Mardirossian, Narbe; Head-Gordon, Teresa; Skylaris, Chris-Kriton; Head-Gordon, Martin
2016-07-01
Recently developed density functionals have good accuracy for both thermochemistry (TC) and non-covalent interactions (NC) if very large atomic orbital basis sets are used. To approach the basis set limit with potentially lower computational cost, a new self-consistent field (SCF) scheme is presented that employs minimal adaptive basis (MAB) functions. The MAB functions are optimized on each atomic site by minimizing a surrogate function. High accuracy is obtained by applying a perturbative correction (PC) to the MAB calculation, similar to dual basis approaches. Compared to exact SCF results, using this MAB-SCF (PC) approach with the same large target basis set produces <0.15 kcal/mol root-mean-square deviations for most of the tested TC datasets, and <0.1 kcal/mol for most of the NC datasets. The performance of density functionals near the basis set limit can be even better reproduced. With further improvement to its implementation, MAB-SCF (PC) is a promising lower-cost substitute for conventional large-basis calculations as a method to approach the basis set limit of modern density functionals.
Basis set expansion for inverse problems in plasma diagnostic analysis
NASA Astrophysics Data System (ADS)
Jones, B.; Ruiz, C. L.
2013-07-01
A basis set expansion method [V. Dribinski, A. Ossadtchi, V. A. Mandelshtam, and H. Reisler, Rev. Sci. Instrum. 73, 2634 (2002)], 10.1063/1.1482156 is applied to recover physical information about plasma radiation sources from instrument data, which has been forward transformed due to the nature of the measurement technique. This method provides a general approach for inverse problems, and we discuss two specific examples relevant to diagnosing fast z pinches on the 20-25 MA Z machine [M. E. Savage, L. F. Bennett, D. E. Bliss, W. T. Clark, R. S. Coats, J. M. Elizondo, K. R. LeChien, H. C. Harjes, J. M. Lehr, J. E. Maenchen, D. H. McDaniel, M. F. Pasik, T. D. Pointon, A. C. Owen, D. B. Seidel, D. L. Smith, B. S. Stoltzfus, K. W. Struve, W. A. Stygar, L. K. Warne, J. R. Woodworth, C. W. Mendel, K. R. Prestwich, R. W. Shoup, D. L. Johnson, J. P. Corley, K. C. Hodge, T. C. Wagoner, and P. E. Wakeland, in Proceedings of the Pulsed Power Plasma Sciences Conference (IEEE, 2007), p. 979]. First, Abel inversion of time-gated, self-emission x-ray images from a wire array implosion is studied. Second, we present an approach for unfolding neutron time-of-flight measurements from a deuterium gas puff z pinch to recover information about emission time history and energy distribution. Through these examples, we discuss how noise in the measured data limits the practical resolution of the inversion, and how the method handles discontinuities in the source function and artifacts in the projected image. We add to the method a propagation of errors calculation for estimating uncertainties in the inverted solution.
Basis set expansion for inverse problems in plasma diagnostic analysis.
Jones, B; Ruiz, C L
2013-07-01
A basis set expansion method [V. Dribinski, A. Ossadtchi, V. A. Mandelshtam, and H. Reisler, Rev. Sci. Instrum. 73, 2634 (2002)] is applied to recover physical information about plasma radiation sources from instrument data, which has been forward transformed due to the nature of the measurement technique. This method provides a general approach for inverse problems, and we discuss two specific examples relevant to diagnosing fast z pinches on the 20-25 MA Z machine [M. E. Savage, L. F. Bennett, D. E. Bliss, W. T. Clark, R. S. Coats, J. M. Elizondo, K. R. LeChien, H. C. Harjes, J. M. Lehr, J. E. Maenchen, D. H. McDaniel, M. F. Pasik, T. D. Pointon, A. C. Owen, D. B. Seidel, D. L. Smith, B. S. Stoltzfus, K. W. Struve, W. A. Stygar, L. K. Warne, J. R. Woodworth, C. W. Mendel, K. R. Prestwich, R. W. Shoup, D. L. Johnson, J. P. Corley, K. C. Hodge, T. C. Wagoner, and P. E. Wakeland, in Proceedings of the Pulsed Power Plasma Sciences Conference (IEEE, 2007), p. 979]. First, Abel inversion of time-gated, self-emission x-ray images from a wire array implosion is studied. Second, we present an approach for unfolding neutron time-of-flight measurements from a deuterium gas puff z pinch to recover information about emission time history and energy distribution. Through these examples, we discuss how noise in the measured data limits the practical resolution of the inversion, and how the method handles discontinuities in the source function and artifacts in the projected image. We add to the method a propagation of errors calculation for estimating uncertainties in the inverted solution. PMID:23902066
Basis set expansion for inverse problems in plasma diagnostic analysis
Jones, B.; Ruiz, C. L.
2013-07-15
A basis set expansion method [V. Dribinski, A. Ossadtchi, V. A. Mandelshtam, and H. Reisler, Rev. Sci. Instrum. 73, 2634 (2002)] is applied to recover physical information about plasma radiation sources from instrument data, which has been forward transformed due to the nature of the measurement technique. This method provides a general approach for inverse problems, and we discuss two specific examples relevant to diagnosing fast z pinches on the 20–25 MA Z machine [M. E. Savage, L. F. Bennett, D. E. Bliss, W. T. Clark, R. S. Coats, J. M. Elizondo, K. R. LeChien, H. C. Harjes, J. M. Lehr, J. E. Maenchen, D. H. McDaniel, M. F. Pasik, T. D. Pointon, A. C. Owen, D. B. Seidel, D. L. Smith, B. S. Stoltzfus, K. W. Struve, W. A. Stygar, L. K. Warne, J. R. Woodworth, C. W. Mendel, K. R. Prestwich, R. W. Shoup, D. L. Johnson, J. P. Corley, K. C. Hodge, T. C. Wagoner, and P. E. Wakeland, in Proceedings of the Pulsed Power Plasma Sciences Conference (IEEE, 2007), p. 979]. First, Abel inversion of time-gated, self-emission x-ray images from a wire array implosion is studied. Second, we present an approach for unfolding neutron time-of-flight measurements from a deuterium gas puff z pinch to recover information about emission time history and energy distribution. Through these examples, we discuss how noise in the measured data limits the practical resolution of the inversion, and how the method handles discontinuities in the source function and artifacts in the projected image. We add to the method a propagation of errors calculation for estimating uncertainties in the inverted solution.
Evaluation of Density Functionals and Basis Sets for Carbohydrates
Technology Transfer Automated Retrieval System (TEKTRAN)
Correlated ab initio wave function calculations using MP2/aug-cc-pVTZ model chemistry have been performed for three test sets of gas phase saccharide conformations to provide reference values for their relative energies. The test sets consist of 15 conformers of alpha and beta-D-allopyranose, 15 of ...
NASA Astrophysics Data System (ADS)
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.
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
On the performance of atomic natural orbital basis sets: A full configuration interaction study
Illas, F. Departament de Quimica Fisica, Grup de Quimica Quantica, Facultat de Quimica, Universitat de Barcelona, C Ricart, J.M. ); Rubio, J. ); Bagus, P.S. )
1990-10-01
The performance of atomic natural orbital (ANO) basis sets has been studied by comparing self-consistant field (SCF) and full configuration interaction (CI) results obtained for the first row atoms and hydrides. The ANO results have been compared with those obtained using a segmented basis set containing the same number of contracted basis functions. The total energies obtained with the ANO basis sets are always lower than the one obtained by using the segmented one. However, for the hydrides, differential electronic correlation energy obtained with the ANO basis set may be smaller than the one recovered with the segmented set. We relate this poorer differential correlation energy for the ANO basis set to the fact that only one contracted {ital d} function is used for the ANO and segmented basis sets.
Is there an optimal basis to maximise optical information transfer?
NASA Astrophysics Data System (ADS)
Chen, Mingzhou; Dholakia, Kishan; Mazilu, Michael
2016-03-01
We establish the concept of the density of the optical degrees of freedom that may be applied to any photonics based system. As a key example of this versatile approach we explore information transfer using optical communication. We demonstrate both experimentally, theoretically and numerically that the use of a basis set with fields containing optical vortices does not increase the telecommunication capacity of an optical system.
Is there an optimal basis to maximise optical information transfer?
Chen, Mingzhou; Dholakia, Kishan; Mazilu, Michael
2016-01-01
We establish the concept of the density of the optical degrees of freedom that may be applied to any photonics based system. As a key example of this versatile approach we explore information transfer using optical communication. We demonstrate both experimentally, theoretically and numerically that the use of a basis set with fields containing optical vortices does not increase the telecommunication capacity of an optical system. PMID:26976626
Time Domain Propagation of Quantum and Classical Systems using a Wavelet Basis Set Method
NASA Astrophysics Data System (ADS)
Lombardini, Richard; Nowara, Ewa; Johnson, Bruce
2015-03-01
The use of an orthogonal wavelet basis set (Optimized Maximum-N Generalized Coiflets) to effectively model physical systems in the time domain, in particular the electromagnetic (EM) pulse and quantum mechanical (QM) wavefunction, is examined in this work. Although past research has demonstrated the benefits of wavelet basis sets to handle computationally expensive problems due to their multiresolution properties, the overlapping supports of neighboring wavelet basis functions poses problems when dealing with boundary conditions, especially with material interfaces in the EM case. Specifically, this talk addresses this issue using the idea of derivative matching creating fictitious grid points (T.A. Driscoll and B. Fornberg), but replaces the latter element with fictitious wavelet projections in conjunction with wavelet reconstruction filters. Two-dimensional (2D) systems are analyzed, EM pulse incident on silver cylinders and the QM electron wave packet circling the proton in a hydrogen atom system (reduced to 2D), and the new wavelet method is compared to the popular finite-difference time-domain technique.
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. PMID:27394093
Optimal Assembly of Tests with Item Sets.
ERIC Educational Resources Information Center
van der Linden, Wim J.
2000-01-01
Presents six computational methods based on mixed-integer programming for assembling tests from a bank with an item-set structure and evaluated these methods using mathematical programming feasibiity and expected solution times. Illustrates these methods with two data sets from the Law School Admission Test. Discusses the best approximations to…
Berry, Rachel H; Hobson, Michael P; Withington, Stafford
2004-05-01
We discuss a general theoretical framework for representing and propagating fully coherent, fully incoherent, and the intermediate regime of partially coherent submillimeter-wave fields by means of general sampled basis functions, which may have any degree of completeness. Partially coherent fields arise when finite-throughput systems induce coherence on incoherent fields. This powerful extension to traditional modal analysis methods by using undercomplete Gaussian-Hermite modes can be employed to analyze and optimize such Gaussian quasi-optical techniques. We focus on one particular basis set, the Gabor basis, which consists of overlapping translated and modulated Gaussian beams. We present high-accuracy numerical results from field reconstructions and propagations. In particular, we perform one-dimensional analyses illustrating the Van Cittert-Zernike theorem and then extend our simulations to two dimensions, including simple models of horn and bolometer arrays. Our methods and results are of practical importance as a method for analyzing terahertz fields, which are often partially coherent and diffraction limited so that ray tracing is inaccurate and physical optics computationally prohibitive. PMID:15139431
Atomic self-consistent-field program by the basis set expansion method: Columbus version
NASA Astrophysics Data System (ADS)
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
The study of basis sets for the calculation of the structure and dynamics of the benzene-Kr complex.
Shirkov, Leonid; Makarewicz, Jan
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. PMID:26026434
The study of basis sets for the calculation of the structure and dynamics of the benzene-Kr complex
Shirkov, Leonid; Makarewicz, Jan
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.
The study of basis sets for the calculation of the structure and dynamics of the benzene-Kr complex
NASA Astrophysics Data System (ADS)
Shirkov, Leonid; Makarewicz, Jan
2015-05-01
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.
Miliordos, Evangelos; Xantheas, Sotiris S.
2015-03-07
We report the variation of the binding energy of the Formic Acid Dimer with the size of the basis set at the Coupled Cluster with iterative Singles, Doubles and perturbatively connected Triple replacements [CCSD(T)] level of theory, estimate the Complete Basis Set (CBS) limit, and examine the validity of the Basis Set Superposition Error (BSSE)-correction for this quantity that was previously challenged by Kalescky, Kraka, and Cremer (KKC) [J. Chem. Phys. 140, 084315 (2014)]. Our results indicate that the BSSE correction, including terms that account for the substantial geometry change of the monomers due to the formation of two strong hydrogen bonds in the dimer, is indeed valid for obtaining accurate estimates for the binding energy of this system as it exhibits the expected decrease with increasing basis set size. We attribute the discrepancy between our current results and those of KKC to their use of a valence basis set in conjunction with the correlation of all electrons (i.e., including the 1s of C and O). We further show that the use of a core-valence set in conjunction with all electron correlation converges faster to the CBS limit as the BSSE correction is less than half than the valence electron/valence basis set case. The uncorrected and BSSE-corrected binding energies were found to produce the same (within 0.1 kcal/mol) CBS limits. We obtain CCSD(T)/CBS best estimates for D{sub e} = − 16.1 ± 0.1 kcal/mol and for D{sub 0} = − 14.3 ± 0.1 kcal/mol, the later in excellent agreement with the experimental value of −14.22 ± 0.12 kcal/mol.
Bi-directional evolutionary level set method for topology optimization
NASA Astrophysics Data System (ADS)
Zhu, Benliang; Zhang, Xianmin; Fatikow, Sergej; Wang, Nianfeng
2015-03-01
A bi-directional evolutionary level set method for solving topology optimization problems is presented in this article. The proposed method has three main advantages over the standard level set method. First, new holes can be automatically generated in the design domain during the optimization process. Second, the dependency of the obtained optimized configurations upon the initial configurations is eliminated. Optimized configurations can be obtained even being started from a minimum possible initial guess. Third, the method can be easily implemented and is computationally more efficient. The validity of the proposed method is tested on the mean compliance minimization problem and the compliant mechanisms topology optimization problem.
Derivation of a formula for the resonance integral for a nonorthogonal basis set
Yim, Yung-Chang; Eyring, Henry
1981-01-01
In a self-consistent field calculation, a formula for the off-diagonal matrix elements of the core Hamiltonian is derived for a nonorthogonal basis set by a polyatomic approach. A set of parameters is then introduced for the repulsion integral formula of Mataga-Nishimoto to fit the experimental data. The matrix elements computed for the nonorthogonal basis set in the π-electron approximation are transformed to those for an orthogonal basis set by the Löwdin symmetrical orthogonalization. PMID:16593009
A Hartree-Fock study of the confined helium atom: Local and global basis set approaches
NASA Astrophysics Data System (ADS)
Young, Toby D.; Vargas, Rubicelia; Garza, Jorge
2016-02-01
Two different basis set methods are used to calculate atomic energy within Hartree-Fock theory. The first is a local basis set approach using high-order real-space finite elements and the second is a global basis set approach using modified Slater-type orbitals. These two approaches are applied to the confined helium atom and are compared by calculating one- and two-electron contributions to the total energy. As a measure of the quality of the electron density, the cusp condition is analyzed.
Dynamical basis sets for algebraic variational calculations in quantum-mechanical scattering theory
NASA Technical Reports Server (NTRS)
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.
Pseudospectral sampling of Gaussian basis sets as a new avenue to high-dimensional quantum dynamics
NASA Astrophysics Data System (ADS)
Heaps, Charles
This thesis presents a novel approach to modeling quantum molecular dynamics (QMD). Theoretical approaches to QMD are essential to understanding and predicting chemical reactivity and spectroscopy. We implement a method based on a trajectory-guided basis set. In this case, the nuclei are propagated in time using classical mechanics. Each nuclear configuration corresponds to a basis function in the quantum mechanical expansion. Using the time-dependent configurations as a basis set, we are able to evolve in time using relatively little information at each time step. We use a basis set of moving frozen (time-independent width) Gaussian functions that are well-known to provide a simple and efficient basis set for nuclear dynamics. We introduce a new perspective to trajectory-guided Gaussian basis sets based on existing numerical methods. The distinction is based on the Galerkin and collocation methods. In the former, the basis set is tested using basis functions, projecting the solution onto the functional space of the problem and requiring integration over all space. In the collocation method, the Dirac delta function tests the basis set, projecting the solution onto discrete points in space. This effectively reduces the integral evaluation to function evaluation, a fundamental characteristic of pseudospectral methods. We adopt this idea for independent trajectory-guided Gaussian basis functions. We investigate a series of anharmonic vibrational models describing dynamics in up to six dimensions. The pseudospectral sampling is found to be as accurate as full integral evaluation, while the former method is fully general and integration is only possible on very particular model potential energy surfaces. Nonadiabatic dynamics are also investigated in models of photodissociation and collinear triatomic vibronic coupling. Using Ehrenfest trajectories to guide the basis set on multiple surfaces, we observe convergence to exact results using hundreds of basis functions
Adapting DFT+U for the Chemically Motivated Correction of Minimal Basis Set Incompleteness.
Kulik, Heather J; Seelam, Natasha; Mar, Brendan D; Martínez, Todd J
2016-07-28
Recent algorithmic and hardware advances have enabled the application of electronic structure methods to the study of large-scale systems such as proteins with O(10(3)) atoms. Most such methods benefit greatly from the use of reduced basis sets to further enhance their speed, but truly minimal basis sets are well-known to suffer from incompleteness error that gives rise to incorrect descriptions of chemical bonding, preventing minimal basis set use in production calculations. We present a strategy for improving these well-known shortcomings in minimal basis sets by selectively tuning the energetics and bonding of nitrogen and oxygen atoms within proteins and small molecules to reproduce polarized double-ζ basis set geometries at minimal basis set cost. We borrow the well-known +U correction from the density functional theory community normally employed for self-interaction errors and demonstrate its power in the context of correcting basis set incompleteness within a formally self-interaction-free Hartree-Fock framework. We tune the Hubbard U parameters for nitrogen and oxygen atoms on small-molecule tautomers (e.g., cytosine), demonstrate the applicability of the approach on a number of amide-containing molecules (e.g., formamide, alanine tripeptide), and test our strategy on a 10 protein test set where anomalous proton transfer events are reduced by 90% from RHF/STO-3G to RHF/STO-3G+U, bringing the latter into quantitative agreement with RHF/6-31G* results. Although developed with the study of biological molecules in mind, this empirically tuned U approach shows promise as an alternative strategy for correction of basis set incompleteness errors. PMID:27383567
What is the most efficient way to reach the canonical MP2 basis set limit?
NASA Astrophysics Data System (ADS)
Liakos, Dimitrios G.; Izsák, Róbert; Valeev, Edward F.; Neese, Frank
2013-09-01
Various ways of reaching the complete basis set limit at the second-order Møller-Plesset perturbation theory (MP2) level are compared with respect to their cost-to-accuracy ratio. These include: (1) traditional MP2 calculations with correlation consistent basis sets of increasing size, with and without the resolution of identity for Coulomb and exchange (RIJK) or the combined RIJ and 'chain of spheres' (RIJCOSX) approximations; (2) basis set extrapolation obtained with the same MP2 variants; and (3) explicitly correlated F12-MP2 methods. The time required to solve the Hartree-Fock equations is part of the evaluation because the overall efficiency is of central interest in this work. Results were obtained for the ISO34, DC9 and S66 test sets and were analysed in terms of efficiency and accuracy for total energies, reaction energies and their effect on the basis set superposition error. Among the methods studied, the RIJK-MP2-F12 and RIJK-MP2-EP1 (where EP1 stands for 'Extrapolation Protocol 1' as explained in the text) methods perform outstandingly well. Although extrapolation is, in general, slightly faster than explicit correlation, it is found that for reaction energies, RIJK-MP2-F12 performs systematically better. This holds especially in combination with a triple zeta basis set, in which case it even outperforms the much more costly extrapolation involving quadruple- and quintuple-zeta correlation consistent basis sets.
On basis set superposition error corrected stabilization energies for large n-body clusters.
Walczak, Katarzyna; Friedrich, Joachim; Dolg, Michael
2011-10-01
In this contribution, we propose an approximate basis set superposition error (BSSE) correction scheme for the site-site function counterpoise and for the Valiron-Mayer function counterpoise correction of second order to account for the basis set superposition error in clusters with a large number of subunits. The accuracy of the proposed scheme has been investigated for a water cluster series at the CCSD(T), CCSD, MP2, and self-consistent field levels of theory using Dunning's correlation consistent basis sets. The BSSE corrected stabilization energies for a series of water clusters are presented. A study regarding the possible savings with respect to computational resources has been carried out as well as a monitoring of the basis set dependence of the approximate BSSE corrections. PMID:21992293
Peterson, Kirk A.; Figgen, Detlev; Goll, Erich; Stoll, Hermann; Dolg, Michael F.
2003-12-01
Series of correlation consistent basis sets have been developed for the post-d group 16-18 elements in conjunction with small-core relativistic pseudopotentials (PPs) of the energy-consistent variety. The latter were adjusted to multiconfiguration Dirac-Hartree-Fock data based on the Dirac-Coulomb-Breit Hamiltonian. The outer-core (n-1)spd shells are explicitly treated together with the nsp valence shell with these PPs. The accompanying cc-pVnZ-PP and aug-cc-pVnZ-PP basis sets range in size from DZ to 5Z quality and yield systematic convergence of both Hartree-Fock and correlated total energies. In addition to the calculation of atomic electron affinities and dipole polarizabilities of the rare gas atoms, numerous molecular benchmark calculations (HBr, HI, HAt, Br2, I2, At2, SiSe, SiTe, SiPo, KrH+, XeH+, and RnH+) are also reported at the coupled cluster level of theory. For the purposes of comparison, all-electron calculations using the Douglas-Kroll-Hess Hamiltonian have also been carried out for the halogen-containing molecules using basis sets of 5Z quality.
NASA Astrophysics Data System (ADS)
Kupka, Teobald
2008-08-01
Based on B3LYP spin-spin coupling constants (SSCC) of several molecules calculated with cc-pV xZ, cc-pCV xZ, cc-pCV xZ-sd and cc-pCV xZ-sd+ t basis sets, a reasonably fit, using the two-parameter formula, to the Kohn-Sham complete basis set limit (CBS) is shown. Improvement in the CBS values going from cc-pV xZ to the most elaborated cc-pCV xZ-sd+ t basis set family is observed: standard deviation for all data drops from 33.7 to 23.1, and from 6.0 to 4.8 Hz after excluding problematic 1J(F,H) and 1J(F,C). Calculation of water's 1J(OH) using B3LYP/cc-pCV xZ and B3LYP/pcJ- n significantly improved the FC term convergence.
Kollmar, Christian Neese, Frank
2014-10-07
The role of the static Kohn-Sham (KS) response function describing the response of the electron density to a change of the local KS potential is discussed in both the theory of the optimized effective potential (OEP) and the so-called inverse Kohn-Sham problem involving the task to find the local KS potential for a given electron density. In a general discussion of the integral equation to be solved in both cases, it is argued that a unique solution of this equation can be found even in case of finite atomic orbital basis sets. It is shown how a matrix representation of the response function can be obtained if the exchange-correlation potential is expanded in terms of a Schmidt-orthogonalized basis comprising orbitals products of occupied and virtual orbitals. The viability of this approach in both OEP theory and the inverse KS problem is illustrated by numerical examples.
An active set algorithm for nonlinear optimization with polyhedral constraints
NASA Astrophysics Data System (ADS)
Hager, William W.; Zhang, Hongchao
2016-08-01
A polyhedral active set algorithm PASA is developed for solving a nonlinear optimization problem whose feasible set is a polyhedron. Phase one of the algorithm is the gradient projection method, while phase two is any algorithm for solving a linearly constrained optimization problem. Rules are provided for branching between the two phases. Global convergence to a stationary point is established, while asymptotically PASA performs only phase two when either a nondegeneracy assumption holds, or the active constraints are linearly independent and a strong second-order sufficient optimality condition holds.
NASA Astrophysics Data System (ADS)
Booth, George H.; Cleland, Deidre; Alavi, Ali; Tew, David P.
2012-10-01
By performing a stochastic dynamic in a space of Slater determinants, the full configuration interaction quantum Monte Carlo (FCIQMC) method has been able to obtain energies which are essentially free from systematic error to the basis set correlation energy, within small and systematically improvable error bars. However, the weakly exponential scaling with basis size makes converging the energy with respect to basis set costly and in larger systems, impossible. To ameliorate these basis set issues, here we use perturbation theory to couple the FCIQMC wavefunction to an explicitly correlated strongly orthogonal basis of geminals, following the { [2]_{{R12}} } approach of Valeev et al. The required one- and two-particle density matrices are computed on-the-fly during the FCIQMC dynamic, using a sampling procedure which incurs relatively little additional computation expense. The F12 energy corrections are shown to converge rapidly as a function of sampling, both in imaginary time and number of walkers. Our pilot calculations on the binding curve for the carbon dimer, which exhibits strong correlation effects as well as substantial basis set dependence, demonstrate that the accuracy of the FCIQMC-F12 method surpasses that of all previous FCIQMC calculations, and that the F12 correction improves results equivalent to increasing the quality of the one-electron basis by two cardinal numbers.
A novel Gaussian-Sinc mixed basis set for electronic structure calculations.
Jerke, Jonathan L; Lee, Young; Tymczak, C J
2015-08-14
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. PMID:26277128
A novel Gaussian-Sinc mixed basis set for electronic structure calculations
Jerke, Jonathan L.; Lee, Young; Tymczak, C. J.
2015-08-14
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 H{sub 2}, O{sub 2}, and N{sub 2}, 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.
Magnetic properties with multiwavelets and DFT: the complete basis set limit achieved.
Jensen, Stig Rune; Flå, Tor; Jonsson, Dan; Monstad, Rune Sørland; Ruud, Kenneth; Frediani, Luca
2016-08-01
Multiwavelets are emerging as an attractive alternative to traditional basis sets such as Gaussian-type orbitals and plane waves. One of their distinctive properties is the ability to reach the basis set limit (often a chimera for traditional approaches) reliably and consistently by fixing the desired precision ε. We present our multiwavelet implementation of the linear response formalism, applied to static magnetic properties, at the self-consistent field level of theory (both for Hartree-Fock and density functional theories). We demonstrate that the multiwavelets consistently improve the accuracy of the results when increasing the desired precision, yielding results that have four to five digits precision, thus providing a very useful benchmark which could otherwise only be estimated by extrapolation methods. Our results show that magnetizabilities obtained with the augmented quadruple-ζ basis (aug-cc-pCVQZ) are practically at the basis set limit, whereas absolute nuclear magnetic resonance shielding tensors are more challenging: even by making use of a standard extrapolation method, the accuracy is not substantially improved. In contrast, our results provide a benchmark that: (1) confirms the validity of the extrapolation ansatz; (2) can be used as a reference to achieve a property-specific extrapolation scheme, thus providing a means to obtain much better extrapolated results; (3) allows us to separate functional-specific errors from basis-set ones and thus to assess the level of cancellation between basis set and functional errors often exploited in density functional theory. PMID:27087397
Atomization Energies of SO and SO2; Basis Set Extrapolation Revisted
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Ricca, Alessandra; Arnold, James (Technical Monitor)
1998-01-01
The addition of tight functions to sulphur and extrapolation to the complete basis set limit are required to obtain accurate atomization energies. Six different extrapolation procedures are tried. The best atomization energies come from the series of basis sets that yield the most consistent results for all extrapolation techniques. In the variable alpha approach, alpha values larger than 4.5 or smaller than 3, appear to suggest that the extrapolation may not be reliable. It does not appear possible to determine a reliable basis set series using only the triple and quadruple zeta based sets. The scalar relativistic effects reduce the atomization of SO and SO2 by 0.34 and 0.81 kcal/mol, respectively, and clearly must be accounted for if a highly accurate atomization energy is to be computed. The magnitude of the core-valence (CV) contribution to the atomization is affected by missing diffuse valence functions. The CV contribution is much more stable if basis set superposition errors are accounted for. A similar study of SF, SF(+), and SF6 shows that the best family of basis sets varies with the nature of the S bonding.
Systematic Study of Locally Dense Basis Sets for NMR Shielding Constants.
Reid, David M; Kobayashi, Rika; Collins, Michael A
2014-01-14
This paper presents a systematic study of partitioning schemes for locally dense basis sets in the context of NMR shielding calculations. The partitionings explored were based exclusively on connectivity and utilized the basis sets from the pcS-n series. Deviations from pcS-4 shieldings were calculated for a set of 28 organic molecules at the HF, B3LYP, and KT3 levels of theory, with the primary goal being the determination of an efficient scheme that achieves maximal deviations of 0.1 ppm for (1)H and 1 ppm for (13)C. Both atom based and group based divisions of basis sets were examined, with the latter providing the most promising results. It is demonstrated that for the systems studied, at least pcS-1 is required for all parts of the molecule. This, coupled with pcS-3 on the group of interest and pcS-2 on the adjacent groups, is sufficient to achieve the desired level of accuracy at a minimal computational expense. In addition, the suitability of the pcS-n basis sets for post-SCF methods was confirmed through a comparison with other standard basis sets at the MP2 level. PMID:26579898
Holden, Zachary C.; Richard, Ryan M.; Herbert, John M.
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.
NASA Astrophysics Data System (ADS)
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.
NASA Astrophysics Data System (ADS)
Park, R.; Jo, D.; Kim, M.; Spracklen, D. V.; Hodzic, A.
2014-12-01
Organic aerosol (OA) constitutes significant mass fractions (20-90%) of total dry fine aerosols in the atmosphere. However, global models of OA have shown large discrepancies when compared to the observations because of the limited capability to simulate secondary OA (SOA). For reducing the discrepancies between observations and models, recent studies have shown that chemical aging reactions in the atmosphere are important because they can lead to decreases in organic volatility, resulting in increase of SOA mass yields. To efficiently simulate chemical aging of SOA in the atmosphere, we implemented the volatility basis set approach in a global 3-D chemical transport model (GEOS-Chem). We present full-year simulations and their comparisons with multiple observations - global aerosol mass spectrometer dataset, the Interagency Monitoring of Protected Visual Environments from the United States, the European Monitoring and Evaluation Programme dataset and water-soluble organic carbon observation data collected over East Asia. Using different input parameters in the model, we also explore the uncertainty of the SOA simulation for which we use an observational constraint to find the optimized values with which the model reduces the discrepancy from the observations. Finally, we estimate the effect of OA on climate using our best simulation results.
Efficient calculation of integrals in mixed ramp-Gaussian basis sets
McKemmish, Laura K.
2015-04-07
Algorithms for the efficient calculation of two-electron integrals in the newly developed mixed ramp-Gaussian basis sets are presented, alongside a Fortran90 implementation of these algorithms, RAMPITUP. These new basis sets have significant potential to (1) give some speed-up (estimated at up to 20% for large molecules in fully optimised code) to general-purpose Hartree-Fock (HF) and density functional theory quantum chemistry calculations, replacing all-Gaussian basis sets, and (2) give very large speed-ups for calculations of core-dependent properties, such as electron density at the nucleus, NMR parameters, relativistic corrections, and total energies, replacing the current use of Slater basis functions or very large specialised all-Gaussian basis sets for these purposes. This initial implementation already demonstrates roughly 10% speed-ups in HF/R-31G calculations compared to HF/6-31G calculations for large linear molecules, demonstrating the promise of this methodology, particularly for the second application. As well as the reduction in the total primitive number in R-31G compared to 6-31G, this timing advantage can be attributed to the significant reduction in the number of mathematically complex intermediate integrals after modelling each ramp-Gaussian basis-function-pair as a sum of ramps on a single atomic centre.
Representability of Bloch states on Projector-augmented-wave (PAW) basis sets
NASA Astrophysics Data System (ADS)
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.
On the performance of large Gaussian basis sets for the computation of total atomization energies
NASA Technical Reports Server (NTRS)
Martin, J. M. L.
1992-01-01
The total atomization energies of a number of molecules have been computed using an augmented coupled-cluster method and (5s4p3d2f1g) and 4s3p2d1f) atomic natural orbital (ANO) basis sets, as well as the correlation consistent valence triple zeta plus polarization (cc-pVTZ) correlation consistent valence quadrupole zeta plus polarization (cc-pVQZ) basis sets. The performance of ANO and correlation consistent basis sets is comparable throughout, although the latter can result in significant CPU time savings. Whereas the inclusion of g functions has significant effects on the computed Sigma D(e) values, chemical accuracy is still not reached for molecules involving multiple bonds. A Gaussian-1 (G) type correction lowers the error, but not much beyond the accuracy of the G1 model itself. Using separate corrections for sigma bonds, pi bonds, and valence pairs brings down the mean absolute error to less than 1 kcal/mol for the spdf basis sets, and about 0.5 kcal/mol for the spdfg basis sets. Some conclusions on the success of the Gaussian-1 and Gaussian-2 models are drawn.
Macedo, Luiz Guilherme M de Borin, Antonio Carlos; Silva, Alberico B.F. da
2007-11-15
Prolapse-free basis sets suitable for four-component relativistic quantum chemical calculations are presented for the superheavy elements up to {sub 118}Uuo ({sub 104}Rf, {sub 105}Db, {sub 106}Sg, {sub 107}Bh, {sub 108}Hs, {sub 109}Mt, {sub 110}Ds, {sub 111}Rg, {sub 112}Uub, {sub 113}Uut, {sub 114}Uuq, {sub 115}Uup, {sub 116}Uuh, {sub 117}Uus, {sub 118}Uuo) and {sub 103}Lr. 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.
Energy-optimal path planning by stochastic dynamically orthogonal level-set optimization
NASA Astrophysics Data System (ADS)
Subramani, Deepak N.; Lermusiaux, Pierre F. J.
2016-04-01
A stochastic optimization methodology is formulated for computing energy-optimal paths from among time-optimal paths of autonomous vehicles navigating in a dynamic flow field. Based on partial differential equations, the methodology rigorously leverages the level-set equation that governs time-optimal reachability fronts for a given relative vehicle-speed function. To set up the energy optimization, the relative vehicle-speed and headings are considered to be stochastic and new stochastic Dynamically Orthogonal (DO) level-set equations are derived. Their solution provides the distribution of time-optimal reachability fronts and corresponding distribution of time-optimal paths. An optimization is then performed on the vehicle's energy-time joint distribution to select the energy-optimal paths for each arrival time, among all stochastic time-optimal paths for that arrival time. Numerical schemes to solve the reduced stochastic DO level-set equations are obtained, and accuracy and efficiency considerations are discussed. These reduced equations are first shown to be efficient at solving the governing stochastic level-sets, in part by comparisons with direct Monte Carlo simulations. To validate the methodology and illustrate its accuracy, comparisons with semi-analytical energy-optimal path solutions are then completed. In particular, we consider the energy-optimal crossing of a canonical steady front and set up its semi-analytical solution using a energy-time nested nonlinear double-optimization scheme. We then showcase the inner workings and nuances of the energy-optimal path planning, considering different mission scenarios. Finally, we study and discuss results of energy-optimal missions in a wind-driven barotropic quasi-geostrophic double-gyre ocean circulation.
Woon, D.E.; Dunning, T.H. Jr. )
1994-02-15
An accurate description of the electrical properties of atoms and molecules is critical for quantitative predictions of the nonlinear properties of molecules and of long-range atomic and molecular interactions between both neutral and charged species. We report a systematic study of the basis sets required to obtain accurate correlated values for the static dipole ([alpha][sub 1]), quadrupole ([alpha][sub 2]), and octopole ([alpha][sub 3]) polarizabilities and the hyperpolarizability ([gamma]) of the rare gas atoms He, Ne, and Ar. Several methods of correlation treatment were examined, including various orders of Moller--Plesset perturbation theory (MP2, MP3, MP4), coupled-cluster theory with and without perturbative treatment of triple excitations [CCSD, CCSD(T)], and singles and doubles configuration interaction (CISD). All of the basis sets considered here were constructed by adding even-tempered sets of diffuse functions to the correlation consistent basis sets of Dunning and co-workers. With multiply-augmented sets we find that the electrical properties of the rare gas atoms converge smoothly to values that are in excellent agreement with the available experimental data and/or previously computed results. As a further test of the basis sets presented here, the dipole polarizabilities of the F[sup [minus
Influence of basis sets and electron correlation on theoretically predicted infrared intensities
Miller, M.D. ); Jensen, F. ); Chapman, O.L.; Houk, K.N. )
1989-06-01
A systematic study of the effects of basis sets and electron correlation on calculated infrared intensities has been performed with ab initio molecular orbital calculations and Moeller-Plesset perturbation theory. Absolute IR intensities of hydrogen fluoride, hydroxy radical, carbon monoxide, hydrogen cyanide, and formaldehyde have been calculated with basis sets ranging from 3-21G to 6-311++G(2dd{prime},2pp{prime}) and with electron correlation corrections up through MP4(SDTQ). A basis set with polarization and diffuse functions is necessary to obtain reasonably accurate intensities. Electron correlation significantly improves the agreement between experimental and calculated values. Except for carbon monoxide, the intensities calculated at the MP4 level compare favorably with experimental intensities, the errors being less than the measured difference between those obtained from inert-gas matrices at low temperature and those reported for the gas phase.
On the basis set convergence of electron-electron entanglement measures: helium-like systems.
Hofer, Thomas S
2013-01-01
A systematic investigation of three different electron-electron entanglement measures, namely the von Neumann, the linear and the occupation number entropy at full configuration interaction level has been performed for the four helium-like systems hydride, helium, Li(+) and Be(2+) using a large number of different basis sets. The convergence behavior of the resulting energies and entropies revealed that the latter do in general not show the expected strictly monotonic increase upon increase of the one-electron basis. Overall, the three different entanglement measures show good agreement among each other, the largest deviations being observed for small basis sets. The data clearly demonstrates that it is important to consider the nature of the chemical system when investigating entanglement phenomena in the framework of Gaussian type basis sets: while in case of hydride the use of augmentation functions is crucial, the application of core functions greatly improves the accuracy in case of cationic systems such as Li(+) and Be(2+). In addition, numerical derivatives of the entanglement measures with respect to the nucleic charge have been determined, which proved to be a very sensitive probe of the convergence leading to qualitatively wrong results (i.e., the wrong sign) if too small basis sets are used. PMID:24790952
On the basis set convergence of electron–electron entanglement measures: helium-like systems
Hofer, Thomas S.
2013-01-01
A systematic investigation of three different electron–electron entanglement measures, namely the von Neumann, the linear and the occupation number entropy at full configuration interaction level has been performed for the four helium-like systems hydride, helium, Li+ and Be2+ using a large number of different basis sets. The convergence behavior of the resulting energies and entropies revealed that the latter do in general not show the expected strictly monotonic increase upon increase of the one–electron basis. Overall, the three different entanglement measures show good agreement among each other, the largest deviations being observed for small basis sets. The data clearly demonstrates that it is important to consider the nature of the chemical system when investigating entanglement phenomena in the framework of Gaussian type basis sets: while in case of hydride the use of augmentation functions is crucial, the application of core functions greatly improves the accuracy in case of cationic systems such as Li+ and Be2+. In addition, numerical derivatives of the entanglement measures with respect to the nucleic charge have been determined, which proved to be a very sensitive probe of the convergence leading to qualitatively wrong results (i.e., the wrong sign) if too small basis sets are used. PMID:24790952
Discrete Optimization with Polynomially Detectable Boundaries and Restricted Level Sets
NASA Astrophysics Data System (ADS)
Zinder, Yakov; Memar, Julia; Singh, Gaurav
The paper describes an optimization procedure for a class of discrete optimization problems which is defined by certain properties of the boundary of the feasible region and level sets of the objective function. It is shown that these properties are possessed, for example, by various scheduling problems, including a number of well-known NP-hard problems which play an important role in scheduling theory. For an important particular case the presented optimization procedure is compared with a version of the branch-and-bound algorithm by means of computational experiments.
NASA Astrophysics Data System (ADS)
Huang, Jingsong; Kertesz, Miklos
2004-05-01
Intermolecular transfer integrals, and associated band-structures of organic molecular materials can be calculated through a dimer approach. Extensive numerical studies are performed on an ethylene π-dimer to investigate the basis sets dependence of transfer integrals. Convergence of calculated transfer integrals is achieved with respect to both Gaussian and plane-wave basis sets, provided the same level of theory is used. Effects of diffuse and polarization Gaussian functions on transfer integrals are identified. Comparison of experimental and theoretical values of transfer integrals of the TTF-TCNQ charge transfer salt is also presented.
Strategies for reducing basis set superposition error (BSSE) in O/AU and O/Ni
NASA Astrophysics Data System (ADS)
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.
NASA Astrophysics Data System (ADS)
Boffi, Nicholas M.; Jain, Manish; Natan, Amir
2016-02-01
A real-space high order finite difference method is used to analyze the effect of spherical domain size on the Hartree-Fock (and density functional theory) virtual eigenstates. We show the domain size dependence of both positive and negative virtual eigenvalues of the Hartree-Fock equations for small molecules. We demonstrate that positive states behave like a particle in spherical well and show how they approach zero. For the negative eigenstates, we show that large domains are needed to get the correct eigenvalues. We compare our results to those of Gaussian basis sets and draw some conclusions for real-space, basis-sets, and plane-waves calculations.
Predicting Pt-195 NMR chemical shift using new relativistic all-electron basis set.
Paschoal, D; Guerra, C Fonseca; de Oliveira, M A L; Ramalho, T C; Dos Santos, H F
2016-10-01
Predicting NMR properties is a valuable tool to assist the experimentalists in the characterization of molecular structure. For heavy metals, such as Pt-195, only a few computational protocols are available. In the present contribution, all-electron Gaussian basis sets, suitable to calculate the Pt-195 NMR chemical shift, are presented for Pt and all elements commonly found as Pt-ligands. The new basis sets identified as NMR-DKH were partially contracted as a triple-zeta doubly polarized scheme with all coefficients obtained from a Douglas-Kroll-Hess (DKH) second-order scalar relativistic calculation. The Pt-195 chemical shift was predicted through empirical models fitted to reproduce experimental data for a set of 183 Pt(II) complexes which NMR sign ranges from -1000 to -6000 ppm. Furthermore, the models were validated using a new set of 75 Pt(II) complexes, not included in the descriptive set. The models were constructed using non-relativistic Hamiltonian at density functional theory (DFT-PBEPBE) level with NMR-DKH basis set for all atoms. For the best model, the mean absolute deviation (MAD) and the mean relative deviation (MRD) were 150 ppm and 6%, respectively, for the validation set (75 Pt-complexes) and 168 ppm (MAD) and 5% (MRD) for all 258 Pt(II) complexes. These results were comparable with relativistic DFT calculation, 200 ppm (MAD) and 6% (MRD). © 2016 Wiley Periodicals, Inc. PMID:27510431
Miliordos, Evangelos; Xantheas, Sotiris S.
2015-03-07
We report the variation of the binding energy of the formic acid dimer at the CCSD(T)/ Complete Basis Set limit and examine the validity of the BSSE-correction, previously challenged by Kalescky, Kraka and Cremer [J. Chem. Phys. 140 (2014) 084315]. Our best estimate of D0=14.3±0.1 kcal/mol is in excellent agreement with the experimental value of 14.22±0.12 kcal/mol. The BSSE correction is indeed valid for this system since it exhibits the expected behavior of decreasing with increasing basis set size and its inclusion produces the same limit (within 0.1 kcal/mol) as the one obtained from extrapolation of the uncorrected binding energy. This work was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated for DOE by Battelle. A portion of this research was performed using the Molecular Science Computing Facility (MSCF) in EMSL, a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research and located at PNNL.
An extended basis set {ital ab} {ital initio} study of alkali metal cation--water clusters
Feller, D.; Glendening, E.D.; Woon, D.E.; Feyereisen, M.W.
1995-09-01
Ionic clusters comprised of a single alkali metal cation and up to eight water molecules were studied at the Hartree--Fock and correlated levels of theory using the correlation consistent sequence of basis sets. Estimates of the degree of convergence in the computed properties with respect to the complete basis set limit were facilitated by the underlying systematic manner in which the correlation consistent sets approach completeness. In favorable cases, improved property values could be obtained by fitting finite basis set results with a simple analytical expression in order to extrapolate to the complete basis set limit. The sensitivity of structures and binding energies were analyzed with regard to the inclusion of valence and core-valence correlation recovery at the MP2, MP4, and CCSD(T) levels of theory. The replacement of metal core electrons and the introduction of relativistic contributions via effective core potentials was compared to corresponding all-electron results. {copyright} {ital 1995} {ital American} {ital Institute} {ital of} {ital Physics}.
Molecular Dipole Moments within the Incremental Scheme Using the Domain-Specific Basis-Set Approach.
Fiedler, Benjamin; Coriani, Sonia; Friedrich, Joachim
2016-07-12
We present the first implementation of the fully automated incremental scheme for CCSD unrelaxed dipole moments using the domain-specific basis-set approach. Truncation parameters are varied, and the accuracy of the method is statistically analyzed for a test set of 20 molecules. The local approximations introduce small errors at second order and negligible ones at third order. For a third-order incremental CCSD expansion with a CC2 error correction, a cc-pVDZ/SV domain-specific basis set (tmain = 3.5 Bohr), and the truncation parameter f = 30 Bohr, we obtain a mean error of 0.00 mau (-0.20 mau) and a standard deviation of 1.95 mau (2.17 mau) for the total dipole moments (Cartesian components of the dipole vectors). By analyzing incremental CCSD energies, we demonstrate that the MP2 and CC2 error correction schemes are an exclusive correction for the domain-specific basis-set error. Our implementation of the incremental scheme provides fully automated computations of highly accurate dipole moments at reduced computational cost and is fully parallelized in terms of the calculation of the increments. Therefore, one can utilize the incremental scheme, on the same hardware, to extend the basis set in comparison to standard CCSD and thus obtain a better total accuracy. PMID:27300371
Zhao, Chunyu; Burge, James H
2013-12-16
Zernike polynomials are an orthonormal set of scalar functions over a circular domain, and are commonly used to represent wavefront phase or surface irregularity. In optical testing, slope or curvature of a surface or wavefront is sometimes measured instead, from which the surface or wavefront map is obtained. Previously we derived an orthonormal set of vector polynomials that fit to slope measurement data and yield the surface or wavefront map represented by Zernike polynomials. Here we define a 3-element curvature vector used to represent the second derivatives of a continuous surface, and derive a set of orthonormal curvature basis functions that are written in terms of Zernike polynomials. We call the new curvature functions the C polynomials. Closed form relations for the complete basis set are provided, and we show how to determine Zernike surface coefficients from the curvature data as represented by the C polynomials. PMID:24514717
How to spoil a good basis set for Rayleigh-Ritz calculations
Pupyshev, Vladimir I.; Montgomery, H. E. Jr.
2013-08-15
For model quantum mechanical systems such as the harmonic oscillator and a particle in an impenetrable box, we consider the set of exact discrete spectrum functions and define the modified basis set by subtraction of the ground state wavefunction from all the other wavefunctions with some real weights. It is demonstrated that the modified set of functions is complete in the space of square integrable functions if and only if the series of the squared weights diverges. A similar, but nonequivalent criterion is derived for convergence of Rayleigh-Ritz ground state energy calculations to the exact ground state energy value with the basis set extension. Some numerical illustrations are provided which demonstrate a wide variety of possible situations for model systems.
A Novel Gaussian-Sinc mixed Basis Set for Electronic Structure calculations
NASA Astrophysics Data System (ADS)
Jerke, Jonathan; Lee, Young; Tymczak, C. J.
2015-03-01
A Gaussian-Sinc mixed basis set for the computation of the electronic structure of atoms and molecules is presented. Excellent bases functions are known for ``core'' and ``valence'' separately, such as Gaussians for the ``core'' wave functions and Plane-waves for ``valance'' wave functions, but as yet no method is known that can accurately deal with both regimes in a single basis. A Gaussian-Sinc mixed basis can do both. This method resolves several issues such as: i) the Sincs basis spans the same space as the plane-waves basis, yet are semi-local enough to define all interaction elements including Exchange; ii) the Gaussians span the spherically symmetric core states and can be mixed with the Sinc functions in a computationally efficient methodology; iii) together, this mixed basis set is a flexible, computationally efficient and a highly accurate method for solving atomic and molecular problems. This methodology has been implemented within the Hartree-Fock level of theory within ultra-strong magnetic fields. To demonstrate the utility of this new method, we calculated the ground state Hartree-Fock energies to five digits accuracy in ultra strong magnetic fields for Helium to Neon, Molecular Hydrogen, Water, Carbon dioxide and Benzene. Welch Foundation (Grant J-1675), the ARO (Grant W911Nf-13-1-0162), the Texas Southern University High Performance Computing Center (http:/hpcc.tsu.edu/; Grant PHY-1126251) and NSF-CREST CRCN project (Grant HRD-1137732).
Training set optimization under population structure in genomic selection
Technology Transfer Automated Retrieval System (TEKTRAN)
The optimization of the training set (TRS) in genomic selection (GS) has received much interest in both animal and plant breeding, because it is critical to the accuracy of the prediction models. In this study, five different TRS sampling algorithms, stratified sampling, mean of the Coefficient of D...
Training set optimization under population structure in genomic selection
Technology Transfer Automated Retrieval System (TEKTRAN)
The optimization of the training set (TRS) in genomic selection has received much interest in both animal and plant breeding, because it is critical to the accuracy of the prediction models. In this study, five different TRS sampling algorithms, stratified sampling, mean of the coefficient of determ...
Correlation consistent basis sets for actinides. I. The Th and U atoms
Peterson, Kirk A.
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
Correlation consistent basis sets for actinides. I. The Th and U atoms.
Peterson, Kirk A
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 ThFn (n = 2 - 4), ThO2, and UFn (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 ThF4, ThF3, ThF2, and ThO2 are all within their experimental uncertainties. Bond dissociation energies of ThF4 and ThF3, as well as UF6 and UF5, 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 ThF4 and ThO2. The DKH3 atomization energy of ThO2 was calculated to be smaller than the DKH2 value by ∼1 kcal/mol. PMID:25702000
Correlation consistent basis sets for actinides. I. The Th and U atoms
NASA Astrophysics Data System (ADS)
Peterson, Kirk A.
2015-02-01
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 ThFn (n = 2 - 4), ThO2, and UFn (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 ThF4, ThF3, ThF2, and ThO2 are all within their experimental uncertainties. Bond dissociation energies of ThF4 and ThF3, as well as UF6 and UF5, 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 ThF4 and ThO2. The DKH3 atomization energy of ThO2 was calculated to be smaller than the DKH2 value by ˜1 kcal/mol.
Level set based structural topology optimization for minimizing frequency response
NASA Astrophysics Data System (ADS)
Shu, Lei; Wang, Michael Yu; Fang, Zongde; Ma, Zhengdong; Wei, Peng
2011-11-01
For the purpose of structure vibration reduction, a structural topology optimization for minimizing frequency response is proposed based on the level set method. The objective of the present study is to minimize the frequency response at the specified points or surfaces on the structure with an excitation frequency or a frequency range, subject to the given amount of the material over the admissible design domain. The sensitivity analysis with respect to the structural boundaries is carried out, while the Extended finite element method (X-FEM) is employed for solving the state equation and the adjoint equation. The optimal structure with smooth boundaries is obtained by the level set evolution with advection velocity, derived from the sensitivity analysis and the optimization algorithm. A number of numerical examples, in the frameworks of two-dimension (2D) and three-dimension (3D), are presented to demonstrate the feasibility and effectiveness of the proposed approach.
NASA Astrophysics Data System (ADS)
Borges, A.; Solomon, G. C.
2016-05-01
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.
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
Sensitivity of the properties of ruthenium ``blue dimer'' to method, basis set, and continuum model
NASA Astrophysics Data System (ADS)
Ozkanlar, Abdullah; Clark, Aurora E.
2012-05-01
The ruthenium "blue dimer" [(bpy)2RuIIIOH2]2O4+ is best known as the first well-defined molecular catalyst for water oxidation. It has been subject to numerous computational studies primarily employing density functional theory. However, those studies have been limited in the functionals, basis sets, and continuum models employed. The controversy in the calculated electronic structure and the reaction energetics of this catalyst highlights the necessity of benchmark calculations that explore the role of density functionals, basis sets, and continuum models upon the essential features of blue-dimer reactivity. In this paper, we report Kohn-Sham complete basis set (KS-CBS) limit extrapolations of the electronic structure of "blue dimer" using GGA (BPW91 and BP86), hybrid-GGA (B3LYP), and meta-GGA (M06-L) density functionals. The dependence of solvation free energy corrections on the different cavity types (UFF, UA0, UAHF, UAKS, Bondi, and Pauling) within polarizable and conductor-like polarizable continuum model has also been investigated. The most common basis sets of double-zeta quality are shown to yield results close to the KS-CBS limit; however, large variations are observed in the reaction energetics as a function of density functional and continuum cavity model employed.
Level-Set Topology Optimization with Aeroelastic Constraints
NASA Technical Reports Server (NTRS)
Dunning, Peter D.; Stanford, Bret K.; Kim, H. Alicia
2015-01-01
Level-set topology optimization is used to design a wing considering skin buckling under static aeroelastic trim loading, as well as dynamic aeroelastic stability (flutter). The level-set function is defined over the entire 3D volume of a transport aircraft wing box. Therefore, the approach is not limited by any predefined structure and can explore novel configurations. The Sequential Linear Programming (SLP) level-set method is used to solve the constrained optimization problems. The proposed method is demonstrated using three problems with mass, linear buckling and flutter objective and/or constraints. A constraint aggregation method is used to handle multiple buckling constraints in the wing skins. A continuous flutter constraint formulation is used to handle difficulties arising from discontinuities in the design space caused by a switching of the critical flutter mode.
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
Parameterization of Model Validating Sets for Uncertainty Bound Optimizations. Revised
NASA Technical Reports Server (NTRS)
Lim, K. B.; Giesy, D. P.
2000-01-01
Given measurement data, a nominal model and a linear fractional transformation uncertainty structure with an allowance on unknown but bounded exogenous disturbances, easily computable tests for the existence of a model validating uncertainty set are given. Under mild conditions, these tests are necessary and sufficient for the case of complex, nonrepeated, block-diagonal structure. For the more general case which includes repeated and/or real scalar uncertainties, the tests are only necessary but become sufficient if a collinearity condition is also satisfied. With the satisfaction of these tests, it is shown that a parameterization of all model validating sets of plant models is possible. The new parameterization is used as a basis for a systematic way to construct or perform uncertainty tradeoff with model validating uncertainty sets which have specific linear fractional transformation structure for use in robust control design and analysis. An illustrative example which includes a comparison of candidate model validating sets is given.
Heats of formation for third-period hydrides: test of an extended basis set
Gordon, M.S.; Heitzinger, J.
1987-04-23
The extended basis set developed by McLean and Chandler for third-period atoms is combined with the -311G hydrogen basis and augmented by polarization functions to predict the heats of formation of third-period hydrides at the full MP4 computational level. The calculated heats of formation are at least as accurate as those predicted for second-row hydrides by using MP4/6-311G(d,p) wave functions, with a root mean square error of approximately 4 kcal/mol.
Truncated spherical-wave basis set for first-principles pseudopotential calculations
NASA Astrophysics Data System (ADS)
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.
Benchmarks of electronically excited states: Basis set effects on CASPT2 results
NASA Astrophysics Data System (ADS)
Silva-Junior, Mario R.; Schreiber, Marko; Sauer, Stephan P. A.; Thiel, Walter
2010-11-01
Vertical excitation energies and one-electron properties are computed for the valence excited states of 28 medium-sized organic benchmark molecules using multistate multiconfigurational second-order perturbation theory (MS-CASPT2) and the augmented correlation-consistent aug-cc-pVTZ basis set. They are compared with previously reported MS-CASPT2 results obtained with the smaller TZVP basis. The basis set extension from TZVP to aug-cc-pVTZ causes rather minor and systematic shifts in the vertical excitation energies that are normally slightly reduced (on average by 0.11 eV for the singlets and by 0.09 eV for the triplets), whereas the changes in the calculated oscillator strengths and dipole moments are somewhat more pronounced on a relative scale. These basis set effects at the MS-CASPT2 level are qualitatively and quantitatively similar to those found at the coupled cluster level for the same set of benchmark molecules. The previously proposed theoretical best estimates (TBE-1) for the vertical excitation energies for 104 singlet and 63 triplet excited states of the benchmark molecules are upgraded by replacing TZVP with aug-cc-pVTZ data that yields a new reference set (TBE-2). Statistical evaluations of the performance of density functional theory (DFT) and semiempirical methods lead to the same ranking and very similar quantitative results for TBE-1 and TBE-2, with slightly better performance measures with respect to TBE-2. DFT/MRCI is most accurate among the investigated DFT-based approaches, while the OMx methods with orthogonalization corrections perform best at the semiempirical level.
Polyatomic molecular Dirac-Hartree-Fock calculations with Gaussian basis sets
NASA Technical Reports Server (NTRS)
Dyall, Kenneth G.; Faegri, Knut, Jr.; Taylor, Peter R.
1990-01-01
Numerical methods have been used successfully in atomic Dirac-Hartree-Fock (DHF) calculations for many years. Some DHF calculations using numerical methods have been done on diatomic molecules, but while these serve a useful purpose for calibration, the computational effort in extending this approach to polyatomic molecules is prohibitive. An alternative more in line with traditional quantum chemistry is to use an analytical basis set expansion of the wave function. This approach fell into disrepute in the early 1980's due to problems with variational collapse and intruder states, but has recently been put on firm theoretical foundations. In particular, the problems of variational collapse are well understood, and prescriptions for avoiding the most serious failures have been developed. Consequently, it is now possible to develop reliable molecular programs using basis set methods. This paper describes such a program and reports results of test calculations to demonstrate the convergence and stability of the method.
Method and Basis Set Analysis of Oxorhenium(V) Complexes for Theoretical Calculations
Demoin, Dustin Wayne; Li, Yawen; Jurisson, Silvia S.; Deakyne, Carol A.
2012-01-01
A variety of method and basis set combinations has been evaluated for monooxorhenium(V) complexes with N, O, P, S, Cl, and Se donor atoms. The geometries and energies obtained are compared to both high-level computations and literature structures. These calculations show that the PBE0 method outperforms the B3LYP method with respect to both structure and energetics. The combination of 6-31G** basis set on the nonmetal atoms and LANL2TZ effective core potential on the rhenium center gives reliable equilibrium structures with minimal computational resources for both model and literature compounds. Single-point energy calculations at the PBE0/LANL2TZ,6-311+G* level of theory are recommended for energetics. PMID:23087847
Elimination of the linearization error and improved basis-set convergence within the FLAPW method
NASA Astrophysics Data System (ADS)
Michalicek, Gregor; Betzinger, Markus; Friedrich, Christoph; Blügel, Stefan
2013-12-01
We analyze in detail the error that arises from the linearization in linearized augmented-plane-wave (LAPW) basis functions around predetermined energies El and show that it can lead to undesirable dependences of the calculated results on method-inherent parameters such as energy parameters El and muffin-tin sphere radii. To overcome these dependences, we evaluate approaches that eliminate the linearization error systematically by adding local orbitals (LOs) to the basis set. We consider two kinds of LOs: (i) constructed from solutions ul(r,E) to the scalar-relativistic approximation of the radial Dirac equation with E>El and (ii) constructed from second energy derivatives ∂2ul(r,E)/∂E2 at E=El. We find that the latter eliminates the error most efficiently and yields the density functional answer to many electronic and materials properties with very high precision. Finally, we demonstrate that the so constructed LAPW +LO basis shows a more favorable convergence behavior than the conventional LAPW basis due to a better decoupling of muffin-tin and interstitial regions, similarly to the related APW +lo approach, which, however, requires an extra set of LOs to reach the same total energy.
Representing the behavior of partially coherent optical systems by using overcomplete basis sets
NASA Astrophysics Data System (ADS)
Withington, Stafford; Hobson, Michael P.; Berry, Rachel H.
2004-02-01
A technique is described for representing the behavior of partially coherent optical systems by using overcomplete basis sets. The scheme is closely related to Gabor function theory. Through singular-value decomposition it is shown that if E is a matrix containing the sampled basis functions, then all of the information needed for optical calculations is contained in S=EE† and R=E†E. For overcomplete sets, S can be inverted to give a dual basis set, E~=S-1E, which can be used to find the correlation matrix elements A of a sampled bimodal expansion of the spatial coherence function. Overcomplete correlation matrices can be scattered easily at optical components. They can be used to determine (i) the natural modes of a field; (ii) the total power in a field, Pt=Tr[RA] (iii) the power coupled between two fields, A and B, that are in different states of coherence, Pc=Tr[RARB] and (iv) the entropy of a field, Q=Tr[ZΣr(I-Z)r/r], where Z=RA/Tr[RA].
Representing the behavior of partially coherent optical systems by using overcomplete basis sets.
Withington, Stafford; Hobson, Michael P; Berry, Rachel H
2004-02-01
A technique is described for representing the behavior of partially coherent optical systems by using overcomplete basis sets. The scheme is closely related to Gabor function theory. Through singular-value decomposition it is shown that if E is a matrix containing the sampled basis functions, then all of the information needed for optical calculations is contained in S = EE(dagger) and R = E(dagger)E. For overcomplete sets, S can be inverted to give a dual basis set, E = S(-1)E, which can be used to find the correlation matrix elements A of a sampled bimodal expansion of the spatial coherence function. Overcomplete correlation matrices can be scattered easily at optical components. They can be used to determine (i) the natural modes of a field; (ii) the total power in a field, Pt = Tr[RA]; (iii) the power coupled between two fields, A and B, that are in different states of coherence, Pc = Tr[RARB]; and (iv) the entropy of a field, Q = Tr[Zsigmar(I-Z)r/r], where Z = RA/Tr[RA]. PMID:14763763
Choi, Sunghwan; Hong, Kwangwoo; Kim, Jaewook; Kim, Woo Youn
2015-03-07
We developed a self-consistent field program based on Kohn-Sham density functional theory using Lagrange-sinc functions as a basis set and examined its numerical accuracy for atoms and molecules through comparison with the results of Gaussian basis sets. The result of the Kohn-Sham inversion formula from the Lagrange-sinc basis set manifests that the pseudopotential method is essential for cost-effective calculations. The Lagrange-sinc basis set shows faster convergence of the kinetic and correlation energies of benzene as its size increases than the finite difference method does, though both share the same uniform grid. Using a scaling factor smaller than or equal to 0.226 bohr and pseudopotentials with nonlinear core correction, its accuracy for the atomization energies of the G2-1 set is comparable to all-electron complete basis set limits (mean absolute deviation ≤1 kcal/mol). The same basis set also shows small mean absolute deviations in the ionization energies, electron affinities, and static polarizabilities of atoms in the G2-1 set. In particular, the Lagrange-sinc basis set shows high accuracy with rapid convergence in describing density or orbital changes by an external electric field. Moreover, the Lagrange-sinc basis set can readily improve its accuracy toward a complete basis set limit by simply decreasing the scaling factor regardless of systems.
On fast computation of finite-time coherent sets using radial basis functions
NASA Astrophysics Data System (ADS)
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.
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
Relativistic correlating basis sets for actinide atoms from 90Th to 103Lr.
Noro, Takeshi; Sekiya, Masahiro; Osanai, You; Koga, Toshikatsu; Matsuyama, Hisashi
2007-12-01
For 14 actinide atoms from (90)Th to (103)Lr, contracted Gaussian-type function sets are developed for the description of correlations of the 5f, 6d, and 7s electrons. Basis sets for the 6d orbitals are also prepared, since the orbitals are important in molecular environments despite their vacancy in the ground state of some actinides. A segmented contraction scheme is employed for the compactness and efficiency. Contraction coefficients and exponents are so determined as to minimize the deviation from accurate natural orbitals of the lowest term arising from the 5f(n-1)6d(1)7s(2) configuration. The spin-free relativistic effects are considered through the third-order Douglas-Kroll approximation. To test the present correlating sets, all-electron calculations are performed on the ground state of (90)ThO molecule. The calculated spectroscopic constants are in excellent agreement with experimental values. PMID:17508413
XFEM schemes for level set based structural optimization
NASA Astrophysics Data System (ADS)
Li, Li; Wang, Michael Yu; Wei, Peng
2012-12-01
In this paper, some elegant extended finite element method (XFEM) schemes for level set method structural optimization are proposed. Firstly, two-dimension (2D) and three-dimension (3D) XFEM schemes with partition integral method are developed and numerical examples are employed to evaluate their accuracy, which indicate that an accurate analysis result can be obtained on the structural boundary. Furthermore, the methods for improving the computational accuracy and efficiency of XFEM are studied, which include the XFEM integral scheme without quadrature sub-cells and higher order element XFEM scheme. Numerical examples show that the XFEM scheme without quadrature sub-cells can yield similar accuracy of structural analysis while prominently reducing the time cost and that higher order XFEM elements can improve the computational accuracy of structural analysis in the boundary elements, but the time cost is increasing. Therefore, the balance of time cost between FE system scale and the order of element needs to be discussed. Finally, the reliability and advantages of the proposed XFEM schemes are illustrated with several 2D and 3D mean compliance minimization examples that are widely used in the recent literature of structural topology optimization. All numerical results demonstrate that the proposed XFEM is a promising structural analysis approach for structural optimization with the level set method.
Ghost transmission: How large basis sets can make electron transport calculations worse
Herrmann, Carmen; Solomon, Gemma C.; Subotnik, Joseph E.; Mujica, Vladimiro; Ratner, Mark A.
2010-01-01
The Landauer approach has proven to be an invaluable tool for calculating the electron transport properties of single molecules, especially when combined with a nonequilibrium Green’s function approach and Kohn–Sham density functional theory. However, when using large nonorthogonal atom-centered basis sets, such as those common in quantum chemistry, one can find erroneous results if the Landauer approach is applied blindly. In fact, basis sets of triple-zeta quality or higher sometimes result in an artificially high transmission and possibly even qualitatively wrong conclusions regarding chemical trends. In these cases, transport persists when molecular atoms are replaced by basis functions alone (“ghost atoms”). The occurrence of such ghost transmission is correlated with low-energy virtual molecular orbitals of the central subsystem and may be interpreted as a biased and thus inaccurate description of vacuum transmission. An approximate practical correction scheme is to calculate the ghost transmission and subtract it from the full transmission. As a further consequence of this study, it is recommended that sensitive molecules be used for parameter studies, in particular those whose transmission functions show antiresonance features such as benzene-based systems connected to the electrodes in meta positions and other low-conducting systems such as alkanes and silanes.
Variational Trajectory Optimization Tool Set: Technical description and user's manual
NASA Technical Reports Server (NTRS)
Bless, Robert R.; Queen, Eric M.; Cavanaugh, Michael D.; Wetzel, Todd A.; Moerder, Daniel D.
1993-01-01
The algorithms that comprise the Variational Trajectory Optimization Tool Set (VTOTS) package are briefly described. The VTOTS is a software package for solving nonlinear constrained optimal control problems from a wide range of engineering and scientific disciplines. The VTOTS package was specifically designed to minimize the amount of user programming; in fact, for problems that may be expressed in terms of analytical functions, the user needs only to define the problem in terms of symbolic variables. This version of the VTOTS does not support tabular data; thus, problems must be expressed in terms of analytical functions. The VTOTS package consists of two methods for solving nonlinear optimal control problems: a time-domain finite-element algorithm and a multiple shooting algorithm. These two algorithms, under the VTOTS package, may be run independently or jointly. The finite-element algorithm generates approximate solutions, whereas the shooting algorithm provides a more accurate solution to the optimization problem. A user's manual, some examples with results, and a brief description of the individual subroutines are included.
Optimizing antibiotic therapy in the intensive care unit setting
Kollef, Marin H
2001-01-01
Antibiotics are one of the most common therapies administered in the intensive care unit setting. In addition to treating infections, antibiotic use contributes to the emergence of resistance among pathogenic microorganisms. Therefore, avoiding unnecessary antibiotic use and optimizing the administration of antimicrobial agents will help to improve patient outcomes while minimizing further pressures for resistance. This review will present several strategies aimed at achieving optimal use of antimicrobial agents. It is important to note that each intensive care unit should have a program in place which monitors antibiotic utilization and its effectiveness. Only in this way can the impact of interventions aimed at improving antibiotic use (e.g. antibiotic rotation, de-escalation therapy) be evaluated at the local level. PMID:11511331
An Optimal Set of Flesh Points on Tongue and Lips for Speech-Movement Classification
Samal, Ashok; Rong, Panying; Green, Jordan R.
2016-01-01
Purpose The authors sought to determine an optimal set of flesh points on the tongue and lips for classifying speech movements. Method The authors used electromagnetic articulographs (Carstens AG500 and NDI Wave) to record tongue and lip movements from 13 healthy talkers who articulated 8 vowels, 11 consonants, a phonetically balanced set of words, and a set of short phrases during the recording. We used a machine-learning classifier (support-vector machine) to classify the speech stimuli on the basis of articulatory movements. We then compared classification accuracies of the flesh-point combinations to determine an optimal set of sensors. Results When data from the 4 sensors (T1: the vicinity between the tongue tip and tongue blade; T4: the tongue-body back; UL: the upper lip; and LL: the lower lip) were combined, phoneme and word classifications were most accurate and were comparable with the full set (including T2: the tongue-body front; and T3: the tongue-body front). Conclusion We identified a 4-sensor set—that is, T1, T4, UL, LL—that yielded a classification accuracy (91%–95%) equivalent to that using all 6 sensors. These findings provide an empirical basis for selecting sensors and their locations for scientific and emerging clinical applications that incorporate articulatory movements. PMID:26564030
Optimization of the generator settings for endobiliary radiofrequency ablation
Barret, Maximilien; Leblanc, Sarah; Vienne, Ariane; Rouquette, Alexandre; Beuvon, Frederic; Chaussade, Stanislas; Prat, Frederic
2015-01-01
AIM: To determine the optimal generator settings for endobiliary radiofrequency ablation. METHODS: Endobiliary radiofrequency ablation was performed in live swine on the ampulla of Vater, the common bile duct and in the hepatic parenchyma. Radiofrequency ablation time, “effect”, and power were allowed to vary. The animals were sacrificed two hours after the procedure. Histopathological assessment of the depth of the thermal lesions was performed. RESULTS: Twenty-five radiofrequency bursts were applied in three swine. In the ampulla of Vater (n = 3), necrosis of the duodenal wall was observed starting with an effect set at 8, power output set at 10 W, and a 30 s shot duration, whereas superficial mucosal damage of up to 350 μm in depth was recorded for an effect set at 8, power output set at 6 W and a 30 s shot duration. In the common bile duct (n = 4), a 1070 μm, safe and efficient ablation was obtained for an effect set at 8, a power output of 8 W, and an ablation time of 30 s. Within the hepatic parenchyma (n = 18), the depth of tissue damage varied from 1620 μm (effect = 8, power = 10 W, ablation time = 15 s) to 4480 μm (effect = 8, power = 8 W, ablation time = 90 s). CONCLUSION: The duration of the catheter application appeared to be the most important parameter influencing the depth of the thermal injury during endobiliary radiofrequency ablation. In healthy swine, the currently recommended settings of the generator may induce severe, supratherapeutic tissue damage in the biliary tree, especially in the high-risk area of the ampulla of Vater. PMID:26566429
Optimal exposure sets for high dynamic range scenes
NASA Astrophysics Data System (ADS)
Valli Kumari, V.; RaviKiran, B.; Raju, K. V. S. V. N.; Shajahan Basha, S. A.
2011-10-01
The dynamic range of many natural scenes is far greater than the dynamic range of the imaging devices. These scenes present a challenge to the consumer digital cameras. The well-known technique to capture the full dynamic range of the scene is by fusing multiple images of the same scene. Usually people combine three or five different exposures to get the full dynamic range of the scene. Some cameras like Pentax K-7, always combines fixed exposures together to produce the output result. However, this should be adaptive to the scene characteristics. We propose an optimal solution for dynamically selecting the exposure sets.
Asynchronous parallel generating set search for linearly-constrained optimization.
Lewis, Robert Michael; Griffin, Joshua D.; Kolda, Tamara Gibson
2006-08-01
Generating set search (GSS) is a family of direct search methods that encompasses generalized pattern search and related methods. We describe an algorithm for asynchronous linearly-constrained GSS, which has some complexities that make it different from both the asynchronous bound-constrained case as well as the synchronous linearly-constrained case. The algorithm has been implemented in the APPSPACK software framework and we present results from an extensive numerical study using CUTEr test problems. We discuss the results, both positive and negative, and conclude that GSS is a reliable method for solving small-to-medium sized linearly-constrained optimization problems without derivatives.
NASA Astrophysics Data System (ADS)
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
Grimme, Stefan Brandenburg, Jan Gerit; Bannwarth, Christoph; Hansen, Andreas
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
Imputing gene expression from optimally reduced probe sets
Donner, Yoni; Feng, Ting; Benoist, Christophe; Koller, Daphne
2012-01-01
Measuring complete gene expression profiles for a large number of experiments is costly. We propose an approach in which a small subset of probes is selected based on a preliminary set of full expression profiles. In subsequent experiments, only the subset is measured, and the missing values are imputed. We develop several algorithms to simultaneously select probes and impute missing values, and demonstrate that these probe selection for imputation (PSI) algorithms can successfully reconstruct missing gene expression values in a wide variety of applications, as evaluated using multiple metrics of biological importance. We analyze the performance of PSI methods under varying conditions, provide guidelines for choosing the optimal method based on the experimental setting, and indicate how to estimate imputation accuracy. Finally, we apply our approach to a large-scale study of immune system variation. PMID:23064520
Correlation consistent basis sets for lanthanides: The atoms La-Lu.
Lu, Qing; Peterson, Kirk A
2016-08-01
Using the 3rd-order Douglas-Kroll-Hess (DKH3) Hamiltonian, all-electron correlation consistent basis sets of double-, triple-, and quadruple-zeta quality have been developed for the lanthanide elements La through Lu. Basis sets designed for the recovery of valence correlation (defined here as 4f5s5p5d6s), cc-pVnZ-DK3, and outer-core correlation (valence + 4s4p4d), cc-pwCVnZ-DK3, are reported (n = D, T, and Q). Systematic convergence of both Hartree-Fock and correlation energies towards their respective complete basis set (CBS) limits are observed. Benchmark calculations of the first three ionization potentials (IPs) of La through Lu are reported at the DKH3 coupled cluster singles and doubles with perturbative triples, CCSD(T), level of theory, including effects of correlation down through the 4s electrons. Spin-orbit coupling is treated at the 2-component HF level. After extrapolation to the CBS limit, the average errors with respect to experiment were just 0.52, 1.14, and 4.24 kcal/mol for the 1st, 2nd, and 3rd IPs, respectively, compared to the average experimental uncertainties of 0.03, 1.78, and 2.65 kcal/mol, respectively. The new basis sets are also used in CCSD(T) benchmark calculations of the equilibrium geometries, atomization energies, and heats of formation for Gd2, GdF, and GdF3. Except for the equilibrium geometry and harmonic frequency of GdF, which are accurately known from experiment, all other calculated quantities represent significant improvements compared to the existing experimental quantities. With estimated uncertainties of about ±3 kcal/mol, the 0 K atomization energies (298 K heats of formation) are calculated to be (all in kcal/mol): 33.2 (160.1) for Gd2, 151.7 (-36.6) for GdF, and 447.1 (-295.2) for GdF3. PMID:27497543
Correlation consistent basis sets for lanthanides: The atoms La-Lu
NASA Astrophysics Data System (ADS)
Lu, Qing; Peterson, Kirk A.
2016-08-01
Using the 3rd-order Douglas-Kroll-Hess (DKH3) Hamiltonian, all-electron correlation consistent basis sets of double-, triple-, and quadruple-zeta quality have been developed for the lanthanide elements La through Lu. Basis sets designed for the recovery of valence correlation (defined here as 4f5s5p5d6s), cc-pVnZ-DK3, and outer-core correlation (valence + 4s4p4d), cc-pwCVnZ-DK3, are reported (n = D, T, and Q). Systematic convergence of both Hartree-Fock and correlation energies towards their respective complete basis set (CBS) limits are observed. Benchmark calculations of the first three ionization potentials (IPs) of La through Lu are reported at the DKH3 coupled cluster singles and doubles with perturbative triples, CCSD(T), level of theory, including effects of correlation down through the 4s electrons. Spin-orbit coupling is treated at the 2-component HF level. After extrapolation to the CBS limit, the average errors with respect to experiment were just 0.52, 1.14, and 4.24 kcal/mol for the 1st, 2nd, and 3rd IPs, respectively, compared to the average experimental uncertainties of 0.03, 1.78, and 2.65 kcal/mol, respectively. The new basis sets are also used in CCSD(T) benchmark calculations of the equilibrium geometries, atomization energies, and heats of formation for Gd2, GdF, and GdF3. Except for the equilibrium geometry and harmonic frequency of GdF, which are accurately known from experiment, all other calculated quantities represent significant improvements compared to the existing experimental quantities. With estimated uncertainties of about ±3 kcal/mol, the 0 K atomization energies (298 K heats of formation) are calculated to be (all in kcal/mol): 33.2 (160.1) for Gd2, 151.7 (-36.6) for GdF, and 447.1 (-295.2) for GdF3.
Using symmetry-adapted optimized sum-of-products basis functions to calculate vibrational spectra
NASA Astrophysics Data System (ADS)
Leclerc, Arnaud; Carrington, Tucker
2016-01-01
Vibrational spectra can be computed without storing full-dimensional vectors by using low-rank sum-of-products (SOP) basis functions. We introduce symmetry constraints in the SOP basis functions to make it possible to separately calculate states in different symmetry subgroups. This is done using a power method to compute eigenvalues and an alternating least squares method to optimize basis functions. Owing to the fact that the power method favours the convergence of the lowest states, one must be careful not to exclude basis functions of some symmetries. Exploiting symmetry facilitates making assignments and improves the accuracy. The method is applied to the acetonitrile molecule.
Optimal sets of measurement data for profile reconstruction in scatterometry
NASA Astrophysics Data System (ADS)
Gross, H.; Rathsfeld, A.; Scholze, F.; Bär, M.; Dersch, U.
2007-06-01
We discuss numerical algorithms for the determination of periodic surface structures from light diffraction patterns. With decreasing feature sizes of lithography masks, increasing demands on metrology techniques arise. Scatterometry as a non-imaging indirect optical method is applied to simple periodic line structures in order to determine parameters like side-wall angles, heights, top and bottom widths and to evaluate the quality of the manufacturing process. The numerical simulation of diffraction is based on the finite element solution of the Helmholtz equation. The inverse problem seeks to reconstruct the grating geometry from measured diffraction patterns. Restricting the class of gratings and the set of measurements, this inverse problem can be reformulated as a non-linear operator equation in Euclidean spaces. The operator maps the grating parameters to special efficiencies of diffracted plane wave modes. We employ a Gauss-Newton type iterative method to solve this operator equation. The reconstruction properties and the convergence of the algorithm, however, is controlled by the local conditioning of the non-linear mapping. To improve reconstruction and convergence, we determine optimal sets of efficiencies optimizing the condition numbers of the corresponding Jacobians. Numerical examples are presented for "chrome on glass" masks under the wavelength 632.8 nm and for EUV masks.
Many-body calculations of molecular electric polarizabilities in asymptotically complete basis sets
NASA Astrophysics Data System (ADS)
Monten, Ruben; Hajgató, Balázs; Deleuze, Michael S.
2011-10-01
The static dipole polarizabilities of Ne, CO, N2, F2, HF, H2O, HCN, and C2H2 (acetylene) have been determined close to the Full-CI limit along with an asymptotically complete basis set (CBS), according to the principles of a Focal Point Analysis. For this purpose the results of Finite Field calculations up to the level of Coupled Cluster theory including Single, Double, Triple, Quadruple and perturbative Pentuple excitations [CCSDTQ(P)] were used, in conjunction with suited extrapolations of energies obtained using augmented and doubly-augmented Dunning's correlation consistent polarized valence basis sets of improving quality. The polarizability characteristics of C2H4 (ethylene) and C2H6 (ethane) have been determined on the same grounds at the CCSDTQ level in the CBS limit. Comparison is made with results obtained using lower levels in electronic correlation, or taking into account the relaxation of the molecular structure due to an adiabatic polarization process. Vibrational corrections to electronic polarizabilities have been empirically estimated according to Born-Oppenheimer Molecular Dynamical simulations employing Density Functional Theory. Confrontation with experiment ultimately indicates relative accuracies of the order of 1 to 2%.
Spectral fitting using basis set modified by measured B0 field distribution.
Li, Ningzhi; An, Li; Shen, Jun
2015-12-01
This study sought to demonstrate and evaluate a novel spectral fitting method to improve quantification accuracy in the presence of large magnetic field distortion, especially with high fields. MRS experiments were performed using a point-resolved spectroscopy (PRESS)-type sequence at 7 T. A double-echo gradient echo (GRE) sequence was used to acquire B0 maps following MRS experiments. The basis set was modified based on the measured B0 distribution within the MRS voxel. Quantification results were obtained after fitting the measured MRS data using the modified basis set. The proposed method was validated using numerical Monte Carlo simulations, phantom measurements, and comparison of occipital lobe MRS measurements under homogeneous and inhomogeneous magnetic field conditions. In vivo results acquired from voxels placed in thalamus and prefrontal cortex regions close to the frontal sinus agreed well with published values. Instead of noise-amplifying complex division, the proposed method treats field variations as part of the signal model, thereby avoiding inherent statistical bias associated with regularization. Simulations and experiments showed that the proposed approach reliably quantified results in the presence of relatively large magnetic field distortion. Published 2015. This article is a U.S. Government work and is in the public domain in the USA. PMID:26503305
Analytic basis set for high-Z atomic QED calculations: Heavy He-like ions
Hylton, D.J.; Snyderman, N.J.
1997-04-01
A relativistic Sturmian analytic basis set representation for the Coulomb-Dirac Green function, previously studied by Zapryagaev, Manakov, and Pal{close_quote}chikov [Opt. Spectrosc. {bold 52}, 248 (1982)], is investigated for application to high-Z atomic QED calculations. This pseudoeigenfunction representation follows from exact identities starting from the Whittaker function representation. It eliminates the radial ordering problem of that representation, and so is particularly useful for numerical calculation of the perturbation theory Feynman diagrams with more than one electron Green function. While the Green function represents discrete bound states, and both positive and negative energy continuum states, the Sturmian (bound-state-like) form for the pseudoeigenfunctions makes it possible to more analytically calculate matrix elements for full photon exchange, reducing numerical problems for high photon frequency. For He-like Fm (Z=100) we calculate the perturbation theory equivalent of the Dirac-Fock-Breit ground-state energy, agreeing well with the Grant code and with the numerical B-spline basis set approach results of Blundell, Mohr, Johnson, and Sapirstein [Phys. Rev. A {bold 48}, 2615 (1993)]. Preliminary results on the relativistic and QED correlation are also reported. {copyright} {ital 1997} {ital The American Physical Society}
Nonlinearly-constrained optimization using asynchronous parallel generating set search.
Griffin, Joshua D.; Kolda, Tamara Gibson
2007-05-01
Many optimization problems in computational science and engineering (CS&E) are characterized by expensive objective and/or constraint function evaluations paired with a lack of derivative information. Direct search methods such as generating set search (GSS) are well understood and efficient for derivative-free optimization of unconstrained and linearly-constrained problems. This paper addresses the more difficult problem of general nonlinear programming where derivatives for objective or constraint functions are unavailable, which is the case for many CS&E applications. We focus on penalty methods that use GSS to solve the linearly-constrained problems, comparing different penalty functions. A classical choice for penalizing constraint violations is {ell}{sub 2}{sup 2}, the squared {ell}{sub 2} norm, which has advantages for derivative-based optimization methods. In our numerical tests, however, we show that exact penalty functions based on the {ell}{sub 1}, {ell}{sub 2}, and {ell}{sub {infinity}} norms converge to good approximate solutions more quickly and thus are attractive alternatives. Unfortunately, exact penalty functions are discontinuous and consequently introduce theoretical problems that degrade the final solution accuracy, so we also consider smoothed variants. Smoothed-exact penalty functions are theoretically attractive because they retain the differentiability of the original problem. Numerically, they are a compromise between exact and {ell}{sub 2}{sup 2}, i.e., they converge to a good solution somewhat quickly without sacrificing much solution accuracy. Moreover, the smoothing is parameterized and can potentially be adjusted to balance the two considerations. Since many CS&E optimization problems are characterized by expensive function evaluations, reducing the number of function evaluations is paramount, and the results of this paper show that exact and smoothed-exact penalty functions are well-suited to this task.
On Optimizing Joint Inversion of Constrained Geophysical Data Sets
NASA Astrophysics Data System (ADS)
Sosa Aguirre, U. A.; Velazquez, L.; Argaez, M.; Velasco, A. A.; Romero, R.
2010-12-01
We implemented a joint inversion least-squares (LSQ) algorithm to characterize 1-D crustal velocity Earth structure using geophysical data sets with two different optimization methods: truncated singular value decomposition (TSVD), and primal-dual interior-point (PDIP). We used receiver function and surface wave dispersion velocity observations, and created a framework to incorporate other data sets. An improvement in the final outcome (model) is expected by providing better physical constraints than using just one single data set. The TSVD and PDIP methods solve a regularized unconstrained and an inherent regularized constrained minimization problems, respectively. Both techniques implement the inclusion of bounds into the layered shear velocities in a different fashion. We conduct a numerical experimentation with synthetic data, and find that the PDID method’s solution was more robust in terms of satisfying geophysical constraints, accuracy, and efficiency than the TSVD approach. Finally, we apply the PDIP method for characterizing material properties of the Rio Grande Rift region using real recorded seismic data with promising numerical results.
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-04
Korringa-Kohn-Rostoker (KKR) Green's function, multiple-scattering theory is an ecient sitecentered, 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 L_{max} = (l,m)_{max}, while scattering matrices, which determine spectral properties, are truncated at L_{tr} = (l,m)_{tr} where phase shifts δl>l_{tr} are negligible. Historically, L_{max} is set equal to L_{tr}, which is correct for large enough L_{max} but not computationally expedient; a better procedure retains higher-order (free-electron and single-site) contributions for L_{max} > L_{tr} with δl>l_{tr} set to zero [Zhang and Butler, Phys. Rev. B 46, 7433]. We present a numerically ecient and accurate augmented-KKR Green's function formalism that solves the KKR equations by exact matrix inversion [R^{3} process with rank N(l_{tr} + 1)^{2}] and includes higher-L contributions via linear algebra [R^{2} process with rank N(l_{max} +1)^{2}]. 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 L1_{0} CoPt, and present the numerical results for accuracy and for the convergence of the total energies, Fermi energies, and magnetic moments versus L_{max} for a given L_{tr}.
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-04
Korringa-Kohn-Rostoker (KKR) Green's function, multiple-scattering theory is an ecient sitecentered, 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,m)max, while scattering matrices, which determine spectral properties, are truncated at Ltr = (l,m)tr where phase shifts δl>ltr are negligible. Historically, Lmax is set equal to Ltr, 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 > Ltr with δl>ltr set to zero [Zhang andmore » Butler, Phys. Rev. B 46, 7433]. We present a numerically ecient 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]. 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 L10 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 Ltr.« less
Green's function multiple-scattering theory with a truncated basis set: An augmented-KKR formalism
NASA Astrophysics Data System (ADS)
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.
TiCl, TiH and TiH+ Bond Energies, a Test of a Correlation Consistent Ti Basis Set
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Arnold, James O. (Technical Monitor)
1999-01-01
Correlation consistent basis sets are developed for Ti atom. The polarization functions are optimized for the average of the 3F and 5F states. One series of correlation consistent basis sets is for 3d and 4s correlation, while the second series includes 3s and 3p correlation as well as 3d and 4s correlation. These basis sets are tested using the Ti 3F-5F separation and the dissociation energies of TiCl X4Phi, TiH X4Phi, and TiH(+) X3Phi. The CCSD(T) complete basis set limit values are determined by extrapolation. The Douglas Kroll approach is used to compute the scalar relativistic effect. Spin-orbit effects are taken from experiment and/or computed at the CASSCF level. The Ti 3F-5F separation is in excellent agreement with experiment, while the TiCl, TiH, and TiH(+) bond energies are in good agreement with experiment. Extrapolation with the valence basis set is consistent with other atoms, while including 3s and 3p correlation appears to make extrapolation.
Crystal orbital Hamilton population (COHP) analysis as projected from plane-wave basis sets.
Deringer, Volker L; Tchougréeff, Andrei L; Dronskowski, Richard
2011-06-01
Simple, yet predictive bonding models are essential achievements of chemistry. In the solid state, in particular, they often appear in the form of visual bonding indicators. Because the latter require the crystal orbitals to be constructed from local basis sets, the application of the most popular density-functional theory codes (namely, those based on plane waves and pseudopotentials) appears as being ill-fitted to retrieve the chemical bonding information. In this paper, we describe a way to re-extract Hamilton-weighted populations from plane-wave electronic-structure calculations to develop a tool analogous to the familiar crystal orbital Hamilton population (COHP) method. We derive the new technique, dubbed "projected COHP" (pCOHP), and demonstrate its viability using examples of covalent, ionic, and metallic crystals (diamond, GaAs, CsCl, and Na). For the first time, this chemical bonding information is directly extracted from the results of plane-wave calculations. PMID:21548594
AgH, Ag/sub 2/, and AgO revisited: Basis set extensions
Martin, R.L.
1987-05-01
An extended basis set has been developed for Ag which significantly improves the agreement between theoretical and experimental spectroscopic parameters for AgH, AgO, and Ag/sub 2/. The major improvement comes about as a result of the improved treatment of electron correlation in the Ag d shell upon the introduction of f functions. Their inclusion produces very slight differences at the SCF level, but significant reductions in r/sub e/ and increases in ..omega../sub e/ and D/sub e/ in the Mo-dash-barller--Plesset perturbation theory expansion. At the MP4(SDTQ) level, typical results are 0.02 A too long for r/sub e/, 4% too low for ..omega../sub e/, and 10 kcal too small for D/sub e/. From a pragmatic standpoint, MP2 give results very similar to this at a much reduced level of effort.
Li, Z.; Pan, Y.K.; Tao, F.M.
1996-01-15
Bond function basis sets combined with the counterpoise procedure are used to calculate the molecular dissociation energies D{sub e} of 24 diatomic molecules and ions. The calculated values of D{sub e} are compared to those without bond functions and/or counterpoise corrections. The equilibrium bond lengths r{sub e}, and harmonic frequencies w{sub e} are also calculated for a few selected molecules. The calculations at the fourth-order-Moller-Plesset approximation (MP4) have consistently recovered about 95-99% of the experimental values for D{sub e}, compared to as low as 75% without use of bond functions. The calculated values of r{sub 3} are typically 0.01 {Angstrom} larger than the experimental values, and the calculated values of w{sub e} are over 95% of the experimental values. 37 refs., 2 tabs.
A two-dimensional volatility basis set - Part 3: Prognostic modeling and NOx dependence
NASA Astrophysics Data System (ADS)
Chuang, W. K.; Donahue, N. M.
2015-06-01
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.
A two-dimensional volatility basis set - Part 3: Prognostic modeling and NOx dependence
NASA Astrophysics Data System (ADS)
Chuang, W. K.; Donahue, N. M.
2016-01-01
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 (2D-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 Aerosol Mass Spectrometer (AMS) measurements. We also discuss the mathematical methods underlying the implementation of the 2D-VBS and provide the complete code in the Supplement. A developer version is available on Bitbucket, an online community repository.
Tsuzuki, Seiji; Tanabe, Kazutoshi )
1991-03-21
Intermolecular interaction energies of methane dimer were calculated by using several basis sets up to 6-311G(3d,4p) with electron correlation energy correction by the Moeller-Plesset perturbation method and basis set superposition error (BSSE) correction by the counterpoise method to evaluate the basis set effect. The calculated interaction energies depended on the basis set considerably. Whereas the interaction energies of repulsive component calculated at HF level were not affected by the change of basis set, the dispersion energy component dependent greatly on the basis set used. The dispersion energies calculated with the Moeller-Plesset second- and third-order perturbation by using 6-311G(2d,2p) basis set were 0-10% and 4-6% smaller than those obtained with the fourth-order (MP4(SDTQ)) perturbation, respectively. The BSSE's calculated by the counterpoise method were still about 30% of the calculated intermolecular interaction energies for the conformers of energy minima event at the MP4(SDTQ)/6-311G(2d,2p) level. The calculated interaction potentials of dimers at the MP4(SDTQ)/6-311G(2d,2p) level were considerably shallower than those obtained by MM2 force fields but were close to the potentials given by the Williams potential and by the recently reported MM3 force field.
Optimal gate-width setting for passive neutrons multiplicity counting
Croft, Stephen; Evans, Louise G; Schear, Melissa A
2010-01-01
When setting up a passive neutron coincidence counter it is natural to ask what coincidence gate settings should be used to optimize the counting precision. If the gate width is too short then signal is lost and the precision is compromised because in a given period only a few coincidence events will be observed. On the other hand if the gate is too large the signal will be maximized but it will also be compromised by the high level of random pile-up or Accidental coincidence events which must be subtracted. In the case of shift register electronics connected to an assay chamber with an exponential dieaway profile operating in the regime where the Accidentals rate dominates the Reals coincidence rate but where dead-time is not a concern, simple arguments allow one to show that the relative precision on the net Reals rate is minimized when the coincidence gate is set to about 1.2 times the lie dieaway time of the system. In this work we show that making the same assumptions it is easy to show that the relative precision on the Triples rates is also at a minimum when the relative precision of the Doubles (or Reals) is at a minimum. Although the analysis is straightforward to our knowledge such a discussion has not been documented in the literature before. Actual measurement systems do not always behave in the ideal we choose to model them. Fortunately however the variation in the relative precision as a function of gate width is rather flat for traditional safeguards counters and so the performance is somewhat forgiving of the exact choice. The derivation further serves to delineate the important parameters which determine the relative counting precision of the Doubles and Triples rates under the regime considered. To illustrate the similarities and differences we consider the relative standard deviation that might be anticipated for a passive correlation count of an axial section of a spent nuclear fuel assembly under practically achievable conditions.
NASA Astrophysics Data System (ADS)
Hsu, Po Jen; Lai, S. K.; Rapallo, Arnaldo
2014-03-01
Improved basis sets for the study of polymer dynamics by means of the diffusion theory, and tests on a melt of cis-1,4-polyisoprene decamers, and a toluene solution of a 71-mer syndiotactic trans-1,2-polypentadiene were presented recently [R. Gaspari and A. Rapallo, J. Chem. Phys. 128, 244109 (2008)]. The proposed hybrid basis approach (HBA) combined two techniques, the long time sorting procedure and the maximum correlation approximation. The HBA takes advantage of the strength of these two techniques, and its basis sets proved to be very effective and computationally convenient in describing both local and global dynamics in cases of flexible synthetic polymers where the repeating unit is a unique type of monomer. The question then arises if the same efficacy continues when the HBA is applied to polymers of different monomers, variable local stiffness along the chain and with longer persistence length, which have different local and global dynamical properties against the above-mentioned systems. Important examples of this kind of molecular chains are the proteins, so that a fragment of the protein transthyretin is chosen as the system of the present study. This peptide corresponds to a sequence that is structured in β-sheets of the protein and is located on the surface of the channel with thyroxin. The protein transthyretin forms amyloid fibrils in vivo, whereas the peptide fragment has been shown [C. P. Jaroniec, C. E. MacPhee, N. S. Astrof, C. M. Dobson, and R. G. Griffin, Proc. Natl. Acad. Sci. U.S.A. 99, 16748 (2002)] to form amyloid fibrils in vitro in extended β-sheet conformations. For these reasons the latter is given considerable attention in the literature and studied also as an isolated fragment in water solution where both experimental and theoretical efforts have indicated the propensity of the system to form β turns or α helices, but is otherwise predominantly unstructured. Differing from previous computational studies that employed implicit
Hsu, Po Jen; Lai, S K; Rapallo, Arnaldo
2014-03-14
Improved basis sets for the study of polymer dynamics by means of the diffusion theory, and tests on a melt of cis-1,4-polyisoprene decamers, and a toluene solution of a 71-mer syndiotactic trans-1,2-polypentadiene were presented recently [R. Gaspari and A. Rapallo, J. Chem. Phys. 128, 244109 (2008)]. The proposed hybrid basis approach (HBA) combined two techniques, the long time sorting procedure and the maximum correlation approximation. The HBA takes advantage of the strength of these two techniques, and its basis sets proved to be very effective and computationally convenient in describing both local and global dynamics in cases of flexible synthetic polymers where the repeating unit is a unique type of monomer. The question then arises if the same efficacy continues when the HBA is applied to polymers of different monomers, variable local stiffness along the chain and with longer persistence length, which have different local and global dynamical properties against the above-mentioned systems. Important examples of this kind of molecular chains are the proteins, so that a fragment of the protein transthyretin is chosen as the system of the present study. This peptide corresponds to a sequence that is structured in β-sheets of the protein and is located on the surface of the channel with thyroxin. The protein transthyretin forms amyloid fibrils in vivo, whereas the peptide fragment has been shown [C. P. Jaroniec, C. E. MacPhee, N. S. Astrof, C. M. Dobson, and R. G. Griffin, Proc. Natl. Acad. Sci. U.S.A. 99, 16748 (2002)] to form amyloid fibrils in vitro in extended β-sheet conformations. For these reasons the latter is given considerable attention in the literature and studied also as an isolated fragment in water solution where both experimental and theoretical efforts have indicated the propensity of the system to form β turns or α helices, but is otherwise predominantly unstructured. Differing from previous computational studies that employed implicit
Hsu, Po Jen; Lai, S. K.; Rapallo, Arnaldo
2014-03-14
Improved basis sets for the study of polymer dynamics by means of the diffusion theory, and tests on a melt of cis-1,4-polyisoprene decamers, and a toluene solution of a 71-mer syndiotactic trans-1,2-polypentadiene were presented recently [R. Gaspari and A. Rapallo, J. Chem. Phys. 128, 244109 (2008)]. The proposed hybrid basis approach (HBA) combined two techniques, the long time sorting procedure and the maximum correlation approximation. The HBA takes advantage of the strength of these two techniques, and its basis sets proved to be very effective and computationally convenient in describing both local and global dynamics in cases of flexible synthetic polymers where the repeating unit is a unique type of monomer. The question then arises if the same efficacy continues when the HBA is applied to polymers of different monomers, variable local stiffness along the chain and with longer persistence length, which have different local and global dynamical properties against the above-mentioned systems. Important examples of this kind of molecular chains are the proteins, so that a fragment of the protein transthyretin is chosen as the system of the present study. This peptide corresponds to a sequence that is structured in β-sheets of the protein and is located on the surface of the channel with thyroxin. The protein transthyretin forms amyloid fibrils in vivo, whereas the peptide fragment has been shown [C. P. Jaroniec, C. E. MacPhee, N. S. Astrof, C. M. Dobson, and R. G. Griffin, Proc. Natl. Acad. Sci. U.S.A. 99, 16748 (2002)] to form amyloid fibrils in vitro in extended β-sheet conformations. For these reasons the latter is given considerable attention in the literature and studied also as an isolated fragment in water solution where both experimental and theoretical efforts have indicated the propensity of the system to form β turns or α helices, but is otherwise predominantly unstructured. Differing from previous computational studies that employed implicit
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
Kaprálová-Žďánská, Petra Ruth; Šmydke, Jan; Department of Radiation and Chemical Physics, Institute of Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, 182 21 Prague 8 ; Civiš, Svatopluk
2013-09-14
Recently optimized exponentially tempered Gaussian basis sets [P. R. Kapralova-Zdanska and J. Smydke, J. Chem. Phys. 138, 024105 (2013)] are employed in quantitative simulations of helium absorption cross-sections and two-photon excitation yields of doubly excited resonances. Linearly polarized half-infinite and Gaussian laser pulses at wavelengths 38–58 nm and large intensities up to 100 TW/cm{sup 2} are considered. The emphasis is laid on convergence of the results with respect to the quality of the Gaussian basis sets (typically limited by a number of partial waves, density, and spatial extent of the basis functions) as well as to the quality of the basis set of field-free states (typically limited by the maximum rotational quantum number and maximum excitation of the lower electron). Particular attention is paid to stability of the results with respect to varying complex scaling parameter. Moreover, the study of the dynamics is preceded by a thorough check of helium energies and oscillator strengths as they are obtained with the exponentially tempered Gaussian basis sets, being also compared with yet unpublished emission wavelengths measured in electric discharge experiments.
A mathematical basis for the design and design optimization of adaptive trusses in precision control
NASA Technical Reports Server (NTRS)
Das, S. K.; Utku, S.; Chen, G.-S.; Wada, B. K.
1991-01-01
A mathematical basis for the optimal design of adaptive trusses to be used in supporting precision equipment is provided. The general theory of adaptive structures is introduced, and the global optimization problem of placing a limited number, q, of actuators, so as to maximally achieve precision control and provide prestress, is stated. Two serialized optimization problems, namely, optimal actuator placement for prestress and optimal actuator placement for precision control, are addressed. In the case of prestressing, the computation of a 'desired' prestress is discussed, the interaction between actuators and redundants in conveying the prestress is shown in its mathematical form, and a methodology for arriving at the optimal placement of actuators and additional redundants is discussed. With regard to precision control, an optimal placement scheme (for q actuators) for maximum 'authority' over the precision points is suggested. The results of the two serialized optimization problems are combined to give a suboptimal solution to the global optimization problem. A method for improving this suboptimal actuator placement scheme by iteration is presented.
Parallel axes gear set optimization in two-parameter space
NASA Astrophysics Data System (ADS)
Theberge, Y.; Cardou, A.; Cloutier, L.
1991-05-01
This paper presents a method for optimal spur and helical gear transmission design that may be used in a computer aided design (CAD) approach. The design objective is generally taken as obtaining the most compact set for a given power input and gear ratio. A mixed design procedure is employed which relies both on heuristic considerations and computer capabilities. Strength and kinematic constraints are considered in order to define the domain of feasible designs. Constraints allowed include: pinion tooth bending strength, gear tooth bending strength, surface stress (resistance to pitting), scoring resistance, pinion involute interference, gear involute interference, minimum pinion tooth thickness, minimum gear tooth thickness, and profile or transverse contact ratio. A computer program was developed which allows the user to input the problem parameters, to select the calculation procedure, to see constraint curves in graphic display, to have an objective function level curve drawn through the design space, to point at a feasible design point and to have constraint values calculated at that point. The user can also modify some of the parameters during the design process.
European air quality modelled by CAMx including the volatility basis set scheme
NASA Astrophysics Data System (ADS)
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
Toward accurate thermochemical models for transition metals : G3large basis sets for atoms Sc-Zn.
Mayhall, N. J.; Raghavachari, K.; Redfern, P. C.; Curtiss, L. A.; Rassolov, V.; Indiana Univ.; Univ. of South Carolina
2008-04-01
An augmented valence triple-zeta basis set, referred to as G3Large, is reported for the first-row transition metal elements Sc through Zn. The basis set is constructed in a manner similar to the G3Large basis set developed previously for other elements (H-Ar, K, Ca, Ga-Kr) and used as a key component in Gaussian-3 theory. It is based on a contraction of a set of 15s13p5d Gaussian primitives to 8s7p3d, and also includes sets of f and g polarization functions, diffuse spd functions, and core df polarization functions. The basis set is evaluated with triples-augmented coupled cluster [CCSD(T)] and Brueckner orbital [BD(T)] methods for a small test set involving energies of atoms, atomic ions, and diatomic hydrides. It performs well for the low-lying s{yields}d excitation energies of atoms, atomic ionization energies, and the dissociation energies of the diatomic hydrides. The Brueckner orbital-based BD(T) method performs substantially better than Hartree-Fock-based CCSD(T) for molecules such as NiH, where the starting unrestricted Hartree-Fock wavefunction suffers from a high degree of spin contamination. Comparison with available data for geometries of transition metal hydrides also shows good agreement. A smaller basis set without core polarization functions, G3MP2Large, is also defined.
Dharmarajan, Venkatasubramanian; Lee, Jeong-Heon; Patel, Anamika; Skalnik, David G.; Cosgrove, Michael S.
2012-01-01
Translocations and amplifications of the mixed lineage leukemia-1 (MLL1) gene are associated with aggressive myeloid and lymphocytic leukemias in humans. MLL1 is a member of the SET1 family of histone H3 lysine 4 (H3K4) methyltransferases, which are required for transcription of genes involved in hematopoiesis and development. MLL1 associates with a subcomplex containing WDR5, RbBP5, Ash2L, and DPY-30 (WRAD), which together form the MLL1 core complex that is required for sequential mono- and dimethylation of H3K4. We previously demonstrated that WDR5 binds the conserved WDR5 interaction (Win) motif of MLL1 in vitro, an interaction that is required for the H3K4 dimethylation activity of the MLL1 core complex. In this investigation, we demonstrate that arginine 3765 of the MLL1 Win motif is required to co-immunoprecipitate WRAD from mammalian cells, suggesting that the WDR5-Win motif interaction is important for the assembly of the MLL1 core complex in vivo. We also demonstrate that peptides that mimic SET1 family Win motif sequences inhibit H3K4 dimethylation by the MLL1 core complex with varying degrees of efficiency. To understand the structural basis for these differences, we determined structures of WDR5 bound to six different naturally occurring Win motif sequences at resolutions ranging from 1.9 to 1.2 Å. Our results reveal that binding energy differences result from interactions between non-conserved residues C-terminal to the Win motif and to a lesser extent from subtle variation of residues within the Win motif. These results highlight a new class of methylation inhibitors that may be useful for the treatment of MLL1-related malignancies. PMID:22665483
Electron affinities of the first-row atoms revisited. Systematic basis sets and wave functions
Kendall, R.A.; Dunning, T.H. Jr. ); Harrison, R.J. )
1992-05-01
The calculation of accurate electron affinities (EAs) of atomic or molecular species is one of the most challenging tasks in quantum chemistry. We describe a reliable procedure for calculating the electron affinity of an atom and present results for hydrogen, boron, carbon, oxygen, and fluorine (hydrogen is included for completeness). This procedure involves the use of the recently proposed correlation-consistent basis sets augmented with functions to describe the more diffuse character of the atomic anion coupled with a straightforward, uniform expansion of the reference space for multireference singles and doubles configuration-interaction (MRSD-CI) calculations. Comparison with previous results and with corresponding full CI calculations are given. The most accurate EAs obtained from the MRSD-CI calculations are (with experimental values in parentheses) hydrogen 0.740 eV (0.754), boron 0.258 (0.277), carbon 1.245 (1.263), oxygen 1.384 (1.461), and fluorine 3.337 (3.401). The EAs obtained from the MR-SDCI calculations differ by less than 0.03 eV from those predicted by the full CI calculations.
Optimization of massive countermeasure design in complex rockfall settings
NASA Astrophysics Data System (ADS)
Agliardi, Federico; Crosta, Giovanni B.
2015-04-01
Rockfall protection is a major need in areas impended by subvertical rockwalls with complex 3D morphology and little or no talus to provide natural rockfall attenuation. The design of massive embankments, usually required to ensure such protection, is particularly difficult in complex rockfall settings, due to: widespread occurrence of rockfall sources; difficult characterization of size distribution and location of unstable volumes; variability of failure mechanisms; spatial scattering of rockfall trajectories; high expected kinetic energies. Moreover, rockwalls in complex lithological and structural settings are often prone to mass falls related to rock mass sector collapses. All these issues may hamper a safe application of classic embankment analysis approaches, using empirical rules or 2D-based height/energy statistics, and point to the need of integrated analyses of rock slope instability and rockfall runout in 3D. We explore the potential of combining advanced rock mass characterisation techniques and 3D rockfall modelling to support challenging countermeasure design at a site near Lecco (Southern Alps, Italy). Here subvertical cliffs up to 600 m high impend on a narrow (< 150 m) strip of flat land along the Como Lake shore. Rock is thickly bedded limestone (Dolomia Principale Fm) involved in a ENE-trending, S-verging kilometre-scale anticline fold. The spatial variability of bedding attitude and fracture intensity is strongly controlled by the geological structure, with individual block sizes varying in the range 0.2-15 m3. This results in spatially variable rockfall susceptibility and mechanisms, from single block falls to mass falls. Several rockfall events occurred between 1981 and 2010 motivated the design of slope benching and a massive embankment. To support reliable design verification and optimization we performed a 3D assessment of both rock slope instability and rockfall runout. We characterised fracture patterns and rock mass quality associated
Tellgren, Erik I; Reine, Simen S; Helgaker, Trygve
2012-07-14
Analytical integral evaluation is a central task of modern quantum chemistry. Here we present a general method for evaluating differentiated integrals over standard Gaussian and mixed Gaussian/plane-wave hybrid orbitals. The main idea is to have a representation of basis sets that is flexible enough to enable differentiated integrals to be reinterpreted as standard integrals over modified basis functions. As an illustration of the method, we report a very simple implementation of Hartree-Fock level geometrical derivatives in finite magnetic fields for gauge-origin independent atomic orbitals, within the London program. As a quantum-chemical application, we optimize the structure of helium clusters and some well-known covalently bound molecules (water, ammonia and benzene) subject to strong magnetic fields. PMID:22653039
The Scientific Basis of Uncertainty Factors Used in Setting Occupational Exposure Limits
Dankovic, D. A.; Naumann, B. D.; Maier, A.; Dourson, M. L.; Levy, L. S.
2015-01-01
The uncertainty factor concept is integrated into health risk assessments for all aspects of public health practice, including by most organizations that derive occupational exposure limits. The use of uncertainty factors is predicated on the assumption that a sufficient reduction in exposure from those at the boundary for the onset of adverse effects will yield a safe exposure level for at least the great majority of the exposed population, including vulnerable subgroups. There are differences in the application of the uncertainty factor approach among groups that conduct occupational assessments; however, there are common areas of uncertainty which are considered by all or nearly all occupational exposure limit-setting organizations. Five key uncertainties that are often examined include interspecies variability in response when extrapolating from animal studies to humans, response variability in humans, uncertainty in estimating a no-effect level from a dose where effects were observed, extrapolation from shorter duration studies to a full life-time exposure, and other insufficiencies in the overall health effects database indicating that the most sensitive adverse effect may not have been evaluated. In addition, a modifying factor is used by some organizations to account for other remaining uncertainties—typically related to exposure scenarios or accounting for the interplay among the five areas noted above. Consideration of uncertainties in occupational exposure limit derivation is a systematic process whereby the factors applied are not arbitrary, although they are mathematically imprecise. As the scientific basis for uncertainty factor application has improved, default uncertainty factors are now used only in the absence of chemical-specific data, and the trend is to replace them with chemical-specific adjustment factors whenever possible. The increased application of scientific data in the development of uncertainty factors for individual chemicals also
The Scientific Basis of Uncertainty Factors Used in Setting Occupational Exposure Limits.
Dankovic, D A; Naumann, B D; Maier, A; Dourson, M L; Levy, L S
2015-01-01
The uncertainty factor concept is integrated into health risk assessments for all aspects of public health practice, including by most organizations that derive occupational exposure limits. The use of uncertainty factors is predicated on the assumption that a sufficient reduction in exposure from those at the boundary for the onset of adverse effects will yield a safe exposure level for at least the great majority of the exposed population, including vulnerable subgroups. There are differences in the application of the uncertainty factor approach among groups that conduct occupational assessments; however, there are common areas of uncertainty which are considered by all or nearly all occupational exposure limit-setting organizations. Five key uncertainties that are often examined include interspecies variability in response when extrapolating from animal studies to humans, response variability in humans, uncertainty in estimating a no-effect level from a dose where effects were observed, extrapolation from shorter duration studies to a full life-time exposure, and other insufficiencies in the overall health effects database indicating that the most sensitive adverse effect may not have been evaluated. In addition, a modifying factor is used by some organizations to account for other remaining uncertainties-typically related to exposure scenarios or accounting for the interplay among the five areas noted above. Consideration of uncertainties in occupational exposure limit derivation is a systematic process whereby the factors applied are not arbitrary, although they are mathematically imprecise. As the scientific basis for uncertainty factor application has improved, default uncertainty factors are now used only in the absence of chemical-specific data, and the trend is to replace them with chemical-specific adjustment factors whenever possible. The increased application of scientific data in the development of uncertainty factors for individual chemicals also has
Optimization of global model composed of radial basis functions using the term-ranking approach
Cai, Peng; Tao, Chao Liu, Xiao-Jun
2014-03-15
A term-ranking method is put forward to optimize the global model composed of radial basis functions to improve the predictability of the model. The effectiveness of the proposed method is examined by numerical simulation and experimental data. Numerical simulations indicate that this method can significantly lengthen the prediction time and decrease the Bayesian information criterion of the model. The application to real voice signal shows that the optimized global model can capture more predictable component in chaos-like voice data and simultaneously reduce the predictable component (periodic pitch) in the residual signal.
Temperature dependence of source specific volatility basis sets for motor vehicle exhaust
NASA Astrophysics Data System (ADS)
Roy, Anirban; Choi, Yunsoo
2015-10-01
Recent work on emissions testing has focused on developing source specific volatility distributions which could be used to improve emissions inventories. One problem about these volatility profiles is that they are evaluated only at one temperature which is usually 298 K. This study uses a simple statistical model to evaluate the temperature dependence of the source-resolved volatility basis set, considering gasoline and diesel vehicle exhaust. The steps involved (a) fitting a distribution to the emissions data (b) evaluating the goodness of fit using a statistical test (c) updating the volatility bins using the Clausius-Clayperon equation; calculating the heats of vaporization of each volatility class using a regression model (d) assessing how the volatility of different VOC classes-Extremely Low Volatile, Low Volatile, Semi-Volatile, Intermediate Volatile and Volatile Organic Compounds - are affected by temperature. The results indicated that there could be significant changes in gas-particle partitioning of these emissions. For diesel exhaust at 298 K, the fractions are 5.4 × 10-4 (ELVOC), 0.074 (LVOC), 0.76 (SVOC), 0.17 (IVOC) and 10-5 (VOC) respectively. Looking at a window of ∓20 K, the partitioning for 278 K is 3 × 10-3 (ELVOC), 0.26 (LVOC), 0.67 (SVOC), 0.07 (IVOC) with no VOC fraction; while at 318 K it is 1.5 × 10-7 (ELVOC), 9 × 10-3 (LVOC), 0.64 (SVOC), 0.35 (IVOC) and 2 × 10-5 (VOC); demonstrating a significant change with temperature. The parameterizations developed in this work could be used to improve motor vehicle emissions inventory models such as MOVES.
Fast Electron Correlation Methods for Molecular Clusters without Basis Set Superposition Errors
Kamiya, Muneaki; Hirata, So; Valiev, Marat
2008-02-19
Two critical extensions to our fast, accurate, and easy-to-implement binary or ternary interaction method for weakly-interacting molecular clusters [Hirata et al. Mol. Phys. 103, 2255 (2005)] have been proposed, implemented, and applied to water hexamers, hydrogen fluoride chains and rings, and neutral and zwitterionic glycine–water clusters with an excellent result for an initial performance assessment. Our original method included up to two- or three-body Coulomb, exchange, and correlation energies exactly and higher-order Coulomb energies in the dipole–dipole approximation. In this work, the dipole moments are replaced by atom-centered point charges determined so that they reproduce the electrostatic potentials of the cluster subunits as closely as possible and also self-consistently with one another in the cluster environment. They have been shown to lead to dramatic improvement in the description of short-range electrostatic potentials not only of large, charge-separated subunits like zwitterionic glycine but also of small subunits. Furthermore, basis set superposition errors (BSSE) known to plague direct evaluation of weak interactions have been eliminated by com-bining the Valiron–Mayer function counterpoise (VMFC) correction with our binary or ternary interaction method in an economical fashion (quadratic scaling n2 with respect to the number of subunits n when n is small and linear scaling when n is large). A new variant of VMFC has also been proposed in which three-body and all higher-order Coulomb effects on BSSE are estimated approximately. The BSSE-corrected ternary interaction method with atom-centered point charges reproduces the VMFC-corrected results of conventional electron correlation calculations within 0.1 kcal/mol. The proposed method is significantly more accurate and also efficient than conventional correlation methods uncorrected of BSSE.
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
NASA Astrophysics Data System (ADS)
Kobus, J.; Moncrieff, D.; Wilson, S.
2004-02-01
In a previous paper, we have made a comparison of the accuracy with which the electric dipole polarizability agrzz and hyperpolarizability bgrzzz can be calculated by using either the finite basis set approach (the algebraic approximation) or the finite difference method in calculations for the ground states of the H2, LiH, BH and FH molecules, at their respective experimental equilibrium geometries, within the Hartree-Fock model. A re-examination of the hyperpolarizability of the BH molecule shows it to be very sensitive both to the choice of grid employed in the finite difference Hartree-Fock calculation and the construction of the basis set used in the matrix Hartree-Fock study. A new comparison of finite difference and finite basis set hyperpolarizabilities for the BH molecule is made, together with new calculations for the LiH and FH ground states.
Simulating the oxygen content of ambient organic aerosol with the 2D volatility basis set
NASA Astrophysics Data System (ADS)
Murphy, B. N.; Donahue, N. M.; Fountoukis, C.; Pandis, S. N.
2011-08-01
A module predicting the oxidation state of organic aerosol (OA) has been developed using the two-dimensional volatility basis set (2D-VBS) framework. This model is an extension of the 1D-VBS framework and tracks saturation concentration and oxygen content of organic species during their atmospheric lifetime. The host model, a one-dimensional Lagrangian transport model, is used to simulate air parcels arriving at Finokalia, Greece during the Finokalia Aerosol Measurement Experiment in May 2008 (FAME-08). Extensive observations were collected during this campaign using an aerosol mass spectrometer (AMS) and a thermodenuder to determine the chemical composition and volatility, respectively, of the ambient OA. Although there are several uncertain model parameters, the consistently high oxygen content of OA measured during FAME-08 (O:C = 0.8) can help constrain these parameters and elucidate OA formation and aging processes that are necessary for achieving the high degree of oxygenation observed. The base-case model reproduces observed OA mass concentrations (measured mean = 3.1 μg m-3, predicted mean = 3.3 μg m-3) and O:C (predicted O:C = 0.78) accurately. A suite of sensitivity studies explore uncertainties due to (1) the anthropogenic secondary OA (SOA) aging rate constant, (2) assumed enthalpies of vaporization, (3) the volatility change and number of oxygen atoms added for each generation of aging, (4) heterogeneous chemistry, (5) the oxidation state of the first generation of compounds formed from SOA precursor oxidation, and (6) biogenic SOA aging. Perturbations in most of these parameters do impact the ability of the model to predict O:C well throughout the simulation period. By comparing measurements of the O:C from FAME-08, several sensitivity cases including a high oxygenation case, a low oxygenation case, and biogenic SOA aging case are found to unreasonably depict OA aging, keeping in mind that this study does not consider possibly important processes
Simulating the oxygen content of ambient organic aerosol with the 2D volatility basis set
NASA Astrophysics Data System (ADS)
Murphy, B. N.; Donahue, N. M.; Fountoukis, C.; Pandis, S. N.
2011-03-01
A module predicting the oxidation state of organic aerosol (OA) has been developed using the two-dimensional volatility basis set (2D-VBS) framework. This model is an extension of the 1D-VBS framework and tracks saturation concentration and oxygen content of organic species during their atmospheric lifetime. The host model, a one-dimensional Lagrangian transport model, is used to simulate air parcels arriving at Finokalia, Greece during the Finokalia Aerosol Measurement Experiment in May 2008 (FAME-08). Extensive observations were collected during this campaign using an aerosol mass spectrometer (AMS) and a thermodenuder to determine the chemical composition and volatility, respectively, of the ambient OA. Although there are several uncertain model parameters, the consistently high oxygen content of OA measured during FAME-08 (O:C = 0.8) can help constrain these parameters and elucidate OA formation and aging processes that are necessary for achieving the high degree of oxygenation observed. The base-case model reproduces observed OA mass concentrations (measured mean = 3.1 μg m-3, predicted mean = 3.3 μg m-3) and O:C ratio (predicted O:C = 0.78) accurately. A suite of sensitivity studies explore uncertainties due to (1) the anthropogenic secondary OA (SOA) aging rate constant, (2) assumed enthalpies of vaporization, (3) the volatility change and number of oxygen atoms added for each generation of aging, (4) heterogeneous chemistry, (5) the oxidation state of the first generation of compounds formed from SOA precursor oxidation, and (6) biogenic SOA aging. Perturbations in most of these parameters do impact the ability of the model to predict O:C ratios well throughout the simulation period. By comparing measurements of the O:C ratio from FAME-08, several sensitivity cases including a high oxygenation case, low oxygenation case, and biogenic SOA aging case are found to unreasonably depict OA aging. However, many of the cases chosen for this study predict average
Finite-basis correction applied to the optimized effective potential within the FLAPW method
NASA Astrophysics Data System (ADS)
Friedrich, Christoph; Betzinger, Markus; Blügel, Stefan
2011-03-01
The optimized-effective-potential (OEP) method is a special technique to construct local exchange-correlation (xc) potentials from general orbital-dependent xc energy functionals for density-functional theory. Recently, we showed that particular care must be taken to construct local potentials within the all-electron full-potential augmented-plane-wave (FLAPW) approach. In fact, we found that the LAPW basis had to be converged to an accuracy that was far beyond that in calculations using conventional functionals, leading to a very high computational cost. This could be traced back to the convergence behavior of the density response function: only a highly converged basis lends the density enough flexibility to react adequately to changes of the potential. In this work we derive a numerical correction for the response function, which vanishes in the limit of an infinite, complete basis. It is constructed in the atomic spheres from the response of the basis functions themselves to changes of the potential. We show that such a finite-basis correction reduces the computational demand of OEP calculations considerably. We also discuss a similar correction scheme for GW calculations.
NASA Astrophysics Data System (ADS)
de Jong, W. A.; Harrison, R. J.; Dixon, D. A.
2001-01-01
A parallel implementation of the spin-free one-electron Douglas-Kroll-Hess (DKH) Hamiltonian in NWChem is discussed. An efficient and accurate method to calculate DKH gradients is introduced. It is shown that the use of a standard (nonrelativistic) contracted basis set can produce erroneous results for elements beyond the first row elements. The generation of DKH contracted cc-pVXZ(X=D,T,Q,5) basis sets for H, He, B-Ne, Al-Ar, and Ga-Br is discussed. The effect of DKH at the Hartree-Fock level on the bond distances, vibrational frequencies, and total dissociation energies for CF4, SiH4, SiF4, and Br2CO is discussed. It is suggested that the predominant effect of the scalar relativistic correction on the total dissociation energy can be calculated at the Hartree-Fock level if an adequate basis set is used.
A Complete Basis Set Estimate of Cation-p Bond Strengths: Na+(ethylene) and Na+(benzene)
Feller, David F.
2000-06-02
Large scale second order perturbation theory and couple cluster theory calculations were performed on the Na+(ethylene) and Na+(benzene) complexes in an effort to estimate binding enthalpies in the complete basis set limit. The resulting best estimates are DH0[Na+(ethylene)] = -13.7 ? 0.2 kcal/mol and DH0[Na+(benzene)] = -23.9 ? 0.3 kcal/mol, which include small corrections for core/valence correlation effects. The former value can be compared to a measurement of -10.3 ? 1.0 kcal/mol obtained from collision induced dissociation, while the latter value is approximately midway between the two existing experimental values which differed by 6.5 kcal/mol. For the basis sets considered in this study, the counterpoise-corrected binding energies were found to be in much worse agreement with the complete basis set limit than the raw values.
NASA Astrophysics Data System (ADS)
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.
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
NASA Astrophysics Data System (ADS)
Renaud, Rémi; Bendahmane, Mounir; Chery, Romain; Martin, Claire; Gurden, Hirac; Pain, Frederic
2012-06-01
Wide field multispectral imaging of light backscattered by brain tissues provides maps of hemodynamics changes (total blood volume and oxygenation) following activation. This technique relies on the fit of the reflectance images obtain at two or more wavelengths using a modified Beer-Lambert law1,2. It has been successfully applied to study the activation of several sensory cortices in the anesthetized rodent using visible light1-5. We have carried out recently the first multispectral imaging in the olfactory bulb6 (OB) of anesthetized rats. However, the optimization of wavelengths choice has not been discussed in terms of cross talk and uniqueness of the estimated parameters (blood volume and saturation maps) although this point was shown to be crucial for similar studies in Diffuse Optical Imaging in humans7-10. We have studied theoretically and experimentally the optimal sets of wavelength for multispectral imaging of rodent brain activation in the visible. Sets of optimal wavelengths have been identified and validated in vivo for multispectral imaging of the OB of rats following odor stimulus. We studied the influence of the wavelengths sets on the magnitude and time courses of the oxy- and deoxyhemoglobin concentration variations as well as on the spatial extent of activated brain areas following stimulation. Beyond the estimation of hemodynamic parameters from multispectral reflectance data, we observed repeatedly and for all wavelengths a decrease of light reflectance. For wavelengths longer than 590 nm, these observations differ from those observed in the somatosensory and barrel cortex and question the basis of the reflectance changes during activation in the OB. To solve this issue, Monte Carlo simulations (MCS) have been carried out to assess the relative contribution of absorption, scattering and anisotropy changes to the intrinsic optical imaging signals in somatosensory cortex (SsC) and OB model.
Brorsen, Kurt R.; Sirjoosingh, Andrew; Pak, Michael V.; Hammes-Schiffer, Sharon
2015-06-07
The nuclear electronic orbital (NEO) reduced explicitly correlated Hartree-Fock (RXCHF) approach couples select electronic orbitals to the nuclear orbital via Gaussian-type geminal functions. This approach is extended to enable the use of a restricted basis set for the explicitly correlated electronic orbitals and an open-shell treatment for the other electronic orbitals. The working equations are derived and the implementation is discussed for both extensions. The RXCHF method with a restricted basis set is applied to HCN and FHF{sup −} and is shown to agree quantitatively with results from RXCHF calculations with a full basis set. The number of many-particle integrals that must be calculated for these two molecules is reduced by over an order of magnitude with essentially no loss in accuracy, and the reduction factor will increase substantially for larger systems. Typically, the computational cost of RXCHF calculations with restricted basis sets will scale in terms of the number of basis functions centered on the quantum nucleus and the covalently bonded neighbor(s). In addition, the RXCHF method with an odd number of electrons that are not explicitly correlated to the nuclear orbital is implemented using a restricted open-shell formalism for these electrons. This method is applied to HCN{sup +}, and the nuclear densities are in qualitative agreement with grid-based calculations. Future work will focus on the significance of nonadiabatic effects in molecular systems and the further enhancement of the NEO-RXCHF approach to accurately describe such effects.
De Jong, Wibe A.; Harrison, Robert J.; Dixon, David A.
2001-01-01
A parallel implementation of the spin-free one-electron Douglas-Kroll(-Hess) Hamiltonian (DKH) in NWChem is discussed. An efficient and accurate method to calculate DKH gradients is introduced. It is shown that the use of standard (non-relativistic) contracted basis set can produce erroneous results for elements beyond the first row elements. The generation of DKH contracted cc-pVXZ (X = D, T, Q, 5) basis sets for H, He, B - Ne, Al - Ar, and Ga - Br will be discussed.
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
Optimal Experience among Campers in a Resident Camp Setting.
ERIC Educational Resources Information Center
Bialeschki, M. Deborah; Henderson, Karla A.
The purpose of this study was to assess optimal experience, also known as "flow" and "quality of experience" in a private, coeducational resident camp program for children. Flow refers to those times in work and leisure when people report feelings of enjoyment, concentration, and deep involvement. Flow theory predicts that an experience will be…
NASA Astrophysics Data System (ADS)
Petković, Dalibor; Gocic, Milan; Shamshirband, Shahaboddin; Qasem, Sultan Noman; Trajkovic, Slavisa
2016-08-01
Accurate estimation of the reference evapotranspiration (ET0) is important for the water resource planning and scheduling of irrigation systems. For this purpose, the radial basis function network with particle swarm optimization (RBFN-PSO) and radial basis function network with back propagation (RBFN-BP) were used in this investigation. The FAO-56 Penman-Monteith equation was used as reference equation to estimate ET0 for Serbia during the period of 1980-2010. The obtained simulation results confirmed the proposed models and were analyzed using the root mean-square error (RMSE), the mean absolute error (MAE), and the coefficient of determination ( R 2). The analysis showed that the RBFN-PSO had better statistical characteristics than RBFN-BP and can be helpful for the ET0 estimation.
NASA Astrophysics Data System (ADS)
Petković, Dalibor; Gocic, Milan; Shamshirband, Shahaboddin; Qasem, Sultan Noman; Trajkovic, Slavisa
2015-06-01
Accurate estimation of the reference evapotranspiration (ET0) is important for the water resource planning and scheduling of irrigation systems. For this purpose, the radial basis function network with particle swarm optimization (RBFN-PSO) and radial basis function network with back propagation (RBFN-BP) were used in this investigation. The FAO-56 Penman-Monteith equation was used as reference equation to estimate ET0 for Serbia during the period of 1980-2010. The obtained simulation results confirmed the proposed models and were analyzed using the root mean-square error (RMSE), the mean absolute error (MAE), and the coefficient of determination (R 2). The analysis showed that the RBFN-PSO had better statistical characteristics than RBFN-BP and can be helpful for the ET0 estimation.
Optimal continuous variable quantum teleportation protocol for realistic settings
NASA Astrophysics Data System (ADS)
Luiz, F. S.; Rigolin, Gustavo
2015-03-01
We show the optimal setup that allows Alice to teleport coherent states | α > to Bob giving the greatest fidelity (efficiency) when one takes into account two realistic assumptions. The first one is the fact that in any actual implementation of the continuous variable teleportation protocol (CVTP) Alice and Bob necessarily share non-maximally entangled states (two-mode finitely squeezed states). The second one assumes that Alice's pool of possible coherent states to be teleported to Bob does not cover the whole complex plane (| α | < ∞). The optimal strategy is achieved by tuning three parameters in the original CVTP, namely, Alice's beam splitter transmittance and Bob's displacements in position and momentum implemented on the teleported state. These slight changes in the protocol are currently easy to be implemented and, as we show, give considerable gain in performance for a variety of possible pool of input states with Alice.
MSE optimal bit-rate allocation in JPEG2000 Part 2 compression applied to a 3D data set
NASA Astrophysics Data System (ADS)
Kosheleva, Olga M.; Cabrera, Sergio D.; Usevitch, Bryan E.; Aguirre, Alberto; Vidal, Edward, Jr.
2004-10-01
A bit rate allocation (BRA) strategy is needed to optimally compress three-dimensional (3-D) data on a per-slice basis, treating it as a collection of two-dimensional (2-D) slices/components. This approach is compatible with the framework of JPEG2000 Part 2 which includes the option of pre-processing the slices with a decorrelation transform in the cross-component direction so that slices of transform coefficients are compressed. In this paper, we illustrate the impact of a recently developed inter-slice rate-distortion optimal bit-rate allocation approach that is applicable to this compression system. The approach exploits the MSE optimality of all JPEG2000 bit streams for all slices when each is produced in the quality progressive mode. Each bit stream can be used to produce a rate-distortion curve (RDC) for each slice that is MSE optimal at each bit rate of interest. The inter-slice allocation approach uses all RDCs for all slices to optimally select an overall optimal set of bit rates for all the slices using a constrained optimization procedure. The optimization is conceptually similar to Post-Compression Rate-Distortion optimization that is used within JPEG2000 to optimize bit rates allocated to codeblocks. Results are presented for two types of data sets: a 3-D computed tomography (CT) medical image, and a 3-D metereological data set derived from a particular modeling program. For comparison purposes, compression results are also illustrated for the traditional log-variance approach and for a uniform allocation strategy. The approach is illustrated using two decorrelation tranforms (the Karhunen Loeve transform, and the discrete wavelet transform) for which the inter-slice allocation scheme has the most impact.
Code of Federal Regulations, 2010 CFR
2010-10-01
... 42 Public Health 3 2010-10-01 2010-10-01 false Conditions for payment on a fee schedule basis for physician services in a teaching setting. 415.170 Section 415.170 Public Health CENTERS FOR MEDICARE & MEDICAID SERVICES, DEPARTMENT OF HEALTH AND HUMAN SERVICES (CONTINUED) MEDICARE PROGRAM SERVICES FURNISHED BY PHYSICIANS IN PROVIDERS,...
Liang, Y; Sun, W Q
2000-11-01
Recalcitrant seed axes were reported to survive to lower water contents under fast-drying conditions. The present study was to examine the hypothesis that drying rate and dehydration duration could interact to determine desiccation tolerance through different physico-chemical mechanisms. The effect of drying rate on desiccation tolerance of Theobroma cacao seed axes at 16 degrees C was examined. Rapid-drying at low relative humidity (RH) and slow-drying at high RH were more harmful to cocoa axes, because electrolyte leakage began to increase and axis viability began to decrease at high water contents. Maximum desiccation tolerance was observed with intermediate drying rates at RH between 88% and 91%, indicating the existence of an optimal drying rate or optimal desiccation duration. This maximum level of desiccation tolerance for cocoa axes (corresponding to a critical water potential of -9 MPa) was also detected using the equilibration method, in which axes were dehydrated over a series of salt solutions or glycerol solutions until the equilibrium. These data confirmed that the physiological basis of the optimal drying rate is related to both mechanical stress during desiccation and the length of desiccation duration during which deleterious reactions may occur. The optimal drying rate represents a situation where combined damages from mechanical and metabolic stresses become minimal. PMID:11113169
Holocene sea level variations on the basis of integration of independent data sets
Sahagian, D.; Berkman, P. . Dept. of Geological Sciences and Byrd Polar Research Center)
1992-01-01
Variations in sea level through earth history have occurred at a wide variety of time scales. Sea level researchers have attacked the problem of measuring these sea level changes through a variety of approaches, each relevant only to the time scale in question, and usually only relevant to the specific locality from which a specific type of data are derived. There is a plethora of different data types that can and have been used (locally) for the measurement of Holocene sea level variations. The problem of merging different data sets for the purpose of constructing a global eustatic sea level curve for the Holocene has not previously been adequately addressed. The authors direct the efforts to that end. Numerous studies have been published regarding Holocene sea level changes. These have involved exposed fossil reef elevations, elevation of tidal deltas, elevation of depth of intertidal peat deposits, caves, tree rings, ice cores, moraines, eolian dune ridges, marine-cut terrace elevations, marine carbonate species, tide gauges, and lake level variations. Each of these data sets is based on particular set of assumptions, and is valid for a specific set of environments. In order to obtain the most accurate possible sea level curve for the Holocene, these data sets must be merged so that local and other influences can be filtered out of each data set. Since each data set involves very different measurements, each is scaled in order to define the sensitivity of the proxy measurement parameter to sea level, including error bounds. This effectively determines the temporal and spatial resolution of each data set. The level of independence of data sets is also quantified, in order to rule out the possibility of a common non-eustatic factor affecting more than one variety of data. The Holocene sea level curve is considered to be independent of other factors affecting the proxy data, and is taken to represent the relation between global ocean water and basin volumes.
NASA Astrophysics Data System (ADS)
Witte, Jonathon; Neaton, Jeffrey B.; Head-Gordon, Martin
2016-05-01
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.
Optimal tooth numbers for compact standard spur gear sets
NASA Technical Reports Server (NTRS)
Savage, M.; Coy, J. J.; Townsend, D. P.
1981-01-01
The design of a standard gear mesh is treated with the objective of minimizing the gear size for a given ratio, pinion torque, and allowable tooth strength. Scoring, pitting fatigue, bending fatigue, and the kinematic limits of contact ratio and interference are considered. A design space is defined in terms of the number of teeth on the pinion and the diametral pitch. This space is then combined with the objective function of minimum center distance to obtain an optimal design region. This region defines the number of pinion teeth for the most compact design. The number is a function of the gear ratio only. A design example illustrating this procedure is also given.
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
Optimizing distance-based methods for large data sets
NASA Astrophysics Data System (ADS)
Scholl, Tobias; Brenner, Thomas
2015-10-01
Distance-based methods for measuring spatial concentration of industries have received an increasing popularity in the spatial econometrics community. However, a limiting factor for using these methods is their computational complexity since both their memory requirements and running times are in {{O}}(n^2). In this paper, we present an algorithm with constant memory requirements and shorter running time, enabling distance-based methods to deal with large data sets. We discuss three recent distance-based methods in spatial econometrics: the D&O-Index by Duranton and Overman (Rev Econ Stud 72(4):1077-1106, 2005), the M-function by Marcon and Puech (J Econ Geogr 10(5):745-762, 2010) and the Cluster-Index by Scholl and Brenner (Reg Stud (ahead-of-print):1-15, 2014). Finally, we present an alternative calculation for the latter index that allows the use of data sets with millions of firms.
Optimal electrocardiographic limb lead set for rapid emphysema screening
Bajaj, Rishi; Chhabra, Lovely; Basheer, Zainab; Spodick, David H
2013-01-01
Background Pulmonary emphysema of any etiology has been shown to be strongly and quasidiagnostically associated with a vertical frontal P wave axis. A vertical P wave axis (>60 degrees) during sinus rhythm can be easily determined by a P wave in lead III greater than the P wave in lead I (bipolar lead set) or a dominantly negative P wave in aVL (unipolar lead set). The purpose of this investigation was to determine which set of limb leads may be better for identifying the vertical P vector of emphysema in adults. Methods Unselected consecutive electrocardiograms from 100 patients with a diagnosis of emphysema were analyzed to determine the P wave axis. Patients aged younger than 45 years, those not in sinus rhythm, and those with poor quality tracings were excluded. The electrocardiographic data were divided into three categories depending on the frontal P wave axis, ie, >60 degrees, 60 degrees, or <60 degrees, by each criterion (P amplitude lead III > lead I and a negative P wave in aVL). Results Sixty-six percent of patients had a P wave axis > 60 degrees based on aVL, and 88% of patients had a P wave axis > 60 degrees based on the P wave in lead III being greater than in lead I. Conclusion A P wave in lead III greater than that in lead I is a more sensitive marker than a negative P wave in aVL for diagnosing emphysema and is recommended for rapid routine screening. PMID:23378754
NASA Astrophysics Data System (ADS)
Ma, Zhonghua; Zhang, Yanli; Tuckerman, Mark E.
2012-07-01
It is generally believed that studies of liquid water using the generalized gradient approximation to density functional theory require dispersion corrections in order to obtain reasonably accurate structural and dynamical properties. Here, we report on an ab initio molecular dynamics study of water in the isothermal-isobaric ensemble using a converged discrete variable representation basis set and an empirical dispersion correction due to Grimme [J. Comp. Chem. 27, 1787 (2006)], 10.1002/jcc.20495. At 300 K and an applied pressure of 1 bar, the density obtained without dispersion corrections is approximately 0.92 g/cm3 while that obtained with dispersion corrections is 1.07 g/cm3, indicating that the empirical dispersion correction overestimates the density by almost as much as it is underestimated without the correction for this converged basis. Radial distribution functions exhibit a loss of structure in the second solvation shell. Comparison of our results with other studies using the same empirical correction suggests the cause of the discrepancy: the Grimme dispersion correction is parameterized for use with a particular basis set; this parameterization is sensitive to this choice and, therefore, is not transferable to other basis sets.
Perturbing engine performance measurements to determine optimal engine control settings
Jiang, Li; Lee, Donghoon; Yilmaz, Hakan; Stefanopoulou, Anna
2014-12-30
Methods and systems for optimizing a performance of a vehicle engine are provided. The method includes determining an initial value for a first engine control parameter based on one or more detected operating conditions of the vehicle engine, determining a value of an engine performance variable, and artificially perturbing the determined value of the engine performance variable. The initial value for the first engine control parameter is then adjusted based on the perturbed engine performance variable causing the engine performance variable to approach a target engine performance variable. Operation of the vehicle engine is controlled based on the adjusted initial value for the first engine control parameter. These acts are repeated until the engine performance variable approaches the target engine performance variable.
Solving the Dirac equation, using the large component only, in a Dirac-type Slater orbital basis set
NASA Astrophysics Data System (ADS)
van Lenthe, E.; Baerends, E. J.; Snijders, J. G.
1995-04-01
We solve the Dirac equation by solving the two-component energy-dependent equation for the large component that results from the elimination of the small component. This requires for every occupied orbital the diagonalization of a Hamiltonian. Advantages are, however, that these Hamiltonians are all bounded from below, unlike the Dirac Hamiltonian, and that only a basis set for the large component is needed. We use Dirac-type Slater orbitals, adapted from solutions to the hydrogen-like atom. This offers the perspective of performing relativistic calculations to the same accuracy as non-relativistic ones, with a comparable number of basis functions.
Non-homogeneous solutions of a Coulomb Schrödinger equation as basis set for scattering problems
Del Punta, J. A.; Ambrosio, M. J.; Gasaneo, G.; Zaytsev, S. A.; Ancarani, L. U.
2014-05-15
We introduce and study two-body Quasi Sturmian functions which are proposed as basis functions for applications in three-body scattering problems. They are solutions of a two-body non-homogeneous Schrödinger equation. We present different analytic expressions, including asymptotic behaviors, for the pure Coulomb potential with a driven term involving either Slater-type or Laguerre-type orbitals. The efficiency of Quasi Sturmian functions as basis set is numerically illustrated through a two-body scattering problem.
A decision support system to determine optimal ventilator settings
2014-01-01
Background Choosing the correct ventilator settings for the treatment of patients with respiratory tract disease is quite an important issue. Since the task of specifying the parameters of ventilation equipment is entirely carried out by a physician, physician’s knowledge and experience in the selection of these settings has a direct effect on the accuracy of his/her decisions. Nowadays, decision support systems have been used for these kinds of operations to eliminate errors. Our goal is to minimize errors in ventilation therapy and prevent deaths caused by incorrect configuration of ventilation devices. The proposed system is designed to assist less experienced physicians working in the facilities without having lung mechanics like cottage hospitals. Methods This article describes a decision support system proposing the ventilator settings required to be applied in the treatment according to the patients’ physiological information. The proposed model has been designed to minimize the possibility of making a mistake and to encourage more efficient use of time in support of the decision making process while the physicians make critical decisions about the patient. Artificial Neural Network (ANN) is implemented in order to calculate frequency, tidal volume, FiO2 outputs, and this classification model has been used for estimation of pressure support / volume support outputs. For the obtainment of the highest performance in both models, different configurations have been tried. Various tests have been realized for training methods, and a number of hidden layers mostly affect factors regarding the performance of ANNs. Results The physiological information of 158 respiratory patients over the age of 60 and were treated in three different hospitals between the years 2010 and 2012 has been used in the training and testing of the system. The diagnosed disease, core body temperature, pulse, arterial systolic pressure, diastolic blood pressure, PEEP, PSO2, pH, pCO2
NASA Astrophysics Data System (ADS)
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.
NASA Astrophysics Data System (ADS)
Zhang, Jun; Dolg, Michael
2014-01-01
The third-order incremental dual-basis set zero-buffer approach was combined with CCSD(T)-F12x (x = a, b) theory to develop a new approach, i.e., the inc3-db-B0-CCSD(T)-F12 method, which can be applied as a black-box procedure to efficiently obtain the near complete basis set (CBS) limit of the CCSD(T) energies also for large systems. We tested this method for several cases of different chemical nature: four complexes taken from the standard benchmark sets S66 and X40, the energy difference between isomers of water hexamer and the rotation barrier of biphenyl. The results show that our method has an error relative to the best estimation of CBS energy of only 0.2 kcal/mol or less. By parallelization, our method can accomplish the CCSD(T)-F12 calculations of about 60 correlated electrons and 800 basis functions in only several days, which by standard implementation are impossible for ordinary hardware. We conclude that the inc3-db-B0-CCSD(T)-F12a/AVTZ method, which is of CCSD(T)/AV5Z quality, is close to the limit of accuracy that one can achieve for large systems currently.
Zhang, Jun Dolg, Michael
2014-01-28
The third-order incremental dual-basis set zero-buffer approach was combined with CCSD(T)-F12x (x = a, b) theory to develop a new approach, i.e., the inc3-db-B0-CCSD(T)-F12 method, which can be applied as a black-box procedure to efficiently obtain the near complete basis set (CBS) limit of the CCSD(T) energies also for large systems. We tested this method for several cases of different chemical nature: four complexes taken from the standard benchmark sets S66 and X40, the energy difference between isomers of water hexamer and the rotation barrier of biphenyl. The results show that our method has an error relative to the best estimation of CBS energy of only 0.2 kcal/mol or less. By parallelization, our method can accomplish the CCSD(T)-F12 calculations of about 60 correlated electrons and 800 basis functions in only several days, which by standard implementation are impossible for ordinary hardware. We conclude that the inc3-db-B0-CCSD(T)-F12a/AVTZ method, which is of CCSD(T)/AV5Z quality, is close to the limit of accuracy that one can achieve for large systems currently.
Optimizing the application of entomopathogenic nematodes: experimental set-up.
Brusselman, E; Steurbaut, W; Sonck, B
2007-01-01
The complex issue concerning the spray application of Entomopathogenic Nematodes (EPNs) with a hydraulic sprayer is still not solved. This research project focuses on the effect of spray application technique on the viability and deposition of EPNs. In this paper the experimental set-up used for this evaluation is described. A modular spray application system has been developed and is currently used to evaluate the effect of different parts of a sprayer on the viability of the EPNs. Based on the results of experiments using this modular spray application system, recommendations regarding pump type, mixing system, nozzle type and filter size will be formulated. Because of the large number of experiments in this research project, an image analysis system for the determination of the viability of the nematodes is developed. This paper describes two experiments comparing the new developed image processing technique with the standard microscopic counting technique. PMID:18399427
NASA Astrophysics Data System (ADS)
Purwanto, Wirawan; Krakauer, Henry; Zhang, Shiwei; Virgus, Yudistira
2011-03-01
Weak H2 physisorption energies present a significant challenge to first-principle theoretical modeling and prediction of materials for H storage. There has been controversy regarding the accuracy of DFT on systems involving Ca cations. We use the auxiliary-field quantum Monte Carlo (AFQMC) method to accurately predict the binding energy of Ca + , - 4{H}2 . AFQMC scales as Nbasis3and has demonstrated accuracy similar to or better than the gold-standard coupled cluster CCSD(T) method. We apply a modified Cholesky decomposition to achieve efficient Hubbard-Stratonovich transformation in AFQMC at large basis sizes. We employ the largest correlation consistent basis sets available, up to Ca/cc-pCV5Z, to extrapolate to the complete basis limit. The calculated potential energy curve exhibits binding with a double-well structure. Supported by DOE and NSF. Calculations were performed at OLCF Jaguar and CPD.
Liquid Water through Density-Functional Molecular Dynamics: Plane-Wave vs Atomic-Orbital Basis Sets.
Miceli, Giacomo; Hutter, Jürg; Pasquarello, Alfredo
2016-08-01
We determine and compare structural, dynamical, and electronic properties of liquid water at near ambient conditions through density-functional molecular dynamics simulations, when using either plane-wave or atomic-orbital basis sets. In both frameworks, the electronic structure and the atomic forces are self-consistently determined within the same theoretical scheme based on a nonlocal density functional accounting for van der Waals interactions. The overall properties of liquid water achieved within the two frameworks are in excellent agreement with each other. Thus, our study supports that implementations with plane-wave or atomic-orbital basis sets yield equivalent results and can be used indiscriminately in study of liquid water or aqueous solutions. PMID:27434607
NASA Technical Reports Server (NTRS)
Almloef, Jan; Taylor, Peter R.
1989-01-01
A recently proposed scheme for using natural orbitals from atomic configuration interaction (CI) wave functions as a basis set for linear combination of atomic orbitals (LCAO) calculations is extended for the calculation of molecular properties. For one-electron properties like multipole moments, which are determined largely by the outermost regions of the molecular wave function, it is necessary to increase the flexibility of the basis in these regions. This is most easily done by uncontracting the outmost Gaussian primitives, and/or by adding diffuse primitives. A similar approach can be employed for the calculation of polarizabilities. Properties which are not dominated by the long-range part of the wave function, such as spectroscopic constants or electric field gradients at the nucleus, can generally be treated satisfactorily with the original atomic natural orbital (ANO) sets.
NASA Technical Reports Server (NTRS)
Almlof, Jan; Taylor, Peter R.
1990-01-01
A recently proposed scheme for using natural orbitals from atomic configuration interaction wave functions as a basis set for linear combination of atomic orbitals (LCAO) calculations is extended for the calculation of molecular properties. For one-electron properties like multipole moments, which are determined largely by the outermost regions of the molecular wave function, it is necessary to increase the flexibility of the basis in these regions. This is most easily done by uncontracting the outermost Gaussian primitives, and/or by adding diffuse primitives. A similar approach can be employed for the calculation of polarizabilities. Properties which are not dominated by the long-range part of the wave function, such as spectroscopic constants or electric field gradients at the nucleus, can generally be treated satisfactorily with the original atomic natural orbital sets.
NASA Astrophysics Data System (ADS)
Jorge, F. E.; Martins, L. S. C.; Franco, M. L.
2016-01-01
Segmented all-electron basis sets of valence double zeta quality plus polarization functions (DZP) for the elements from Ce to Lu are generated to be used with the non-relativistic and Douglas-Kroll-Hess (DKH) Hamiltonians. At the B3LYP level, the DZP-DKH atomic ionization energies and equilibrium bond lengths and atomization energies of the lanthanide trifluorides are evaluated and compared with benchmark theoretical and experimental data reported in the literature. In general, this compact size set shows to have a regular, efficient, and reliable performance. It can be particularly useful in molecular property calculations that require explicit treatment of the core electrons.
NASA Astrophysics Data System (ADS)
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.
Selection of optimal muscle set for 16-channel standing neuroprosthesis
Gartman, Steven J.; Audu, Musa L.; Kirsch, Robert F.; Triolo, Ronald J.
2009-01-01
The Case Western Reserve University/Department of Veterans Affairs 8-channel lower-limb neuroprosthesis can restore standing to selected individuals with paraplegia by application of functional electrical stimulation. The second generation of this system will include 16 channels of stimulation and a closed-loop control scheme to provide automatic postural corrections. This study used a musculoskeletal model of the legs and trunk to determine which muscles to target with the new system in order to maximize the range of postures that can be statically maintained, which should increase the system’s ability to provide adequate support to maintain standing when the user’s posture moves away from a neutral stance, either by an external disturbance or a volitional change in posture by the user. The results show that the prime muscle targets should be the medial gastrocnemius, tibialis anterior, vastus lateralis, semimembranosus, gluteus maximus, gluteus medius, adductor magnus, and erector spinae. This set of 16 muscles supports 42 percent of the standing postures that are attainable by the nondisabled model. Coactivation of the lateral gastrocnemius and peroneus longus with the medial gastrocnemius and of the peroneus tertius with the tibialis anterior increased the percentage of feasible postures to 71 percent. PMID:16847793
NASA Astrophysics Data System (ADS)
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.
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. PMID:20136303
Řezáč, Jan; de la Lande, Aurélien
2015-02-10
Separation of the energetic contribution of charge transfer to interaction energy in noncovalent complexes would provide important insight into the mechanisms of the interaction. However, the calculation of charge-transfer energy is not an easy task. It is not a physically well-defined term, and the results might depend on how it is described in practice. Commonly, the charge transfer is defined in terms of molecular orbitals; in this framework, however, the charge transfer vanishes as the basis set size increases toward the complete basis set limit. This can be avoided by defining the charge transfer in terms of the spatial extent of the electron densities of the interacting molecules, but the schemes used so far do not reflect the actual electronic structure of each particular system and thus are not reliable. We propose a spatial partitioning of the system, which is based on a charge transfer-free reference state, namely superimposition of electron densities of the noninteracting fragments. We show that this method, employing constrained DFT for the calculation of the charge-transfer energy, yields reliable results and is robust with respect to the strength of the charge transfer, the basis set size, and the DFT functional used. Because it is based on DFT, the method is applicable to rather large systems. PMID:26580910
Matrix-product-state method with local basis optimization for nonequilibrium electron-phonon systems
NASA Astrophysics Data System (ADS)
Heidrich-Meisner, Fabian; Brockt, Christoph; Dorfner, Florian; Vidmar, Lev; Jeckelmann, Eric
We present a method for simulating the time evolution of quasi-one-dimensional correlated systems with strongly fluctuating bosonic degrees of freedom (e.g., phonons) using matrix product states. For this purpose we combine the time-evolving block decimation (TEBD) algorithm with a local basis optimization (LBO) approach. We discuss the performance of our approach in comparison to TEBD with a bare boson basis, exact diagonalization, and diagonalization in a limited functional space. TEBD with LBO can reduce the computational cost by orders of magnitude when boson fluctuations are large and thus it allows one to investigate problems that are out of reach of other approaches. First, we test our method on the non-equilibrium dynamics of a Holstein polaron and show that it allows us to study the regime of strong electron-phonon coupling. Second, the method is applied to the scattering of an electronic wave packet off a region with electron-phonon coupling. Our study reveals a rich physics including transient self-trapping and dissipation. Supported by Deutsche Forschungsgemeinschaft (DFG) via FOR 1807.
[The BASYS observation system for the analysis of aggressive behavior in classroom-settings].
Wettstein, Alexander
2012-01-01
Educational or therapeutic measures of aggressive student behavior are often based on the judgments of teachers. However, empirical studies show that the objectivity of these judgments is generally low. In order to assess aggressive behavior in classroom settings, we developed a context-sensitive observational system. The observation system exists in a version for teachers in action as well as a version for the uninvolved observer. The teacher version allows categorizing aggressive behavior while teaching. The aim is to differentiate the perception and the judgments of teachers, so that the judgments can serve as trustable diagnostic information. The version for an independent observer, in addition, contains categories to collect information about the context in which aggressions take place. The behavior observation system was tested in four field-studies in regular and special classes. The empirical results show that, after training, teachers were able to make objective observations, and that aggressive behavior depends to a large extent on situational factors. The system allows identification of problematic people-environment relationships and the derivation of intervention measures. PMID:22748725
Ermler, Walter V.; Tilson, Jeffrey L.
2012-12-15
A procedure for structuring generally contracted valence-core/valence basis sets of Gaussian-type functions for use with relativistic effective core potentials (gcv-c/v-RECP basis sets) is presented. Large valence basis sets are enhanced using a compact basis set derived for outer core electrons in the presence of small-core RECPs. When core electrons are represented by relativistic effective core potentials (RECPs), and appropriate levels of theory, these basis sets are shown to provide accurate representations of atomic and molecular valence and outer-core electrons. Core/valence polarization and correlation effects can be calculated using these basis sets through standard methods for treating electron correlation. Calculations of energies and spectra for Ru, Os, Ir, In and Cs are reported. Spectroscopic constants for RuO2+, OsO2+, Cs2 and InH are calculated and compared with experiment.
The Study of the Optimal Parameter Settings in a Hospital Supply Chain System in Taiwan
Liao, Hung-Chang; Chen, Meng-Hao; Wang, Ya-huei
2014-01-01
This study proposed the optimal parameter settings for the hospital supply chain system (HSCS) when either the total system cost (TSC) or patient safety level (PSL) (or both simultaneously) was considered as the measure of the HSCS's performance. Four parameters were considered in the HSCS: safety stock, maximum inventory level, transportation capacity, and the reliability of the HSCS. A full-factor experimental design was used to simulate an HSCS for the purpose of collecting data. The response surface method (RSM) was used to construct the regression model, and a genetic algorithm (GA) was applied to obtain the optimal parameter settings for the HSCS. The results show that the best method of obtaining the optimal parameter settings for the HSCS is the simultaneous consideration of both the TSC and the PSL to measure performance. Also, the results of sensitivity analysis based on the optimal parameter settings were used to derive adjustable strategies for the decision-makers. PMID:25250397
Kruse, Holger; Grimme, Stefan
2012-04-21
chemistry yields MAD=0.68 kcal/mol, which represents a huge improvement over plain B3LYP/6-31G* (MAD=2.3 kcal/mol). Application of gCP-corrected B97-D3 and HF-D3 on a set of large protein-ligand complexes prove the robustness of the method. Analytical gCP gradients make optimizations of large systems feasible with small basis sets, as demonstrated for the inter-ring distances of 9-helicene and most of the complexes in Hobza's S22 test set. The method is implemented in a freely available FORTRAN program obtainable from the author's website. PMID:22519309
Full Glowworm Swarm Optimization Algorithm for Whole-Set Orders Scheduling in Single Machine
Yu, Zhang; Yang, Xiaomei
2013-01-01
By analyzing the characteristics of whole-set orders problem and combining the theory of glowworm swarm optimization, a new glowworm swarm optimization algorithm for scheduling is proposed. A new hybrid-encoding schema combining with two-dimensional encoding and random-key encoding is given. In order to enhance the capability of optimal searching and speed up the convergence rate, the dynamical changed step strategy is integrated into this algorithm. Furthermore, experimental results prove its feasibility and efficiency. PMID:24294135
Analysis of the optimal laminated target made up of discrete set of materials
NASA Technical Reports Server (NTRS)
Aptukov, Valery N.; Belousov, Valentin L.
1991-01-01
A new class of problems was analyzed to estimate an optimal structure of laminated targets fabricated from the specified set of homogeneous materials. An approximate description of the perforation process is based on the model of radial hole extension. The problem is solved by using the needle-type variation technique. The desired optimization conditions and quantitative/qualitative estimations of optimal targets were obtained and are discussed using specific examples.
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
Structural basis for inhibition of the histone chaperone activity of SET/TAF-Iβ by cytochrome c
González-Arzola, Katiuska; Díaz-Moreno, Irene; Cano-González, Ana; Díaz-Quintana, Antonio; Velázquez-Campoy, Adrián; Moreno-Beltrán, Blas; López-Rivas, Abelardo; De la Rosa, Miguel A.
2015-01-01
Chromatin is pivotal for regulation of the DNA damage process insofar as it influences access to DNA and serves as a DNA repair docking site. Recent works identify histone chaperones as key regulators of damaged chromatin’s transcriptional activity. However, understanding how chaperones are modulated during DNA damage response is still challenging. This study reveals that the histone chaperone SET/TAF-Iβ interacts with cytochrome c following DNA damage. Specifically, cytochrome c is shown to be translocated into cell nuclei upon induction of DNA damage, but not upon stimulation of the death receptor or stress-induced pathways. Cytochrome c was found to competitively hinder binding of SET/TAF-Iβ to core histones, thereby locking its histone-binding domains and inhibiting its nucleosome assembly activity. In addition, we have used NMR spectroscopy, calorimetry, mutagenesis, and molecular docking to provide an insight into the structural features of the formation of the complex between cytochrome c and SET/TAF-Iβ. Overall, these findings establish a framework for understanding the molecular basis of cytochrome c-mediated blocking of SET/TAF-Iβ, which subsequently may facilitate the development of new drugs to silence the oncogenic effect of SET/TAF-Iβ’s histone chaperone activity. PMID:26216969
Structural basis for inhibition of the histone chaperone activity of SET/TAF-Iβ by cytochrome c.
González-Arzola, Katiuska; Díaz-Moreno, Irene; Cano-González, Ana; Díaz-Quintana, Antonio; Velázquez-Campoy, Adrián; Moreno-Beltrán, Blas; López-Rivas, Abelardo; De la Rosa, Miguel A
2015-08-11
Chromatin is pivotal for regulation of the DNA damage process insofar as it influences access to DNA and serves as a DNA repair docking site. Recent works identify histone chaperones as key regulators of damaged chromatin's transcriptional activity. However, understanding how chaperones are modulated during DNA damage response is still challenging. This study reveals that the histone chaperone SET/TAF-Iβ interacts with cytochrome c following DNA damage. Specifically, cytochrome c is shown to be translocated into cell nuclei upon induction of DNA damage, but not upon stimulation of the death receptor or stress-induced pathways. Cytochrome c was found to competitively hinder binding of SET/TAF-Iβ to core histones, thereby locking its histone-binding domains and inhibiting its nucleosome assembly activity. In addition, we have used NMR spectroscopy, calorimetry, mutagenesis, and molecular docking to provide an insight into the structural features of the formation of the complex between cytochrome c and SET/TAF-Iβ. Overall, these findings establish a framework for understanding the molecular basis of cytochrome c-mediated blocking of SET/TAF-Iβ, which subsequently may facilitate the development of new drugs to silence the oncogenic effect of SET/TAF-Iβ's histone chaperone activity. PMID:26216969
Sakhanenko, Nikita A.; Skupin, Alexander; Ignac, Tomasz
2014-01-01
Abstract Context dependence is central to the description of complexity. Keying on the pairwise definition of “set complexity,” we use an information theory approach to formulate general measures of systems complexity. We examine the properties of multivariable dependency starting with the concept of interaction information. We then present a new measure for unbiased detection of multivariable dependency, “differential interaction information.” This quantity for two variables reduces to the pairwise “set complexity” previously proposed as a context-dependent measure of information in biological systems. We generalize it here to an arbitrary number of variables. Critical limiting properties of the “differential interaction information” are key to the generalization. This measure extends previous ideas about biological information and provides a more sophisticated basis for the study of complexity. The properties of “differential interaction information” also suggest new approaches to data analysis. Given a data set of system measurements, differential interaction information can provide a measure of collective dependence, which can be represented in hypergraphs describing complex system interaction patterns. We investigate this kind of analysis using simulated data sets. The conjoining of a generalized set complexity measure, multivariable dependency analysis, and hypergraphs is our central result. While our focus is on complex biological systems, our results are applicable to any complex system. PMID:24377753
Zhao, Bin; Wang, Shuxiao; Donahue, Neil M; Chuang, Wayne; Hildebrandt Ruiz, Lea; Ng, Nga L; Wang, Yangjun; Hao, Jiming
2015-02-17
We evaluate the one-dimensional volatility basis set (1D-VBS) and two-dimensional volatility basis set (2D-VBS) in simulating the aging of SOA derived from toluene and α-pinene against smog-chamber experiments. If we simulate the first-generation products with empirical chamber fits and the subsequent aging chemistry with a 1D-VBS or a 2D-VBS, the models mostly overestimate the SOA concentrations in the toluene oxidation experiments. This is because the empirical chamber fits include both first-generation oxidation and aging; simulating aging in addition to this results in double counting of the initial aging effects. If the first-generation oxidation is treated explicitly, the base-case 2D-VBS underestimates the SOA concentrations and O:C increase of the toluene oxidation experiments; it generally underestimates the SOA concentrations and overestimates the O:C increase of the α-pinene experiments. With the first-generation oxidation treated explicitly, we could modify the 2D-VBS configuration individually for toluene and α-pinene to achieve good model-measurement agreement. However, we are unable to simulate the oxidation of both toluene and α-pinene with the same 2D-VBS configuration. We suggest that future models should implement parallel layers for anthropogenic (aromatic) and biogenic precursors, and that more modeling studies and laboratory research be done to optimize the "best-guess" parameters for each layer. PMID:25581402
Shcherbinin, Konstantin
2011-10-01
The design requirements for physical protection systems currently adopted at Ukrainian NPPs were established in the middle of the nineties on the basis of deterministic analyses and expert evaluations of the vulnerabilities of reactor facilities. At the present time the in-depth assessment of the nuclear safety of Ukrainian NPPs with VVER-1000/B320 reactors has been completed using Probabilistic Safety Assessment (PSA). PSA has established and provided a qualitative assessment of the significance of the equipment for maintaining the integrity of a reactor core and preventing an abnormal radioactive release. The availability of qualitative assessments of the importance of equipment for nuclear safety allows one to assess the existing physical protection system using (1) comparative analysis: to determine whether all equipment and zones that may be affected, as established by the nuclear safety assessment, are actually included into the vital zones protected by the existing physical protection system; (2) specific analysis of dominant contributors: since nuclear safety analyses provide qualitative assessment of the equipment’s importance for safety, it is easy to select a limited group of essential equipment that makes the major contribution to safety, and (3) specific equipment analysis included in dominant emergency sets: part of the components might not be included in the essential equipment group, but it might be included in the dominant emergency sets. If some equipment is found not to be covered by the physical protection system, it is possible, using qualitative assessment of the importance of this equipment for safety, to assess the required degree of enhancement of physical protection. Such analysis will permit an assessment of the sufficiency of the existing physical protection system, to define the "tight" system areas and, therefore, to develop a justified optimization of a PPS with the objective of implementing priority upgrades to enhance safety, and
An Optimal Set of Flesh Points on Tongue and Lips for Speech-Movement Classification
ERIC Educational Resources Information Center
Wang, Jun; Samal, Ashok; Rong, Panying; Green, Jordan R.
2016-01-01
Purpose: The authors sought to determine an optimal set of flesh points on the tongue and lips for classifying speech movements. Method: The authors used electromagnetic articulographs (Carstens AG500 and NDI Wave) to record tongue and lip movements from 13 healthy talkers who articulated 8 vowels, 11 consonants, a phonetically balanced set of…
NASA Astrophysics Data System (ADS)
Feller, David; Peterson, Kirk A.
2013-08-01
The effectiveness of the recently developed, explicitly correlated coupled cluster method CCSD(T)-F12b is examined in terms of its ability to reproduce atomization energies derived from complete basis set extrapolations of standard CCSD(T). Most of the standard method findings were obtained with aug-cc-pV7Z or aug-cc-pV8Z basis sets. For a few homonuclear diatomic molecules it was possible to push the basis set to the aug-cc-pV9Z level. F12b calculations were performed with the cc-pVnZ-F12 (n = D, T, Q) basis set sequence and were also extrapolated to the basis set limit using a Schwenke-style, parameterized formula. A systematic bias was observed in the F12b method with the (VTZ-F12/VQZ-F12) basis set combination. This bias resulted in the underestimation of reference values associated with small molecules (valence correlation energies <0.5 Eh) and an even larger overestimation of atomization energies for bigger systems. Consequently, caution should be exercised in the use of F12b for high accuracy studies. Root mean square and mean absolute deviation error metrics for this basis set combination were comparable to complete basis set values obtained with standard CCSD(T) and the aug-cc-pVDZ through aug-cc-pVQZ basis set sequence. However, the mean signed deviation was an order of magnitude larger. Problems partially due to basis set superposition error were identified with second row compounds which resulted in a weak performance for the smaller VDZ-F12/VTZ-F12 combination of basis sets.
Aerostructural Level Set Topology Optimization for a Common Research Model Wing
NASA Technical Reports Server (NTRS)
Dunning, Peter D.; Stanford, Bret K.; Kim, H. Alicia
2014-01-01
The purpose of this work is to use level set topology optimization to improve the design of a representative wing box structure for the NASA common research model. The objective is to minimize the total compliance of the structure under aerodynamic and body force loading, where the aerodynamic loading is coupled to the structural deformation. A taxi bump case was also considered, where only body force loads were applied. The trim condition that aerodynamic lift must balance the total weight of the aircraft is enforced by allowing the root angle of attack to change. The level set optimization method is implemented on an unstructured three-dimensional grid, so that the method can optimize a wing box with arbitrary geometry. Fast matching and upwind schemes are developed for an unstructured grid, which make the level set method robust and efficient. The adjoint method is used to obtain the coupled shape sensitivities required to perform aerostructural optimization of the wing box structure.
Generation of pareto optimal ensembles of calibrated parameter sets for climate models.
Dalbey, Keith R.; Levy, Michael Nathan
2010-12-01
Climate models have a large number of inputs and outputs. In addition, diverse parameters sets can match observations similarly well. These factors make calibrating the models difficult. But as the Earth enters a new climate regime, parameters sets may cease to match observations. History matching is necessary but not sufficient for good predictions. We seek a 'Pareto optimal' ensemble of calibrated parameter sets for the CCSM climate model, in which no individual criteria can be improved without worsening another. One Multi Objective Genetic Algorithm (MOGA) optimization typically requires thousands of simulations but produces an ensemble of Pareto optimal solutions. Our simulation budget of 500-1000 runs allows us to perform the MOGA optimization once, but with far fewer evaluations than normal. We devised an analytic test problem to aid in the selection MOGA settings. The test problem's Pareto set is the surface of a 6 dimensional hypersphere with radius 1 centered at the origin, or rather the portion of it in the [0,1] octant. We also explore starting MOGA from a space-filling Latin Hypercube sample design, specifically Binning Optimal Symmetric Latin Hypercube Sampling (BOSLHS), instead of Monte Carlo (MC). We compare the Pareto sets based on: their number of points, N, larger is better; their RMS distance, d, to the ensemble's center, 0.5553 is optimal; their average radius, {mu}(r), 1 is optimal; their radius standard deviation, {sigma}(r), 0 is optimal. The estimated distributions for these metrics when starting from MC and BOSLHS are shown in Figs. 1 and 2.
Yao, Y. X.; Wang, C. Z.; Ho, K. M.
2010-06-16
A chemical bonding scheme is presented for the analysis of solid-state systems. The scheme is based on the intrinsic oriented quasiatomic minimal-basis-set orbitals (IO-QUAMBOs) previously developed by Ivanic and Ruedenberg for molecular systems. In the solid-state scheme, IO-QUAMBOs are generated by a unitary transformation of the quasiatomic orbitals located at each site of the system with the criteria of maximizing the sum of the fourth power of interatomic orbital bond order. Possible bonding and antibonding characters are indicated by the single particle matrix elements, and can be further examined by the projected density of states. We demonstrate the method by applications to graphene and (6,0) zigzag carbon nanotube. The oriented-orbital scheme automatically describes the system in terms of sp{sup 2} hybridization. The effect of curvature on the electronic structure of the zigzag carbon nanotube is also manifested in the deformation of the intrinsic oriented orbitals as well as a breaking of symmetry leading to nonzero single particle density matrix elements. In an additional study, the analysis is performed on the Al{sub 3}V compound. The main covalent bonding characters are identified in a straightforward way without resorting to the symmetry analysis. Our method provides a general way for chemical bonding analysis of ab initio electronic structure calculations with any type of basis sets.
Cost optimization for series-parallel execution of a collection of intersecting operation sets
NASA Astrophysics Data System (ADS)
Dolgui, Alexandre; Levin, Genrikh; Rozin, Boris; Kasabutski, Igor
2016-05-01
A collection of intersecting sets of operations is considered. These sets of operations are performed successively. The operations of each set are activated simultaneously. Operation durations can be modified. The cost of each operation decreases with the increase in operation duration. In contrast, the additional expenses for each set of operations are proportional to its time. The problem of selecting the durations of all operations that minimize the total cost under constraint on completion time for the whole collection of operation sets is studied. The mathematical model and method to solve this problem are presented. The proposed method is based on a combination of Lagrangian relaxation and dynamic programming. The results of numerical experiments that illustrate the performance of the proposed method are presented. This approach was used for optimization multi-spindle machines and machining lines, but the problem is common in engineering optimization and thus the techniques developed could be useful for other applications.
Yoo, Sung-Hoon; Oh, Sung-Kwun; Pedrycz, Witold
2015-09-01
In this study, we propose a hybrid method of face recognition by using face region information extracted from the detected face region. In the preprocessing part, we develop a hybrid approach based on the Active Shape Model (ASM) and the Principal Component Analysis (PCA) algorithm. At this step, we use a CCD (Charge Coupled Device) camera to acquire a facial image by using AdaBoost and then Histogram Equalization (HE) is employed to improve the quality of the image. ASM extracts the face contour and image shape to produce a personal profile. Then we use a PCA method to reduce dimensionality of face images. In the recognition part, we consider the improved Radial Basis Function Neural Networks (RBF NNs) to identify a unique pattern associated with each person. The proposed RBF NN architecture consists of three functional modules realizing the condition phase, the conclusion phase, and the inference phase completed with the help of fuzzy rules coming in the standard 'if-then' format. In the formation of the condition part of the fuzzy rules, the input space is partitioned with the use of Fuzzy C-Means (FCM) clustering. In the conclusion part of the fuzzy rules, the connections (weights) of the RBF NNs are represented by four kinds of polynomials such as constant, linear, quadratic, and reduced quadratic. The values of the coefficients are determined by running a gradient descent method. The output of the RBF NNs model is obtained by running a fuzzy inference method. The essential design parameters of the network (including learning rate, momentum coefficient and fuzzification coefficient used by the FCM) are optimized by means of Differential Evolution (DE). The proposed P-RBF NNs (Polynomial based RBF NNs) are applied to facial recognition and its performance is quantified from the viewpoint of the output performance and recognition rate. PMID:26163042
Kim, D.S.; Seong, P.H. . Dept. of Nuclear Engineering)
1994-02-01
This paper describes the optimal testing input sets required for the fault diagnosis of the nuclear power plant digital electronic circuits. With the complicated systems such as very large scale integration (VLSI), nuclear power plant (NPP), and aircraft, testing is the major factor of the maintenance of the system. Particularly, diagnosis time grows quickly with the complexity of the component. In this research, for reduce diagnosis time the authors derived the optimal testing sets that are the minimal testing sets required for detecting the failure and for locating of the failed component. For reduced diagnosis time, the technique presented by Hayes fits best for the approach to testing sets generation among many conventional methods. However, this method has the following disadvantages: (a) it considers only the simple network (b) it concerns only whether the system is in failed state or not and does not provide the way to locate the failed component. Therefore the authors have derived the optimal testing input sets that resolve these problems by Hayes while preserving its advantages. When they applied the optimal testing sets to the automatic fault diagnosis system (AFDS) which incorporates the advanced fault diagnosis method of artificial intelligence technique, they found that the fault diagnosis using the optimal testing sets makes testing the digital electronic circuits much faster than that using exhaustive testing input sets; when they applied them to test the Universal (UV) Card which is a nuclear power plant digital input/output solid state protection system card, they reduced the testing time up to about 100 times.
Seijo, Luis; Barandiarán, Zoila
2004-10-01
We present a linear scaling method for the energy minimization step of semiempirical and first-principles Hartree-Fock and Kohn-Sham calculations. It is based on the self-consistent calculation of the optimum localized orbitals of any localization method of choice and on the use of orbital-specific basis sets. The full set of localized orbitals of a large molecule is seen as an orbital mosaic where each tessera is made of only a few of them. The orbital tesserae are computed out of a set of embedded cluster pseudoeigenvalue coupled equations which are solved in a building-block self-consistent fashion. In each iteration, the embedded cluster equations are solved independently of each other and, as a result, the method is parallel at a high level of the calculation. In addition to full system calculations, the method enables to perform simpler, much less demanding embedded cluster calculations, where only a fraction of the localized molecular orbitals are variational while the rest is frozen, taking advantage of the transferability of the localized orbitals of a given localization method between similar molecules. Monitoring single point energy calculations of large poly(ethylene oxide) molecules and three dimensional carbon monoxide clusters using an extended Huckel Hamiltonian are presented. PMID:15473725
NASA Astrophysics Data System (ADS)
Rach, N.; Müller, M. M.; Calarco, T.; Montangero, S.
2015-12-01
In quantum optimal control theory the success of an optimization algorithm is highly influenced by how the figure of merit to be optimized behaves as a function of the control field, i.e., by the control landscape. Constraints on the control field introduce local minima in the landscape—false traps—which might prevent an efficient solution of the optimal control problem. Rabitz et al. [Science 303, 1998 (2004), 10.1126/science.1093649] showed that local minima occur only rarely for unconstrained optimization. Here, we extend this result to the case of bandwidth-limited control pulses showing that in this case one can eliminate the false traps arising from the constraint. Based on this theoretical understanding, we modify the chopped-random-basis (CRAB) optimal control algorithm and show that this development exploits the advantages of both (unconstrained) gradient algorithms and of truncated basis methods, allowing one to always follow the gradient of the unconstrained landscape by bandwidth-limited control functions. We study the effects of additional constraints and show that for reasonable constraints the convergence properties are still maintained. Finally, we numerically show that this approach saturates the theoretical bound on the minimal bandwidth of the control needed to optimally drive the system.
NASA Astrophysics Data System (ADS)
Guan, Qingze; Blume, Doerte
2016-05-01
The explicit correlated Gaussian (ECG) basis set expansion approach is a variational approach that has been used in various areas, including molecular, nuclear, atomic, and chemical physics. In the world of cold atoms, e.g., the ECG approach has been used to calculate the eigenenergies and eigenstates of few-body systems governed by Efimov physics. Since the first experimental realization of synthesized gauge fields, few-body systems with spin-orbit coupling have attracted a great deal of attention. Here, the ECG approach is customized to few-body systems with both short-range interactions and spin-orbit couplings. Benchmark tests and a performance analysis will be presented. Support by the NSF is gratefully acknowledged.