Isodesmic reaction for accurate theoretical pKa calculations of amino acids and peptides.
Sastre, S; Casasnovas, R; Muñoz, F; Frau, J
2016-04-20
Theoretical and quantitative prediction of pKa values at low computational cost is a current challenge in computational chemistry. We report that the isodesmic reaction scheme provides semi-quantitative predictions (i.e. mean absolute errors of 0.5-1.0 pKa unit) for the pKa1 (α-carboxyl), pKa2 (α-amino) and pKa3 (sidechain groups) of a broad set of amino acids and peptides. This method fills the gaps of thermodynamic cycles for the computational pKa calculation of molecules that are unstable in the gas phase or undergo proton transfer reactions or large conformational changes from solution to the gas phase. We also report the key criteria to choose a reference species to make accurate predictions. This method is computationally inexpensive and makes use of standard density functional theory (DFT) and continuum solvent models. It is also conceptually simple and easy to use for researchers not specialized in theoretical chemistry methods. PMID:27052591
Accurate quantum chemical calculations
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.; Taylor, Peter R.
1989-01-01
An important goal of quantum chemical calculations is to provide an understanding of chemical bonding and molecular electronic structure. A second goal, the prediction of energy differences to chemical accuracy, has been much harder to attain. First, the computational resources required to achieve such accuracy are very large, and second, it is not straightforward to demonstrate that an apparently accurate result, in terms of agreement with experiment, does not result from a cancellation of errors. Recent advances in electronic structure methodology, coupled with the power of vector supercomputers, have made it possible to solve a number of electronic structure problems exactly using the full configuration interaction (FCI) method within a subspace of the complete Hilbert space. These exact results can be used to benchmark approximate techniques that are applicable to a wider range of chemical and physical problems. The methodology of many-electron quantum chemistry is reviewed. Methods are considered in detail for performing FCI calculations. The application of FCI methods to several three-electron problems in molecular physics are discussed. A number of benchmark applications of FCI wave functions are described. Atomic basis sets and the development of improved methods for handling very large basis sets are discussed: these are then applied to a number of chemical and spectroscopic problems; to transition metals; and to problems involving potential energy surfaces. Although the experiences described give considerable grounds for optimism about the general ability to perform accurate calculations, there are several problems that have proved less tractable, at least with current computer resources, and these and possible solutions are discussed.
How Accurately can we Calculate Thermal Systems?
Cullen, D; Blomquist, R N; Dean, C; Heinrichs, D; Kalugin, M A; Lee, M; Lee, Y; MacFarlan, R; Nagaya, Y; Trkov, A
2004-04-20
I would like to determine how accurately a variety of neutron transport code packages (code and cross section libraries) can calculate simple integral parameters, such as K{sub eff}, for systems that are sensitive to thermal neutron scattering. Since we will only consider theoretical systems, we cannot really determine absolute accuracy compared to any real system. Therefore rather than accuracy, it would be more precise to say that I would like to determine the spread in answers that we obtain from a variety of code packages. This spread should serve as an excellent indicator of how accurately we can really model and calculate such systems today. Hopefully, eventually this will lead to improvements in both our codes and the thermal scattering models that they use in the future. In order to accomplish this I propose a number of extremely simple systems that involve thermal neutron scattering that can be easily modeled and calculated by a variety of neutron transport codes. These are theoretical systems designed to emphasize the effects of thermal scattering, since that is what we are interested in studying. I have attempted to keep these systems very simple, and yet at the same time they include most, if not all, of the important thermal scattering effects encountered in a large, water-moderated, uranium fueled thermal system, i.e., our typical thermal reactors.
Accurate theoretical chemistry with coupled pair models.
Neese, Frank; Hansen, Andreas; Wennmohs, Frank; Grimme, Stefan
2009-05-19
Quantum chemistry has found its way into the everyday work of many experimental chemists. Calculations can predict the outcome of chemical reactions, afford insight into reaction mechanisms, and be used to interpret structure and bonding in molecules. Thus, contemporary theory offers tremendous opportunities in experimental chemical research. However, even with present-day computers and algorithms, we cannot solve the many particle Schrodinger equation exactly; inevitably some error is introduced in approximating the solutions of this equation. Thus, the accuracy of quantum chemical calculations is of critical importance. The affordable accuracy depends on molecular size and particularly on the total number of atoms: for orientation, ethanol has 9 atoms, aspirin 21 atoms, morphine 40 atoms, sildenafil 63 atoms, paclitaxel 113 atoms, insulin nearly 800 atoms, and quaternary hemoglobin almost 12,000 atoms. Currently, molecules with up to approximately 10 atoms can be very accurately studied by coupled cluster (CC) theory, approximately 100 atoms with second-order Møller-Plesset perturbation theory (MP2), approximately 1000 atoms with density functional theory (DFT), and beyond that number with semiempirical quantum chemistry and force-field methods. The overwhelming majority of present-day calculations in the 100-atom range use DFT. Although these methods have been very successful in quantum chemistry, they do not offer a well-defined hierarchy of calculations that allows one to systematically converge to the correct answer. Recently a number of rather spectacular failures of DFT methods have been found-even for seemingly simple systems such as hydrocarbons, fueling renewed interest in wave function-based methods that incorporate the relevant physics of electron correlation in a more systematic way. Thus, it would be highly desirable to fill the gap between 10 and 100 atoms with highly correlated ab initio methods. We have found that one of the earliest (and now
Accurate radiative transfer calculations for layered media.
Selden, Adrian C
2016-07-01
Simple yet accurate results for radiative transfer in layered media with discontinuous refractive index are obtained by the method of K-integrals. These are certain weighted integrals applied to the angular intensity distribution at the refracting boundaries. The radiative intensity is expressed as the sum of the asymptotic angular intensity distribution valid in the depth of the scattering medium and a transient term valid near the boundary. Integrated boundary equations are obtained, yielding simple linear equations for the intensity coefficients, enabling the angular emission intensity and the diffuse reflectance (albedo) and transmittance of the scattering layer to be calculated without solving the radiative transfer equation directly. Examples are given of half-space, slab, interface, and double-layer calculations, and extensions to multilayer systems are indicated. The K-integral method is orders of magnitude more accurate than diffusion theory and can be applied to layered scattering media with a wide range of scattering albedos, with potential applications to biomedical and ocean optics. PMID:27409700
Accurate calculation of diffraction-limited encircled and ensquared energy.
Andersen, Torben B
2015-09-01
Mathematical properties of the encircled and ensquared energy functions for the diffraction-limited point-spread function (PSF) are presented. These include power series and a set of linear differential equations that facilitate the accurate calculation of these functions. Asymptotic expressions are derived that provide very accurate estimates for the relative amount of energy in the diffraction PSF that fall outside a square or rectangular large detector. Tables with accurate values of the encircled and ensquared energy functions are also presented. PMID:26368873
Strategy Guideline. Accurate Heating and Cooling Load Calculations
Burdick, Arlan
2011-06-01
This guide presents the key criteria required to create accurate heating and cooling load calculations and offers examples of the implications when inaccurate adjustments are applied to the HVAC design process. The guide shows, through realistic examples, how various defaults and arbitrary safety factors can lead to significant increases in the load estimate. Emphasis is placed on the risks incurred from inaccurate adjustments or ignoring critical inputs of the load calculation.
Strategy Guideline: Accurate Heating and Cooling Load Calculations
Burdick, A.
2011-06-01
This guide presents the key criteria required to create accurate heating and cooling load calculations and offers examples of the implications when inaccurate adjustments are applied to the HVAC design process. The guide shows, through realistic examples, how various defaults and arbitrary safety factors can lead to significant increases in the load estimate. Emphasis is placed on the risks incurred from inaccurate adjustments or ignoring critical inputs of the load calculation.
Theoretical Calculations of Atomic Data for Spectroscopy
NASA Technical Reports Server (NTRS)
Bautista, Manuel A.
2000-01-01
Several different approximations and techniques have been developed for the calculation of atomic structure, ionization, and excitation of atoms and ions. These techniques have been used to compute large amounts of spectroscopic data of various levels of accuracy. This paper presents a review of these theoretical methods to help non-experts in atomic physics to better understand the qualities and limitations of various data sources and assess how reliable are spectral models based on those data.
A general, accurate procedure for calculating molecular interaction force.
Yang, Pinghai; Qian, Xiaoping
2009-09-15
The determination of molecular interaction forces, e.g., van der Waals force, between macroscopic bodies is of fundamental importance for understanding sintering, adhesion and fracture processes. In this paper, we develop an accurate, general procedure for van der Waals force calculation. This approach extends a surface formulation that converts a six-dimensional (6D) volume integral into a 4D surface integral for the force calculation. It uses non-uniform rational B-spline (NURBS) surfaces to represent object surfaces. Surface integrals are then done on the parametric domain of the NURBS surfaces. It has combined advantages of NURBS surface representation and surface formulation, including (1) molecular interactions between arbitrary-shaped objects can be represented and evaluated by the NURBS model further common geometries such as spheres, cones, planes can be represented exactly and interaction forces are thus calculated accurately; (2) calculation efficiency is improved by converting the volume integral to the surface integral. This approach is implemented and validated via its comparison with analytical solutions for simple geometries. Calculation of van der Waals force between complex geometries with surface roughness is also demonstrated. A tutorial on the NURBS approach is given in Appendix A. PMID:19596335
Fast and accurate Coulomb calculation with Gaussian functions.
Füsti-Molnár, László; Kong, Jing
2005-02-15
Coulomb interaction is one of the major time-consuming components in a density functional theory (DFT) calculation. In the last decade, dramatic progresses have been made to improve the efficiency of Coulomb calculation, including continuous fast multipole method (CFMM) and J-engine method, all developed first inside Q-Chem. The most recent development is the advent of Fourier transform Coulomb method developed by Fusti-Molnar and Pulay, and an improved version of the method has been recently implemented in Q-Chem. It replaces the least efficient part of the previous Coulomb methods with an accurate numerical integration scheme that scales in O(N2) instead of O(N4) with the basis size. The result is a much smaller slope in the linear scaling with respect to the molecular size and we will demonstrate through a series of benchmark calculations that it speeds up the calculation of Coulomb energy by several folds over the efficient existing code, i.e., the combination of CFMM and J-engine, without loss of accuracy. Furthermore, we will show that it is complementary to the latter and together the three methods offer the best performance for Coulomb part of DFT calculations, making the DFT calculations affordable for very large systems involving thousands of basis functions. PMID:15743222
Accurate pressure gradient calculations in hydrostatic atmospheric models
NASA Technical Reports Server (NTRS)
Carroll, John J.; Mendez-Nunez, Luis R.; Tanrikulu, Saffet
1987-01-01
A method for the accurate calculation of the horizontal pressure gradient acceleration in hydrostatic atmospheric models is presented which is especially useful in situations where the isothermal surfaces are not parallel to the vertical coordinate surfaces. The present method is shown to be exact if the potential temperature lapse rate is constant between the vertical pressure integration limits. The technique is applied to both the integration of the hydrostatic equation and the computation of the slope correction term in the horizontal pressure gradient. A fixed vertical grid and a dynamic grid defined by the significant levels in the vertical temperature distribution are employed.
Accurate and efficient linear scaling DFT calculations with universal applicability.
Mohr, Stephan; Ratcliff, Laura E; Genovese, Luigi; Caliste, Damien; Boulanger, Paul; Goedecker, Stefan; Deutsch, Thierry
2015-12-21
Density functional theory calculations are computationally extremely expensive for systems containing many atoms due to their intrinsic cubic scaling. This fact has led to the development of so-called linear scaling algorithms during the last few decades. In this way it becomes possible to perform ab initio calculations for several tens of thousands of atoms within reasonable walltimes. However, even though the use of linear scaling algorithms is physically well justified, their implementation often introduces some small errors. Consequently most implementations offering such a linear complexity either yield only a limited accuracy or, if one wants to go beyond this restriction, require a tedious fine tuning of many parameters. In our linear scaling approach within the BigDFT package, we were able to overcome this restriction. Using an ansatz based on localized support functions expressed in an underlying Daubechies wavelet basis - which offers ideal properties for accurate linear scaling calculations - we obtain an amazingly high accuracy and a universal applicability while still keeping the possibility of simulating large system with linear scaling walltimes requiring only a moderate demand of computing resources. We prove the effectiveness of our method on a wide variety of systems with different boundary conditions, for single-point calculations as well as for geometry optimizations and molecular dynamics. PMID:25958954
Accurate calculations of bound rovibrational states for argon trimer
Brandon, Drew; Poirier, Bill
2014-07-21
This work presents a comprehensive quantum dynamics calculation of the bound rovibrational eigenstates of argon trimer (Ar{sub 3}), using the ScalIT suite of parallel codes. The Ar{sub 3} rovibrational energy levels are computed to a very high level of accuracy (10{sup −3} cm{sup −1} or better), and up to the highest rotational and vibrational excitations for which bound states exist. For many of these rovibrational states, wavefunctions are also computed. Rare gas clusters such as Ar{sub 3} are interesting because the interatomic interactions manifest through long-range van der Waals forces, rather than through covalent chemical bonding. As a consequence, they exhibit strong Coriolis coupling between the rotational and vibrational degrees of freedom, as well as highly delocalized states, all of which renders accurate quantum dynamical calculation difficult. Moreover, with its (comparatively) deep potential well and heavy masses, Ar{sub 3} is an especially challenging rare gas trimer case. There are a great many rovibrational eigenstates to compute, and a very high density of states. Consequently, very few previous rovibrational state calculations for Ar{sub 3} may be found in the current literature—and only for the lowest-lying rotational excitations.
Accurate calculations of bound rovibrational states for argon trimer.
Brandon, Drew; Poirier, Bill
2014-07-21
This work presents a comprehensive quantum dynamics calculation of the bound rovibrational eigenstates of argon trimer (Ar3), using the ScalIT suite of parallel codes. The Ar3 rovibrational energy levels are computed to a very high level of accuracy (10(-3) cm(-1) or better), and up to the highest rotational and vibrational excitations for which bound states exist. For many of these rovibrational states, wavefunctions are also computed. Rare gas clusters such as Ar3 are interesting because the interatomic interactions manifest through long-range van der Waals forces, rather than through covalent chemical bonding. As a consequence, they exhibit strong Coriolis coupling between the rotational and vibrational degrees of freedom, as well as highly delocalized states, all of which renders accurate quantum dynamical calculation difficult. Moreover, with its (comparatively) deep potential well and heavy masses, Ar3 is an especially challenging rare gas trimer case. There are a great many rovibrational eigenstates to compute, and a very high density of states. Consequently, very few previous rovibrational state calculations for Ar3 may be found in the current literature-and only for the lowest-lying rotational excitations. PMID:25053315
Efficient determination of accurate atomic polarizabilities for polarizeable embedding calculations.
Schröder, Heiner; Schwabe, Tobias
2016-08-15
We evaluate embedding potentials, obtained via various methods, used for polarizable embedding computations of excitation energies of para-nitroaniline in water and organic solvents as well as of the green fluorescent protein. We found that isotropic polarizabilities derived from DFTD3 dispersion coefficients correlate well with those obtained via the LoProp method. We show that these polarizabilities in conjunction with appropriately derived point charges are in good agreement with calculations employing static multipole moments up to quadrupoles and anisotropic polarizabilities for both computed systems. The (partial) use of these easily-accessible parameters drastically reduces the computational effort to obtain accurate embedding potentials especially for proteins. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc. PMID:27317509
Highly Accurate Calculations of the Phase Diagram of Cold Lithium
NASA Astrophysics Data System (ADS)
Shulenburger, Luke; Baczewski, Andrew
The phase diagram of lithium is particularly complicated, exhibiting many different solid phases under the modest application of pressure. Experimental efforts to identify these phases using diamond anvil cells have been complemented by ab initio theory, primarily using density functional theory (DFT). Due to the multiplicity of crystal structures whose enthalpy is nearly degenerate and the uncertainty introduced by density functional approximations, we apply the highly accurate many-body diffusion Monte Carlo (DMC) method to the study of the solid phases at low temperature. These calculations span many different phases, including several with low symmetry, demonstrating the viability of DMC as a method for calculating phase diagrams for complex solids. Our results can be used as a benchmark to test the accuracy of various density functionals. This can strengthen confidence in DFT based predictions of more complex phenomena such as the anomalous melting behavior predicted for lithium at high pressures. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. DOE's National Nuclear Security Administration under contract DE-AC04-94AL85000.
Accurate theoretical and experimental characterization of optical grating coupler.
Fesharaki, Faezeh; Hossain, Nadir; Vigne, Sebastien; Chaker, Mohamed; Wu, Ke
2016-09-01
Periodic structures, acting as reflectors, filters, and couplers, are a fundamental building block section in many optical devices. In this paper, a three-dimensional simulation of a grating coupler, a well-known periodic structure, is conducted. Guided waves and leakage characteristics of an out-of-plane grating coupler are studied in detail, and its coupling efficiency is examined. Furthermore, a numerical calibration analysis is applied through a commercial software package on the basis of a full-wave finite-element method to calculate the complex propagation constant of the structure and to evaluate the radiation pattern. For experimental evaluation, an optimized grating coupler is fabricated using electron-beam lithography technique and plasma etching. An excellent agreement between simulations and measurements is observed, thereby validating the demonstrated method. PMID:27607706
Accurate ionization potential of semiconductors from efficient density functional calculations
NASA Astrophysics Data System (ADS)
Ye, Lin-Hui
2016-07-01
Despite its huge successes in total-energy-related applications, the Kohn-Sham scheme of density functional theory cannot get reliable single-particle excitation energies for solids. In particular, it has not been able to calculate the ionization potential (IP), one of the most important material parameters, for semiconductors. We illustrate that an approximate exact-exchange optimized effective potential (EXX-OEP), the Becke-Johnson exchange, can be used to largely solve this long-standing problem. For a group of 17 semiconductors, we have obtained the IPs to an accuracy similar to that of the much more sophisticated G W approximation (GWA), with the computational cost of only local-density approximation/generalized gradient approximation. The EXX-OEP, therefore, is likely as useful for solids as for finite systems. For solid surfaces, the asymptotic behavior of the vx c has effects similar to those of finite systems which, when neglected, typically cause the semiconductor IPs to be underestimated. This may partially explain why standard GWA systematically underestimates the IPs and why using the same GWA procedures has not been able to get an accurate IP and band gap at the same time.
Accurate calculation of field and carrier distributions in doped semiconductors
NASA Astrophysics Data System (ADS)
Yang, Wenji; Tang, Jianping; Yu, Hongchun; Wang, Yanguo
2012-06-01
We use the numerical squeezing algorithm(NSA) combined with the shooting method to accurately calculate the built-in fields and carrier distributions in doped silicon films (SFs) in the micron and sub-micron thickness range and results are presented in graphical form for variety of doping profiles under different boundary conditions. As a complementary approach, we also present the methods and the results of the inverse problem (IVP) - finding out the doping profile in the SFs for given field distribution. The solution of the IVP provides us the approach to arbitrarily design field distribution in SFs - which is very important for low dimensional (LD) systems and device designing. Further more, the solution of the IVP is both direct and much easy for all the one-, two-, and three-dimensional semiconductor systems. With current efforts focused on the LD physics, knowing of the field and carrier distribution details in the LD systems will facilitate further researches on other aspects and hence the current work provides a platform for those researches.
Accurate calculation of (31)P NMR chemical shifts in polyoxometalates.
Pascual-Borràs, Magda; López, Xavier; Poblet, Josep M
2015-04-14
We search for the best density functional theory strategy for the determination of (31)P nuclear magnetic resonance (NMR) chemical shifts, δ((31)P), in polyoxometalates. Among the variables governing the quality of the quantum modelling, we tackle herein the influence of the functional and the basis set. The spin-orbit and solvent effects were routinely included. To do so we analysed the family of structures α-[P2W18-xMxO62](n-) with M = Mo(VI), V(V) or Nb(V); [P2W17O62(M'R)](n-) with M' = Sn(IV), Ge(IV) and Ru(II) and [PW12-xMxO40](n-) with M = Pd(IV), Nb(V) and Ti(IV). The main results suggest that, to date, the best procedure for the accurate calculation of δ((31)P) in polyoxometalates is the combination of TZP/PBE//TZ2P/OPBE (for NMR//optimization step). The hybrid functionals (PBE0, B3LYP) tested herein were applied to the NMR step, besides being more CPU-consuming, do not outperform pure GGA functionals. Although previous studies on (183)W NMR suggested that the use of very large basis sets like QZ4P were needed for geometry optimization, the present results indicate that TZ2P suffices if the functional is optimal. Moreover, scaling corrections were applied to the results providing low mean absolute errors below 1 ppm for δ((31)P), which is a step forward in order to confirm or predict chemical shifts in polyoxometalates. Finally, via a simplified molecular model, we establish how the small variations in δ((31)P) arise from energy changes in the occupied and virtual orbitals of the PO4 group. PMID:25738630
Spectroscopically Accurate Calculations of the Rovibrational Energies of Diatomic Hydrogen
NASA Astrophysics Data System (ADS)
Perry, Jason
2005-05-01
The Born-Oppenheimer approximation has been used to calculate the rotational and vibrational states of diatomic hydrogen. Because it is an approximation, our group now wants to use a Born-Oppenheimer potential to calculate the electronic energy that has been corrected to match closely with spectroscopic results. We are using a code that has corrections for adiabatic, relativistic, radiative, and non-adiabatic effects. The rovibrational energies have now been calculated for both bound and quasi-bound states. We also want to compute quadrupole transition probabilities for diatomic hydrogen. These calculations aspire to investigate diatomic hydrogen in astrophysical environments.
Theoretical Calculations of Equations of State
NASA Astrophysics Data System (ADS)
Prakash, Madappa
2016-07-01
The modeling of core-collapse supernovae, neutron stars from their birth to old age, and binary mergers of compact stars requires a detailed knowledge of the equation of state (EOS) of matter at finite temperature. Thermodynamic state variables such as the free energy, energy per baryon, pressure, entropy per baryon, specific heats, chemical potentials of the various species and their derivatives with respect to number densities, thermal and adiabatic indices, etc., all play distinct roles in large-scale computer simulations involving compact objects. In this talk, recent developments in the calculation of the thermal properties of dense matter will be reviewed. Recent results from beyond relativistic mean field theory will be discussed. Highlights will include the role of non-nucleonic degrees of freedom at finite temperature, and possible avenues for future investigations.
Time-accurate Navier-Stokes calculations with multigrid acceleration
NASA Technical Reports Server (NTRS)
Melson, N. D.; Sanetrik, Mark D.; Atkins, Harold L.
1993-01-01
An efficient method for calculating unsteady flows is presented, with emphasis on a modified version of the thin-layer Navier-Stokes equations. Fourier stability analysis is used to illustrate the effect of treating the source term implicitly instead of explicity, as well as to illustrate other algorithmic choices. A 2D circular cylinder (with a Reynolds number of 1200 and a Mach number of 0.3) is calculated. The present scheme requires only about 10 percent of the computer time required by global minimum time stepping.
General formula for accurate calculation of halofullerenes polarizability
NASA Astrophysics Data System (ADS)
Sabirov, D. Sh.; Garipova, R. R.; Bulgakov, R. G.
2012-01-01
At the first time mean polarizabilities of fluoro-, chloro-, bromo[60]fullerenes have been calculated. The dependences of halofullerenes polarizabilities on the structure and the number of halogen atoms in a molecule have been analyzed. The phenomenon of polarizability depression is typical for all compounds under study. General formula for calculation of all classes of halofullerenes mean polarizabilities has been derived based on polarizability depression. Its applicability to related compounds (C70 fullerene choro-derivatives and hypothetical iodo[60]fullerenes) has been shown.
Machine learning of parameters for accurate semiempirical quantum chemical calculations
Dral, Pavlo O.; von Lilienfeld, O. Anatole; Thiel, Walter
2015-04-14
We investigate possible improvements in the accuracy of semiempirical quantum chemistry (SQC) methods through the use of machine learning (ML) models for the parameters. For a given class of compounds, ML techniques require sufficiently large training sets to develop ML models that can be used for adapting SQC parameters to reflect changes in molecular composition and geometry. The ML-SQC approach allows the automatic tuning of SQC parameters for individual molecules, thereby improving the accuracy without deteriorating transferability to molecules with molecular descriptors very different from those in the training set. The performance of this approach is demonstrated for the semiempiricalmore » OM2 method using a set of 6095 constitutional isomers C7H10O2, for which accurate ab initio atomization enthalpies are available. The ML-OM2 results show improved average accuracy and a much reduced error range compared with those of standard OM2 results, with mean absolute errors in atomization enthalpies dropping from 6.3 to 1.7 kcal/mol. They are also found to be superior to the results from specific OM2 reparameterizations (rOM2) for the same set of isomers. The ML-SQC approach thus holds promise for fast and reasonably accurate high-throughput screening of materials and molecules.« less
Machine learning of parameters for accurate semiempirical quantum chemical calculations
Dral, Pavlo O.; von Lilienfeld, O. Anatole; Thiel, Walter
2015-04-14
We investigate possible improvements in the accuracy of semiempirical quantum chemistry (SQC) methods through the use of machine learning (ML) models for the parameters. For a given class of compounds, ML techniques require sufficiently large training sets to develop ML models that can be used for adapting SQC parameters to reflect changes in molecular composition and geometry. The ML-SQC approach allows the automatic tuning of SQC parameters for individual molecules, thereby improving the accuracy without deteriorating transferability to molecules with molecular descriptors very different from those in the training set. The performance of this approach is demonstrated for the semiempirical OM2 method using a set of 6095 constitutional isomers C_{7}H_{10}O_{2}, for which accurate ab initio atomization enthalpies are available. The ML-OM2 results show improved average accuracy and a much reduced error range compared with those of standard OM2 results, with mean absolute errors in atomization enthalpies dropping from 6.3 to 1.7 kcal/mol. They are also found to be superior to the results from specific OM2 reparameterizations (rOM2) for the same set of isomers. The ML-SQC approach thus holds promise for fast and reasonably accurate high-throughput screening of materials and molecules.
The Calculation of Accurate Metal-Ligand Bond Energies
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W.; Partridge, Harry, III; Ricca, Alessandra; Arnold, James O. (Technical Monitor)
1997-01-01
The optimization of the geometry and calculation of zero-point energies are carried out at the B3LYP level of theory. The bond energies are determined at this level, as well as at the CCSD(T) level using very large basis sets. The successive OH bond energies to the first row transition metal cations are reported. For most systems there has been an experimental determination of the first OH. In general, the CCSD(T) values are in good agreement with experiment. The bonding changes from mostly covalent for the early metals to mostly electrostatic for the late transition metal systems.
How Accurately Can We Calculate Neutrons Slowing Down In Water ?
Cullen, D E; Blomquist, R; Greene, M; Lent, E; MacFarlane, R; McKinley, S; Plechaty, E; Sublet, J C
2006-03-30
We have compared the results produced by a variety of currently available Monte Carlo neutron transport codes for the relatively simple problem of a fast source of neutrons slowing down and thermalizing in water. Initial comparisons showed rather large differences in the calculated flux; up to 80% differences. By working together we iterated to improve the results by: (1) insuring that all codes were using the same data, (2) improving the models used by the codes, and (3) correcting errors in the codes; no code is perfect. Even after a number of iterations we still found differences, demonstrating that our Monte Carlo and supporting codes are far from perfect; in particularly we found that the often overlooked nuclear data processing codes can be the weakest link in our systems of codes. The results presented here represent the today's state-of-the-art, in the sense that all of the Monte Carlo codes are modern, widely available and used codes. They all use the most up-to-date nuclear data, and the results are very recent, weeks or at most a few months old; these are the results that current users of these codes should expect to obtain from them. As such, the accuracy and limitations of the codes presented here should serve as guidelines to code users in interpreting their results for similar problems. We avoid crystal ball gazing, in the sense that we limit the scope of this report to what is available to code users today, and we avoid predicting future improvements that may or may not actual come to pass. An exception that we make is in presenting results for an improved thermal scattering model currently being testing using advanced versions of NJOY and MCNP that are not currently available to users, but are planned for release in the not too distant future. The other exception is to show comparisons between experimentally measured water cross sections and preliminary ENDF/B-VII thermal scattering law, S({alpha},{beta}) data; although these data are strictly
Theoretical Calculation and Validation of the Water Vapor Continuum Absorption
NASA Technical Reports Server (NTRS)
Ma, Qiancheng; Tipping, Richard H.
1998-01-01
The primary objective of this investigation is the development of an improved parameterization of the water vapor continuum absorption through the refinement and validation of our existing theoretical formalism. The chief advantage of our approach is the self-consistent, first principles, basis of the formalism which allows us to predict the frequency, temperature and pressure dependence of the continuum absorption as well as provide insights into the physical mechanisms responsible for the continuum absorption. Moreover, our approach is such that the calculated continuum absorption can be easily incorporated into satellite retrieval algorithms and climate models. Accurate determination of the water vapor continuum is essential for the next generation of retrieval algorithms which propose to use the combined constraints of multi-spectral measurements such as those under development for EOS data analysis (e.g., retrieval algorithms based on MODIS and AIRS measurements); current Pathfinder activities which seek to use the combined constraints of infrared and microwave (e.g., HIRS and MSU) measurements to improve temperature and water profile retrievals, and field campaigns which seek to reconcile spectrally-resolved and broad-band measurements such as those obtained as part of FIRE. Current widely used continuum treatments have been shown to produce spectrally dependent errors, with the magnitude of the error dependent on temperature and abundance which produces errors with a seasonal and latitude dependence. Translated into flux, current water vapor continuum parameterizations produce flux errors of order 10 W/ml, which compared to the 4 W/m' magnitude of the greenhouse gas forcing and the 1-2 W/m' estimated aerosol forcing is certainly climatologically significant and unacceptably large. While it is possible to tune the empirical formalisms, the paucity of laboratory measurements, especially at temperatures of interest for atmospheric applications, preclude tuning
Theoretical Calculation and Validation of the Water Vapor Continuum Absorption
NASA Technical Reports Server (NTRS)
Ma, Qiancheng; Tipping, Richard H.
1998-01-01
The primary objective of this investigation is the development of an improved parameterization of the water vapor continuum absorption through the refinement and validation of our existing theoretical formalism. The chief advantage of our approach is the self-consistent, first principles, basis of the formalism which allows us to predict the frequency, temperature and pressure dependence of the continuum absorption as well as provide insights into the physical mechanisms responsible for the continuum absorption. Moreover, our approach is such that the calculated continuum absorption can be easily incorporated into satellite retrieval algorithms and climate models. Accurate determination of the water vapor continuum is essential for the next generation of retrieval algorithms which propose to use the combined constraints of multispectral measurements such as those under development for EOS data analysis (e.g., retrieval algorithms based on MODIS and AIRS measurements); current Pathfinder activities which seek to use the combined constraints of infrared and microwave (e.g., HIRS and MSU) measurements to improve temperature and water profile retrievals, and field campaigns which seek to reconcile spectrally-resolved and broad-band measurements such as those obtained as part of FIRE. Current widely used continuum treatments have been shown to produce spectrally dependent errors, with the magnitude of the error dependent on temperature and abundance which produces errors with a seasonal and latitude dependence. Translated into flux, current water vapor continuum parameterizations produce flux errors of order 10 W/sq m, which compared to the 4 W/sq m magnitude of the greenhouse gas forcing and the 1-2 W/sq m estimated aerosol forcing is certainly climatologically significant and unacceptably large. While it is possible to tune the empirical formalisms, the paucity of laboratory measurements, especially at temperatures of interest for atmospheric applications, preclude
NASA Technical Reports Server (NTRS)
Liu, Yen; Vinokur, Marcel
1989-01-01
This paper treats the accurate and efficient calculation of thermodynamic properties of arbitrary gas mixtures for equilibrium flow computations. New improvements in the Stupochenko-Jaffe model for the calculation of thermodynamic properties of diatomic molecules are presented. A unified formulation of equilibrium calculations for gas mixtures in terms of irreversible entropy is given. Using a highly accurate thermo-chemical data base, a new, efficient and vectorizable search algorithm is used to construct piecewise interpolation procedures with generate accurate thermodynamic variable and their derivatives required by modern computational algorithms. Results are presented for equilibrium air, and compared with those given by the Srinivasan program.
SMARTIES: User-friendly codes for fast and accurate calculations of light scattering by spheroids
NASA Astrophysics Data System (ADS)
Somerville, W. R. C.; Auguié, B.; Le Ru, E. C.
2016-05-01
We provide a detailed user guide for SMARTIES, a suite of MATLAB codes for the calculation of the optical properties of oblate and prolate spheroidal particles, with comparable capabilities and ease-of-use as Mie theory for spheres. SMARTIES is a MATLAB implementation of an improved T-matrix algorithm for the theoretical modelling of electromagnetic scattering by particles of spheroidal shape. The theory behind the improvements in numerical accuracy and convergence is briefly summarized, with reference to the original publications. Instructions of use, and a detailed description of the code structure, its range of applicability, as well as guidelines for further developments by advanced users are discussed in separate sections of this user guide. The code may be useful to researchers seeking a fast, accurate and reliable tool to simulate the near-field and far-field optical properties of elongated particles, but will also appeal to other developers of light-scattering software seeking a reliable benchmark for non-spherical particles with a challenging aspect ratio and/or refractive index contrast.
Discretely disordered photonic bandgap structures: a more accurate invariant measure calculation
NASA Astrophysics Data System (ADS)
Kissel, Glen J.
2009-02-01
In the one-dimensional optical analog to Anderson localization, a periodically layered medium has one or more parameters randomly disordered. Such a randomized system can be modeled by an infinite product of 2x2 random transfer matrices with the upper Lyapunov exponent of the matrix product identified as the localization factor (inverse localization length) for the model. The theorem of Furstenberg allows us, at least theoretically, to calculate this upper Lyapunov exponent. In Furstenberg's formula we not only integrate with respect to the probability measure of the random matrices, but also with respect to the invariant probability measure of the direction of the vector propagated by the random matrices. This invariant measure is difficult to find analytically, and, as a result, the most successful approach is to determine the invariant measure numerically. A Monte Carlo simulation which uses accumulated bin counts to track the direction of the propagated vector through a long chain of random matrices does a good job of estimating the invariant probability measure, but with a level of uncertainty. A potentially more accurate numerical technique by Froyland and Aihara obtains the invariant measure as a left eigenvector of a large sparse matrix containing probability values determined by the action of the random matrices on input vectors. We first apply these two techniques to a random Fibonacci sequence whose Lyapunov exponent was determined by Viswanath. We then demonstrate these techniques on a quarter-wave stack model with binary discrete disorder in layer thickness, and compare results to the continuously disordered counterpart.
An Effective Method to Accurately Calculate the Phase Space Factors for β - β - Decay
Neacsu, Andrei; Horoi, Mihai
2016-01-01
Accurate calculations of the electron phase space factors are necessary for reliable predictions of double-beta decay rates and for the analysis of the associated electron angular and energy distributions. We present an effective method to calculate these phase space factors that takes into account the distorted Coulomb field of the daughter nucleus, yet it allows one to easily calculate the phase space factors with good accuracy relative to the most exact methods available in the recent literature.
Potential theoretic methods for far field sound radiation calculations
NASA Technical Reports Server (NTRS)
Hariharan, S. I.; Stenger, Edward J.; Scott, J. R.
1995-01-01
In the area of computational acoustics, procedures which accurately predict the far-field sound radiation are much sought after. A systematic development of such procedures are found in a sequence of papers by Atassi. The method presented here is an alternate approach to predicting far field sound based on simple layer potential theoretic methods. The main advantages of this method are: it requires only a simple free space Green's function, it can accommodate arbitrary shapes of Kirchoff surfaces, and is readily extendable to three-dimensional problems. Moreover, the procedure presented here, though tested for unsteady lifting airfoil problems, can easily be adapted to other areas of interest, such as jet noise radiation problems. Results are presented for lifting airfoil problems and comparisons are made with the results reported by Atassi. Direct comparisons are also made for the flat plate case.
Tung, Wei-Cheng; Adamowicz, Ludwik
2014-03-28
Very accurate calculations of the ground-state potential energy curve (PEC) of the LiH(+) ion performed with all-electron explicitly correlated Gaussian functions with shifted centers are presented. The variational method is employed. The calculations involve optimization of nonlinear exponential parameters of the Gaussians performed with the aid of the analytical first derivatives of the energy determined with respect to the parameters. The diagonal adiabatic correction is also calculated for each PEC point. The PEC is then used to calculate the vibrational energies of the system. In that calculation, the non-adiabatic effects are accounted for by using an effective vibrational mass obtained by the minimization of the difference between the vibrational energies obtained from the calculations where the Born-Oppenheimer approximation was not assumed and the results of the present calculations. PMID:24697449
NASA Astrophysics Data System (ADS)
Tung, Wei-Cheng; Adamowicz, Ludwik
2014-03-01
Very accurate calculations of the ground-state potential energy curve (PEC) of the LiH+ ion performed with all-electron explicitly correlated Gaussian functions with shifted centers are presented. The variational method is employed. The calculations involve optimization of nonlinear exponential parameters of the Gaussians performed with the aid of the analytical first derivatives of the energy determined with respect to the parameters. The diagonal adiabatic correction is also calculated for each PEC point. The PEC is then used to calculate the vibrational energies of the system. In that calculation, the non-adiabatic effects are accounted for by using an effective vibrational mass obtained by the minimization of the difference between the vibrational energies obtained from the calculations where the Born-Oppenheimer approximation was not assumed and the results of the present calculations.
Accurate calculation of conductive conductances in complex geometries for spacecrafts thermal models
NASA Astrophysics Data System (ADS)
Garmendia, Iñaki; Anglada, Eva; Vallejo, Haritz; Seco, Miguel
2016-02-01
The thermal subsystem of spacecrafts and payloads is always designed with the help of Thermal Mathematical Models. In the case of the Thermal Lumped Parameter (TLP) method, the non-linear system of equations that is created is solved to calculate the temperature distribution and the heat power that goes between nodes. The accuracy of the results depends largely on the appropriate calculation of the conductive and radiative conductances. Several established methods for the determination of conductive conductances exist but they present some limitations for complex geometries. Two new methods are proposed in this paper to calculate accurately these conductive conductances: The Extended Far Field method and the Mid-Section method. Both are based on a finite element calculation but while the Extended Far Field method uses the calculation of node mean temperatures, the Mid-Section method is based on assuming specific temperature values. They are compared with traditionally used methods showing the advantages of these two new methods.
Ab Initio Potential Energy Surfaces and the Calculation of Accurate Vibrational Frequencies
NASA Technical Reports Server (NTRS)
Lee, Timothy J.; Dateo, Christopher E.; Martin, Jan M. L.; Taylor, Peter R.; Langhoff, Stephen R. (Technical Monitor)
1995-01-01
Due to advances in quantum mechanical methods over the last few years, it is now possible to determine ab initio potential energy surfaces in which fundamental vibrational frequencies are accurate to within plus or minus 8 cm(exp -1) on average, and molecular bond distances are accurate to within plus or minus 0.001-0.003 Angstroms, depending on the nature of the bond. That is, the potential energy surfaces have not been scaled or empirically adjusted in any way, showing that theoretical methods have progressed to the point of being useful in analyzing spectra that are not from a tightly controlled laboratory environment, such as vibrational spectra from the interstellar medium. Some recent examples demonstrating this accuracy will be presented and discussed. These include the HNO, CH4, C2H4, and ClCN molecules. The HNO molecule is interesting due to the very large H-N anharmonicity, while ClCN has a very large Fermi resonance. The ab initio studies for the CH4 and C2H4 molecules present the first accurate full quartic force fields of any kind (i.e., whether theoretical or empirical) for a five-atom and six-atom system, respectively.
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.
NASA Technical Reports Server (NTRS)
Avrett, E. H.
1984-01-01
Models and spectra of sunspots were studied, because they are important to energy balance and variability discussions. Sunspot observations in the ultraviolet region 140 to 168 nn was obtained by the NRL High Resolution Telescope and Spectrograph. Extensive photometric observations of sunspot umbrae and prenumbrae in 10 chanels covering the wavelength region 387 to 3800 nm were made. Cool star opacities and model atmospheres were computed. The Sun is the first testcase, both to check the opacity calculations against the observed solar spectrum, and to check the purely theoretical model calculation against the observed solar energy distribution. Line lists were finally completed for all the molecules that are important in computing statistical opacities for energy balance and for radiative rate calculations in the Sun (except perhaps for sunspots). Because many of these bands are incompletely analyzed in the laboratory, the energy levels are not well enough known to predict wavelengths accurately for spectrum synthesis and for detailed comparison with the observations.
NASA Astrophysics Data System (ADS)
Skone, Jonathan; Govoni, Marco; Galli, Giulia
Dielectric-dependent hybrid [DDH] functionals have recently been shown to yield highly accurate energy gaps and dielectric constants for a wide variety of solids, at a computational cost considerably less than standard GW calculations. The fraction of exact exchange included in the definition of DDH functionals depends (self-consistently) on the dielectric constant of the material. In the present talk we introduce a range-separated (RS) version of DDH functionals where short and long-range components are matched using material dependent, non-empirical parameters. Comparing with state of the art GW calculations and experiment, we show that such RS hybrids yield accurate electronic properties of both molecules and solids, including energy gaps, photoelectron spectra and absolute ionization potentials. This work was supported by NSF-CCI Grant Number NSF-CHE-0802907 and DOE-BES.
Accurate near-field calculation in the rigorous coupled-wave analysis method
NASA Astrophysics Data System (ADS)
Weismann, Martin; Gallagher, Dominic F. G.; Panoiu, Nicolae C.
2015-12-01
The rigorous coupled-wave analysis (RCWA) is one of the most successful and widely used methods for modeling periodic optical structures. It yields fast convergence of the electromagnetic far-field and has been adapted to model various optical devices and wave configurations. In this article, we investigate the accuracy with which the electromagnetic near-field can be calculated by using RCWA and explain the observed slow convergence and numerical artifacts from which it suffers, namely unphysical oscillations at material boundaries due to the Gibbs phenomenon. In order to alleviate these shortcomings, we also introduce a mathematical formulation for accurate near-field calculation in RCWA, for one- and two-dimensional straight and slanted diffraction gratings. This accurate near-field computational approach is tested and evaluated for several representative test-structures and configurations in order to illustrate the advantages provided by the proposed modified formulation of the RCWA.
Accurate calculation of phase shifts for electron collisions with positive ions
NASA Astrophysics Data System (ADS)
Gien, T. T.
2003-06-01
The Harris-Nesbet variational method was considered for the calculation of phase shifts of electron collisions with hydrogen-like ions (Li2+, Be3+, and B4+). Calculations were carried out for both singlet and triplet scattering. Very accurate results of phase shift of electron collisions with these ionic targets were obtained for the first time for partial waves of L up to six. The phase shifts that we obtained for low partial wave (S, P, and D) scattering were compared with those available in the literature by a few other research groups employing different numerical methods.
A localized basis that allows fast and accurate second order Moller-Plesset calculations
Subotnik, Joseph E.; Head-Gordon, Martin
2004-10-27
We present a method for computing a basis of localized orthonormal orbitals (both occupied and virtual), in whose representation the Fock matrix is extremely diagonal-dominant. The existence of these orbitals is shown empirically to be sufficient for achieving highly accurate MP@ energies, calculated according to Kapuy's method. This method (which we abbreviate KMP2), which involves a different partitioning of the n-electron Hamiltonian, scales at most quadratically with potential for linearity in the number of electrons. As such, we believe the KMP2 algorithm presented here could be the basis of a viable approach to local correlation calculations.
Theoretical model for calculation of helicity in solar active regions
NASA Astrophysics Data System (ADS)
Chatterjee, P.
We (Choudhuri, Chatterjee and Nandy, 2005) calculate helicities of solar active regions based on the idea of Choudhuri (2003) that poloidal flux lines get wrapped around a toroidal flux tube rising through the convection zone, thereby giving rise to the helicity. Rough estimates based on this idea compare favourably with the observed magnitude of helicity. We use our solar dynamo model based on the Babcock--Leighton α-effect to study how helicity varies with latitude and time. At the time of solar maximum, our theoretical model gives negative helicity in the northern hemisphere and positive helicity in the south, in accordance with observed hemispheric trends. However, we find that, during a short interval at the beginning of a cycle, helicities tend to be opposite of the preferred hemispheric trends. Next we (Chatterjee, Choudhuri and Petrovay 2006) use the above idea along with the sunspot decay model of Petrovay and Moreno-Insertis, (1997) to estimate the distribution of helicity inside a flux tube as it keeps collecting more azimuthal flux during its rise through the convection zone and as turbulent diffusion keeps acting on it. By varying parameters over reasonable ranges in our simple 1-d model, we find that the azimuthal flux penetrates the flux tube to some extent instead of being confined to a narrow sheath outside.
Accurate variational calculations and analysis of the HOCl vibrational energy spectrum
Skokov, S.; Qi, J.; Bowman, J.M.; Yang, C.; Gray, S.K.; Peterson, K.A. |; Mandelshtam, V.A.
1998-12-01
Large scale variational calculations for the vibrational states of HOCl are performed using a recently developed, accurate {ital ab initio} potential energy surface. Three different approaches for obtaining vibrational states are employed and contrasted; a truncation/recoupling scheme with direct diagonalization, the Lanczos method, and Chebyshev iteration with filter diagonalization. The complete spectrum of bound states for nonrotating HOCl is computed and analyzed within a random matrix theory framework. This analysis indicates almost entirely regular dynamics with only a small degree of chaos. The nearly regular spectral structure allows us to make assignments for the most significant part of the spectrum, based on analysis of coordinate expectation values and eigenfunctions. Ground state dipole moments and dipole transition probabilities are also calculated using accurate {ital ab initio} data. Computed values are in good agreement with available experimental data. Some exact rovibrational calculations for J=1, including Coriolis coupling, are performed. The exact results are nearly identical with those obtained from the adiabatic rotation approximation and very close to those from the centrifugal sudden approximation, thus indicating a very small degree of asymmetry and Coriolis coupling for the HOCl molecule. {copyright} {ital 1998 American Institute of Physics.}
Accurate calculation of computer-generated holograms using angular-spectrum layer-oriented method.
Zhao, Yan; Cao, Liangcai; Zhang, Hao; Kong, Dezhao; Jin, Guofan
2015-10-01
Fast calculation and correct depth cue are crucial issues in the calculation of computer-generated hologram (CGH) for high quality three-dimensional (3-D) display. An angular-spectrum based algorithm for layer-oriented CGH is proposed. Angular spectra from each layer are synthesized as a layer-corresponded sub-hologram based on the fast Fourier transform without paraxial approximation. The proposed method can avoid the huge computational cost of the point-oriented method and yield accurate predictions of the whole diffracted field compared with other layer-oriented methods. CGHs of versatile formats of 3-D digital scenes, including computed tomography and 3-D digital models, are demonstrated with precise depth performance and advanced image quality. PMID:26480062
Rapid and accurate calculation of protein 1H, 13C and 15N chemical shifts.
Neal, Stephen; Nip, Alex M; Zhang, Haiyan; Wishart, David S
2003-07-01
A computer program (SHIFTX) is described which rapidly and accurately calculates the diamagnetic 1H, 13C and 15N chemical shifts of both backbone and sidechain atoms in proteins. The program uses a hybrid predictive approach that employs pre-calculated, empirically derived chemical shift hypersurfaces in combination with classical or semi-classical equations (for ring current, electric field, hydrogen bond and solvent effects) to calculate 1H, 13C and 15N chemical shifts from atomic coordinates. The chemical shift hypersurfaces capture dihedral angle, sidechain orientation, secondary structure and nearest neighbor effects that cannot easily be translated to analytical formulae or predicted via classical means. The chemical shift hypersurfaces were generated using a database of IUPAC-referenced protein chemical shifts--RefDB (Zhang et al., 2003), and a corresponding set of high resolution (<2.1 A) X-ray structures. Data mining techniques were used to extract the largest pairwise contributors (from a list of approximately 20 derived geometric, sequential and structural parameters) to generate the necessary hypersurfaces. SHIFTX is rapid (<1 CPU second for a complete shift calculation of 100 residues) and accurate. Overall, the program was able to attain a correlation coefficient (r) between observed and calculated shifts of 0.911 (1Halpha), 0.980 (13Calpha), 0.996 (13Cbeta), 0.863 (13CO), 0.909 (15N), 0.741 (1HN), and 0.907 (sidechain 1H) with RMS errors of 0.23, 0.98, 1.10, 1.16, 2.43, 0.49, and 0.30 ppm, respectively on test data sets. We further show that the agreement between observed and SHIFTX calculated chemical shifts can be an extremely sensitive measure of the quality of protein structures. Our results suggest that if NMR-derived structures could be refined using heteronuclear chemical shifts calculated by SHIFTX, their precision could approach that of the highest resolution X-ray structures. SHIFTX is freely available as a web server at http
On the accurate theoretical determination of the static hyperpolarizability of trans-butadiene
NASA Astrophysics Data System (ADS)
Maroulis, George
1999-07-01
Finite-field many-body perturbation theory and coupled cluster calculations are reported for the static second dipole hyperpolarizability γαβγδ of trans-butadiene. A very large basis set of [9s6p4d1f/6s3p1d] size (336 contracted Gaussian-type functions) should lead to self-consistent field (SCF) values of near-Hartree-Fock quality. We report γxxxx=6.19, γxxxz=-0.44, γxxyy=3.42, γzzxx=2.07, γxyyz=-0.50, γxzzz=1.73, γyyyy=14.72, γyyzz=8.46, γzzzz=24.10 and γ¯=14.58 for 10-3×γαβγδ/e4a04Eh-3 at the experimental geometry (molecule on the xz plane with z as the main axis). γ¯=(14.6±0.4)×103e4a04Eh-3 should be a very reliable estimate of the Hartree-Fock limit of the mean hyperpolarizability. Keeping all other molecular geometry parameters constant, we find that near the Hartree-Fock limit the mean hyperpolarizability varies with the C=C bond length as 10-3×γ¯(RC=C)/e4a04Eh-3=14.93+31.78ΔR+30.88ΔR2-2.96ΔR3 and with the C-C bond length as 10-3×γ¯(RC-C)/e4a04Eh-3=14.93-7.20ΔR+3.04ΔR2, where ΔR/a0 is the displacement from the respective experimental value. The dependence of the components of γαβγδ on the molecular geometry parameters is not uniform. Electron correlation corrections have been calculated at various molecular geometries at the coupled-cluster single, double and perturbatively linked triple excitations level of theory for all independent components of γαβγδ. In absolute terms, electron correlation affects strongly the γzzzz, less strongly the γxxxx, and even less strongly the out-of-plane component γyyyy. The present analysis suggests a conservative estimate of (3.0±0.6)×103e4a04Eh-3 for the electron correlation correction to γ¯ at the experimental molecular geometry. Most of this value is appropriate to γzzzz. A static limit of γ¯=(17.6±1.0)×103e4a04Eh-3 is advanced (neglecting vibrational averaging). Even if a crude theoretical estimate of the dispersion of γ¯ at 1064 nm is added to this value, the
Abate-Pella, Daniel; Freund, Dana M; Ma, Yan; Simón-Manso, Yamil; Hollender, Juliane; Broeckling, Corey D; Huhman, David V; Krokhin, Oleg V; Stoll, Dwight R; Hegeman, Adrian D; Kind, Tobias; Fiehn, Oliver; Schymanski, Emma L; Prenni, Jessica E; Sumner, Lloyd W; Boswell, Paul G
2015-09-18
Identification of small molecules by liquid chromatography-mass spectrometry (LC-MS) can be greatly improved if the chromatographic retention information is used along with mass spectral information to narrow down the lists of candidates. Linear retention indexing remains the standard for sharing retention data across labs, but it is unreliable because it cannot properly account for differences in the experimental conditions used by various labs, even when the differences are relatively small and unintentional. On the other hand, an approach called "retention projection" properly accounts for many intentional differences in experimental conditions, and when combined with a "back-calculation" methodology described recently, it also accounts for unintentional differences. In this study, the accuracy of this methodology is compared with linear retention indexing across eight different labs. When each lab ran a test mixture under a range of multi-segment gradients and flow rates they selected independently, retention projections averaged 22-fold more accurate for uncharged compounds because they properly accounted for these intentional differences, which were more pronounced in steep gradients. When each lab ran the test mixture under nominally the same conditions, which is the ideal situation to reproduce linear retention indices, retention projections still averaged 2-fold more accurate because they properly accounted for many unintentional differences between the LC systems. To the best of our knowledge, this is the most successful study to date aiming to calculate (or even just to reproduce) LC gradient retention across labs, and it is the only study in which retention was reliably calculated under various multi-segment gradients and flow rates chosen independently by labs. PMID:26292625
Efficient yet accurate approximations for ab initio calculations of alcohol cluster thermochemistry
NASA Astrophysics Data System (ADS)
Umer, Muhammad; Kopp, Wassja A.; Leonhard, Kai
2015-12-01
We have calculated the binding enthalpies and entropies of gas phase alcohol clusters from ethanol to 1-decanol. In addition to the monomers, we have investigated dimers, tetramers, and pentamers. Geometries have been obtained at the B3LYP/TZVP level and single point energy calculations have been performed with the Resolution of the Identity-MP2 (RIMP2) method and basis set limit extrapolation using aug-cc-pVTZ and aug-cc-pVQZ basis sets. Thermochemistry is calculated with decoupled hindered rotor treatment for large amplitude motions. The results show three points: First, it is more accurate to transfer the rigid-rotor harmonic oscillator entropies from propanol to longer alcohols than to compute them with an ultra-fine grid and tight geometry convergence criteria. Second, the computational effort can be reduced considerably by using dimerization energies of longer alcohols at density functional theory (B3LYP) level plus a RIMP2 correction obtained from 1-propanol. This approximation yields results almost with the same accuracy as RIMP2 — both methods differ for 1-decanol only 0.4 kJ/mol. Third, the entropy of dimerization including the hindered rotation contribution is converged at 1-propanol with respect to chain length. This allows for a transfer of hindered rotation contributions from smaller alcohols to longer ones which reduces the required computational and man power considerably.
Efficient yet accurate approximations for ab initio calculations of alcohol cluster thermochemistry.
Umer, Muhammad; Kopp, Wassja A; Leonhard, Kai
2015-12-01
We have calculated the binding enthalpies and entropies of gas phase alcohol clusters from ethanol to 1-decanol. In addition to the monomers, we have investigated dimers, tetramers, and pentamers. Geometries have been obtained at the B3LYP/TZVP level and single point energy calculations have been performed with the Resolution of the Identity-MP2 (RIMP2) method and basis set limit extrapolation using aug-cc-pVTZ and aug-cc-pVQZ basis sets. Thermochemistry is calculated with decoupled hindered rotor treatment for large amplitude motions. The results show three points: First, it is more accurate to transfer the rigid-rotor harmonic oscillator entropies from propanol to longer alcohols than to compute them with an ultra-fine grid and tight geometry convergence criteria. Second, the computational effort can be reduced considerably by using dimerization energies of longer alcohols at density functional theory (B3LYP) level plus a RIMP2 correction obtained from 1-propanol. This approximation yields results almost with the same accuracy as RIMP2 - both methods differ for 1-decanol only 0.4 kJ/mol. Third, the entropy of dimerization including the hindered rotation contribution is converged at 1-propanol with respect to chain length. This allows for a transfer of hindered rotation contributions from smaller alcohols to longer ones which reduces the required computational and man power considerably. PMID:26646881
Takahashi, F; Endo, A
2007-01-01
A system utilising radiation transport codes has been developed to derive accurate dose distributions in a human body for radiological accidents. A suitable model is quite essential for a numerical analysis. Therefore, two tools were developed to setup a 'problem-dependent' input file, defining a radiation source and an exposed person to simulate the radiation transport in an accident with the Monte Carlo calculation codes-MCNP and MCNPX. Necessary resources are defined by a dialogue method with a generally used personal computer for both the tools. The tools prepare human body and source models described in the input file format of the employed Monte Carlo codes. The tools were validated for dose assessment in comparison with a past criticality accident and a hypothesized exposure. PMID:17510203
Temperature dependent effective potential method for accurate free energy calculations of solids
NASA Astrophysics Data System (ADS)
Hellman, Olle; Steneteg, Peter; Abrikosov, I. A.; Simak, S. I.
2013-03-01
We have developed a thorough and accurate method of determining anharmonic free energies, the temperature dependent effective potential technique (TDEP). It is based on ab initio molecular dynamics followed by a mapping onto a model Hamiltonian that describes the lattice dynamics. The formalism and the numerical aspects of the technique are described in detail. A number of practical examples are given, and results are presented, which confirm the usefulness of TDEP within ab initio and classical molecular dynamics frameworks. In particular, we examine from first principles the behavior of force constants upon the dynamical stabilization of the body centered phase of Zr, and show that they become more localized. We also calculate the phase diagram for 4He modeled with the Aziz potential and obtain results which are in favorable agreement both with respect to experiment and established techniques.
Accurate calculation of Coulomb sums: Efficacy of Pade-like methods
Sarkar, B. ); Bhattacharyya, K. )
1993-09-01
The adequacy of numerical sequence accelerative transforms in providing accurate estimates of Coulomb sums is considered, referring particularly to distorted lattices. Performance of diagonal Pade approximants (DPA) in this context is critically assessed. Failure in the case of lattice vacancies is also demonstrated. The method of multiple-point Pade approximants (MPA) has been introduced for slowly convergent sequences and is shown to work well for both regular and distorted lattices, the latter being due either to impurities or vacancies. Viability of the two methods is also compared. In divergent situations with distortions owing to vacancies, a strategy of obtaining reliable results by separate applications of both DPA and MPA at appropriate places is also sketched. Representative calculations involve two basic cubic-lattice sums, one slowly convergent and the other divergent, from which very good quality estimates of Madelung constants for a number of common lattices follow.
Accurate band gaps of semiconductors and insulators from Quantum Monte Carlo calculations
NASA Astrophysics Data System (ADS)
Nazarov, Roman; Hood, Randolph; Morales, Miguel
2015-03-01
Ab initio calculations are useful tools in developing materials with targeted band gaps for semiconductor industry. Unfortunately, the main workhorse of ab initio calculations - density functional theory (DFT) in local density approximation (LDA) or generalized gradient approximation (GGA) underestimates band gaps. Several approaches have been proposed starting from empirical corrections to more elaborate exchange-correlation functionals to deal with this problem. But none of these work well for the entire range of semiconductors and insulators. Deficiencies of DFT as a mean field method can be overcome using many-body techniques. Quantum Monte Carlo (QMC) methods can obtain a nearly exact numerical solutions of both total energies and spectral properties. Diffusion Monte Carlo (DMC), the most widely used QMC method, has been shown to provide gold standard results for different material properties, including spectroscopic constants of dimers and clusters, equation of state for solids, accurate descriptions of defects in metals and insulators. To test DMC's accuracy in a wider range of semiconductors and insulators we have computed band gaps of several semiconductors and insulators. We show that DMC can provide superior agreement with experiment compared with more traditional DFT approaches including high level exchange-correlation functionals (e.g. HSE).
Accurate calculation of the dissociation energy of the highly anharmonic system ClHCl(-).
Stein, Christopher; Oswald, Rainer; Botschwina, Peter; Peterson, Kirk A
2015-05-28
Accurate bond dissociation energies (D0) are reported for different isotopologues of the highly anharmonic system ClHCl(-). The mass-independent equilibrium dissociation energy De was obtained by a composite method with frozen-core (fc) CCSD(T) as the basic contribution. Basis sets as large as aug-cc-pV8(+d)Z were employed, and extrapolation to the complete basis set (CBS) limit was carried out. Explicitly correlated calculations with the CCSD(T)-F12b method were also performed to support the conventionally calculated values. Core-core and core-valence correlation, scalar relativity, and higher-order correlation were considered as well. Two mass-dependent contributions, namely, the diagonal Born-Oppenheimer correction and the difference in zero-point energies between the complex and the HCl fragment, were then added in order to arrive at precise D0 values. Results for (35)ClH(35)Cl(-) and (35)ClD(35)Cl(-) are 23.81 and 23.63 kcal/mol, respectively, with estimated uncertainties of 0.05 kcal/mol. In contrast to FHF(-) ( Stein , C. ; Oswald , R. ; Sebald , P. ; Botschwina , P. ; Stoll , H. , Peterson , K. A. Mol. Phys. 2013 , 111 , 2647 - 2652 ), the D0 values of the bichloride species are larger than their De counterparts, which is an unusual situation in hydrogen-bonded systems. PMID:25405989
Accurate 2d finite element calculations for hydrogen in magnetic fields of arbitrary strength
NASA Astrophysics Data System (ADS)
Schimeczek, C.; Wunner, G.
2014-02-01
Recent observations of hundreds of hydrogen-rich magnetic white dwarf stars with magnetic fields up to 105 T (103 MG) have called for more comprehensive and accurate databases for wavelengths and oscillator strengths of the H atom in strong magnetic fields for all states evolving from the field-free levels with principal quantum numbers n≤10. We present a code to calculate the energy eigenvalues and wave functions of such states which is capable of covering the entire regime of field strengths B=0 T to B˜109 T. We achieve this high flexibility by using a two-dimensional finite element expansion of the wave functions in terms of B-splines in the directions parallel and perpendicular to the magnetic field, instead of using asymptotically valid basis expansions in terms of spherical harmonics or Landau orbitals. We have paid special attention to the automation of the program such that the data points for the magnetic field strengths at which the energy of a given state are calculated can be selected automatically. Furthermore, an elaborate method for varying the basis parameters is applied to ensure that the results reach a pre-selected precision, which also can be adjusted freely. Energies and wave functions are stored in a convenient format for further analysis, e.g. for the calculation of transition energies and oscillator strengths. The code has been tested to work for 300 states with an accuracy of better than 10-6 Rydberg across several symmetry subspaces over the entire regime of magnetic field strengths.
NASA Astrophysics Data System (ADS)
Kopparla, P.; Natraj, V.; Spurr, R. J. D.; Shia, R. L.; Yung, Y. L.
2014-12-01
Radiative transfer (RT) computations are an essential component of energy budget calculations in climate models. However, full treatment of RT processes is computationally expensive, prompting usage of 2-stream approximations in operational climate models. This simplification introduces errors of the order of 10% in the top of the atmosphere (TOA) fluxes [Randles et al., 2013]. Natraj et al. [2005, 2010] and Spurr and Natraj [2013] demonstrated the ability of a technique using principal component analysis (PCA) to speed up RT simulations. In the PCA method for RT performance enhancement, empirical orthogonal functions are developed for binned sets of inherent optical properties that possess some redundancy; costly multiple-scattering RT calculations are only done for those (few) optical states corresponding to the most important principal components, and correction factors are applied to approximate radiation fields. Here, we extend the PCA method to a broadband spectral region from the ultraviolet to the shortwave infrared (0.3-3 micron), accounting for major gas absorptions in this region. Comparisons between the new model, called Universal Principal Component Analysis model for Radiative Transfer (UPCART), 2-stream models (such as those used in climate applications) and line-by-line RT models are performed, in order for spectral radiances, spectral fluxes and broadband fluxes. Each of these are calculated at the TOA for several scenarios with varying aerosol types, extinction and scattering optical depth profiles, and solar and viewing geometries. We demonstrate that very accurate radiative forcing estimates can be obtained, with better than 1% accuracy in all spectral regions and better than 0.1% in most cases as compared to an exact line-by-line RT model. The model is comparable in speeds to 2-stream models, potentially rendering UPCART useful for operational General Circulation Models (GCMs). The operational speed and accuracy of UPCART can be further
Morales, Giovanni; Martínez, Ramiro
2009-07-30
This research's main goals were to analyze ketene dimers' relative stability and expand group additivity value (GAV) methodology for estimating the thermochemical properties of high-weight ketene polymers (up to tetramers). The CBS-Q multilevel procedure and statistical thermodynamics were used for calculating the thermochemical properties of 20 cyclic structures, such as diketenes, cyclobutane-1,3-diones, cyclobut-2-enones and pyran-4-ones, as well as 57 acyclic base compounds organized into five groups. According to theoretical heat of formation predictions, diketene was found to be thermodynamically favored over cyclobutane-1,3-dione and its enol-tautomeric form (3-hydroxycyclobut-2-enone). This result did not agree with old combustion experiments. 3-Hydroxycyclobut-2-enone was found to be the least stable dimer and its reported experimental detection in solution may have been due to solvent effects. Substituted diketenes had lower stability than substituted cyclobutane-1,3-diones with an increased number of methyl substituents, suggesting that cyclobutane-1,3-dione type dimers are the major products because of thermodynamic control of alkylketene dimerization. Missing GAVs for the ketene dimers and related structures were calculated through linear regression on the 57 acyclic base compounds. Corrections for non next neighbor interactions (such as gauche, eclipses, and internal hydrogen bond) were needed for obtaining a highly accurate and precise regression model. To the best of our knowledge, the hydrogen bond correction for GAV methodology is the first reported in the literature; this correction was correlated to MP2/6-31Gdagger and HF/6-31Gdagger derived geometries to facilitate its application. GAVs assessed by the linear regression model were able to reproduce acyclic compounds' theoretical thermochemical properties and experimental heat of formation for acetylacetone. Ring formation and substituent position corrections were calculated by consecutively
The focus of light - theoretical calculation and experimental tomographic reconstruction
NASA Astrophysics Data System (ADS)
Quabis, S.; Dorn, R.; Eberler, M.; Glöckl, O.; Leuchs, G.
2001-01-01
We present numerical calculations on the field distribution in the focus of an optical system with high numerical aperture. The diffraction integrals which are based on the Debye approximation are derived and evaluated for a radially polarized input field with a doughnut-shaped intensity distribution. It is shown that this mode focusses down to a spot size significantly smaller as compared to the case of linear polarization. An experimental setup to measure the three-dimensional intensity distribution in the focal region is presented, which is based on the knife-edge method and on tomographic reconstruction.
Calorimetric determinations and theoretical calculations of polymorphs of thalidomide
NASA Astrophysics Data System (ADS)
Lara-Ochoa, F.; Pérez, G. Espinosa; Mijangos-Santiago, F.
2007-09-01
The analysis of the thermograms of thalidomide obtained for the two reported polymorphs α and β by differential scanning calorimetry (DSC) shows some inconsistencies that are discussed in the present work. The conception of a new polymorph form, named β ∗, allowed us to explain the observed thermal behavior more satisfactorily. This new polymorph shows enantiotropy with both α and β polymorphs, reflected in the unique endotherm obtained in the DSC-thermograms, when a heating rate of 10 °C/min is applied. Several additional experiments, such as re-melting of both polymorph forms, showed that there is indeed a new polymorph with an endotherm located between the endotherms of α and β. IR, Raman, and powder X-ray permit us to characterize the isolated compound, resulting from the re-melting of both polymorph forms. Mechanical calculations were performed to elucidate the conformations of each polymorph, and ab initio quantum chemical calculations were performed to determine the energy of the more stable conformers and the spatial cell energy for both polymorphs α and β. These results suggested a possible conformation for the newly discovered polymorph β ∗.
Theoretical calculation of spectra of dibutyl phthalate and dioctyl phthalate
NASA Astrophysics Data System (ADS)
Du, Jian-Bin; Tang, Yan-Lin; Long, Zheng-Wen; Hu, Shuang-Hui; Li, Tao
2014-05-01
Dibutyl phthalate DBP and dioctyl phthalate DOP are the main components of the plasticizers. In order to investigate their molecular structure, chemical bond and spectrum, the geometrical parameters of the ground state and infrared (IR) spectrum are calculated using the density functional theory B3LYP method at the level of 6-311++G( d, p). On this basis, the first twenty-six excited states and the UV-Vis absorption spectra of DBP and DOP are studied using the time-dependent density functional theory (TDDFT) in the same fundamental group and compared with the ultraviolet absorption peak of the molecules measured with UNICO UV-Vis spectrophotometer. The two kinds of molecular spectra are then classified and compared with that in reference. The results show that the strong absorption of IR spectra of DOP and DBP are produced by C-H bending in-plane vibration and C=O telescopic vibration producing. The most absorption of UV-Vis absorption spectra appears in the end absorption belt from n to σ* transition, and the stronger absorption in the E belt of benzene electronic transition from π to π*. There are blue shift for DOP end absorption belt from n to σ* transition and red shift for DOP E absorption belt from π to π* transition relative to that of DBP. This calculation results are better in accord with the spectral data measured by UNICO ultraviolet and visible spectrophotometer.
Plasmon spectroscopy: Theoretical and numerical calculations, and optimization techniques
NASA Astrophysics Data System (ADS)
Rodríguez-Oliveros, Rogelio; Paniagua-Domínguez, Ramón; Sánchez-Gil, José A.; Macías, Demetrio
2016-02-01
We present an overview of recent advances in plasmonics, mainly concerning theoretical and numerical tools required for the rigorous determination of the spectral properties of complex-shape nanoparticles exhibiting strong localized surface plasmon resonances (LSPRs). Both quasistatic approaches and full electrodynamic methods are described, providing a thorough comparison of their numerical implementations. Special attention is paid to surface integral equation formulations, giving examples of their performance in complicated nanoparticle shapes of interest for their LSPR spectra. In this regard, complex (single) nanoparticle configurations (nanocrosses and nanorods) yield a hierarchy of multiple-order LSPR s with evidence of a rich symmetric or asymmetric (Fano-like) LSPR line shapes. In addition, means to address the design of complex geometries to retrieve LSPR spectra are commented on, with special interest in biologically inspired algorithms. Thewealth of LSPRbased applications are discussed in two choice examples, single-nanoparticle surface-enhanced Raman scattering (SERS) and optical heating, and multifrequency nanoantennas for fluorescence and nonlinear optics.
NASA Astrophysics Data System (ADS)
Kopparla, P.; Natraj, V.; Shia, R. L.; Spurr, R. J. D.; Crisp, D.; Yung, Y. L.
2015-12-01
Radiative transfer (RT) computations form the engine of atmospheric retrieval codes. However, full treatment of RT processes is computationally expensive, prompting usage of two-stream approximations in current exoplanetary atmospheric retrieval codes [Line et al., 2013]. Natraj et al. [2005, 2010] and Spurr and Natraj [2013] demonstrated the ability of a technique using principal component analysis (PCA) to speed up RT computations. In the PCA method for RT performance enhancement, empirical orthogonal functions are developed for binned sets of inherent optical properties that possess some redundancy; costly multiple-scattering RT calculations are only done for those few optical states corresponding to the most important principal components, and correction factors are applied to approximate radiation fields. Kopparla et al. [2015, in preparation] extended the PCA method to a broadband spectral region from the ultraviolet to the shortwave infrared (0.3-3 micron), accounting for major gas absorptions in this region. Here, we apply the PCA method to a some typical (exo-)planetary retrieval problems. Comparisons between the new model, called Universal Principal Component Analysis Radiative Transfer (UPCART) model, two-stream models and line-by-line RT models are performed, for spectral radiances, spectral fluxes and broadband fluxes. Each of these are calculated at the top of the atmosphere for several scenarios with varying aerosol types, extinction and scattering optical depth profiles, and stellar and viewing geometries. We demonstrate that very accurate radiance and flux estimates can be obtained, with better than 1% accuracy in all spectral regions and better than 0.1% in most cases, as compared to a numerically exact line-by-line RT model. The accuracy is enhanced when the results are convolved to typical instrument resolutions. The operational speed and accuracy of UPCART can be further improved by optimizing binning schemes and parallelizing the codes, work
NASA Technical Reports Server (NTRS)
Du, Ping
1993-01-01
As a theoretical component of the joint effort with the laboratory of Dr. Lou Allamandola to search for potential candidates for interstellar organic carbon compound that are responsible for the visible diffuse interstellar absorption bands (DIB's), quantum mechanical calculations were performed on the electron absorption spectra of selected polycyclic aromatic hydrocarbons (PAH) and derivatives. In the completed project, 15 different species of naphthalene, its hydrogen abstraction and addition derivatives, and corresponding cations and anions were studied. Using semiempirical quantum mechanical method INDO/S, the ground electronic state of each species was evaluated with restricted Hartree-Fock scheme and limited configuration interaction. The lowest energy spin state for each species was used for electron absorption calculations. Results indicate that these calculations are accurate enough to reproduce the spectra of naphthalene cation and anion observed in neon matrix. The spectral pattern of the hydrogen abstraction and addition derivatives predicted based on these results indicate that the electron configuration of the pi orbitals of these species is the dominant determinant. A combined list of 19 absorptions calculated from 4500 A to 10,400 A were compiled and suggested as potential candidates that are relevant for the DIB's absorptions. Continued studies on pyrene and derivatives revealed the ground state symmetries and multiplicities of its neutral, anionic, and cationic species. Spectral calculations show that the cation (B(sub 3g)-2) and the anion (A(sub u)-2) are more likely to have low energy absorptions in the regions between 10 kK and 20 kK, similar to naphthalene. These absorptions, together with those to be determined from the hydrogen abstraction and addition derivatives of pyrene, can be used to provide additional candidates and suggest experimental work in the search for interstellar compounds that are responsible for DIB's.
Properties of Solar Thermal Fuels by Accurate Quantum Monte Carlo Calculations
NASA Astrophysics Data System (ADS)
Saritas, Kayahan; Ataca, Can; Grossman, Jeffrey C.
2014-03-01
Efficient utilization of the sun as a renewable and clean energy source is one of the major goals of this century due to increasing energy demand and environmental impact. Solar thermal fuels are materials that capture and store the sun's energy in the form of chemical bonds, which can then be released as heat on demand and charged again. Previous work on solar thermal fuels faced challenges related to the cyclability of the fuel over time, as well as the need for higher energy densities. Recently, it was shown that by templating photoswitches onto carbon nanostructures, both high energy density as well as high stability can be achieved. In this work, we explore alternative molecules to azobenzene in such a nano-templated system. We employ the highly accurate quantum Monte Carlo (QMC) method to predict the energy storage potential for each molecule. Our calculations show that in many cases the level of accuracy provided by density functional theory (DFT) is sufficient. However, in some cases, such as dihydroazulene, the drastic change in conjugation upon light absorption causes the DFT predictions to be inconsistent and incorrect. For this case, we compare our QMC results for the geometric structure, band gap and reaction enthalpy with different DFT functionals.
Aryl sulfoxide radical cations. Generation, spectral properties, and theoretical calculations.
Baciocchi, Enrico; Del Giacco, Tiziana; Gerini, Maria Francesca; Lanzalunga, Osvaldo
2006-08-17
Aromatic sulfoxide radical cations have been generated by pulse radiolysis and laser flash photolysis techniques. In water (pulse radiolysis) the radical cations showed an intense absorption band in the UV region (ca. 300 nm) and a broad less intense band in the visible region (from 500 to 1000 nm) whose position depends on the nature of the ring substituent. At very low pulse energy, the radical cations decayed by first-order kinetics, the decay rate increasing as the pH increases. It is suggested that the decay involves a nucleophilic attack of H(2)O or OH(-) (in basic solutions) to the positively charged sulfur atom to give the radical ArSO(OH)CH(3)(*). By sensitized [N-methylquinolinium tetrafluoborate (NMQ(+))] laser flash photolysis (LFP) the aromatic sulfoxide radical cations were generated in acetonitrile. In these experiments, however, only the band of the radical cation in the visible region could be observed, the UV band being covered by the UV absorption of NMQ(+). The lambda(max) values of the bands in the visible region resulted almost identical to those observed in water for the same radical cations. In the LFP experiments the sulfoxide radical cations decayed by second-order kinetics at a diffusion-controlled rate, and the decay is attributed to the back electron transfer between the radical cation and NMQ(*). DFT calculations were also carried out for a number of 4-X ring substituted (X = H, Me, Br, OMe, CN) aromatic sulfoxide radical cations (and their neutral parents). In all radical cations, the conformation with the S-O bond almost coplanar with the aromatic ring is the only one corresponding to the energy minimum. The maximum of energy corresponds to the conformation where the S-O bond is perpendicular to the aromatic ring. The rotational energy barriers are not very high, ranging from 3.9 to 6.9 kcal/mol. In all radical cations, the major fraction of charge and spin density is localized on the SOMe group. However, a substantial delocalization
Kouznetsov, Alexei; Tambasco, Mauro
2011-03-15
Purpose: To develop and validate a fast and accurate method that uses computed tomography (CT) voxel data to estimate absorbed radiation dose at a point of interest (POI) or series of POIs from a kilovoltage (kV) imaging procedure. Methods: The authors developed an approach that computes absorbed radiation dose at a POI by numerically evaluating the linear Boltzmann transport equation (LBTE) using a combination of deterministic and Monte Carlo (MC) techniques. This hybrid approach accounts for material heterogeneity with a level of accuracy comparable to the general MC algorithms. Also, the dose at a POI is computed within seconds using the Intel Core i7 CPU 920 2.67 GHz quad core architecture, and the calculations are performed using CT voxel data, making it flexible and feasible for clinical applications. To validate the method, the authors constructed and acquired a CT scan of a heterogeneous block phantom consisting of a succession of slab densities: Tissue (1.29 cm), bone (2.42 cm), lung (4.84 cm), bone (1.37 cm), and tissue (4.84 cm). Using the hybrid transport method, the authors computed the absorbed doses at a set of points along the central axis and x direction of the phantom for an isotropic 125 kVp photon spectral point source located along the central axis 92.7 cm above the phantom surface. The accuracy of the results was compared to those computed with MCNP, which was cross-validated with EGSnrc, and served as the benchmark for validation. Results: The error in the depth dose ranged from -1.45% to +1.39% with a mean and standard deviation of -0.12% and 0.66%, respectively. The error in the x profile ranged from -1.3% to +0.9%, with standard deviations of -0.3% and 0.5%, respectively. The number of photons required to achieve these results was 1x10{sup 6}. Conclusions: The voxel-based hybrid method evaluates the LBTE rapidly and accurately to estimate the absorbed x-ray dose at any POI or series of POIs from a kV imaging procedure.
NASA Astrophysics Data System (ADS)
Kripal, Ram; Yadav, Awadhesh Kumar
2015-06-01
Zero field splitting parameters (ZFSPs) D and E of Cr3+ ion doped ammonium oxalate monohydrate (AOM) are calculated with formula using the superposition model. The theoretically calculated ZFSPs for Cr3+ in AOM crystal are compared with the experimental value obtained by electron paramagnetic resonance (EPR). Theoretical ZFSPs are in good agreement with the experimental ones. The energy band positions of optical absorption spectra of Cr3+ in AOM crystal calculated with CFA package are in good match with the experimental values.
Accurate calculation and modeling of the adiabatic connection in density functional theory
NASA Astrophysics Data System (ADS)
Teale, A. M.; Coriani, S.; Helgaker, T.
2010-04-01
AC. When parametrized in terms of the same input data, the AC-CI model offers improved performance over the corresponding AC-D model, which is shown to be the lowest-order contribution to the AC-CI model. The utility of the accurately calculated AC curves for the analysis of standard density functionals is demonstrated for the BLYP exchange-correlation functional and the interaction-strength-interpolation (ISI) model AC integrand. From the results of this analysis, we investigate the performance of our proposed two-parameter AC-D and AC-CI models when a simple density functional for the AC at infinite interaction strength is employed in place of information at the fully interacting point. The resulting two-parameter correlation functionals offer a qualitatively correct behavior of the AC integrand with much improved accuracy over previous attempts. The AC integrands in the present work are recommended as a basis for further work, generating functionals that avoid spurious error cancellations between exchange and correlation energies and give good accuracy for the range of densities and types of correlation contained in the systems studied here.
Accurate and efficient calculation of discrete correlation functions and power spectra
NASA Astrophysics Data System (ADS)
Xu, Y. F.; Liu, J. M.; Zhu, W. D.
2015-07-01
Operational modal analysis (OMA), or output-only modal analysis, has been widely conducted especially when excitation applied on a structure is unknown or difficult to measure. Discrete cross-correlation functions and cross-power spectra between a reference data series and measured response data series are bases for OMA to identify modal properties of a structure. Such functions and spectra can be efficiently transformed from each other using the discrete Fourier transform (DFT) and inverse DFT (IDFT) based on the cross-correlation theorem. However, a direct application of the theorem and transforms, including the DFT and IDFT, can yield physically erroneous results due to periodic extension of the DFT on a function of a finite length to be transformed, which is false most of the time. Padding zero series to ends of data series before applying the theorem and transforms can reduce the errors, but the results are still physically erroneous. A new methodology is developed in this work to calculate discrete cross-correlation functions of non-negative time delays and associated cross-power spectra, referred to as half spectra, for OMA. The methodology can be extended to cross-correlation functions of any time delays and associated cross-power spectra, referred to as full spectra. The new methodology is computationally efficient due to use of the transforms. Data series are properly processed to avoid the errors caused by the periodic extension, and the resulting cross-correlation functions and associated cross-power spectra perfectly comply with their definitions. A coherence function, a convergence function, and a convergence index are introduced to evaluate qualities of measured cross-correlation functions and associated cross-power spectra. The new methodology was numerically and experimentally applied to an ideal two-degree-of-freedom (2-DOF) mass-spring-damper system and a damaged aluminum beam, respectively, and OMA was conducted using half spectra to estimate
NASA Astrophysics Data System (ADS)
Shi, De-Heng; Liu, Yu-Fang; Sun, Jin-Feng; Zhu, Zun-Lue; Yang, Xiang-Dong
2006-12-01
The reasonable dissociation limit of the second excited singlet state B1Π of 7LiH molecule is obtained. The accurate dissociation energy and equilibrium geometry of the B1Π state are calculated using a symmetry-adapted-cluster configuration-interaction method in full active space. The whole potential energy curve for the B1Π state is obtained over the internuclear distance ranging from about 0.10 nm to 0.54 nm, and has a least-square fit to the analytic Murrell-Sorbie function form. The vertical excitation energy is calculated from the ground state to the B1Π state and compared with previous theoretical results. The equilibrium internuclear distance obtained by geometry optimization is found to be quite different from that obtained by single-point energy scanning under the same calculation condition. Based on the analytic potential energy function, the harmonic frequency value of the B1Π state is estimated. A comparison of the theoretical calculations of dissociation energies, equilibrium interatomic distances and the analytic potential energy function with those obtained by previous theoretical results clearly shows that the present work is more comprehensive and in better agreement with experiments than previous theories, thus it is an improvement on previous theories.
Tung, Wei-Cheng; Adamowicz, Ludwik
2014-03-28
Very accurate calculations of the ground-state potential energy curve (PEC) of the LiH{sup +} ion performed with all-electron explicitly correlated Gaussian functions with shifted centers are presented. The variational method is employed. The calculations involve optimization of nonlinear exponential parameters of the Gaussians performed with the aid of the analytical first derivatives of the energy determined with respect to the parameters. The diagonal adiabatic correction is also calculated for each PEC point. The PEC is then used to calculate the vibrational energies of the system. In that calculation, the non-adiabatic effects are accounted for by using an effective vibrational mass obtained by the minimization of the difference between the vibrational energies obtained from the calculations where the Born-Oppenheimer approximation was not assumed and the results of the present calculations.
NASA Astrophysics Data System (ADS)
Moreira, Rafael Y. O.; Brasil, Davi S. B.; Alves, Cláudio N.; Guilhon, Giselle M. S. P.; Santos, Lourivaldo S.; Arruda, Mara S. P.; Müller, Adolfo H.; Barbosa, Patrícia S.; Abreu, Alcicley S.; Silva, Edilene O.; Rumjanek, Victor M.; Souza, Jaime, Jr.; da Silva, Albérico B. F.; Santos, Regina H. De A.
Julocrotine, N-(2,6-dioxo-1-phenethyl-piperidin-3-yl)-2-methyl-butyramide, is a potent antiproliferative agent against the promastigote and amastigote forms of Leishmania amazonensis (L.). In this work, the crystal structure of Julocrotine was solved by X-ray diffraction, and its geometrical parameters were compared with theoretical calculations at the B3LYP and HF level of theory. IR and NMR spectra also have been obtained and compared with theoretical calculations. IR absorptions calculated with the B3LYP level of theory employed together with the 6-311G+(d,p) basis set, are close to those observed experimentally. Theoretical NMR calculations show little deviation from experimental results. The results show that the theory is in accordance with the experimental data.0
DICOM organ dose does not accurately represent calculated dose in mammography
NASA Astrophysics Data System (ADS)
Suleiman, Moayyad E.; Brennan, Patrick C.; McEntee, Mark F.
2016-03-01
This study aims to analyze the agreement between the mean glandular dose estimated by the mammography unit (organ dose) and mean glandular dose calculated using Dance et al published method (calculated dose). Anonymised digital mammograms from 50 BreastScreen NSW centers were downloaded and exposure information required for the calculation of dose was extracted from the DICOM header along with the organ dose estimated by the system. Data from quality assurance annual tests for the included centers were collected and used to calculate the mean glandular dose for each mammogram. Bland-Altman analysis and a two-tailed paired t-test were used to study the agreement between calculated and organ dose and the significance of any differences. A total of 27,869 dose points from 40 centers were included in the study, mean calculated dose and mean organ dose (+/- standard deviation) were 1.47 (+/-0.66) and 1.38 (+/-0.56) mGy respectively. A statistically significant 0.09 mGy bias (t = 69.25; p<0.0001) with 95% limits of agreement between calculated and organ doses ranging from -0.34 and 0.52 were shown by Bland-Altman analysis, which indicates a small yet highly significant difference between the two means. The use of organ dose for dose audits is done at the risk of over or underestimating the calculated dose, hence, further work is needed to identify the causal agents for differences between organ and calculated doses and to generate a correction factor for organ dose.
Baer, M.R.; Hobbs, M.L.; McGee, B.C.
1998-11-03
Exponential-13,6 (EXP-13,6) potential pammeters for 750 gases composed of 48 elements were determined and assembled in a database, referred to as the JCZS database, for use with the Jacobs Cowperthwaite Zwisler equation of state (JCZ3-EOS)~l) The EXP- 13,6 force constants were obtained by using literature values of Lennard-Jones (LJ) potential functions, by using corresponding states (CS) theory, by matching pure liquid shock Hugoniot data, and by using molecular volume to determine the approach radii with the well depth estimated from high-pressure isen- tropes. The JCZS database was used to accurately predict detonation velocity, pressure, and temperature for 50 dif- 3 Accurate predictions were also ferent explosives with initial densities ranging from 0.25 glcm3 to 1.97 g/cm . obtained for pure liquid shock Hugoniots, static properties of nitrogen, and gas detonations at high initial pressures.
Myint, P. C.; Hao, Y.; Firoozabadi, A.
2015-03-27
Thermodynamic property calculations of mixtures containing carbon dioxide (CO_{2}) and water, including brines, are essential in theoretical models of many natural and industrial processes. The properties of greatest practical interest are density, solubility, and enthalpy. Many models for density and solubility calculations have been presented in the literature, but there exists only one study, by Spycher and Pruess, that has compared theoretical molar enthalpy predictions with experimental data [1]. In this report, we recommend two different models for enthalpy calculations: the CPA equation of state by Li and Firoozabadi [2], and the CO_{2} activity coefficient model by Duan and Sun [3]. We show that the CPA equation of state, which has been demonstrated to provide good agreement with density and solubility data, also accurately calculates molar enthalpies of pure CO_{2}, pure water, and both CO_{2}-rich and aqueous (H_{2}O-rich) mixtures of the two species. It is applicable to a wider range of conditions than the Spycher and Pruess model. In aqueous sodium chloride (NaCl) mixtures, we show that Duan and Sun’s model yields accurate results for the partial molar enthalpy of CO_{2}. It can be combined with another model for the brine enthalpy to calculate the molar enthalpy of H_{2}O-CO_{2}-NaCl mixtures. We conclude by explaining how the CPA equation of state may be modified to further improve agreement with experiments. This generalized CPA is the basis of our future work on this topic.
Wills, John M; Mattsson, Ann E
2012-06-06
Brooks, Johansson, and Skriver, using the LMTO-ASA method and considerable insight, were able to explain many of the ground state properties of the actinides. In the many years since this work was done, electronic structure calculations of increasing sophistication have been applied to actinide elements and compounds, attempting to quantify the applicability of DFT to actinides and actinide compounds and to try to incorporate other methodologies (i.e. DMFT) into DFT calculations. Through these calculations, the limits of both available density functionals and ad hoc methodologies are starting to become clear. However, it has also become clear that approximations used to incorporate relativity are not adequate to provide rigorous tests of the underlying equations of DFT, not to mention ad hoc additions. In this talk, we describe the result of full-potential LMTO calculations for the elemental actinides, comparing results obtained with a full Dirac basis with those obtained from scalar-relativistic bases, with and without variational spin-orbit. This comparison shows that the scalar relativistic treatment of actinides does not have sufficient accuracy to provide a rigorous test of theory and that variational spin-orbit introduces uncontrolled errors in the results of electronic structure calculations on actinide elements.
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Langhoff, Stephen R.; Arnold, James O. (Technical Monitor)
1996-01-01
The vibrational frequencies and infrared intensities of naphthalene neutral and cation are studied at the self-consistent-field (SCF), second-order Moller-Plesset (MP2), and density functional theory (DFT) levels using a variety of one-particle basis sets. Very accurate frequencies can be obtained at the DFT level in conjunction with large basis sets if they are scaled with two factors, one for the C-H stretches and a second for all other modes. We also find remarkably good agreement at the B3LYP/4-31G level using only one scale factor. Unlike the neutral PAHs where all methods do reasonably well for the intensities, only the DFT results are accurate for the PAH cations. The failure of the SCF and MP2 methods is caused by symmetry breaking and an inability to describe charge delocalization. We present several interesting cases of symmetry breaking in this study. An assessment is made as to whether an ensemble of PAH neutrals or cations could account for the unidentified infrared bands observed in many astronomical sources.
Do Bond Functions Help for the Calculation of Accurate Bond Energies?
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Arnold, James (Technical Monitor)
1998-01-01
The bond energies of 8 chemically bound diatomics are computed using several basis sets with and without bond functions (BF). The bond energies obtained using the aug-pVnZ+BF basis sets (with a correction for basis set superposition error, BSSE) tend to be slightly smaller that the results obtained using the aug-pV(n+I)Z basis sets, but slightly larger than the BSSE corrected aug-pV(n+I)Z results. The aug-cc-pVDZ+BF and aug-cc-pVTZ+BF basis sets yield reasonable estimates of bond energies, but, in most cases, these results cannot be considered highly accurate. Extrapolation of the results obtained with basis sets including bond functions appears to be inferior to the results obtained by extrapolation using atom-centered basis sets. Therefore bond functions do not appear to offer a path for obtaining highly accurate results for chemically bound systems at a lower computational cost than atom centered basis sets.
PyVCI: A flexible open-source code for calculating accurate molecular infrared spectra
NASA Astrophysics Data System (ADS)
Sibaev, Marat; Crittenden, Deborah L.
2016-06-01
The PyVCI program package is a general purpose open-source code for simulating accurate molecular spectra, based upon force field expansions of the potential energy surface in normal mode coordinates. It includes harmonic normal coordinate analysis and vibrational configuration interaction (VCI) algorithms, implemented primarily in Python for accessibility but with time-consuming routines written in C. Coriolis coupling terms may be optionally included in the vibrational Hamiltonian. Non-negligible VCI matrix elements are stored in sparse matrix format to alleviate the diagonalization problem. CPU and memory requirements may be further controlled by algorithmic choices and/or numerical screening procedures, and recommended values are established by benchmarking using a test set of 44 molecules for which accurate analytical potential energy surfaces are available. Force fields in normal mode coordinates are obtained from the PyPES library of high quality analytical potential energy surfaces (to 6th order) or by numerical differentiation of analytic second derivatives generated using the GAMESS quantum chemical program package (to 4th order).
Validation of the new code package APOLLO2.8 for accurate PWR neutronics calculations
Santamarina, A.; Bernard, D.; Blaise, P.; Leconte, P.; Palau, J. M.; Roque, B.; Vaglio, C.; Vidal, J. F.
2013-07-01
This paper summarizes the Qualification work performed to demonstrate the accuracy of the new APOLLO2.S/SHEM-MOC package based on JEFF3.1.1 nuclear data file for the prediction of PWR neutronics parameters. This experimental validation is based on PWR mock-up critical experiments performed in the EOLE/MINERVE zero-power reactors and on P.I. Es on spent fuel assemblies from the French PWRs. The Calculation-Experiment comparison for the main design parameters is presented: reactivity of UOX and MOX lattices, depletion calculation and fuel inventory, reactivity loss with burnup, pin-by-pin power maps, Doppler coefficient, Moderator Temperature Coefficient, Void coefficient, UO{sub 2}-Gd{sub 2}O{sub 3} poisoning worth, Efficiency of Ag-In-Cd and B4C control rods, Reflector Saving for both standard 2-cm baffle and GEN3 advanced thick SS reflector. From this qualification process, calculation biases and associated uncertainties are derived. This code package APOLLO2.8 is already implemented in the ARCADIA new AREVA calculation chain for core physics and is currently under implementation in the future neutronics package of the French utility Electricite de France. (authors)
Frolov, Andrey I
2015-05-12
Accurate calculation of solvation free energies (SFEs) is a fundamental problem of theoretical chemistry. In this work we perform a careful validation of the theory of solutions in energy representation (ER method) developed by Matubayasi et al. [J. Chem. Phys. 2000, 113, 6070-6081] for SFE calculations in supercritical solvents. This method can be seen as a bridge between the molecular simulations and the classical (not quantum) density functional theory (DFT) formulated in energy representation. We performed extensive calculations of SFEs of organic molecules of different chemical natures in pure supercritical CO2 (sc-CO2) and in sc-CO2 with addition of 6 mol % of ethanol, acetone, and n-hexane as cosolvents. We show that the ER method reproduces SFE data calculated by a method free of theoretical approximations (the Bennett's acceptance ratio) with the mean absolute error of only 0.05 kcal/mol. However, the ER method requires by an order less computational resources. Also, we show that the quality of ER calculations should be carefully monitored since the lack of sampling can result into a considerable bias in predictions. The present calculations reproduce the trends in the cosolvent-induced solubility enhancement factors observed in experimental data. Thus, we think that molecular simulations coupled with the ER method can be used for quick calculations of the effect of variation of temperature, pressure, and cosolvent concentration on SFE and hence solubility of bioactive compounds in supercritical fluids. This should dramatically reduce the burden of experimental work on optimizing solvency of supercritical solvents. PMID:26574423
Theoretical calculation of light-induced forces and torques on complex microrotors
NASA Astrophysics Data System (ADS)
Liu, Yuxiang; Zhu, Anding; Huang, Wenhao
2004-12-01
In this letter, we propose the new theoretical investigation on the optical forces and torques on complex microrotors. On the basis of R. C. Gauthier"s hybrid ray-wave model, the optical forces and toques on two complex asymmetric micro-objects, the conical microrotor and the helical microrotor, are analyzed. The viscous drag torque is estimated by Stokes flow to obtain the rotational speed. The results of our computation agree well with the previously published experimental results, which indicates that our approach of the optical torque calculation is suitable for other complex microrotors and that the theoretical calculation is very helpful to optimum design of light-driven microrotors.
Kuwada, Takeshi; Fukui, Miyako; Hata, Toshiyuki; Choshi, Tominari; Nobuhiro, Junko; Ono, Yukio; Hibino, Satoshi
2003-01-01
The specific optical rotations of (R)-oxopropaline D calculated by two ab initio MO methods were +52+/-31 degrees and +61+/-29 degrees, respectively, and (+)-oxopropaline D (3) was presumed to have an R-configuration. On the basis of this theoretical result, the reaction of 1-litio-beta-carboline with (R)-glyceraldehyde acetonide followed by oxidation with MnO(2) gave (R)-oxopropaline D acetonide (4a), which was consistent with the previously synthesized (+)-oxopropaline D acetonide (4) in all respects. From the results of theoretical calculations and the experimental synthesis, we determined that natural (+)-oxopropaline D (3) has an R-configuration. PMID:12520122
NASA Technical Reports Server (NTRS)
Boughner, Robert E.
1986-01-01
A method for calculating the photodissociation rates needed for photochemical modeling of the stratosphere, which includes the effects of molecular scattering, is described. The procedure is based on Sokolov's method of averaging functional correction. The radiation model and approximations used to calculate the radiation field are examined. The approximated diffuse fields and photolysis rates are compared with exact data. It is observed that the approximate solutions differ from the exact result by 10 percent or less at altitudes above 15 km; the photolysis rates differ from the exact rates by less than 5 percent for altitudes above 10 km and all zenith angles, and by less than 1 percent for altitudes above 15 km.
Santiago, Régis Tadeu; Haiduke, Roberto Luiz Andrade
2015-10-30
This research provides a performance investigation of density functional theory and also proposes new functional parameterizations to deal with electric field gradient (EFG) calculations at nuclear positions. The entire procedure is conducted within the four-component formalism. First, we noticed that traditional hybrid and long-range corrected functionals are more efficient in the description of EFG variations for a set of elements (indium, antimony, iodine, lutetium, and hafnium) among linear molecules. Thus, we selected the PBE0, B3LYP, and CAM-B3LYP functionals and promoted a reoptimization of their parameters for a better description of these EFG changes. The PBE0q variant developed here showed an overall promising performance in a validation test conducted with potassium, iodine, copper, and gold. In general, the correlation coefficients found in linear regressions between experimental nuclear quadrupole coupling constants and calculated EFGs are improved while the systematic EFG errors also decrease as a result of this reparameterization. PMID:26284820
Accurate Calculation of Oscillator Strengths for CI II Lines Using Non-orthogonal Wavefunctions
NASA Technical Reports Server (NTRS)
Tayal, S. S.
2004-01-01
Non-orthogonal orbitals technique in the multiconfiguration Hartree-Fock approach is used to calculate oscillator strengths and transition probabilities for allowed and intercombination lines in Cl II. The relativistic corrections are included through the Breit-Pauli Hamiltonian. The Cl II wave functions show strong term dependence. The non-orthogonal orbitals are used to describe the term dependence of radial functions. Large sets of spectroscopic and correlation functions are chosen to describe adequately strong interactions in the 3s(sup 2)3p(sup 3)nl (sup 3)Po, (sup 1)Po and (sup 3)Do Rydberg series and to properly account for the important correlation and relaxation effects. The length and velocity forms of oscillator strength show good agreement for most transitions. The calculated radiative lifetime for the 3s3p(sup 5) (sup 3)Po state is in good agreement with experiment.
Accurate calculations of the high-pressure elastic constants based on the first-principles
NASA Astrophysics Data System (ADS)
Wang, Chen-Ju; Gu, Jian-Bing; Kuang, Xiao-Yu; Yang, Xiang-Dong
2015-08-01
The energy term corresponding to the first order of the strain in Taylor series expansion of the energy with respect to strain is always ignored when high-pressure elastic constants are calculated. Whether the modus operandi would affect the results of the high-pressure elastic constants is still unsolved. To clarify this query, we calculate the high-pressure elastic constants of tantalum and rhenium when the energy term mentioned above is considered and neglected, respectively. Results show that the neglect of the energy term corresponding to the first order of the strain indeed would influence the veracity of the high-pressure elastic constants, and this influence becomes larger with pressure increasing. Therefore, the energy term corresponding to the first-order of the strain should be considered when the high-pressure elastic constants are calculated. Project supported by the National Natural Science Foundation of China (Grant No. 11274235), the Young Scientist Fund of the National Natural Science Foundation of China (Grant No. 11104190), and the Doctoral Education Fund of Education Ministry of China (Grant Nos. 20100181110086 and 20110181120112).
Ab-initio Calculations of Accurate Electronic Properties of ZnS
NASA Astrophysics Data System (ADS)
Khamala, Bethuel; Franklin, Loushanda; Malozovski, Yuriy; Stewart, Anthony; Bagayoko, Diola; Bagayoko Research Group Team
2014-03-01
We present the results from ab-initio, self consistent, local density approximation (LDA) calculations of the electronic and related properties of zinc-blende zinc sulphide (zb-ZnS). We employed the Ceperley and Alder LDA potential and the linear combination of atomic orbital (LCAO) formalism in our non-relativistic computations. The implementation of the LCAO formalism followed the Bagayoko, Zhao, and Williams method as enhanced by Ekuma and Franklin (BZW-EF). The BZW-EF method includes a methodical search for the optimal basis set that yields the minima of the occupied energies. This search entails increasing the size of the basis set and related modifications of angular symmetry and of radial orbitals. Our calculated, direct gap of 3.725 eV, at the Γ point, is in excellent agreement with experiment. We have also calculated the total (DOS) and partial (pDOS) densities of states, electron and hole effective masses and total energies that agree very well with available, corresponding experimental results. Acknowledgement: This research is funded in part by the National Science Foundation (NSF) and the Louisiana Board of Regents, through LASiGMA [Award Nos. EPS- 1003897, NSF (2010-15)-RII-SUBR] and NSF HRD-1002541, the US Department of Energy - National, Nuclear Security Administration (NNSA) (Award No. DE-NA0001861), LaSPACE, and LONI-SUBR.
Ab-initio Calculations of Accurate Electronic Properties of Wurzite AlN
NASA Astrophysics Data System (ADS)
Nwigboji, Ifeanyi; Malozovsky, Yuriy; Bagayoko, Diola; Bagayoko Research Group Team
2014-03-01
We present results from ab-initio, self consistent local density approximation (LDA) calculations of electronic and related properties of wurtzite Aluminum Nitride (w-AlN). Our non-relativistic computations employed the Ceperley and Alder LDA potential and the linear combination of atomic orbital (LCAO) formalism. The implementation of the LCAO formalism followed the Bagayoko, Zhao, and Williams' method as enhanced by Ekuma and Franklin (BZW-EF). The BZW-EF method verifiably obtains the minima of the occupied energies; these minima provide the most variationally and physically valid density functional theory (DFT) description of the ground states of materials under study. Our preliminary results for w-AlN show that w-AlN has a direct band gap of 5.82 eV at the Γ point. The preliminary energy bands were obtained with a basis set comprising 48 functions. None of the several, larger basis sets tested to date led to occupied energies lower than those obtained with the above 48. While most previous LDA calculations are 2 eV smaller or more than the experimental value of 5.9 eV that is in excellent agreement with our finding, considering the typical experimental uncertainty of 0.2 eV for absorption measurements on AlN. We also discuss our calculated density of states (DOS) and partial densities of states (pDOS).
Accurate calculation of Stokes drag for point-particle tracking in two-way coupled flows
NASA Astrophysics Data System (ADS)
Horwitz, J. A. K.; Mani, A.
2016-08-01
In this work, we propose and test a method for calculating Stokes drag applicable to particle-laden fluid flows where two-way momentum coupling is important. In the point-particle formulation, particle dynamics are coupled to fluid dynamics via a source term that appears in the respective momentum equations. When the particle Reynolds number is small and the particle diameter is smaller than the fluid scales, it is common to approximate the momentum coupling source term as the Stokes drag. The Stokes drag force depends on the difference between the undisturbed fluid velocity evaluated at the particle location, and the particle velocity. However, owing to two-way coupling, the fluid velocity is modified in the neighborhood of a particle, relative to its undisturbed value. This causes the computed Stokes drag force to be underestimated in two-way coupled point-particle simulations. We develop estimates for the drag force error as function of the particle size relative to the grid size. Because the disturbance field created by the particle contaminates the surrounding fluid, correctly calculating the drag force cannot be done solely by direct interpolation of the fluid velocity. Instead, we develop a correction method that calculates the undisturbed fluid velocity from the computed disturbed velocity field by adding an estimate of the velocity disturbance created by the particle. The correction scheme is tested for a particle settling in an otherwise quiescent fluid and is found to reduce the error in computed settling velocity by an order of magnitude compared with common interpolation schemes.
Shaughnessy, M C; Jones, R E
2016-02-01
We develop and demonstrate a method to efficiently use density functional calculations to drive classical dynamics of complex atomic and molecular systems. The method has the potential to scale to systems and time scales unreachable with current ab initio molecular dynamics schemes. It relies on an adapting dataset of independently computed Hellmann-Feynman forces for atomic configurations endowed with a distance metric. The metric on configurations enables fast database lookup and robust interpolation of the stored forces. We discuss mechanisms for the database to adapt to the needs of the evolving dynamics, while maintaining accuracy, and other extensions of the basic algorithm. PMID:26669825
A new class of atomic basis functions for accurate electronic structure calculations of molecules
NASA Astrophysics Data System (ADS)
Laikov, Dimitri N.
2005-11-01
A new general approach is developed for obtaining systematic sequences of atomic single-particle basis sets for use in correlated electronic structure calculations of molecules. All the constituent functions are defined as the solutions of variational problems and are of three types: a minimal Hartree-Fock set, additional functions to represent low-lying excited configurations, and general functions for describing electron correlation. The latter are determined to minimize a functional derived from the closed-shell second-order correlation energy expression. Generally-contracted Gaussian expansions are developed to approximate these general functions in the non-relativistic case and within a scalar-relativistic approximation.
NASA Astrophysics Data System (ADS)
Orwat, Justyna; Mielimaka, Ryszard
2016-06-01
The article demonstrates the comparison of prognosed values of deformation indicators with their values obtained from geodetic measurements. The value of correlation coefficient R between its theoretical and practical values was calculated for each indicator. Following the measured values as basic ones, the evaluation of efficiency of prognosis via the use of Knothe`s theory (a mathematical model, in which lowering of point is calculated as a double integral from Gauss error function with adequate parameters).
Accurate calculation of binding energies for molecular clusters - Assessment of different models
NASA Astrophysics Data System (ADS)
Friedrich, Joachim; Fiedler, Benjamin
2016-06-01
In this work we test different strategies to compute high-level benchmark energies for medium-sized molecular clusters. We use the incremental scheme to obtain CCSD(T)/CBS energies for our test set and carefully validate the accuracy for binding energies by statistical measures. The local errors of the incremental scheme are <1 kJ/mol. Since they are smaller than the basis set errors, we obtain higher total accuracy due to the applicability of larger basis sets. The final CCSD(T)/CBS benchmark values are ΔE = - 278.01 kJ/mol for (H2O)10, ΔE = - 221.64 kJ/mol for (HF)10, ΔE = - 45.63 kJ/mol for (CH4)10, ΔE = - 19.52 kJ/mol for (H2)20 and ΔE = - 7.38 kJ/mol for (H2)10 . Furthermore we test state-of-the-art wave-function-based and DFT methods. Our benchmark data will be very useful for critical validations of new methods. We find focal-point-methods for estimating CCSD(T)/CBS energies to be highly accurate and efficient. For foQ-i3CCSD(T)-MP2/TZ we get a mean error of 0.34 kJ/mol and a standard deviation of 0.39 kJ/mol.
NASA Astrophysics Data System (ADS)
Komsa, Hannu-Pekka; Berseneva, Natalia; Krasheninnikov, Arkady V.; Nieminen, Risto M.
2014-07-01
Impurities and defects frequently govern materials properties, with the most prominent example being the doping of bulk semiconductors where a minute amount of foreign atoms can be responsible for the operation of the electronic devices. Several computational schemes based on a supercell approach have been developed to get insights into types and equilibrium concentrations of point defects, which successfully work in bulk materials. Here, we show that many of these schemes cannot directly be applied to two-dimensional (2D) systems, as formation energies of charged point defects are dominated by large spurious electrostatic interactions between defects in inhomogeneous environments. We suggest two approaches that solve this problem and give accurate formation energies of charged defects in 2D systems in the dilute limit. Our methods, which are applicable to all kinds of charged defects in any 2D system, are benchmarked for impurities in technologically important h-BN and MoS2 2D materials, and they are found to perform equally well for substitutional and adatom impurities.
NASA Technical Reports Server (NTRS)
Mcgrath, W. R.; Richards, P. L.; Face, D. W.; Prober, D. E.; Lloyd, F. L.
1988-01-01
A systematic study of the gain and noise in superconductor-insulator-superconductor mixers employing Ta based, Nb based, and Pb-alloy based tunnel junctions was made. These junctions displayed both weak and strong quantum effects at a signal frequency of 33 GHz. The effects of energy gap sharpness and subgap current were investigated and are quantitatively related to mixer performance. Detailed comparisons are made of the mixing results with the predictions of a three-port model approximation to the Tucker theory. Mixer performance was measured with a novel test apparatus which is accurate enough to allow for the first quantitative tests of theoretical noise predictions. It is found that the three-port model of the Tucker theory underestimates the mixer noise temperature by a factor of about 2 for all of the mixers. In addition, predicted values of available mixer gain are in reasonable agreement with experiment when quantum effects are weak. However, as quantum effects become strong, the predicted available gain diverges to infinity, which is in sharp contrast to the experimental results. Predictions of coupled gain do not always show such divergences.
Rucker, Robert; Oelschlaeger, Peter; Warshel, Arieh
2010-01-01
DNA polymerase β (pol β) is a small eukaryotic enzyme with the ability to repair short single-stranded DNA gaps that has found use as a model system for larger replicative DNA polymerases. For all DNA polymerases, the factors determining their catalytic power and fidelity are the interactions between the bases of the base pair, amino acids near the active site, and the two magnesium ions. In this report, we study effects of all three aspects on human pol β transition state (TS) binding free energies by reproducing a consistent set of experimentally determined data for different structures. Our calculations comprise the combination of four different base pairs (incoming pyrimidine nucleotides incorporated opposite both matched and mismatched purines) with four different pol β structures (wild type and three separate mutations of ionized residues to alanine). We decompose the incoming deoxynucleoside 5′-triphosphate-TS, and run separate calculations for the neutral base part and the highly charged triphosphate part, using different dielectric constants in order to account for the specific electrostatic environments. This new approach improves our ability to predict the effect of matched and mismatched base pairing and of mutations in DNA polymerases on fidelity and may be a useful tool in studying the potential of DNA polymerase mutations in the development of cancer. It also supports our point of view with regards to the origin of the structural control of fidelity, allowing for a quantified description of the fidelity of DNA polymerases. PMID:19842163
Halverson, Thomas; Poirier, Bill
2012-12-14
In a series of earlier articles [B. Poirier, J. Theor. Comput. Chem. 2, 65 (2003); B. Poirier and A. Salam, J. Chem. Phys. 121, 1690 (2004); and ibid. 121, 1704 (2004)], a new method was introduced for performing exact quantum dynamics calculations. The method uses a 'weylet' basis set (orthogonalized Weyl-Heisenberg wavelets) combined with phase space truncation, to defeat the exponential scaling of CPU effort with system dimensionality-the first method ever able to achieve this long-standing goal. Here, we develop another such method, which uses a much more convenient basis of momentum-symmetrized Gaussians. Despite being non-orthogonal, symmetrized Gaussians are collectively local, allowing for effective phase space truncation. A dimension-independent code for computing energy eigenstates of both coupled and uncoupled systems has been created, exploiting massively parallel algorithms. Results are presented for model isotropic uncoupled harmonic oscillators and coupled anharmonic oscillators up to 27 dimensions. These are compared with the previous weylet calculations (uncoupled harmonic oscillators up to 15 dimensions), and found to be essentially just as efficient. Coupled system results are also compared to corresponding exact results obtained using a harmonic oscillator basis, and also to approximate results obtained using first-order perturbation theory up to the maximum dimensionality for which the latter may be feasibly obtained (four dimensions).
Halverson, Thomas; Poirier, Bill
2012-12-14
In a series of earlier articles [B. Poirier, J. Theor. Comput. Chem. 2, 65 (2003); B. Poirier and A. Salam, J. Chem. Phys. 121, 1690 (2004); and ibid. 121, 1704 (2004)], a new method was introduced for performing exact quantum dynamics calculations. The method uses a "weylet" basis set (orthogonalized Weyl-Heisenberg wavelets) combined with phase space truncation, to defeat the exponential scaling of CPU effort with system dimensionality--the first method ever able to achieve this long-standing goal. Here, we develop another such method, which uses a much more convenient basis of momentum-symmetrized Gaussians. Despite being non-orthogonal, symmetrized Gaussians are collectively local, allowing for effective phase space truncation. A dimension-independent code for computing energy eigenstates of both coupled and uncoupled systems has been created, exploiting massively parallel algorithms. Results are presented for model isotropic uncoupled harmonic oscillators and coupled anharmonic oscillators up to 27 dimensions. These are compared with the previous weylet calculations (uncoupled harmonic oscillators up to 15 dimensions), and found to be essentially just as efficient. Coupled system results are also compared to corresponding exact results obtained using a harmonic oscillator basis, and also to approximate results obtained using first-order perturbation theory up to the maximum dimensionality for which the latter may be feasibly obtained (four dimensions). PMID:23248981
NASA Astrophysics Data System (ADS)
Halverson, Thomas; Poirier, Bill
2012-12-01
In a series of earlier articles [B. Poirier, J. Theor. Comput. Chem. 2, 65 (2003);, 10.1142/S0219633603000380 B. Poirier and A. Salam, J. Chem. Phys. 121, 1690 (2004);, 10.1063/1.1767511 B. Poirier and A. Salam, J. Chem. Phys. 121, 1704 (2004), 10.1063/1.1767512], a new method was introduced for performing exact quantum dynamics calculations. The method uses a "weylet" basis set (orthogonalized Weyl-Heisenberg wavelets) combined with phase space truncation, to defeat the exponential scaling of CPU effort with system dimensionality—the first method ever able to achieve this long-standing goal. Here, we develop another such method, which uses a much more convenient basis of momentum-symmetrized Gaussians. Despite being non-orthogonal, symmetrized Gaussians are collectively local, allowing for effective phase space truncation. A dimension-independent code for computing energy eigenstates of both coupled and uncoupled systems has been created, exploiting massively parallel algorithms. Results are presented for model isotropic uncoupled harmonic oscillators and coupled anharmonic oscillators up to 27 dimensions. These are compared with the previous weylet calculations (uncoupled harmonic oscillators up to 15 dimensions), and found to be essentially just as efficient. Coupled system results are also compared to corresponding exact results obtained using a harmonic oscillator basis, and also to approximate results obtained using first-order perturbation theory up to the maximum dimensionality for which the latter may be feasibly obtained (four dimensions).
NASA Astrophysics Data System (ADS)
Wang, Xinxin; Shi, Deheng; Sun, Jinfeng; Zhu, Zunlue
2016-08-01
The potential energy curves were calculated for the 21 states (X2Π, A2Π, 32Π, 42Π, 52Π, 12Σ+, 22Σ+, 32Σ+, 12Σ-, 22Σ-, 32Σ-, 12Δ, 22Δ, 32Δ, 12Φ, 14Σ+, a4Σ-, 24Σ-, 14Π, 24Π and 14Δ), which originated from the two lowest dissociation channels of ClO radical. The calculations were done for internuclear separations approximately from 0.08 to 1.10 nm using the CASSCF method, which was followed by the icMRCI approach with the aug-cc-pV5Z basis set. Of these 21 states, the 14Π, 24Π, 32Δ, 42Π, 52Π, 12Φ, 32Σ+, 14Δ and 24Σ- states are repulsive. The 12Δ, 12Σ-, 14Σ+, 22Σ-, 12Σ+, 22Σ+, 22Δ and 32Σ- states are very weakly bound. Only the A2Π state has one barrier. The avoided crossing exists between the A2Π and the 32Π state. However, the avoided crossing does not generate any double wells. Core- valence correlation correction was accounted for at the level of an aug-cc-pCVQZ basis set. Scalar relativistic correction was included by the third-order Douglas-Kroll Hamiltonian approximation at the level of an aug-cc-pVQZ basis set. All the potential energy curves were extrapolated to the complete basis set limit. The spectroscopic parameters were determined. The 12Σ-, 22Σ-, 32Σ- and 14Σ+ states may be very difficult to be detected in an experiment, since each of these Λ-S states has only one or two vibrational states. The Franck-Condon factors and radiative lifetimes were calculated for several low vibrational levels of the A2Π - X2Π, 32Π - a4Σ-, 22Δ - a4Σ- and 32Σ- - 12Σ- transitions. The spin-orbit coupling effect on the spectroscopic parameters of the X2Π, A2Π, 32Π, a4Σ- and 22Σ+ states were discussed. The spectroscopic properties reported here can be expected to be reliably predicted ones.
NASA Astrophysics Data System (ADS)
Bačić, Z.
1991-09-01
We show that the triatomic adiabatic vibrational eigenstates (AVES) provide a convenient basis for accurate discrete variable representation (DVR) calculation and automatic assignment of highly excited, large amplitude motion vibrational states of floppy triatomic molecules. The DVR-AVES states are eigenvectors of the diagonal (in the stretch states) blocks of the adiabatically rearranged triatomic DVR-ray eigenvector (DVR-REV) Hamiltonian [J. C. Light and Z. Bačić, J. Chem. Phys. 87, 4008 (1987)]. The transformation of the full triatomic vibrational Hamiltonian from the DVR-REV basis to the new DVR-AVES basis is simple, and does not involve calculation of any new matrix elements. No dynamical approximation is made in the energy level calculation by the DVR-AVES approach; its accuracy and efficiency are identical to those of the DVR-REV method. The DVR-AVES states, as the adiabatic approximation to the vibrational states of a triatomic molecule, are labeled by three vibrational quantum numbers. Consequently, accurate large amplitude motion vibrational levels obtained by diagonalizing the full vibrational Hamiltonian transformed to the DVR-AVES basis, can be assigned automatically by the code, with the three quantum numbers of the dominant DVR-AVES state associated with the largest (by modulus) eigenvector element in the DVR-AVES basis. The DVR-AVES approach is used to calculate accurate highly excited localized and delocalized vibrational levels of HCN/HNC and LiCN/LiNC. A significant fraction of localized states of both systems, below and above the isomerization barrier, is assigned automatically, without inspection of wave function plots or separate approximate calculations.
Highly accurate incremental CCSD(T) calculations on aqua- and amine-complexes
NASA Astrophysics Data System (ADS)
Anacker, Tony; Friedrich, Joachim
2013-07-01
In this work, the accuracy of the second-order incremental expansion using the domain-specific basis set approach is tested for 20 cationic metal-aqua and 25 cationic metal-amine complexes. The accuracy of the approach is analysed by the statistical measures range, arithmetic mean, mean absolute deviation, root mean square deviation and standard deviation. Using these measures we find that the error due to the local approximations decreases with increasing basis set. Next we construct a local virtual space using projected atomic orbitals (PAOs). The accuracy of the incremental series in combination with a distance-based truncation of the PAO space is analysed and compared to the convergence of the incremental series within the domain-specific basis set approach. Furthermore, we establish the recently proposed incremental CCSD(T)|MP2 method as a benchmark method to obtain highly accurate CCSD(T) energies. In combination with a basis set of quintuple-ζ quality we establish benchmarks for the binding energies of the investigated complexes. Finally, we use the inc-CCSD(T)|MP2/aV5Z' binding energies of 45 complexes and 34 dissociation reactions to compute the accuracy of several state of the art density functional theory (DFT) functionals like BP86, B3LYP, CAM-B3LYP, M06, PBE0 and TPSSh. With our implementation of the incremental scheme it was possible to compute the inc-CCSD(T)|MP2/aV5Z' energy for Al(H2O)3+ 25 (6106 AOs).
Santra, Robin; Christ, Kevin V.; Greene, Chris H.
2004-04-01
The first three electronically excited states in the alkaline-earth-metal atoms magnesium, calcium, and strontium comprise the (nsnp){sup 3}P{sub J}{sup o}(J=0,1,2) fine-structure manifold. All three states are metastable and are of interest for optical atomic clocks as well as for cold-collision physics. An efficient technique--based on a physically motivated potential that models the presence of the ionic core--is employed to solve the Schroedinger equation for the two-electron valence shell. In this way, radiative lifetimes, laser-induced clock shifts, and long-range interaction parameters are calculated for metastable Mg, Ca, and Sr.
Structure and properties of electronic and hole centers in CsBr from theoretical calculations
Halliday, Matthew T.; Hess, Wayne P.; Shluger, Alexander L.
2015-06-24
The electronic structure, geometry, diffusion barriers and optical properties of fundamental defects of CsBr are calculated using hybrid functional DFT and TD- DFT methods. The B3LYP functional with a modified exchange contribution has been used in an embedded cluster scheme to model the structure and spectroscopic properties of self-trapped triplet exciton, interstitial Br atoms and ions, self-trapped holes and Br vacancies. The calculated migration barriers and positions of maxima of optical absorption bands are in good agreement with experiment, justifying the obtained defect geometries. The o*-center triplet exciton luminescence energy is also accurately calculated.
Accurate ab initio calculations which demonstrate a 3 Pi u ground state for Al2
NASA Technical Reports Server (NTRS)
Bauschlicher, Charles W., Jr.; Partridge, Harry; Langhoff, Stephen R.; Taylor, Peter R.; Walch, Stephen P.
1986-01-01
The spectroscopic parameters and separations between the three low-lying X 3 Pi u, A 3 Sigma g -, and a 1 Sigma g + states of Al2 are studied as a function of both the one-particle and n-particle basis set. Approximate correlation treatments are calibrated against full Cl calculations correlating the six valence electrons in a double-zeta plus two d-function basis set. Since the CASSCF/MRCI 3 Pi u to 3 Sigma g - separation is in excellent agreement wtih the FCI value, the MRCI calculations were carried out in an extended (20s13p6d4f)/(6s5p3d2f) gaussian basis. Including a small correction for relativistic effects, the best estimate is that 3 Sigma g - state lies 174/cm above the 3 Pi u ground state. The 1 Sigma g + state lies at least 2000/cm higher in energy. At the CPF level, inclusion of 2s and 2p correlation has little effect on D sub e, reduces T sub e by only 26/cm, and shortens the bond lengths by about 0.02 a sub o. Further strong support for a 3 Pi u ground state comes from the experimental absorption spectra, since both observed transitions can be convincingly assigned as 3 Pi u yields 3 Pi g. The (2) 3 Pi g state is observed to be sensitive to the level of correlation treatment, and to have its minimum shifted to shorter rho values, such that the strongest experimental absorption peak probably corresponds to the 0 yields 2 transition.
Technology Transfer Automated Retrieval System (TEKTRAN)
Time-dependent density functional theory (TDDFT) was employed for theoretical calculation of electronic circular dichroism (ECD) of a hexahydroxydiphenoyl (HHDP)-containing flavanone glycoside, mattucinol-7-O-[4'',6''-O-(aS)-hexahydroxydiphenoyl]-ß-d-glucopyranoside (2). It identified the roles of t...
Microcomputer Calculation of Theoretical Pre-Exponential Factors for Bimolecular Reactions.
ERIC Educational Resources Information Center
Venugopalan, Mundiyath
1991-01-01
Described is the application of microcomputers to predict reaction rates based on theoretical atomic and molecular properties taught in undergraduate physical chemistry. Listed is the BASIC program which computes the partition functions for any specific bimolecular reactants. These functions are then used to calculate the pre-exponential factor of…
Invariant-theoretic method for calculating Clebsch-Gordan coefficients for space groups
Aizenberg, A.Ya.; Gufan, Yu.M.
1995-03-01
A new invariant-theoretic method to directly calculate Clebsch-Gordan coefficients for space and point groups representations is proposed. The method is exemplified with the space groups O{sub h}{sup 5} and D{sub 6h}{sup 1}. 34 refs.
Efficiency and power loss in d. c. chopper circuits. [Theoretical calculation
Beck, M.O.
1981-01-01
The object of this paper was to investigate the efficiency and source of power losses of various classes of dc chopper circuits. The study involved a theoretical calculation of the power losses, supported by a considerable amount of practical work on full power-rated traction motor test bed. 3 refs.
Wang, Xinxin; Shi, Deheng; Sun, Jinfeng; Zhu, Zunlue
2016-08-01
The potential energy curves were calculated for the 21 states (X(2)Π, A(2)Π, 3(2)Π, 4(2)Π, 5(2)Π, 1(2)Σ(+), 2(2)Σ(+), 3(2)Σ(+), 1(2)Σ(-), 2(2)Σ(-), 3(2)Σ(-), 1(2)Δ, 2(2)Δ, 3(2)Δ, 1(2)Φ, 1(4)Σ(+), a(4)Σ(-), 2(4)Σ(-), 1(4)Π, 2(4)Π and 1(4)Δ), which originated from the two lowest dissociation channels of ClO radical. The calculations were done for internuclear separations approximately from 0.08 to 1.10nm using the CASSCF method, which was followed by the icMRCI approach with the aug-cc-pV5Z basis set. Of these 21 states, the 1(4)Π, 2(4)Π, 3(2)Δ, 4(2)Π, 5(2)Π, 1(2)Φ, 3(2)Σ(+), 1(4)Δ and 2(4)Σ(-) states are repulsive. The 1(2)Δ, 1(2)Σ(-), 1(4)Σ(+), 2(2)Σ(-), 1(2)Σ(+), 2(2)Σ(+), 2(2)Δ and 3(2)Σ(-) states are very weakly bound. Only the A(2)Π state has one barrier. The avoided crossing exists between the A(2)Π and the 3(2)Π state. However, the avoided crossing does not generate any double wells. Core- valence correlation correction was accounted for at the level of an aug-cc-pCVQZ basis set. Scalar relativistic correction was included by the third-order Douglas-Kroll Hamiltonian approximation at the level of an aug-cc-pVQZ basis set. All the potential energy curves were extrapolated to the complete basis set limit. The spectroscopic parameters were determined. The 1(2)Σ(-), 2(2)Σ(-), 3(2)Σ(-) and 1(4)Σ(+) states may be very difficult to be detected in an experiment, since each of these Λ-S states has only one or two vibrational states. The Franck-Condon factors and radiative lifetimes were calculated for several low vibrational levels of the A(2)Π - X(2)Π, 3(2)Π - a(4)Σ(-), 2(2)Δ - a(4)Σ(-) and 3(2)Σ(-) - 1(2)Σ(-) transitions. The spin-orbit coupling effect on the spectroscopic parameters of the X(2)Π, A(2)Π, 3(2)Π, a(4)Σ(-) and 2(2)Σ(+) states were discussed. The spectroscopic properties reported here can be expected to be reliably predicted ones. PMID:27111157
Theoretical basis for two-phase disengagement and vapor fraction calculations
Sheppard, C.M.
1995-12-31
This paper reviews the theoretical basis for disengagement calculations. That is the assumptions of pseudo-steady-state, volumetric gas production and applicable drift-flux correlations. Using these an axial void fraction profile can be calculated. The use of graphical, analytical, and simplified-equation-based correlations for vertical cylinders, horizontal cylinders, and spheres are covered. Before disengagement occurs the vent stream void fraction can be estimated using the coupling equation. The implications of the converging cross-sectional area of horizontal cylinders and spheres are discussed for both disengagement and two-phase flow calculations. 22 refs., 9 figs., 5 tabs.
NASA Astrophysics Data System (ADS)
Roberts, B. M.; Dzuba, V. A.; Flambaum, V. V.; Pospelov, M.; Stadnik, Y. V.
2016-06-01
We revisit the WIMP-type dark matter scattering on electrons that results in atomic ionization and can manifest itself in a variety of existing direct-detection experiments. Unlike the WIMP-nucleon scattering, where current experiments probe typical interaction strengths much smaller than the Fermi constant, the scattering on electrons requires a much stronger interaction to be detectable, which in turn requires new light force carriers. We account for such new forces explicitly, by introducing a mediator particle with scalar or vector couplings to dark matter and to electrons. We then perform state-of-the-art numerical calculations of atomic ionization relevant to the existing experiments. Our goals are to consistently take into account the atomic physics aspect of the problem (e.g., the relativistic effects, which can be quite significant) and to scan the parameter space—the dark matter mass, the mediator mass, and the effective coupling strength—to see if there is any part of the parameter space that could potentially explain the DAMA modulation signal. While we find that the modulation fraction of all events with energy deposition above 2 keV in NaI can be quite significant, reaching ˜50 %, the relevant parts of the parameter space are excluded by the XENON10 and XENON100 experiments.
Bubin, Sergiy; Stanke, Monika; Adamowicz, Ludwik
2011-08-21
In this work we report very accurate variational calculations of the complete pure vibrational spectrum of the D(2) molecule performed within the framework where the Born-Oppenheimer (BO) approximation is not assumed. After the elimination of the center-of-mass motion, D(2) becomes a three-particle problem in this framework. As the considered states correspond to the zero total angular momentum, their wave functions are expanded in terms of all-particle, one-center, spherically symmetric explicitly correlated Gaussian functions multiplied by even non-negative powers of the internuclear distance. The nonrelativistic energies of the states obtained in the non-BO calculations are corrected for the relativistic effects of the order of α(2) (where α = 1/c is the fine structure constant) calculated as expectation values of the operators representing these effects. PMID:21861559
Derenzo, Stephen E.; Klintenberg, Mattias K.; Weber, Marvin J.
2000-02-01
In performing atomic cluster calculations of local electronic structure defects in ionic crystals, the crystal is often modeled as a central cluster of 5-50 ions embedded in an array of point charges. For most crystals, however, a finite three-dimensional repeated array of unit cells generates electrostatic potentials that are in significant disagreement with the Madelung (infinite crystal) potentials computed by the Ewald method. This is illustrated for the cubic crystal CaF{sub 2}. We present a novel algorithm for solving this problem for any crystal whose unit cell information is known: (1) the unit cell is used to generate a neutral array containing typically 10 000 point charges at their normal crystallographic positions; (2) the array is divided into zone 1 (a volume defined by the atomic cluster of interest), zone 2 (several hundred additional point charges that together with zone 1 fill a spherical volume), and zone 3 (all other point charges); (3) the Ewald formula is used to compute the site potentials at all point charges in zones 1 and 2; (4) a system of simultaneous linear equations is solved to find the zone 3 charge values that make the zone 1 and zone 2 site potentials exactly equal to their Ewald values and the total charge and dipole moments equal to zero, and (5) the solution is checked at 1000 additional points randomly chosen in zone 1. The method is applied to 33 different crystal types with 50-71 ions in zone 1. In all cases the accuracy determined in step 5 steadily improves as the sizes of zones 2 and 3 are increased, reaching a typical rms error of 1 {mu}V in zone 1 for 500 point charges in zone 2 and 10 000 in zone 3. (c) 2000 American Institute of Physics.
NASA Astrophysics Data System (ADS)
Majumder, Moumita; Dawes, Richard; Wang, Xiao-Gang; Carrington, Tucker; Li, Jun; Guo, Hua; Manzhos, Sergei
2014-06-01
New potential energy surfaces for methane were constructed, represented as analytic fits to about 100,000 individual high-level ab initio data. Explicitly-correlated multireference data (MRCI-F12(AE)/CVQZ-F12) were computed using Molpro [1] and fit using multiple strategies. Fits with small to negligible errors were obtained using adaptations of the permutation-invariant-polynomials (PIP) approach [2,3] based on neural-networks (PIP-NN) [4,5] and the interpolative moving least squares (IMLS) fitting method [6] (PIP-IMLS). The PESs were used in full-dimensional vibrational calculations with an exact kinetic energy operator by representing the Hamiltonian in a basis of products of contracted bend and stretch functions and using a symmetry adapted Lanczos method to obtain eigenvalues and eigenvectors. Very close agreement with experiment was produced from the purely ab initio PESs. References 1- H.-J. Werner, P. J. Knowles, G. Knizia, 2012.1 ed. 2012, MOLPRO, a package of ab initio programs. see http://www.molpro.net. 2- Z. Xie and J. M. Bowman, J. Chem. Theory Comput 6, 26, 2010. 3- B. J. Braams and J. M. Bowman, Int. Rev. Phys. Chem. 28, 577, 2009. 4- J. Li, B. Jiang and Hua Guo, J. Chem. Phys. 139, 204103 (2013). 5- S Manzhos, X Wang, R Dawes and T Carrington, JPC A 110, 5295 (2006). 6- R. Dawes, X-G Wang, A.W. Jasper and T. Carrington Jr., J. Chem. Phys. 133, 134304 (2010).
NASA Astrophysics Data System (ADS)
Luo, Ye; Sorella, Sandro
2014-03-01
We introduce a general and efficient method for the calculation of vibrational frequencies of electronic systems, ranging from molecules to solids. By performing damped molecular dynamics with ab initio forces, we show that quantum vibrational frequencies can be evaluated by diagonalizing the time averaged position-position or force-force correlation matrices, although the ionic motion is treated on the classical level within the Born-Oppenheimer approximation. The novelty of our approach is to evaluate atomic forces with QMC by means of a highly accurate and correlated variational wave function which is optimized simultaneously during the dynamics. QMC is an accurate and promising many-body technique for electronic structure calculation thanks to massively parallel computers. However, since infinite statistics is not feasible, property evaluation may be affected by large noise that is difficult to harness. Our approach controls the QMC stochastic bias systematically and gives very accurate results with moderate computational effort, namely even with noisy forces. We prove the accuracy and efficiency of our method on the water monomer[A. Zen et al., JCTC 9 (2013) 4332] and dimer. We are currently working on the challenging problem of simulating liquid water at ambient conditions.
Bai, D.; Levine, S.L. ); Luoma, J.; Mahgerefteh, M. )
1992-01-01
The Three Mile Island unit 1 core reloads have been designed using fast but accurate scoping codes, PSUI-LEOPARD and ADMARC. PSUI-LEOPARD has been normalized to EPRI-CPM2 results and used to calculate the two-group constants, whereas ADMARC is a modern two-dimensional, two-group diffusion theory nodal code. Problems in accuracy were encountered for cycles 8 and higher as the core lifetime was increased beyond 500 effective full-power days. This is because the heavier loaded cores in both {sup 235}U and {sup 10}B have harder neutron spectra, which produces a change in the transport effect in the baffle reflector region, and the burnable poison (BP) simulations were not accurate enough for the cores containing the increased amount of {sup 10}B required in the BP rods. In the authors study, a technique has been developed to take into account the change in the transport effect in the baffle region by modifying the fast neutron diffusion coefficient as a function of cycle length and core exposure or burnup. A more accurate BP simulation method is also developed, using integral transport theory and CPM2 data, to calculate the BP contribution to the equivalent fuel assembly (supercell) two-group constants. The net result is that the accuracy of the scoping codes is as good as that produced by CASMO/SIMULATE or CPM2/SIMULATE when comparing with measured data.
Calculations of steady and transient channel flows with a time-accurate L-U factorization scheme
NASA Technical Reports Server (NTRS)
Kim, S.-W.
1991-01-01
Calculations of steady and unsteady, transonic, turbulent channel flows with a time accurate, lower-upper (L-U) factorization scheme are presented. The L-U factorization scheme is formally second-order accurate in time and space, and it is an extension of the steady state flow solver (RPLUS) used extensively to solve compressible flows. A time discretization method and the implementation of a consistent boundary condition specific to the L-U factorization scheme are also presented. The turbulence is described by the Baldwin-Lomax algebraic turbulence model. The present L-U scheme yields stable numerical results with the use of much smaller artificial dissipations than those used in the previous steady flow solver for steady and unsteady channel flows. The capability to solve time dependent flows is shown by solving very weakly excited and strongly excited, forced oscillatory, channel flows.
NASA Astrophysics Data System (ADS)
Santos, Marcela F.; Braga, Carolyne B.; Rozada, Thiago C.; Basso, Ernani A.; Fiorin, Barbara C.
2014-08-01
The geometries involved in the conformational equilibria of 2,2-dichloro-N-cyclohexyl-N-methyl-acetamide (DCCMA) and 2-chloro-N,N-dicyclohexylacetamide (CDCA) were investigated. Theoretical calculations at the B3LYP/cc-pVDZ level of theory showed that gauche forms (Clsbnd Csbnd Cdbnd O) are the most stable and the predominant conformers in isolated phase. Both compounds had the conformational behavior in solvents of different polarities estimated from theoretical calculations with the PCM (Polarizable Continuum Model), at the same level of theory, using infrared data from deconvolution of the carbonyl absorption bands and 13C NMR spectra. Their IR spectra showed two carbonyl absorptions and that the conformer with the highest dipole moment had its population increased when the most polar solvents were used, in accordance with the theoretical calculation in solution. 1JCH coupling constants were obtained from their NMR spectra, and revealed that there was population variation of conformers with solvent exchange. Experimental data (NMR and IR) as well as calculations including the solvent effects followed the same trend.
van Wyk, Marnus J; Bingle, Marianne; Meyer, Frans J C
2005-09-01
International bodies such as International Commission on Non-Ionizing Radiation Protection (ICNIRP) and the Institute for Electrical and Electronic Engineering (IEEE) make provision for human exposure assessment based on SAR calculations (or measurements) and basic restrictions. In the case of base station exposure this is mostly applicable to occupational exposure scenarios in the very near field of these antennas where the conservative reference level criteria could be unnecessarily restrictive. This study presents a variety of critical aspects that need to be considered when calculating SAR in a human body close to a mobile phone base station antenna. A hybrid FEM/MoM technique is proposed as a suitable numerical method to obtain accurate results. The verification of the FEM/MoM implementation has been presented in a previous publication; the focus of this study is an investigation into the detail that must be included in a numerical model of the antenna, to accurately represent the real-world scenario. This is accomplished by comparing numerical results to measurements for a generic GSM base station antenna and appropriate, representative canonical and human phantoms. The results show that it is critical to take the disturbance effect of the human phantom (a large conductive body) on the base station antenna into account when the antenna-phantom spacing is less than 300 mm. For these small spacings, the antenna structure must be modeled in detail. The conclusion is that it is feasible to calculate, using the proposed techniques and methodology, accurate occupational compliance zones around base station antennas based on a SAR profile and basic restriction guidelines. PMID:15931680
NASA Astrophysics Data System (ADS)
Teale, Andrew M.; Lutnæs, Ola B.; Helgaker, Trygve; Tozer, David J.; Gauss, Jürgen
2013-01-01
Accurate sets of benchmark nuclear-magnetic-resonance shielding constants and spin-rotation constants are calculated using coupled-cluster singles-doubles (CCSD) theory and coupled-cluster singles-doubles-perturbative-triples [CCSD(T)] theory, in a variety of basis sets consisting of (rotational) London atomic orbitals. The accuracy of the calculated coupled-cluster constants is established by a careful comparison with experimental data, taking into account zero-point vibrational corrections. Coupled-cluster basis-set convergence is analyzed and extrapolation techniques are employed to estimate basis-set-limit quantities, thereby establishing an accurate benchmark data set. Together with the set provided for rotational g-tensors and magnetizabilities in our previous work [O. B. Lutnæs, A. M. Teale, T. Helgaker, D. J. Tozer, K. Ruud, and J. Gauss, J. Chem. Phys. 131, 144104 (2009)], 10.1063/1.3242081, it provides a substantial source of consistently calculated high-accuracy data on second-order magnetic response properties. The utility of this benchmark data set is demonstrated by examining a wide variety of Kohn-Sham exchange-correlation functionals for the calculation of these properties. None of the existing approximate functionals provide an accuracy competitive with that provided by CCSD or CCSD(T) theory. The need for a careful consideration of vibrational effects is clearly illustrated. Finally, the pure coupled-cluster results are compared with the results of Kohn-Sham calculations constrained to give the same electronic density. Routes to future improvements are discussed in light of this comparison.
NASA Technical Reports Server (NTRS)
Srivastava, R. C.; Coen, J. L.
1992-01-01
The traditional explicit growth equation has been widely used to calculate the growth and evaporation of hydrometeors by the diffusion of water vapor. This paper reexamines the assumptions underlying the traditional equation and shows that large errors (10-30 percent in some cases) result if it is used carelessly. More accurate explicit equations are derived by approximating the saturation vapor-density difference as a quadratic rather than a linear function of the temperature difference between the particle and ambient air. These new equations, which reduce the error to less than a few percent, merit inclusion in a broad range of atmospheric models.
NASA Technical Reports Server (NTRS)
Gamache, Robert R.; Pollack, James B.
1995-01-01
Halfwidths were calculated for H2O with H2 as a broadening gas and were estimated for He as the broadening species. The calculations used the model of Robert and Bonamy with parabolic trajectories and all relevant terms in the interaction potential. The calculations investigated the dependence of the halfwidth on the order of the atom-atom expansion, the rotational states, and the temperature in the range 200 to 400K. Finally, calculations were performed for many transitions of interest in the 5 micrometer window region of the spectrum. The resulting data will be supplied to Dr. R. Freedman for extracting accurate water mixing ratios from the analysis of the thermal channels for the Net Flux experiment on the Galileo probe.
Takahashi, F; Shigemori, Y; Seki, A
2009-01-01
A system has been developed to assess radiation dose distribution inside the body of exposed persons in a radiological accident by utilising radiation transport calculation codes-MCNP and MCNPX. The system consists mainly of two parts, pre-processor and post-processor of the radiation transport calculation. Programs for the pre-processor are used to set up a 'problem-dependent' input file, which defines the accident condition and dosimetric quantities to be estimated. The program developed for the post-processor part can effectively indicate dose information based upon the output file of the code. All of the programs in the dosimetry system can be executed with a generally used personal computer and accurately give the dose profile to an exposed person in a radiological accident without complicated procedures. An experiment using a physical phantom was carried out to verify the availability of the dosimetry system with the developed programs in a gamma ray irradiation field. PMID:19181661
NASA Astrophysics Data System (ADS)
Fedorov, V. M.
2016-05-01
Based on the astronomical ephemerides DE-406, theoretical calculations have been performed of the interannual variability of the Earth's insolation related to celestial-mechanical processes for 365 points of a tropical year in the time period from 1900 to 2050. It has been determined that the average amplitude of variations of the interannual insolation is 0.310 W/m2 (0.023% of the solar constant). The calculated variations are characterized by strict periodicity that corresponds with the length of a synodic month. Connection between the extreme values of the calculated insolation variability and syzygies has been defined. The average amplitude of the calculated variability exceeds by 1.7 times (0.01% of the solar constant) the amplitude of the interannual variability in the 11-year variation of the total Earth's insolation.
Nuclear radii calculations in various theoretical approaches for nucleus-nucleus interactions
Merino, C.; Novikov, I. S.; Shabelski, Yu.
2009-12-15
The information about sizes and nuclear density distributions in unstable (radioactive) nuclei is usually extracted from the data on interaction of radioactive nuclear beams with a nuclear target. We show that in the case of nucleus-nucleus collisions the values of the parameters depend somewhat strongly on the considered theoretical approach and on the assumption about the parametrization of the nuclear density distribution. The obtained values of root-mean-square radii (R{sub rms}) for stable nuclei with atomic weights A=12-40 vary by approximately 0.1 fm when calculated in the optical approximation, in the rigid target approximation, and using the exact expression of the Glauber theory. We present several examples of R{sub rms} radii calculations using these three theoretical approaches and compare these results with the data obtained from electron-nucleus scattering.
NASA Astrophysics Data System (ADS)
Sievänen, Elina; Toušek, Jaromír; Lunerová, Kamila; Marek, Jaromír; Jankovská, Dagmar; Dvorská, Margita; Marek, Radek
2010-08-01
In this article we present a detailed structural investigation for five homoisoflavonoids, molecules important from the pharmacological point of view. For studying the electron distribution as well as its influence on the physicochemical properties, NMR spectroscopy, X-ray diffraction, and theoretical calculations have been used. Nuclear magnetic shieldings obtained by using DFT calculations for optimized molecular geometries are correlated with the experimentally determined chemical shifts. The theoretical data are well in agreement with the experimental values. The single crystal X-ray structures of homoisoflavonoid derivatives 1, 3, and 4 have been solved. The molecular geometries and crystal packing determined by X-ray diffraction are used for characterizing the intermolecular interactions. Electron distribution is crucial for the stability of radicals and hence the antioxidant efficiency of flavonoid structures. The hydrogen bonding governs the formation of complexes of homoisoflavonoids with biological targets.
Theoretical calculation of plane wave speeds for alkali metals under pressure.
NASA Technical Reports Server (NTRS)
Eftis, J.; Macdonald, D. E.; Arkilic, G. M.
1971-01-01
Theoretical calculations of the variation with pressure of small amplitude plane wave speeds are performed for sodium and potassium at zero temperature. The results obtained for wave speeds associated with volume dependent second-order elastic coefficients show better agreement with experimental data than for wave speeds associated with shear dependent coefficients. This result is believed to be due to omission of the band structure correction to the strain energy density.
Tian, Lian; Henningsen, Joseph; Salick, Max R; Crone, Wendy C; Gunderson, McLean; Dailey, Seth H; Chesler, Naomi C
2015-07-01
The mechanical properties of vascular tissues affect hemodynamics and can alter disease progression. The uniaxial tensile test is a simple and effective method for determining the stress-strain relationship in arterial tissue ex vivo. To enable calculation of strain, stretch can be measured directly with image tracking of markers on the tissue or indirectly from the distance between the grips used to hold the specimen. While the imaging technique is generally considered more accurate, it also requires more analysis, and the grip distance method is more widely used. The purpose of this study is to compare the stretch of the testing specimen calculated from the grip distance method to that obtained from the imaging method for canine descending aortas and large proximal pulmonary arteries. Our results showed a significant difference in stretch between the two methods; however, this difference was consistently less than 2%. Therefore, the grip distance method is an accurate approximation of the stretch in large elastic arteries in the uniaxial tensile test. PMID:25881308
Tian, Lian; Henningsen, Joseph; Salick, Max R.; Crone, Wendy C.; Gunderson, McLean; Dailey, Seth H.; Chesler, Naomi C.
2015-01-01
The mechanical properties of vascular tissues affect hemodynamics and can alter disease progression. The uniaxial tensile test is a simple and effective method for determining the stress-strain relationship in arterial tissue ex vivo. To enable calculation of strain, stretch can be measured directly with image tracking of markers on the tissue or indirectly from the distance between the grips used to hold the specimen. While the imaging technique is generally considered more accurate, it also requires more analysis, and the grip distance method is more widely used. The purpose of this study is to compare the stretch of the testing specimen calculated from the grip distance method to that obtained from the imaging method for canine descending aortas and large proximal pulmonary arteries. Our results showed a significant difference in stretch between the two methods; however, this difference was consistently less than 2%. Therefore, the grip distance method is an accurate approximation of the stretch in large elastic arteries in the uniaxial tensile test. PMID:25881308
Wu, Baolin; Guan, Weihua; Pankow, James S
2016-03-01
The objective of this paper is to discuss and develop alternative computational methods to accurately and efficiently calculate significance P-values for the commonly used sequence kernel association test (SKAT) and adaptive sum of SKAT and burden test (SKAT-O) for variant set association. We show that the existing software can lead to either conservative or inflated type I errors. We develop alternative and efficient computational algorithms that quickly compute the SKAT P-value and have well-controlled type I errors. In addition, we derive an alternative and simplified formula for calculating the significance P-value of SKAT-O, which sheds light on the development of efficient and accurate numerical algorithms. We implement the proposed methods in the publicly available R package that can be readily used or adapted to large-scale sequencing studies. Given that more and more large-scale exome and whole genome sequencing or re-sequencing studies are being conducted, the proposed methods are practically very important. We conduct extensive numerical studies to investigate the performance of the proposed methods. We further illustrate their usefulness with application to associations between rare exonic variants and fasting glucose levels in the Atherosclerosis Risk in Communities (ARIC) study. PMID:26757198
NASA Astrophysics Data System (ADS)
Moore, Keith; McLaughlin, Brendan M.; Lane, Ian C.
2016-04-01
BaH (and its isotopomers) is an attractive molecular candidate for laser cooling to ultracold temperatures and a potential precursor for the production of ultracold gases of hydrogen and deuterium. The theoretical challenge is to simulate the laser cooling cycle as reliably as possible and this paper addresses the generation of a highly accurate ab initio 2Σ+ potential for such studies. The performance of various basis sets within the multi-reference configuration-interaction (MRCI) approximation with the Davidson correction is tested and taken to the Complete Basis Set (CBS) limit. It is shown that the calculated molecular constants using a 46 electron effective core-potential and even-tempered augmented polarized core-valence basis sets (aug-pCVnZ-PP, n = 4 and 5) but only including three active electrons in the MRCI calculation are in excellent agreement with the available experimental values. The predicted dissociation energy De for the X2Σ+ state (extrapolated to the CBS limit) is 16 895.12 cm-1 (2.094 eV), which agrees within 0.1% of a revised experimental value of <16 910.6 cm-1, while the calculated re is within 0.03 pm of the experimental result.
NASA Technical Reports Server (NTRS)
Richards, P. G.; Torr, D. G.
1981-01-01
A simplified method for the evaluation of theoretical photoelectron fluxes in the upper atmosphere resulting from the solar radiation at 304 A is presented. The calculation is based on considerations of primary and cascade (secondary) photoelectron production in the two-stream model, where photoelectron transport is described by two electron streams, one moving up and one moving down, and of loss rates due to collisions with neutral gases and thermal electrons. The calculation is illustrated for the case of photoelectrons at an energy of 24.5 eV, and it is noted that the 24.5-eV photoelectron flux may be used to monitor variations in the solar 304 A flux. Theoretical calculations based on various ionization and excitation cross sections of Banks et al. (1974) are shown to be in generally good agreement with AE-E measurements taken between 200 and 235 km, however the use of more recent, larger cross sections leads to photoelectron values a factor of two smaller than observations but in agreement with previous calculations. It is concluded that a final resolution of the photoelectron problem may depend on a reevaluation of the inelastic electron collision cross sections.
NASA Astrophysics Data System (ADS)
He, C.; Liou, K. N.; Takano, Y.; Li, Q.; Yang, P.; Zhang, R.
2014-12-01
The optical properties of black carbon (BC) are significantly affected by its aging process in the atmosphere. We have built a conceptual model defining three BC aging stages, including freshly emitted BC aggregates, coating by soluble material and hygroscopic growth. We apply an improved geometric-optics surface-wave approach (Liou et al., 2011; Takano et al., 2013) to calculate the absorption and scattering properties of BC at each stage and compare the theoretical results with those obtained from laboratory experiments (Zhang et al., 2008; Khalizov et al., 2009). Preliminary results show a general agreement between calculated and measured BC absorption cross sections (bias < 10%) and scattering cross sections (bias < 30%) for BC aerosols with mobility diameters of 155, 245 and 320 nm at Stages 1 and 2, where BC is coated by sulfuric acid and its water solution, respectively. We find that the calculated scattering and absorption cross sections for fresh BC aggregates (Stage 0) with different sizes are invariably larger than experimental results partly because of the uncertainty in theoretical calculations for BC with size parameters less than 1. It appears that the uncertainty in the experiment could also contribute to the discrepancy, considering that the measuring instrument missed some scattering in certain angles (0-7° and 170-180°). Finally, we will apply the conceptual model and the single-scattering results to assess the effects of BC aging processes on direct radiative forcing using observed BC vertical profiles.
Lundahl, Gunnel
2007-01-01
When calculating of the physical F121.1 degrees c-value by the equation F121.1 degrees C = t x 10(T-121.1/z the temperature (T), in combination with the z-value, influences the F121.1 degrees c-value exponentially. Because the z-value for spores of Geobacillus stearothermophilus often varies between 6 and 9, the biological F-value (F(Bio) will not always correspond to the F0-value based on temperature records from the sterilization process calculated with a z-value of 10, even if the calibration of both of them are correct. Consequently an error in calibration of thermocouples and difference in z-values influences the F121.1 degrees c-values logarithmically. The paper describes how results from measurements with different z-values can be compared. The first part describes the mathematics of a calculation program, which makes it easily possible to compare F0-values based on temperature records with the F(BIO)-value based on analysis of bioindicators such as glycerin-water-suspension sensors. For biological measurements, a suitable bioindicator with a high D121-value can be used (such a bioindicator can be manufactured as described in the article "A Method of Increasing Test Range and Accuracy of Bioindicators-Geobacillus stearothermophilus Spores"). By the mathematics and calculations described in this macro program it is possible to calculate for every position the theoretical temperature difference (deltaT(th)) needed to explain the difference in results between the thermocouple and the biointegrator. Since the temperature difference is a linear function and constant all over the process this value is an indication of the magnitude of an error. A graph and table from these calculations gives a picture of the run. The second part deals with product characteristics, the sterilization processes, loading patterns. Appropriate safety margins have to be chosen in the development phase of a sterilization process to achieve acceptable safety limits. Case studies are
NASA Astrophysics Data System (ADS)
McCann, Kathleen; Laane, Jaan
2008-11-01
The Raman and infrared spectra of dipicolinic acid (DPA) and dinicotinic acid (DNic) and their salts (CaDPA, Na 2DPA, and CaDNic) have been recorded and the spectra have been assigned. Ab initio and DFT calculations were carried out to predict the structures and vibrational spectra and were compared to the experimental results. Because of extensive intermolecular hydrogen bonding in the crystals of these molecules, the calculated structures and spectra for the individual molecules agree only moderately well with the experimental values. Theoretical calculations were also carried out for DPA dimers and DPA·2H 2O to better understand the intermolecular interactions. The spectra do show that DPA and its calcium salt, which are present in anthrax spores, can be distinguished from the very similar DNic and CaDNic.
NASA Technical Reports Server (NTRS)
Binienda, W. K.; Arnold, S. M.; Tan, H. Q.
1992-01-01
An essential part of describing the damage state and predicting the damage growth in a multicracked plate is the accurate calculation of stress intensity factors (SIF's). Here, a methodology and rigorous solution formulation for SIF's of a multicracked plate, with fully interacting cracks, subjected to a far-field arbitrary stress state is presented. The fundamental perturbation problem is derived, and the steps needed to formulate the system of singular integral equations whose solution gives rise to the evaluation of the SIF's are identified. This analytical derivation and numerical solution are obtained by using intelligent application of symbolic computations and automatic FORTRAN generation capabilities (described in the second part of this paper). As a result, a symbolic/FORTRAN package, named SYMFRAC, that is capable of providing accurate SIF's at each crack tip was developed and validated.
Esque, Jeremy; Cecchini, Marco
2015-04-23
The calculation of the free energy of conformation is key to understanding the function of biomolecules and has attracted significant interest in recent years. Here, we present an improvement of the confinement method that was designed for use in the context of explicit solvent MD simulations. The development involves an additional step in which the solvation free energy of the harmonically restrained conformers is accurately determined by multistage free energy perturbation simulations. As a test-case application, the newly introduced confinement/solvation free energy (CSF) approach was used to compute differences in free energy between conformers of the alanine dipeptide in explicit water. The results are in excellent agreement with reference calculations based on both converged molecular dynamics and umbrella sampling. To illustrate the general applicability of the method, conformational equilibria of met-enkephalin (5 aa) and deca-alanine (10 aa) in solution were also analyzed. In both cases, smoothly converged free-energy results were obtained in agreement with equilibrium sampling or literature calculations. These results demonstrate that the CSF method may provide conformational free-energy differences of biomolecules with small statistical errors (below 0.5 kcal/mol) and at a moderate computational cost even with a full representation of the solvent. PMID:25807150
Sun, Y. Y.; Kim, Y. H.; Lee, K.; Zhang, S. B.
2008-01-01
Density functional theory (DFT) in the commonly used local density or generalized gradient approximation fails to describe van der Waals (vdW) interactions that are vital to organic, biological, and other molecular systems. Here, we propose a simple, efficient, yet accurate local atomic potential (LAP) approach, named DFT+LAP, for including vdW interactions in the framework of DFT. The LAPs for H, C, N, and O are generated by fitting the DFT+LAP potential energy curves of small molecule dimers to those obtained from coupled cluster calculations with single, double, and perturbatively treated triple excitations, CCSD(T). Excellent transferability of the LAPs is demonstrated by remarkable agreement with the JSCH-2005 benchmark database [P. Jurecka et al. Phys. Chem. Chem. Phys. 8, 1985 (2006)], which provides the interaction energies of CCSD(T) quality for 165 vdW and hydrogen-bonded complexes. For over 100 vdW dominant complexes in this database, our DFT+LAP calculations give a mean absolute deviation from the benchmark results less than 0.5 kcal/mol. The DFT+LAP approach involves no extra computational cost other than standard DFT calculations and no modification of existing DFT codes, which enables straightforward quantum simulations, such as ab initio molecular dynamics, on biomolecular systems, as well as on other organic systems.
NASA Astrophysics Data System (ADS)
Shi, De-Heng; Liu, Yu-Fang; Sun, Jin-Feng; Yang, Xiang-Dong; Zhu, Zun-Lue
2006-05-01
The reasonable dissociation limit of the A1Σ+ state 7LiH molecule is obtained. The accurate dissociation energy and the equilibrium geometry of this state are calculated using a symmetry-adapted-cluster configuration-interaction method in complete active space for the first time. The whole potential energy curve and the dipole moment function for the A1Σ+ state are calculated over a wide internuclear separation range from about 0.1 to 1.4 nm. The calculated equilibrium geometry and dissociation energy of this potential energy curve are of Re=0.2487 nm and De=1.064 eV, respectively. The unusual negative values of the anharmonicity constant and the vibration-rotational coupling constant are of ωeχe=-4.7158cm-1 and αe=-0.08649cm-1, respectively. The vertical excitation energy from the ground to the A1Σ+ state is calculated and the value is of 3.613 eV at 0.15875 nm (the equilibrium position of the ground state). The highly anomalous shape of this potential energy curve, which is exceptionally flat over a wide radial range around the equilibrium position, is discussed in detail. The harmonic frequency value of 502.47cm-1 about this state is approximately estimated. Careful comparison of the theoretical determinations with those obtained by previous theories about the A1Σ+ state dissociation energy clearly shows that the present calculations are much closer to the experiments than previous theories, thus represents an improvement.
Mayhall, Nicholas J; Raghavachari, Krishnan
2011-05-10
We present a new extrapolated fragment-based approach, termed molecules-in-molecules (MIM), for accurate energy calculations on large molecules. In this method, we use a multilevel partitioning approach coupled with electronic structure studies at multiple levels of theory to provide a hierarchical strategy for systematically improving the computed results. In particular, we use a generalized hybrid energy expression, similar in spirit to that in the popular ONIOM methodology, that can be combined easily with any fragmentation procedure. In the current work, we explore a MIM scheme which first partitions a molecule into nonoverlapping fragments and then recombines the interacting fragments to form overlapping subsystems. By including all interactions with a cheaper level of theory, the MIM approach is shown to significantly reduce the errors arising from a single level fragmentation procedure. We report the implementation of energies and gradients and the initial assessment of the MIM method using both biological and materials systems as test cases. PMID:26610128
Son, Sang-Kil
2011-03-01
We introduce a new numerical grid-based method on unstructured grids in the three-dimensional real-space to investigate the electronic structure of polyatomic molecules. The Voronoi-cell finite difference (VFD) method realizes a discrete Laplacian operator based on Voronoi cells and their natural neighbors, featuring high adaptivity and simplicity. To resolve multicenter Coulomb singularity in all-electron calculations of polyatomic molecules, this method utilizes highly adaptive molecular grids which consist of spherical atomic grids. It provides accurate and efficient solutions for the Schroedinger equation and the Poisson equation with the all-electron Coulomb potentials regardless of the coordinate system and the molecular symmetry. For numerical examples, we assess accuracy of the VFD method for electronic structures of one-electron polyatomic systems, and apply the method to the density-functional theory for many-electron polyatomic molecules.
Rodgers, R P; Tannen, R
1983-01-01
Since its appearance in 1960, the method of Barnhard and associates for the determination of total lung capacity (TLC) from routine chest radiograms has been widely studied in normal and diseased subjects. The method appears to be as accurate as the current definitive procedure, total body plethysmography. The method is in routine use in major medical institutions where the procedure has been automated, but the method does not seem to have gained the wide use it deserves. This is likely due to the tedium of the technique when performed manually--a single determination can require 30 min. We present here an implementation of the Barnhard method for the HP41-C hand-held programmable calculator. In conjunction with the use of a transparent reticle used for obtaining the required measurements, the program allows a single measurement to be made in under 12 minutes. We hope this technique will make radiographic TLC measurements more broadly accessible to the medical profession. PMID:6872526
Han, Min Cheol; Yeom, Yeon Soo; Kim, Chan Hyeong; Kim, Seonghoon; Sohn, Jason W
2015-02-21
In the present study, to achieve accurate 4D Monte Carlo dose calculation in radiation therapy, we devised a new approach that combines (1) modeling of the patient body using tetrahedral-mesh geometry based on the patient's 4D CT data, (2) continuous movement/deformation of the tetrahedral patient model by interpolation of deformation vector fields acquired through deformable image registration, and (3) direct transportation of radiation particles during the movement and deformation of the tetrahedral patient model. The results of our feasibility study show that it is certainly possible to construct 4D patient models (= phantoms) with sufficient accuracy using the tetrahedral-mesh geometry and to directly transport radiation particles during continuous movement and deformation of the tetrahedral patient model. This new approach not only produces more accurate dose distribution in the patient but also replaces the current practice of using multiple 3D voxel phantoms and combining multiple dose distributions after Monte Carlo simulations. For routine clinical application of our new approach, the use of fast automatic segmentation algorithms is a must. In order to achieve, simultaneously, both dose accuracy and computation speed, the number of tetrahedrons for the lungs should be optimized. Although the current computation speed of our new 4D Monte Carlo simulation approach is slow (i.e. ~40 times slower than that of the conventional dose accumulation approach), this problem is resolvable by developing, in Geant4, a dedicated navigation class optimized for particle transportation in tetrahedral-mesh geometry. PMID:25615567
NASA Astrophysics Data System (ADS)
Lonardoni, D.; Pederiva, F.; Gandolfi, S.
2014-01-01
Background: An accurate assessment of the hyperon-nucleon interaction is of great interest in view of recent observations of very massive neutron stars. The challenge is to build a realistic interaction that can be used over a wide range of masses and in infinite matter starting from the available experimental data on the binding energy of light hypernuclei. To this end, accurate calculations of the hyperon binding energy in a hypernucleus are necessary. Purpose: We present a quantum Monte Carlo study of Λ and ΛΛ hypernuclei up to A =91. We investigate the contribution of two- and three-body Λ-nucleon forces to the Λ binding energy. Method: Ground state energies are computed solving the Schrödinger equation for nonrelativistic baryons by means of the auxiliary field diffusion Monte Carlo algorithm extended to the hypernuclear sector. Results: We show that a simple adjustment of the parameters of the ΛNN three-body force yields a very good agreement with available experimental data over a wide range of hypernuclear masses. In some cases no experiments have been performed yet, and we give new predictions. Conclusions: The newly fitted ΛNN force properly describes the physics of medium-heavy Λ hypernuclei, correctly reproducing the saturation property of the hyperon separation energy.
Theoretical calculation of the melting curve of Cu-Zr binary alloys
Gunawardana, K. G.S.H.; Wilson, S. R.; Mendelev, M. I.; Song, Xueyu
2014-11-14
Helmholtz free energies of the dominant binary crystalline solids found in the Cu-Zr system at high temperatures close to the melting curve are calculated. This theoretical approach combines fundamental measure density functional theory (applied to the hard-sphere reference system) and a perturbative approach to include the attractive interactions. The studied crystalline solids are Cu(fcc), Cu51Zr14(β), CuZr(B2), CuZr2(C11b), Zr(hcp), and Zr(bcc). The calculated Helmholtz free energies of crystalline solids are in good agreement with results from molecular-dynamics (MD) simulations. Using the same perturbation approach, the liquid phase free energies are calculated as a function of composition and temperature, from which themore » melting curve of the entire composition range of this system can be obtained. Phase diagrams are determined in this way for two leading embedded atom method potentials, and the results are compared with experimental data. Furthermore, theoretical melting temperatures are compared both with experimental values and with values obtained directly from MD simulations at several compositions.« less
Theoretical calculation of the melting curve of Cu-Zr binary alloys
Gunawardana, K. G.S.H.; Wilson, S. R.; Mendelev, M. I.; Song, Xueyu
2014-11-14
Helmholtz free energies of the dominant binary crystalline solids found in the Cu-Zr system at high temperatures close to the melting curve are calculated. This theoretical approach combines fundamental measure density functional theory (applied to the hard-sphere reference system) and a perturbative approach to include the attractive interactions. The studied crystalline solids are Cu(fcc), Cu_{51}Zr_{14}(β), CuZr(B_{2}), CuZr_{2}(C11b), Zr(hcp), and Zr(bcc). The calculated Helmholtz free energies of crystalline solids are in good agreement with results from molecular-dynamics (MD) simulations. Using the same perturbation approach, the liquid phase free energies are calculated as a function of composition and temperature, from which the melting curve of the entire composition range of this system can be obtained. Phase diagrams are determined in this way for two leading embedded atom method potentials, and the results are compared with experimental data. Furthermore, theoretical melting temperatures are compared both with experimental values and with values obtained directly from MD simulations at several compositions.
Theoretical study of the electronic structure with dipole moment calculations of barium monofluoride
NASA Astrophysics Data System (ADS)
Tohme, Samir N.; Korek, Mahmoud
2015-12-01
The potential energy curves have been investigated for the 41 lowest doublet and quartet electronic states in the 2s+1Λ± representation below 55,000 cm-1 of the molecule BaF via CASSCF and MRCI (single and double excitations with Davidson correction) calculations. Twenty-five electronic states have been studied here theoretically for the first time. The crossing and avoided crossing of 20 doublet electronic states have been studied in the region 30,000-50,000 cm-1. The harmonic frequency ωe, the internuclear distance Re, the rotational constant Be, the electronic energy with respect to the ground state Te, and the permanent and transition dipole moments have been calculated in addition to static dipole polarizability of the ground state. By using the canonical functions approach, the eigenvalue Ev, the rotational constant Bv, and the abscissas of the turning points Rmin and Rmax have been calculated for the electronic states up to the vibrational level v=98. The comparison of these values with the theoretical results available in the literature shows a very good agreement.
Hoo, Zhe Hui; Curley, Rachael; Campbell, Michael J; Walters, Stephen J; Hind, Daniel; Wildman, Martin J
2016-01-01
Background Preventative inhaled treatments in cystic fibrosis will only be effective in maintaining lung health if used appropriately. An accurate adherence index should therefore reflect treatment effectiveness, but the standard method of reporting adherence, that is, as a percentage of the agreed regimen between clinicians and people with cystic fibrosis, does not account for the appropriateness of the treatment regimen. We describe two different indices of inhaled therapy adherence for adults with cystic fibrosis which take into account effectiveness, that is, “simple” and “sophisticated” normative adherence. Methods to calculate normative adherence Denominator adjustment involves fixing a minimum appropriate value based on the recommended therapy given a person’s characteristics. For simple normative adherence, the denominator is determined by the person’s Pseudomonas status. For sophisticated normative adherence, the denominator is determined by the person’s Pseudomonas status and history of pulmonary exacerbations over the previous year. Numerator adjustment involves capping the daily maximum inhaled therapy use at 100% so that medication overuse does not artificially inflate the adherence level. Three illustrative cases Case A is an example of inhaled therapy under prescription based on Pseudomonas status resulting in lower simple normative adherence compared to unadjusted adherence. Case B is an example of inhaled therapy under-prescription based on previous exacerbation history resulting in lower sophisticated normative adherence compared to unadjusted adherence and simple normative adherence. Case C is an example of nebulizer overuse exaggerating the magnitude of unadjusted adherence. Conclusion Different methods of reporting adherence can result in different magnitudes of adherence. We have proposed two methods of standardizing the calculation of adherence which should better reflect treatment effectiveness. The value of these indices can
Bulut, N; Castillo, J F; Jambrina, P G; Kłos, J; Roncero, O; Aoiz, F J; Bañares, L
2015-12-17
Accurate quantum reactive scattering time-dependent wave packet close-coupling calculations have been carried out to determine total reaction probabilities and integral cross sections for the O(+) + H2 → OH(+) + H reaction in a range of collision energies from 10(-3) eV up to 1.0 eV for the H2 rovibrational states (v = 0; j = 0, 1, 2) and (v = 1; j = 0) using the potential energy surface (PES) by Martı́nez et al. As expected for a barrierless reaction, the reaction cross section decays rapidly with collision energy, Ec, following a behavior that nearly corresponds to that predicted by the Langevin model. Rotational excitation of H2 into j = 1, 2 has a very moderate effect on reactivity, similarly to what happens with vibrational excitation below Ec ≈ 0.3 eV. However, at higher collision energies the cross section increases notably when H2 is promoted to v = 1. This effect is explained by resorting to the effective potentials in the entrance channel. The integral cross sections have been used to calculate rate constants in the temperature range 200-1000 K. A good overall agreement has been found with the available experimental data on integral cross sections and rate constants. In addition, time-independent quantum mechanical and quasi-classical trajectory (QCT) calculations have been performed on the same PES aimed to compare the various methodologies and to discern the detailed mechanism of the title reaction. In particular, the analysis of individual trajectories has made it possible to explain, in terms of the coupling between reagent relative velocity and the topography of the PES, the presence of a series of alternating maxima and minima in the collision energy dependence of the QCT reaction probabilities for the reactions with H2(v=0,1,j=0), which are absent in the quantum mechanical calculations. PMID:25822338
Matsuzawa, Nobuyuki; Ishitani, Akihiko; Dixon, David A.; Uda, Tsuyoshi
2001-06-13
In order to aid in the design of transparent materials for use as photoresists for F2 lithography (157 nm), we have performed time-dependent density functional theory (TD-DFT) calculations of the photoabsorption of molecules in the vacuum ultraviolet region. The application of this TD-DFT method to the prediction of photoabsorption was benchmarked using model molecules such as formaldehyde, and an empirical equation for correcting the calculated transition energy was obtained. The TD-DFT method with the empirical correction equation provides dramatically more accurate results than those obtained with the CIS (single-excitation configuration interaction) method, which we employed in previous studies. We used it to predict the photoabsorption of various molecules such as methanol, t-butylalcohol, acetic acid, methyl acetate, cycloalkane, norbonane, tricyclodecane, tetrahydropyrane, adamantane, maleic anhydride and their fluorinated derivatives.
Experimental investigation and theoretical calculation for 3He induced nuclear reactions on vanadium
NASA Astrophysics Data System (ADS)
Ali, B. M.; Al-Abyad, M.; Seddik, U.; El-Kameesy, S. U.; Ditrói, F.; Takács, S.; Tárkányi, F.
2016-04-01
Using stacked-foil activation technique and gamma-ray spectrometry, excitation functions for 3He induced nuclear reactions on natV were measured. Cross-sections for natV(3He, xn)52m,gMn and natV(3He, pxn)51Cr nuclear reactions were measured up to 27 MeV utilizing the MGC-20E cyclotron of ATOMKI. The measurements establish for the first time consistent excitation curves. Comparisons with results for values derived from different theoretical codes were included. Integral yield were calculated.
Pavanello, Michele; Van Voorhis, Troy; Visscher, Lucas; Neugebauer, Johannes
2013-02-07
Quantum-mechanical methods that are both computationally fast and accurate are not yet available for electronic excitations having charge transfer character. In this work, we present a significant step forward towards this goal for those charge transfer excitations that take place between non-covalently bound molecules. In particular, we present a method that scales linearly with the number of non-covalently bound molecules in the system and is based on a two-pronged approach: The molecular electronic structure of broken-symmetry charge-localized states is obtained with the frozen density embedding formulation of subsystem density-functional theory; subsequently, in a post-SCF calculation, the full-electron Hamiltonian and overlap matrix elements among the charge-localized states are evaluated with an algorithm which takes full advantage of the subsystem DFT density partitioning technique. The method is benchmarked against coupled-cluster calculations and achieves chemical accuracy for the systems considered for intermolecular separations ranging from hydrogen-bond distances to tens of Angstroms. Numerical examples are provided for molecular clusters comprised of up to 56 non-covalently bound molecules.
NASA Astrophysics Data System (ADS)
Mohr, Stephan; Genovese, Luigi; Ratcliff, Laura; Masella, Michel
The quantum mechanics/molecular mechanis (QM/MM) method is a popular approach that allows to perform atomistic simulations using different levels of accuracy. Since only the essential part of the simulation domain is treated using a highly precise (but also expensive) QM method, whereas the remaining parts are handled using a less accurate level of theory, this approach allows to considerably extend the total system size that can be simulated without a notable loss of accuracy. In order to couple the QM and MM regions we use an approximation of the electrostatic potential based on a multipole expansion. The multipoles of the QM region are determined based on the results of a linear scaling Density Functional Theory (DFT) calculation using a set of adaptive, localized basis functions, as implemented within the BigDFT software package. As this determination comes at virtually no extra cost compared to the QM calculation, the coupling between QM and MM region can be done very efficiently. In this presentation I will demonstrate the accuracy of both the linear scaling DFT approach itself as well as of the approximation of the electrostatic potential based on the multipole expansion, and show some first QM/MM applications using the aforementioned approach.
Liu, Hui; Shi, Deheng; Sun, Jinfeng; Zhu, Zunlue
2016-11-01
The potential energy curves were calculated for the 24 Λ-S states correlating with the lowest four dissociation channels of the BO(+) cation. The potential energy curves were also computed for the 60 Ω states generated from the 24 Λ-S states. Calculations were made for internuclear separations from 0.08 to 1.05nm using the CASSCF method, which was followed by the icMRCI approach with the correlation-consistent basis sets. Core-valence correlation, scalar relativistic and basis extrapolation were accounted for. Of the 24 Λ-S states, only three states (2(5)Π, 1(5)Σ(-), and 2(5)Σ(-)) were found to be repulsive; only the 1(5)Δ state was found to be a very weakly-bound state; and the E(1)Π, 2(3)Π, and 1(5)Π states were found to be very strong bound. In addition, the B(1)Σ(+) and 3(1)Σ(+) states have double wells by the avoided crossing between the two states. The a(3)Π, 1(3)Σ(-), and 2(3)Σ(-) states are inverted with the spin-orbit coupling effect included. The spectroscopic parameters were determined and the vibrational properties of several Λ-S states were predicted. Comparison with available experimental data shows that the methodology employed is highly accurate for this system. PMID:27289351
NASA Astrophysics Data System (ADS)
Michalska, Katarzyna; Mizera, Mikołaj; Lewandowska, Kornelia; Cielecka-Piontek, Judyta
2016-07-01
Tedizolid is the newest antibacterial agent from the oxazolidinone class. For its identification, FT-IR (2000-400 cm-1) and Raman (2000-400 cm-1) analyses were proposed. Studies of the enantiomeric purity of tedizolid were conducted based on ultraviolet-circular dichroism (UV-CD) analysis. Density functional theory (DFT) with the B3LYP hybrid functional and 6-311G(2df,2pd) basis set was used for support of the analysis of the FT-IR and Raman spectra. Theoretical methods made it possible to conduct HOMO and LUMO analysis, which was used to determine the charge transfer for two tedizolid enantiomers. Molecular electrostatic potential maps were calculated with the DFT method for both tedizolid enantiomers. The relationship between the results of ab initio calculations and knowledge about the chemical-biological properties of R- and S-tedizolid enantiomers is also discussed.
Calculation of membrane bending rigidity using field-theoretic umbrella sampling.
Smirnova, Y G; Müller, M
2015-12-28
The free-energy change of membrane shape transformations can be small, e.g., as in the case of membrane bending. Therefore, the calculation of the free-energy difference between different membrane morphologies is a challenge. Here, we discuss a computational method - field-theoretic umbrella sampling - to compute the local chemical potential of a non-equilibrium configuration and illustrate how one can apply it to study free-energy changes of membrane transformations using simulations. Specifically, the chemical potential profile of the bent membrane and the bending rigidity of membrane are calculated for a soft, coarse-grained amphiphile model and the MARTINI model of a dioleoylphosphatidylcholine (DOPC) membrane. PMID:26723640
Synthesis, characterization, optical properties and theoretical calculations of 6-fluoro coumarin
NASA Astrophysics Data System (ADS)
Bai, Yihui; Du, Jinyan; Weng, Xuexiang
6-Fluoro coumarin is synthesized and characterized by 1H NMR and 13C NMR. The optical properties of the title compound are investigated by UV-vis absorption and fluorescence emission spectra, the results show the title compound can absorb UV-vis light at 319, 269 and 215 nm, moreover it exhibits blue-purple fluorescence emission at 416 nm. Theoretical studies on molecular structure, infrared spectra (IR), nuclear magnetic resonance (1H NMR, 13C NMR) chemical shifts, UV-vis absorption and fluorescence emission of the synthesized compound have been worked out. Most chemical calculations were performed by density functional theory (DFT) method at the B3LYP/6-311G(d,p) level (NMR at B3LYP/Aug-CC-Pvdz level) using Gaussian 09 program. The compared results reveal that the scaled theoretical vibrational frequencies are in good accordance with the observed spectra; computational chemical shifts are consistent with the experimental values in most parts, except for some minor deviations; the UV-vis absorption calculated matches the experimental one very well, and the fluorescence emission spectrum is in good agreement with the experimental one when the solute-solvent hydrogen-bonding interaction is considered. These good coincidences prove that the computational methods selected can be used to predict these properties of other similar materials where it is difficult to arrive at experimental results.
Investigation on the Gas-Phase Decomposition of Trichlorfon by GC-MS and Theoretical Calculation
Jiang, Kezhi; Zhang, Ningwen; Zhang, Hu; Wang, Jianmei; Qian, Mingrong
2015-01-01
The gas phase pyrolysis of trichlorfon was investigated by the on-line gas chromatography – mass spectrometry (GC-MS) pyrolysis and theoretical calculations. Two reaction channels were proposed in the pyrolytic reaction, by analyzing the detected pyrolytic products in the total ion chromatography, including 2,2,2-trichloroacetaldehyde, dimethyl phosphite, and dichlorvos. Theoretical calculations showed that there is an intramolecular hydrogen bond between the hydroxyl group and the phosphate O atom in trichlorfon, through which the hydroxyl H atom can be easily transferred to phosphate O atom to trigger two pyrolytic channels. In path-a, migration of H atom results in direct decomposition of trichlorfon to give 2,2,2-trichloroacetaldehyde and dimethyl phosphite in one step. In path-b, migration of H atom in trichlorfon is combined with formation of the O-P bond to give an intermediate, followed by HCl elimination to afford dichlorvos. Path-a is kinetically more favorable than path-b, which is consistent with the GC-MS results. PMID:25856549
NASA Technical Reports Server (NTRS)
Schwenke, David W.
1990-01-01
The dissociation and recombination of H2 over the temperature range 1000-5000 K are calculated in a nonempirical manner. The computation procedure involves the calculation of the state-to-state energy transfer rate coefficients, the solution of the 349 coupled equations which form the master equation, and the determination of the phenomenological rate coefficients. The nonempirical results presented here are in good agreement with experimental data at 1000 and 3000 K.
NASA Astrophysics Data System (ADS)
Lu, Hua; Zhang, Shushu; Liu, Hanzhuang; Wang, Yanwei; Shen, Zhen; Liu, Chungen; You, Xiaozeng
2009-12-01
A boron-dipyrromethene (BODIPY)-based fluorescence probe with a N,N'-(pyridine-2, 6-diylbis(methylene))-dianiline substituent (1) has been prepared by condensation of 2,6-pyridinedicarboxaldehyde with 8-(4-amino)-4,4-difluoro-1,3,5,7-tetramethyl-4-bora-3a,4a-diaza-s-indacene and reduction by NaBH4. The sensing properties of compound 1 toward various metal ions are investigated via fluorometric titration in methanol, which show highly selective fluorescent turn-on response in the presence of Hg2+ over the other metal ions, such as Li+, Na+, K+, Ca2+, Mg2+, Pb2+, Fe2+, Co2+, Ni2+, Cu2+, Zn2+, Cd2+, Ag+, and Mn2+. Computational approach has been carried out to investigate the mechanism why compound 1 provides different fluorescent signal for Hg2+ and other ions. Theoretic calculations of the energy levels show that the quenching of the bright green fluorescence of boradiazaindacene fluorophore is due to the reductive photoinduced electron transfer (PET) from the aniline subunit to the excited state of BODIPY fluorophore. In metal complexes, the frontier molecular orbital energy levels changes greatly. Binding Zn2+ or Cd2+ ion leads to significant decreasing of both the HOMO and LUMO energy levels of the receptor, thus inhibit the reductive PET process, whereas an oxidative PET from the excited state fluorophore to the receptor occurs, vice versa, which also quenches the fluorescence. However, for 1-Hg2+ complex, both the reductive and oxidative PETs are prohibited; therefore, strong fluorescence emission from the fluorophore can be observed experimentally. The agreement of the experimental results and theoretic calculations suggests that our calculation method can be applicable as guidance for the design of new chemosensors for other metal ions.
NASA Astrophysics Data System (ADS)
Infantino, Angelo; Marengo, Mario; Baschetti, Serafina; Cicoria, Gianfranco; Longo Vaschetto, Vittorio; Lucconi, Giulia; Massucci, Piera; Vichi, Sara; Zagni, Federico; Mostacci, Domiziano
2015-11-01
Biomedical cyclotrons for production of Positron Emission Tomography (PET) radionuclides and radiotherapy with hadrons or ions are widely diffused and established in hospitals as well as in industrial facilities and research sites. Guidelines for site planning and installation, as well as for radiation protection assessment, are given in a number of international documents; however, these well-established guides typically offer analytic methods of calculation of both shielding and materials activation, in approximate or idealized geometry set up. The availability of Monte Carlo codes with accurate and up-to-date libraries for transport and interactions of neutrons and charged particles at energies below 250 MeV, together with the continuously increasing power of nowadays computers, makes systematic use of simulations with realistic geometries possible, yielding equipment and site specific evaluation of the source terms, shielding requirements and all quantities relevant to radiation protection. In this work, the well-known Monte Carlo code FLUKA was used to simulate two representative models of cyclotron for PET radionuclides production, including their targetry; and one type of proton therapy cyclotron including the energy selection system. Simulations yield estimates of various quantities of radiological interest, including the effective dose distribution around the equipment, the effective number of neutron produced per incident proton and the activation of target materials, the structure of the cyclotron, the energy degrader, the vault walls and the soil. The model was validated against experimental measurements and comparison with well-established reference data. Neutron ambient dose equivalent H*(10) was measured around a GE PETtrace cyclotron: an average ratio between experimental measurement and simulations of 0.99±0.07 was found. Saturation yield of 18F, produced by the well-known 18O(p,n)18F reaction, was calculated and compared with the IAEA recommended
Li, Yongqing; Yuan, Jiuchuang; Chen, Maodu; Ma, Fengcai; Sun, Mengtao
2013-07-15
An accurate single-sheeted double many-body expansion potential energy surface is reported for the title system. A switching function formalism has been used to warrant the correct behavior at the H2(X1Σg+)+N(2D) and NH (X3Σ-)+H(2S) dissociation channels involving nitrogen in the ground N(4S) and first excited N(2D) states. The topographical features of the novel global potential energy surface are examined in detail, and found to be in good agreement with those calculated directly from the raw ab initio energies, as well as previous calculations available in the literature. The novel surface can be using to treat well the Renner-Teller degeneracy of the 12A″ and 12A' states of NH 2. Such a work can both be recommended for dynamics studies of the N(2D)+H2 reaction and as building blocks for constructing the double many-body expansion potential energy surface of larger nitrogen/hydrogen-containing systems. In turn, a test theoretical study of the reaction N(2D)+H2(X1Σg+)(ν=0,j=0)→NH (X3Σ-)+H(2S) has been carried out with the method of quantum wave packet on the new potential energy surface. Reaction probabilities, integral cross sections, and differential cross sections have been calculated. Threshold exists because of the energy barrier (68.5 meV) along the minimum energy path. On the curve of reaction probability for total angular momentum J = 0, there are two sharp peaks just above threshold. The value of integral cross section increases quickly from zero to maximum with the increase of collision energy, and then stays stable with small oscillations. The differential cross section result shows that the reaction is a typical forward and backward scatter in agreement with experimental measurement result. PMID:23666848
Theoretical calculations and vibrational potential energy surface of 4-silaspiro(3,3)heptane
Ocola, Esther J.; Medders, Cross; Laane, Jaan; Meinander, Niklas
2014-04-28
Theoretical computations have been carried out on 4-silaspiro(3,3)heptane (SSH) in order to calculate its molecular structure and conformational energies. The molecule has two puckered four-membered rings with dihedral angles of 34.2° and a tilt angle of 9.4° between the two rings. Energy calculations were carried out for different conformations of SSH. These results allowed the generation of a two-dimensional ring-puckering potential energy surface (PES) of the form V = a(x{sub 1}{sup 4} + x{sub 2}{sup 4}) – b(x{sub 1}{sup 2} + x{sub 2}{sup 2}) + cx{sub 1}{sup 2}x{sub 2}{sup 2}, where x{sub 1} and x{sub 2} are the ring-puckering coordinates for the two rings. The presence of sufficiently high potential energy barriers prevents the molecule from undergoing pseudorotation. The quantum states, wave functions, and predicted spectra resulting from the PESs were calculated.
Skokov, S.; Peterson, K.A.; Bowman, J.M.
1998-08-01
Accurate {ital ab initio} multireference configuration interaction (CI) calculations with large correlation-consistent basis sets are performed for HOCl. After extrapolation to the complete basis set limit, the {ital ab initio} data are precisely fit to give a semiglobal three-dimensional potential energy surface to describe HOCl{r_arrow}Cl+OH from high overtone excitation of the OH-stretch. The average absolute deviation between the {ital ab initio} and fitted energies is 4.2thinspcm{sup {minus}1} for energies up to 60 kcal/mol relative to the HOCl minimum. Vibrational energies of HOCl including the six overtones of the OH-stretch are computed using a vibrational-Cl method on the fitted potential and also on a slightly adjusted potential. Near-spectroscopic accuracy is obtained using the adjusted potential; the average absolute deviation between theory and experiment for 19 experimentally reported states is 4.8thinspcm{sup {minus}1}. Very good agreement with experiment is also obtained for numerous rotational energies for the ground vibrational state, the ClO-stretch fundamental, and the fifth overtone of the OH-stretch. {copyright} {ital 1998 American Institute of Physics.}
NASA Astrophysics Data System (ADS)
Deguchi, Daiki; Sato, Kazunori; Kino, Hiori; Kotani, Takao
2016-05-01
We have recently implemented a new version of the quasiparticle self-consistent GW (QSGW) method in the ecalj package released at http://github.com/tkotani/ecalj. Since the new version of the ecalj package is numerically stable and more accurate than the previous versions, we can perform calculations easily without being bothered with tuning input parameters. Here we examine its ability to describe energy band properties, e.g., band-gap energy, eigenvalues at special points, and effective mass, for a variety of semiconductors and insulators. We treat C, Si, Ge, Sn, SiC (in 2H, 3C, and 4H structures), (Al, Ga, In) × (N, P, As, Sb), (Zn, Cd, Mg) × (O, S, Se, Te), SiO2, HfO2, ZrO2, SrTiO3, PbS, PbTe, MnO, NiO, and HgO. We propose that a hybrid QSGW method, where we mix 80% of QSGW and 20% of LDA, gives universally good agreement with experiments for these materials.
Den, Takuya S.; Frey, Hans-Martin; Leutwyler, Samuel
2014-11-21
The gas-phase rotational motion of hexafluorobenzene has been measured in real time using femtosecond (fs) time-resolved rotational Raman coherence spectroscopy (RR-RCS) at T = 100 and 295 K. This four-wave mixing method allows to probe the rotation of non-polar gas-phase molecules with fs time resolution over times up to ∼5 ns. The ground state rotational constant of hexafluorobenzene is determined as B{sub 0} = 1029.740(28) MHz (2σ uncertainty) from RR-RCS transients measured in a pulsed seeded supersonic jet, where essentially only the v = 0 state is populated. Using this B{sub 0} value, RR-RCS measurements in a room temperature gas cell give the rotational constants B{sub v} of the five lowest-lying thermally populated vibrationally excited states ν{sub 7/8}, ν{sub 9}, ν{sub 11/12}, ν{sub 13}, and ν{sub 14/15}. Their B{sub v} constants differ from B{sub 0} by between −1.02 MHz and +2.23 MHz. Combining the B{sub 0} with the results of all-electron coupled-cluster CCSD(T) calculations of Demaison et al. [Mol. Phys. 111, 1539 (2013)] and of our own allow to determine the C-C and C-F semi-experimental equilibrium bond lengths r{sub e}(C-C) = 1.3866(3) Å and r{sub e}(C-F) = 1.3244(4) Å. These agree with the CCSD(T)/wCVQZ r{sub e} bond lengths calculated by Demaison et al. within ±0.0005 Å. We also calculate the semi-experimental thermally averaged bond lengths r{sub g}(C-C)=1.3907(3) Å and r{sub g}(C-F)=1.3250(4) Å. These are at least ten times more accurate than two sets of experimental gas-phase electron diffraction r{sub g} bond lengths measured in the 1960s.
NASA Astrophysics Data System (ADS)
Kashinski, D. O.; Nelson, R. G.; Chase, G. M.; di Nallo, O. E.; Byrd, E. F. C.
2016-05-01
We are investigating the accuracy of theoretical models used to predict the visible, ultraviolet, and infrared spectra, as well as other properties, of product materials ejected from the muzzle of currently fielded systems. Recent advances in solid propellants has made the management of muzzle signature (flash) a principle issue in weapons development across the calibers. A priori prediction of the electromagnetic spectra of formulations will allow researchers to tailor blends that yield desired signatures and determine spectrographic detection ranges. Quantum chemistry methods at various levels of sophistication have been employed to optimize molecular geometries, compute unscaled harmonic frequencies, and determine the optical spectra of specific gas-phase species. Electronic excitations are being computed using Time Dependent Density Functional Theory (TD-DFT). Calculation of approximate global harmonic frequency scaling factors for specific DFT functionals is also in progress. A full statistical analysis and reliability assessment of computational results is currently underway. Work supported by the ARL, DoD-HPCMP, and USMA.
NASA Astrophysics Data System (ADS)
da Silva, Julio Cesar A.; Ducati, Lucas C.; Rittner, Roberto
2012-05-01
NMR solvent effects and theoretical calculations showed muscarinic agonists present a large stability for their near synclinal conformations, indicating the presence of significant stabilization factors. Analysis of the results clearly indicated that this stability is not determined by the dihedral around the substituted C-C ethane bond, as stated by some authors, but a consequence of the geometry adopted in order to maximize N+/O interactions in this type of molecules. It can be assumed that acetylcholine and its muscarinic agonists exhibit their biologic activity when the positively charged nitrogen and the oxygen atoms are in the same side of the molecule within an interatomic distance ranging from 3.0 to 6.0 Å.
Theoretical calculations of a compound formed by Fe(+3) and tris(catechol).
Kara, İzzet; Kara, Yeşim; Öztürk Kiraz, Aslı; Mammadov, Ramazan
2015-10-01
Phenolic compounds generally have special smell, easily soluble in water, organic solvents (alcohols, esters, chloroform, ethyl acetate), in aqueous solutions of bases, colorless or colorful, crystalline and amorphous materials. Phenols form colorful complexes when they form compounds with heavy metals. In this study, the structural properties of a compound formed by catechol and Fe(+3) are investigated theoretically. The electronic and thermodynamic properties of the complex were also investigated in gas phase and organic solvents at B3LYP/6-31+G(d,p) and B3LYP/6-311++G(d,p) basis set. The formation of Fe-tris(catechol) complex compound is exothermic, and it is difficult to obtain the complex as the temperature increases. The observed and calculated FT-IR and geometric parameters spectra are in good agreement with empirical. PMID:25983060
Rokhina, Ekaterina V; Lahtinen, Manu; Makarova, Katerina; Jegatheesan, Veeriah; Virkutyte, Jurate
2012-06-01
The nitric acid-functionalized commercial carbon nanofibers (CNFs) were comprehensively studied by instrumental (XRD, BET, SEM, TGA) and theoretical (DFT calculations) methods. The detailed surface study revealed the variation in the characteristics of functionalized CNFs, such as a decreased (up to 34%) surface area and impacted structural, electronic and chemical properties. The effects of functional groups were studied by comparison with pristine nanofibers. The results showed that the C-C bond lengths of the modified CNFs varied significantly. Chemical functionalization altered the frontier orbitals of the pristine material, and therefore altered the nature of their interactions with other substances. Moreover, the pristine and modified CNFs were tested for the removal of phenol from aqueous solutions. It was observed that surface modification tuned the adsorption capacity of carbon nanofibers (up to 0.35 mmol g(-1)), whereas original fibers did not demonstrate any adsorption capacity of phenol. PMID:22209137
Bellazzini, R.; Brez, A.; Massai, M.M.; Torquati, M.R.
1985-02-01
A fast and accurate algorithm to calculate the charge and current induced on all the electrodes of wire chambers (MWPC, 'pad chambers', TPC. . . .) is presented. The algorithm is completely three dimensional so that it is possible to calculate the induced charge on anode wires and cathode strips or 'pads' regardless of their orientation in space.
A 3-D Theoretical Model for Calculating Plasma Effects in Germanium Detectors
NASA Astrophysics Data System (ADS)
Wei, Wenzhao; Liu, Jing; Mei, Dongming; Cubed Collaboration
2015-04-01
In the detection of WIMP-induced nuclear recoil with Ge detectors, the main background source is the electron recoil produced by natural radioactivity. The capability of discriminating nuclear recoil (n) from electron recoil (γ) is crucial to WIMP searches. Digital pulse shape analysis is an encouraging approach to the discrimination of nuclear recoil from electron recoil since nucleus is much heavier than electron and heavier particle generates ionization more densely along its path, which forms a plasma-like cloud of charge that shields the interior from the influence of the electric field. The time needed for total disintegration of this plasma region is called plasma time. The plasma time depends on the initial density and radius of the plasma-like cloud, diffusion constant for charge carriers, and the strength of electric field. In this work, we developed a 3-D theoretical model for calculating the plasma time in Ge detectors. Using this model, we calculated the plasma time for both nuclear recoils and electron recoils to study the possibility for Ge detectors to realize n/ γ discrimination and improve detector sensitivity in detecting low-mass WIMPs. This work is supported by NSF in part by the NSF PHY-0758120, DOE Grant DE-FG02-10ER46709, and the State of South Dakota.
Accuracy of Theoretical Calculations for Electron-Impact Ionization of atoms and Molecules
NASA Astrophysics Data System (ADS)
Madison, Don
2015-09-01
In the last two decades, there have been several close-coupling approaches developed which can accurately calculate the triply differential cross sections for electron impact ionization of effective one and two electron atoms. The agreement between experiment and theory is not particularly good for more complicated atoms and molecules. Very recently, a B-spline R-matrix with pseudostates (BSRPS) approach was used to investigate low energy electron impact ionization of neon and very good agreement with experiment was found. The perturbative 3-body distorted wave (3DW) approach which includes the exact final state electron-electron interaction (post collision interaction - PCI) gave comparably good agreement with experiment. For ionization of molecules, there have been numerous studies of high-energy electron impact. These studies are called EMS (Electron Momentum Spectroscopy) and they were very valuable in determining the accuracy of molecular wavefunctions since the measured cross sections were proportional to the momentum space molecular wavefunction. More recently, lower energy collisions have started to be measured and these cross sections are much more difficult for theory since the detailed kinematics of the experiment become important. So far, the only close coupling calculation reported for ionization of molecules is the time-dependent close-coupling calculation (TDCC) which has been developed for ionization of H2 and it yields relative good agreement with experiment. Again the molecular 3-body distorted wave (M3DW) gave equally good agreement with experiment. For polyatomic molecules, the only theory available is the M3DW. In this talk, I will show the current status of agreement between experiment and theory for low and intermediate energy single ionization of atoms and molecules. Work supported by the NSF and XSEDE.
NASA Astrophysics Data System (ADS)
Wei, Guangfei; Li, Xiongyao; Wang, Shijie
2015-02-01
Terrestrial radiation is another possible source of heat in lunar thermal environment at its nearside besides the solar illumination. On the basis of Clouds and the Earth's Radiant Energy System (CERES) data products, the effect of terrestrial radiation on the brightness temperature (TBe) of the lunar nearside has been theoretically calculated. It shows that the mafic lunar mare with high TBe is more sensitive to terrestrial radiation than the feldspathic highland with low TBe value. According to the synchronous rotation of the Moon, we extract TBe on lunar nearside using the microwave radiometer data from the first Chinese lunar probe Chang'E-1 (CE-1). Consistently, the average TBe at Mare Serenitatis is about 1.2 K while the highland around the Geber crater (19.4°S, 13.9°E) is relatively small at ∼0.4 K. Our results indicate that there is no significant effect of terrestrial radiation on TBe at the lunar nearside. However, to extract TBe accurately, effects of heat flow, rock abundance and subsurface rock fragments which are more significant should be considered in the future work.
NASA Astrophysics Data System (ADS)
Predoi-Cross, Adriana; Holladay, Christopher; Heung, Henry; Bouanich, Jean-Pierre; Mellau, Georg Ch.; Keller, Reimund; Hurtmans, Daniel R.
2008-09-01
We report measurements for N 2-broadening, pressure-shift and line mixing coefficients for 55 oxygen transitions in the A-band retrieved using a multispectrum fitting technique. Nineteen laboratory absorption spectra were recorded at 0.02 cm -1 resolution using a multi-pass absorption cell with path length of 1636.9 cm and the IFS 120 Fourier transform spectrometer located at Justus-Liebig-University in Giessen, Germany. The total sample pressures ranged from 8.8 to 3004.5 Torr with oxygen volume mixing ratios in nitrogen ranging between 0.057 and 0.62. An Exponential Power Gap (EPG) scaling law was used to calculate the N 2-broadening and N 2-line mixing coefficients. The line broadening and shift coefficients for the A-band of oxygen self-perturbed and perturbed by N 2 are modeled using semiclassical calculations based on the Robert-Bonamy formalism and two intermolecular potentials. These potentials involve electrostatic contributions including the hexadecapole moment of the molecules and (a) a simple dispersion contribution with one adjustable parameter to fit the broadening coefficients or (b) the atom-atom Lennard-Jones model without such adjustable parameters. The first potential leads to very weak broadening coefficients for high J transitions whereas the second potential gives much more improved results at medium and large J values, in reasonable agreement with the experimental data. For the line shifts which mainly arise in our calculation from the electronic state dependence of the isotropic potential, their general trends with increasing J values can be well predicted, especially from the first potential. From the theoretical results, we have derived air-broadening and air-induced shift coefficients with an agreement comparable to that obtained for O 2-O 2 and O 2-N 2.
Duan, Yuhua
2012-11-02
Since current technologies for capturing CO{sub 2} to fight global climate change are still too energy intensive, there is a critical need for development of new materials that can capture CO{sub 2} reversibly with acceptable energy costs. Accordingly, solid sorbents have been proposed to be used for CO{sub 2} capture applications through a reversible chemical transformation. By combining thermodynamic database mining with first principles density functional theory and phonon lattice dynamics calculations, a theoretical screening methodology to identify the most promising CO{sub 2} sorbent candidates from the vast array of possible solid materials has been proposed and validated. The calculated thermodynamic properties of different classes of solid materials versus temperature and pressure changes were further used to evaluate the equilibrium properties for the CO{sub 2} adsorption/desorption cycles. According to the requirements imposed by the pre- and post- combustion technologies and based on our calculated thermodynamic properties for the CO{sub 2} capture reactions by the solids of interest, we were able to screen only those solid materials for which lower capture energy costs are expected at the desired pressure and temperature conditions. Only those selected CO{sub 2} sorbent candidates were further considered for experimental validations. The ab initio thermodynamic technique has the advantage of identifying thermodynamic properties of CO{sub 2} capture reactions without any experimental input beyond crystallographic structural information of the solid phases involved. Such methodology not only can be used to search for good candidates from existing database of solid materials, but also can provide some guidelines for synthesis new materials. In this presentation, we first introduce our screening methodology and the results on a testing set of solids with known thermodynamic properties to validate our methodology. Then, by applying our computational method
NASA Astrophysics Data System (ADS)
Pavlov, A. V.; Sitnov, Iu. S.
1985-10-01
Pavlov's (1984) method is used to determine the relative errors (due to errors in measuring the input parameters of the model) in theoretical calculations of the main parameters of the daytime F2-layer under quiet conditions. The parameters calculated are the height of the F2-layer maximum and the electron density.
Habegger, Maria L; Motta, Philip J; Huber, Daniel R; Dean, Mason N
2012-12-01
Evaluations of bite force, either measured directly or calculated theoretically, have been used to investigate the maximum feeding performance of a wide variety of vertebrates. However, bite force studies of fishes have focused primarily on small species due to the intractable nature of large apex predators. More massive muscles can generate higher forces and many of these fishes attain immense sizes; it is unclear how much of their biting performance is driven purely by dramatic ontogenetic increases in body size versus size-specific selection for enhanced feeding performance. In this study, we investigated biting performance and feeding biomechanics of immature and mature individuals from an ontogenetic series of an apex predator, the bull shark, Carcharhinus leucas (73-285cm total length). Theoretical bite force ranged from 36 to 2128N at the most anterior bite point, and 170 to 5914N at the most posterior bite point over the ontogenetic series. Scaling patterns differed among the two age groups investigated; immature bull shark bite force scaled with positive allometry, whereas adult bite force scaled isometrically. When the bite force of C. leucas was compared to those of 12 other cartilaginous fishes, bull sharks presented the highest mass-specific bite force, greater than that of the white shark or the great hammerhead shark. A phylogenetic independent contrast analysis of anatomical and dietary variables as determinants of bite force in these 13 species indicated that the evolution of large adult bite forces in cartilaginous fishes is linked predominantly to the evolution of large body size. Multiple regressions based on mass-specific standardized contrasts suggest that the evolution of high bite forces in Chondrichthyes is further correlated with hypertrophication of the jaw adductors, increased leverage for anterior biting, and widening of the head. Lastly, we discuss the ecological significance of positive allometry in bite force as a possible
Theoretical calculations of nonlinear refraction and absorption coefficients of doped graphene
NASA Astrophysics Data System (ADS)
Margulis, Vl A.; Muryumin, E. E.; Gaiduk, E. A.
2014-12-01
In this study, we present the first theoretical predictions concerning the nonlinear refractive and absorptive properties of the doped graphene in which the Fermi energy {{E}F} of charge carriers (noninteracting massless Dirac fermions) is controlled by an external gate voltage. We base our study on the original perturbation theory technique developed by Genkin and Mednis (1968 Sov. Phys. JETP 27 609) for calculating the nonlinear-optical (NLO) response coefficients of bulk crystalline semiconductors with partially filled bands. Using a simple tight-binding model for the π-electron energy bands of graphene, we obtain analytic expressions for the nonlinear refractive index {{n}2}(ω ) and the nonlinear absorption coefficient {{α }2}(ω ) of the doped graphene at photon energies above twice the value of the Fermi energy (\\hbar ω \\gt 2{{E}F}). We show that in this spectral region, both the nonlinear refraction ant the nonlinear absorption are determined predominantly by the combined processes which simultaneously involve intraband and interband motion of π-electrons. Our calculations indicate that extremely large negative values of n2 (of the order of -{{10}-6} cm2 W-1) can be achieved in the graphene at a relatively low doping level (of about 1012 cm-2) provided that the excitation frequency slightly exceeds the threshold frequency corresponding to the onset of interband transitions. With a further increase of the radiation frequency, the {{n}2}(ω ) becomes positive and begins to decrease in its absolute magnitude. The peculiar frequency dispersion of n2 and a negative sign of the {{α }2} (absorption bleaching), as predicted by our theory, suggest that the doped graphene is a prospective NLO material to be used in practical optical switching applications.
NASA Astrophysics Data System (ADS)
Rusakova, I. L.; Rusakov, Yu Yu; Krivdin, L. B.
2016-04-01
The theoretical grounds of the modern relativistic methods for quantum chemical calculation of spin–spin coupling constants in nuclear magnetic resonance spectra are considered. Examples and prospects of application of relativistic calculations of these constants in the structural studies of organic and heteroorganic compounds are discussed. Practical recommendations on relativistic calculations of spin–spin coupling constants using the available software are given. The bibliography includes 622 references.
Fu, Haohao; Shao, Xueguang; Chipot, Christophe; Cai, Wensheng
2016-08-01
Proper use of the adaptive biasing force (ABF) algorithm in free-energy calculations needs certain prerequisites to be met, namely, that the Jacobian for the metric transformation and its first derivative be available and the coarse variables be independent and fully decoupled from any holonomic constraint or geometric restraint, thereby limiting singularly the field of application of the approach. The extended ABF (eABF) algorithm circumvents these intrinsic limitations by applying the time-dependent bias onto a fictitious particle coupled to the coarse variable of interest by means of a stiff spring. However, with the current implementation of eABF in the popular molecular dynamics engine NAMD, a trajectory-based post-treatment is necessary to derive the underlying free-energy change. Usually, such a posthoc analysis leads to a decrease in the reliability of the free-energy estimates due to the inevitable loss of information, as well as to a drop in efficiency, which stems from substantial read-write accesses to file systems. We have developed a user-friendly, on-the-fly code for performing eABF simulations within NAMD. In the present contribution, this code is probed in eight illustrative examples. The performance of the algorithm is compared with traditional ABF, on the one hand, and the original eABF implementation combined with a posthoc analysis, on the other hand. Our results indicate that the on-the-fly eABF algorithm (i) supplies the correct free-energy landscape in those critical cases where the coarse variables at play are coupled to either each other or to geometric restraints or holonomic constraints, (ii) greatly improves the reliability of the free-energy change, compared to the outcome of a posthoc analysis, and (iii) represents a negligible additional computational effort compared to regular ABF. Moreover, in the proposed implementation, guidelines for choosing two parameters of the eABF algorithm, namely the stiffness of the spring and the mass
NASA Technical Reports Server (NTRS)
Cohen, S. C.
1979-01-01
A model of viscoelastic deformations associated with earthquakes is presented. A strike-slip fault is represented by a rectangular dislocation in a viscoelastic layer (lithosphere) lying over a viscoelastic half-space (asthenosphere). Deformations occur on three time scales. The initial response is governed by the instantaneous elastic properties of the earth. A slower response is associated with viscoelastic relaxation of the lithosphere and a yet slower response is due to viscoelastic relaxation of the asthenosphere. The major conceptual contribution is the inclusion of lithospheric viscoelastic properties into a dislocation model of earthquake related deformations and stresses. Numerical calculations using typical fault parameters reveal that the postseismic displacements and strains are small compared to the coseismic ones near the fault, but become significant further away. Moreover, the directional sense of the deformations attributable to the elastic response, the lithospheric viscoelastic softening, and the asthenospheric viscoelastic flow may differ and depend on location and model details. The results and theoretical arguments suggest that the stress changes accompanying lithospheric relaxation may also be in a different sense than and be larger than the strain changes.
Wu, Yanlin; Shi, Jin; Chen, Hongche; Zhao, Jianfu; Dong, Wenbo
2016-10-01
4-tert-butylphenol (4-t-BP), an endocrine disrupting chemical, is widely distributed in natural bodies of water but is difficult to biodegrade. In this study, we focused on the transformation of 4-t-BP in photo-initiated degradation processes. The steady-state photolysis and laser flash photolysis (LFP) experiments were conducted in order to elucidate its degradation mechanism. Identification of products was performed using the GC-MS, LC-MS and theoretical calculation techniques. The oxidation pathway of 4-t-BP by hydroxyl radical (HO) was also studied and H2O2 was added to produce HO. 4-tert-butylcatechol and 4-tert-butylphenol dimer were produced in 4-t-BP direct photolysis. 4-tert-butylcatechol and hydroquinone were produced by the oxidation of HO. But the formation mechanism of 4-tert-butylcatechol in the two processes was different. The benzene ring was fractured in 4-t-BP oxidation process and 29% of TOC was degraded after 16h irradiation. PMID:27213674
NASA Astrophysics Data System (ADS)
Takeshima, Tsuguhide; Takeuchi, Hiroshi; Egawa, Toru; Konaka, Shigehiro
2005-01-01
The molecular structure of arecoline (methyl 1,2,5,6-tetrahydro-1-methylnicotinate, ? has been determined by gas electron diffraction. Diffraction patterns were taken at about 370 K. Structural constraints for the data analysis were obtained from MP2/6-31G** calculations. Vibrational mean amplitudes and shrinkage corrections were calculated from the force constants obtained from the gas-phase vibrational frequencies and the B3LYP/6-31G** calculations. The electron diffraction data were well reproduced by assuming the mixture of four conformers. The determined structural parameters ( rg (Å) and ∠ (°)) for the main conformer with 3 σ in parentheses are as follows: < rg(N-C ring)>=1.456(4); rg(N-C methyl)=1.451 (d.p.); rg(C dbnd6 C)=1.339(9); < rg(C-C)>=1.512(3); rg(O-C methyl)=1.434(5); rg(C(O)-O)=1.355 (d.p.); rg(C dbnd6 O)=1.209(4); the out-of-plane angle of the methyl group=50.3(23); ∠C ringN ringC ring=112.8(30); ∠N ringC ringC ring(H 2)=110.5(16); <∠C ringC ringC ring>=118.4(5); ∠C dbnd6 CC(O)=116.8(7); ∠CC dbnd6 O=127.6(9); ∠CC-O=109.8(8), where the angle brackets denote averaged values and d.p. denotes dependent parameters. Fixing the abundances of the minor conformers, Ax-s- cis and Ax-s- trans, at the theoretical values (13% in total), those of the Eq-s- cis and Eq-s- trans conformers were determined to be 46(16) and 41(16)%, respectively. Here Ax and Eq denote the axial and equatorial directions of the N-CH 3 bond and s- cis and s- trans show the orientation of the methoxycarbonyl group expressed by the configuration of the C dbnd6 O and C dbnd6 C bonds. The N⋯O carbonyl distances of the Eq-s- cis and Ax-s- cis conformers are 4.832(13) and 4.874(16) Å, respectively. They are close to the N⋯N distance of the most abundant conformer of nicotine, 4.885(6) Å, suggesting that the Eq-s- cis and Ax-s- cis conformers have nicotinic activity.
Chen, Xiao-bo; Wang, Ce; Li, Song; Naruhito, Sawanobori; Kang, Dong-guo
2010-07-01
A photonic phenomenon of fluorescence intensity reverse between red and green fluorescence was studied theoretically and experimentally in the present article. It was found by experiment that Er(0.5) Yb(9.5) : FOV oxyfluoride vitroceramics exhibits strong fluorescence intensity reverse phenomenon. The range of the intensity reverse of Er(0.5) Yb(9.5) : FOV was measured to be 877. Moreover, all basic spectroscopic parameters were calculated. The theoretical basis of numerical calculation for dynamics processes of all levels was established. PMID:20827959
NASA Astrophysics Data System (ADS)
Sun, Yuansheng; Periasamy, Ammasi
2010-03-01
Förster resonance energy transfer (FRET) microscopy is commonly used to monitor protein interactions with filter-based imaging systems, which require spectral bleedthrough (or cross talk) correction to accurately measure energy transfer efficiency (E). The double-label (donor+acceptor) specimen is excited with the donor wavelength, the acceptor emission provided the uncorrected FRET signal and the donor emission (the donor channel) represents the quenched donor (qD), the basis for the E calculation. Our results indicate this is not the most accurate determination of the quenched donor signal as it fails to consider the donor spectral bleedthrough (DSBT) signals in the qD for the E calculation, which our new model addresses, leading to a more accurate E result. This refinement improves E comparisons made with lifetime and spectral FRET imaging microscopy as shown here using several genetic (FRET standard) constructs, where cerulean and venus fluorescent proteins are tethered by different amino acid linkers.
Theoretical methods for the calculation of Bragg curves and 3D distributions of proton beams
NASA Astrophysics Data System (ADS)
Ulmer, W.; Matsinos, E.
2010-12-01
buildup is explained by a theoretical analysis of impinging proton beamlets. In order to obtain a complete dose calculation model for proton treatment planning, some further aspects have to be accounted for: the decrease of the fluence of the primary protons due to nuclear interactions, the transport of released secondary protons, the dose contribution of heavy recoil nuclei, the inclusion of lateral scatter of the primary and secondary protons based on Molière's multiple-scatter theory, and the scatter contributions of collimators. This study also presents some results which go beyond proton dose calculation models; namely, the application of the relativistic generalization of the Bragg-Kleeman rule to electrons and, in an appendix, a method to determine inelastic cross-sections of therapeutic protons in media of therapeutic interest.
NASA Astrophysics Data System (ADS)
Umari, P.; Fabris, S.
2012-05-01
The quasi-particle energy levels of the Zn-Phthalocyanine (ZnPc) molecule calculated with the GW approximation are shown to depend sensitively on the explicit description of the metal-center semicore states. We find that the calculated GW energy levels are in good agreement with the measured experimental photoemission spectra only when explicitly including the Zn 3s and 3p semicore states in the valence. The main origin of this effect is traced back to the exchange term in the self-energy GW approximation. Based on this finding, we propose a simplified approach for correcting GW calculations of metal phthalocyanine molecules that avoids the time-consuming explicit treatment of the metal semicore states. Our method allows for speeding up the calculations without compromising the accuracy of the computed spectra.
Communication: An accurate calculation of the S1 C2H2 cis-trans isomerization barrier height
NASA Astrophysics Data System (ADS)
Baraban, Joshua H.; Matthews, Devin A.; Stanton, John F.
2016-03-01
A high level ab initio calculation of the cis-trans isomerization barrier height in the first excited singlet electronic state of acetylene is found to agree very well with a recent experimental determination.
Peterson, K.A.; Skokov, S.; Bowman, J.M.
1999-10-01
A new, global analytical potential energy surface is constructed for the X&hthinsp;{sup 1}A{sup {prime}} electronic ground state of HOCl that accurately includes the HClO isomer. The potential is obtained by using accurate {ital ab initio} data from a previously published surface [Skokov {ital et al.}, J. Chem. Phys. {bold 109}, 2662 (1998)], as well as a significant number of new data for the HClO region of the surface at the same multireference configuration interaction, complete basis set limit level of theory. Vibrational energy levels and intensities are computed for both HOCl and HClO up to the OH+Cl dissociation limit and above the isomerization barrier. After making only minor adjustments to the {ital ab initio} surface, the errors with respect to experiment for HOCl are generally within a few cm{sup {minus}1} for 22 vibrational levels with the largest error being 26 cm{sup {minus}1}. A total of 813 bound vibrational states are calculated for HOCl. The HClO potential well supports 57 localized states, of which only the first 3 are bound. The strongest dipole transitions for HClO were computed for the fundamentals{emdash}33, 2.9, and 25 km/mol for {nu}{sub 1}, {nu}{sub 2}, and {nu}{sub 3}, respectively. From exact J=1 ro-vibrational calculations, state dependent rotational constants have been calculated for HClO. Lastly, resonance calculations with the new potential demonstrate that the presence of the HClO minimum has a negligible effect on the resonance states of HOCl near the dissociation threshold due to the relatively high and wide isomerization barrier. {copyright} {ital 1999 American Institute of Physics.}
NASA Astrophysics Data System (ADS)
Takeshima, Tsuguhide; Takeuchi, Hiroshi; Egawa, Toru; Konaka, Shigehiro
2007-09-01
The molecular structure of cotinine (( S)-1-methyl-5-(3-pyridinyl)-2-pyrrolidinone), the major metabolite of nicotine, has been determined at about 182 °C by gas electron diffraction combined with MP2 and DFT calculations. The diffraction data are consistent with the existence of the (ax, sc), (ax, ap), (eq, sp) and (eq, ap) conformers, where ax and eq indicate the configuration of the pyrrolidinone ring by means of the position (axial and equatorial) of the pyridine ring, and sc, sp and ap distinguish the isomers arising from the internal rotation around the bond connecting the two rings. The (CH 3)NCCC(N) dihedral angles, ϕ, of the (ax, sc) and (eq, sp) conformers were determined independently to be 158(12)° and 129(13)°, respectively, where the numbers in parentheses are three times the standard errors, 3 σ. According to the MP2 calculations, the corresponding dihedral angles for the (ax, ap) and (eq, ap) conformers were assumed to differ by 180° from their syn counterparts. The ratios x(ax, sc)/ x(ax, ap) and x(eq, sp)/ x(eq, ap) were taken from the theoretically estimated free energy differences, Δ G, where x is the abundance of the conformer. The resultant abundances of (ax, sc), (ax, ap), (eq, sp) and (eq, ap) conformers are 34(6)%, 21% (d.p.), 28% (d.p.), and 17% (d.p.), respectively, where d.p. represents dependent parameters. The determined structural parameters ( rg (Å) and ∠ α (°)) of the most abundant conformer, (ax, sc), are as follows: r(N sbnd C) pyrrol = 1.463(5); r(N sbnd C methyl) = 1.457(←); r(N sbnd C( dbnd O)) = 1.384(12); r(C dbnd O) = 1.219(5); < r(C sbnd C) pyrrol> = 1.541(3); r(C pyrrolsbnd C pyrid) = 1.521(←); < r(C sbnd C) pyrid> = 1.396(2); < r(C sbnd N) pyrid> = 1.343(←); ∠(CNC) pyrrol = 113.9(11); ∠CCC pyrrol(-C pyrid) = 103.6(←); ∠NCO = 124.1(13); ∠NC pyrrolC pyrid = 113.1(12); ∠C pyrrolC pyrrolC pyrid = 113.3(←); ∠(CNC) pyrid = 117.1(2); <∠(NCC) pyrid> = 124.4(←); ∠C methylNC( dbnd O) =
NASA Astrophysics Data System (ADS)
Brawand, Nicholas; Vörös, Márton; Govoni, Marco; Galli, Giulia
The accurate prediction of optoelectronic properties of molecules and solids is a persisting challenge for current density functional theory (DFT) based methods. We propose a hybrid functional where the mixing fraction of exact and local exchange is determined by a non-empirical, system dependent function. This functional yields ionization potentials, fundamental and optical gaps of many, diverse systems in excellent agreement with experiments, including organic and inorganic molecules and nanocrystals. We further demonstrate that the newly defined hybrid functional gives the correct alignment between the energy level of the exemplary TTF-TCNQ donor-acceptor system. DOE-BES: DE-FG02-06ER46262.
Ishida, Toyokazu
2008-09-17
To further understand the catalytic role of the protein environment in the enzymatic process, the author has analyzed the reaction mechanism of the Claisen rearrangement of Bacillus subtilis chorismate mutase (BsCM). By introducing a new computational strategy that combines all-electron QM calculations with ab initio QM/MM modelings, it was possible to simulate the molecular interactions between the substrate and the protein environment. The electrostatic nature of the transition state stabilization was characterized by performing all-electron QM calculations based on the fragment molecular orbital technique for the entire enzyme.
Beare, Richard; Brown, Michael J. I.; Pimbblet, Kevin
2014-12-20
We describe an accurate new method for determining absolute magnitudes, and hence also K-corrections, that is simpler than most previous methods, being based on a quadratic function of just one suitably chosen observed color. The method relies on the extensive and accurate new set of 129 empirical galaxy template spectral energy distributions from Brown et al. A key advantage of our method is that we can reliably estimate random errors in computed absolute magnitudes due to galaxy diversity, photometric error and redshift error. We derive K-corrections for the five Sloan Digital Sky Survey filters and provide parameter tables for use by the astronomical community. Using the New York Value-Added Galaxy Catalog, we compare our K-corrections with those from kcorrect. Our K-corrections produce absolute magnitudes that are generally in good agreement with kcorrect. Absolute griz magnitudes differ by less than 0.02 mag and those in the u band by ∼0.04 mag. The evolution of rest-frame colors as a function of redshift is better behaved using our method, with relatively few galaxies being assigned anomalously red colors and a tight red sequence being observed across the whole 0.0 < z < 0.5 redshift range.
Nalli, Nicholas R; Minnett, Peter J; van Delst, Paul
2008-07-20
Although published sea surface infrared (IR) emissivity models have gained widespread acceptance for remote sensing applications, discrepancies have been identified against field observations obtained from IR Fourier transform spectrometers at view angles approximately > 40 degrees. We therefore propose, in this two-part paper, an alternative approach for calculating surface-leaving IR radiance that treats both emissivity and atmospheric reflection in a systematic yet practical manner. This first part presents the theoretical basis, development, and computations of the proposed model. PMID:18641735
Accurate calculation of the x-ray absorption spectrum of water via the GW/Bethe-Salpeter equation
NASA Astrophysics Data System (ADS)
Gilmore, Keith; Vinson, John; Kas, Josh; Vila, Fernando; Rehr, John
2014-03-01
We calculate x-ray absorption spectra (XAS) of water within the OCEAN code, which combines plane-wave, pseudopotential electronic structure, PAW transition elements, GW self-energy corrections, and the NIST BSE solver. Due to the computational demands of this approach, our initial XAS calculations were limited to 17 molecule super cells. This lead to unphysical, size dependent effects in the calculated spectra. To treat larger systems, we extended the OCEAN interface to support well-parallelized codes such as QuantumESPRESSO. We also implemented an efficient interpolation scheme of Shirley. We applied this large-scale GW/BSE approach to 64 molecule unit cell structures of water obtained from classical DFT/MD and PIMD simulations. In concurrence with previous work, we find the calculated spectrum both qualitatively and quantitatively reproduces the experimental features. The agreement implies that structures based on PIMD, which are similar to the traditional distorted tetrahedral view, are consistent with experimental observations. Supported by the DOE CMCSN through DOE award DE-SC0005180 (Princeton University) and in part by DOE Grant No. DE-FG03-97ER45623 (JJR) with computer support from NERSC.
NASA Astrophysics Data System (ADS)
Xu, Zhongnan; Joshi, Yogesh V.; Raman, Sumathy; Kitchin, John R.
2015-04-01
We validate the usage of the calculated, linear response Hubbard U for evaluating accurate electronic and chemical properties of bulk 3d transition metal oxides. We find calculated values of U lead to improved band gaps. For the evaluation of accurate reaction energies, we first identify and eliminate contributions to the reaction energies of bulk systems due only to changes in U and construct a thermodynamic cycle that references the total energies of unique U systems to a common point using a DFT + U(V ) method, which we recast from a recently introduced DFT + U(R) method for molecular systems. We then introduce a semi-empirical method based on weighted DFT/DFT + U cohesive energies to calculate bulk oxidation energies of transition metal oxides using density functional theory and linear response calculated U values. We validate this method by calculating 14 reactions energies involving V, Cr, Mn, Fe, and Co oxides. We find up to an 85% reduction of the mean average error (MAE) compared to energies calculated with the Perdew-Burke-Ernzerhof functional. When our method is compared with DFT + U with empirically derived U values and the HSE06 hybrid functional, we find up to 65% and 39% reductions in the MAE, respectively.
Xu, Zhongnan; Kitchin, John R.; Joshi, Yogesh V.; Raman, Sumathy
2015-04-14
We validate the usage of the calculated, linear response Hubbard U for evaluating accurate electronic and chemical properties of bulk 3d transition metal oxides. We find calculated values of U lead to improved band gaps. For the evaluation of accurate reaction energies, we first identify and eliminate contributions to the reaction energies of bulk systems due only to changes in U and construct a thermodynamic cycle that references the total energies of unique U systems to a common point using a DFT + U(V ) method, which we recast from a recently introduced DFT + U(R) method for molecular systems. We then introduce a semi-empirical method based on weighted DFT/DFT + U cohesive energies to calculate bulk oxidation energies of transition metal oxides using density functional theory and linear response calculated U values. We validate this method by calculating 14 reactions energies involving V, Cr, Mn, Fe, and Co oxides. We find up to an 85% reduction of the mean average error (MAE) compared to energies calculated with the Perdew-Burke-Ernzerhof functional. When our method is compared with DFT + U with empirically derived U values and the HSE06 hybrid functional, we find up to 65% and 39% reductions in the MAE, respectively.
NASA Technical Reports Server (NTRS)
Burr, D. M.; Emery, J. P.; Lorenz, R. D.
2005-01-01
The Cassini Imaging Science System (ISS) has been returning images of Titan, along with other Saturnian satellites. Images taken through the 938 nm methane window see down to Titan's surface. One of the purposes of the Cassini mission is to investigate possible fluid cycling on Titan. Lemniscate features shown recently and radar evidence of surface flow prompted us to consider theoretically the creation by methane fluid flow of streamlined forms on Titan. This follows work by other groups in theoretical consideration of fluid motion on Titan's surface.
Chikayama, Eisuke; Shimbo, Yudai; Komatsu, Keiko; Kikuchi, Jun
2016-04-14
NMR spectroscopy is a powerful method for analyzing metabolic mixtures. The information obtained from an NMR spectrum is in the form of physical parameters, such as chemical shifts, and construction of databases for many metabolites will be useful for data interpretation. To increase the accuracy of theoretical chemical shifts for development of a database for a variety of metabolites, the effects of sets of conformations (structural ensembles) and the levels of theory on computations of theoretical chemical shifts were systematically investigated for a set of 29 small molecules in the present study. For each of the 29 compounds, 101 structures were generated by classical molecular dynamics at 298.15 K, and then theoretical chemical shifts for 164 (1)H and 123 (13)C atoms were calculated by ab initio quantum chemical methods. Six levels of theory were used by pairing Hartree-Fock, B3LYP (density functional theory), or second order Møller-Plesset perturbation with 6-31G or aug-cc-pVDZ basis set. The six average fluctuations in the (1)H chemical shift were ±0.63, ± 0.59, ± 0.70, ± 0.62, ± 0.75, and ±0.66 ppm for the structural ensembles, and the six average errors were ±0.34, ± 0.27, ± 0.32, ± 0.25, ± 0.32, and ±0.25 ppm. The results showed that chemical shift fluctuations with changes in the conformation because of molecular motion were larger than the differences between computed and experimental chemical shifts for all six levels of theory. In conclusion, selection of an appropriate structural ensemble should be performed before theoretical chemical shift calculations for development of an accurate database for a variety of metabolites. PMID:26963288
NASA Astrophysics Data System (ADS)
Hsieh, Tiane-Jye; Su, Chia-Ching; Chen, Chung-Yi; Liou, Chyong-Huey; Lu, Li-Hwa
2005-05-01
Three natural products, Coumarin ( 1), p-hydroxybenzoic acid ( 2), trans-cinnamic acid ( 3) were isolated from the natural plant of indigenous cinnamon and the structures including relative stereochemistry were elucidated on the basis of spectroscopic data and theoretical calculations. Their sterochemical structures were determined by NMR spectroscopy, mass spectroscopy, and X-ray crystallography. The p-hydroxybenzoic acid complex with water is reported to show the existence of two hydrogen bonds. The two hydrogen bonds are formed in the water molecule of two hydrogen-accepting oxygen of carbonyl group of the p-hydroxybenzoic acid. The intermolecular interaction two hydrogen bond of the model system of the water- p-hydroxybenzoic acid was investigated. An experimental study and a theoretical analysis using the B3LYP/6-31G* method in the GAUSSIAN-03 package program were conducted on the three natural products. The theoretical results are supplemented by experimental data. Optimal geometric structures of three compounds were also determined. The calculated molecular mechanics compared quite well with those obtained from the experimental data. The ionization potentials, highest occupied molecular orbital energy, lowest unoccupied molecular orbital energy, energy gaps, heat of formation, atomization energies, and vibration frequencies of the compounds were also calculated. The results of the calculations show that three natural products are stable molecules with high reactive and various other physical properties. The study also provided an explicit understanding of the sterochemical structure and thermodynamic properties of the three natural products.
Karton, Amir; O'Reilly, Robert J; Radom, Leo
2012-04-26
Accurate electronic barrier heights are obtained for a set of nine proton-transfer tautomerization reactions, which are either (i) uncatalyzed, (ii) catalyzed by one water molecule, or (iii) catalyzed by two water molecules. The barrier heights for reactions (i) and (ii) are obtained by means of the high-level ab initio W2.2 thermochemical protocol, while those for reaction (iii) are obtained using the W1 protocol. These three sets of benchmark barrier heights allow an assessment of the performance of more approximate theoretical procedures for the calculation of barrier heights of uncatalyzed and water-catalyzed reactions. We evaluate initially the performance of the composite G4 procedure and variants thereof (e.g., G4(MP2) and G4(MP2)-6X), as well as that of standard ab initio procedures (e.g., MP2, SCS-MP2, and MP4). We find that the performance of the G4(MP2)-type thermochemical procedures deteriorates with the number of water molecules involved in the catalysis. This behavior is linked to deficiencies in the MP2-based basis-set-correction term in the G4(MP2)-type procedures. This is remedied in the MP4-based G4 procedure, which shows good performance for both the uncatalyzed and the water-catalyzed reactions, with mean absolute deviations (MADs) from the benchmark values lying below the threshold of "chemical accuracy" (arbitrarily defined as 1 kcal mol(-1) ≈ 4.2 kJ mol(-1)). We also examine the performance of a large number of density functional theory (DFT) and double-hybrid DFT (DHDFT) procedures. We find that, with few exceptions (most notably PW6-B95 and B97-2), the performance of the DFT procedures that give good results for the uncatalyzed reactions deteriorates with the number of water molecules involved in the catalysis. The DHDFT procedures, on the other hand, show excellent performance for both the uncatalyzed and catalyzed reactions. Specifically, almost all of them afford MADs below the "chemical accuracy" threshold, with ROB2-PLYP and B2K
Kanai, Y; Takeuchi, N
2009-10-14
We revisit the molecular line growth mechanism of styrene on the hydrogenated Si(001) 2x1 surface. In particular, we investigate the energetics of the radical chain reaction mechanism by means of diffusion quantum Monte Carlo (QMC) and density functional theory (DFT) calculations. For the exchange correlation (XC) functional we use the non-empirical generalized-gradient approximation (GGA) and meta-GGA. We find that the QMC result also predicts the intra dimer-row growth of the molecular line over the inter dimer-row growth, supporting the conclusion based on DFT results. However, the absolute magnitudes of the adsorption and reaction energies, and the heights of the energy barriers differ considerably between the QMC and DFT with the GGA/meta-GGA XC functionals.
Yamamoto, Atsushi; Matsui, Yasunori; Asada, Toshio; Kumeda, Motoki; Takagi, Kenichiro; Suenaga, Yu; Nagae, Kunihiko; Ohta, Eisuke; Sato, Hiroyasu; Koseki, Shiro; Naito, Hiroyoshi; Ikeda, Hiroshi
2016-04-15
The p-type organic semiconductor (OSC) material tetrathieno[2,3-a:3',2'-c:2″,3″-f:3‴,2‴-h]naphthalene (2TTN) and its alkyl-substituted derivatives Cn-2TTNs (n = 6, 8, and 10) have been developed based on the results of theoretical calculation-inspired investigation. A hole mobility for amorphous Cn-2TTNs (10(-2)-10(-3) cm(2) V(-1) s(-1)) was accurately predicted by using a novel statistical method in which the geometric mean of the mobilities for many individual small molecular flocks in an amorphous solid was obtained by using molecular mechanical molecular dynamics simulations and quantum chemical calculations. The simulation also suggests that upon increasing the length of alkyl chains in Cn-2TTNs the mobilities become smaller as a consequence of a decrease in transfer integral values. Cn-2TTNs are synthesized in a microflow reactor through photoreactions of the corresponding precursors. Cn-2TTNs are then utilized in the fabrication of organic field-effect transistors (OFETs). Although spin-coated thin films of Cn-2TTNs are crystalline, the hole mobilities (10(-2)-10(-3) cm(2) V(-1) s(-1)) of trial OFETs decrease upon elongation of the alkyl chains. This finding parallels the results of theoretical simulation. The simulation method for amorphous solids developed in this effort should become a useful tool in studies aimed at designing new OSC materials. PMID:27010327
Wood, R.F.; Wilson, T.M.
1981-01-01
The structure of the Hartree-Fock one-electron equations for simple point defects in ionic crystals are discussed. The importance of polarization effects due to the diffuse nature of the wavefunctions in the relaxed excited states are emphasized, and the usefulness of an effective mass approximation indicated. Several approaches to the calculation of the electronic structure are discussed and evaluated. The connection between electronic structure calculations and phonon perturbations are pointed out through a brief discussion of localized perturbation theory.
NASA Technical Reports Server (NTRS)
Mccluney, W. R.
1974-01-01
The development is considered of procedures for measuring a number of subsurface oceanographic parameters using remotely sensed ocean color data. It is proposed that the first step in this effort should be the development of adequate theoretical models relating the desired oceanographic parameters to the upwelling radiances to be observed. A portion of a contributory theoretical model is shown to be described by a modified single scattering approach based upon a simple treatment of multiple scattering. The resulting quasi-single scattering model can be used to predict the upwelling distribution of spectral radiance emerging from the sea. The shape of the radiance spectrum predicted by this model for clear ocean water shows encouraging agreement with measurments made at the edge of the Sargasso Sea off Cape Hatteras.
Theoretical calculations of emission of wolframite and scheelite-type tungstate crystals
Nikolaenko, T.; Hizhnyi, Y.; Nedilko, S.
2009-01-21
Tungstate crystals AWO{sub 4} (A = Zn,Cd,Pb) are well-known scintillation materials for various applications in science and technology. In recent years the optical properties of these crystals were intensively studied experimentally and theoretically. However, the origin of luminescence in lead, cadmium and zinc tungstates is still the subject of discussion. According to generally accepted view, the centers of luminescence in AWO{sub 4} crystals are in some or other way related to the tungstate anionic groups. We developed a cluster approach in theoretical investigation of the electronic structure of AWO{sub 4} tungstate crystals based on the configuration interaction (CI) computation in which the lattice vibrations were taken into account.
NASA Astrophysics Data System (ADS)
Badnell, N. R.; Spruck, K.; Krantz, C.; Novotný, O.; Becker, A.; Bernhardt, D.; Grieser, M.; Hahn, M.; Repnow, R.; Savin, D. W.; Wolf, A.; Müller, A.; Schippers, S.
2016-05-01
Experimentally measured and theoretically calculated rate coefficients for the recombination of W19 +([Kr ] 4 d10 4 f9 ) ions with free electrons (forming W18 +) are presented. At low electron-ion collision energies, the merged-beam rate coefficient is dominated by strong, mutually overlapping, recombination resonances as already found previously for the neighboring charge-state ions W18 + and W20 +. In the temperature range where W19 + is expected to form in a collisionally ionized plasma, the experimentally derived recombination rate coefficient deviates by up to a factor of about 20 from the theoretical rate coefficient obtained from the Atomic Data and Analysis Structure database. The present calculations, which employ a Breit-Wigner redistributive partitioning of autoionizing widths for dielectronic recombination via multi-electron resonances, reproduce the experimental findings over the entire temperature range.
Johnston, Ryne C.; Zhou, Jing; Smith, Jeremy C.; Parks, Jerry M.
2016-07-08
In redox processes in complex transition metal-containing species are often intimately associated with changes in ligand protonation states and metal coordination number. Moreover, a major challenge is therefore to develop consistent computational approaches for computing pH-dependent redox and ligand dissociation properties of organometallic species. Reduction of the Co center in the vitamin B12 derivative aquacobalamin can be accompanied by ligand dissociation, protonation, or both, making these properties difficult to compute accurately. We examine this challenge here by using density functional theory and continuum solvation to compute Co ligand binding equilibrium constants (Kon/off), pKas and reduction potentials for models of aquacobalaminmore » in aqueous solution. We consider two models for cobalamin ligand coordination: the first follows the hexa, penta, tetra coordination scheme for CoIII, CoII, and CoI species, respectively, and the second model features saturation of each vacant axial coordination site on CoII and CoI species with a single, explicit water molecule to maintain six directly interacting ligands or water molecules in each oxidation state. Comparing these two coordination schemes in combination with five dispersion-corrected density functionals, we find that the accuracy of the computed properties is largely independent of the scheme used, but including only a continuum representation of the solvent yields marginally better results than saturating the first solvation shell around Co throughout. PBE performs best, displaying balanced accuracy and superior performance overall, with RMS errors of 80 mV for seven reduction potentials, 2.0 log units for five pKas and 2.3 log units for two log Kon/off values for the aquacobalamin system. Furthermore, we find that the BP86 functional commonly used in corrinoid studies suffers from erratic behavior and inaccurate descriptions of Co axial ligand binding, leading to substantial errors in predicted
Zhang, Xinyue; Lourenco, Daniela; Aguilar, Ignacio; Legarra, Andres; Misztal, Ignacy
2016-01-01
former. Manhattan plots had higher resolution with 5 and 100 QTL. Using a common weight for a window of 20 SNP that sums or averages the SNP variance enhances accuracy of predicting GEBV and provides accurate estimation of marker effects. PMID:27594861
Zhang, Xinyue; Lourenco, Daniela; Aguilar, Ignacio; Legarra, Andres; Misztal, Ignacy
2016-01-01
. Manhattan plots had higher resolution with 5 and 100 QTL. Using a common weight for a window of 20 SNP that sums or averages the SNP variance enhances accuracy of predicting GEBV and provides accurate estimation of marker effects. PMID:27594861
Johnston, Ryne C; Zhou, Jing; Smith, Jeremy C; Parks, Jerry M
2016-08-01
Redox processes in complex transition metal-containing species are often intimately associated with changes in ligand protonation states and metal coordination number. A major challenge is therefore to develop consistent computational approaches for computing pH-dependent redox and ligand dissociation properties of organometallic species. Reduction of the Co center in the vitamin B12 derivative aquacobalamin can be accompanied by ligand dissociation, protonation, or both, making these properties difficult to compute accurately. We examine this challenge here by using density functional theory and continuum solvation to compute Co-ligand binding equilibrium constants (Kon/off), pKas, and reduction potentials for models of aquacobalamin in aqueous solution. We consider two models for cobalamin ligand coordination: the first follows the hexa, penta, tetra coordination scheme for Co(III), Co(II), and Co(I) species, respectively, and the second model features saturation of each vacant axial coordination site on Co(II) and Co(I) species with a single, explicit water molecule to maintain six directly interacting ligands or water molecules in each oxidation state. Comparing these two coordination schemes in combination with five dispersion-corrected density functionals, we find that the accuracy of the computed properties is largely independent of the scheme used, but including only a continuum representation of the solvent yields marginally better results than saturating the first solvation shell around Co throughout. PBE performs best, displaying balanced accuracy and superior performance overall, with RMS errors of 80 mV for seven reduction potentials, 2.0 log units for five pKas and 2.3 log units for two log Kon/off values for the aquacobalamin system. Furthermore, we find that the BP86 functional commonly used in corrinoid studies suffers from erratic behavior and inaccurate descriptions of Co-axial ligand binding, leading to substantial errors in predicted pKas and
Hodder, Joanne N; La Delfa, Nicholas J; Potvin, Jim R
2016-08-01
To predict shoulder strength, most current ergonomics software assume independence of the strengths about each of the orthopedic axes. Using this independent axis approach (IAA), the shoulder can be predicted to have strengths as high as the resultant of the maximum moment about any two or three axes. We propose that shoulder strength is not independent between axes, and propose an approach that calculates the weighted average (WAA) between the strengths of the axes involved in the demand. Fifteen female participants performed maximum isometric shoulder exertions with their right arm placed in a rigid adjustable brace affixed to a tri-axial load cell. Maximum exertions were performed in 24 directions, including four primary directions, horizontal flexion-extension, abduction-adduction, and at 15° increments in between those axes. Moments were computed and comparisons made between the experimentally collected strengths and those predicted by the IAA and WAA methods. The IAA over-predicted strength in 14 of 20 non-primary exertions directions, while the WAA underpredicted strength in only 2 of these directions. Therefore, it is not valid to assume that shoulder axes are independent when predicting shoulder strengths between two orthopedic axes, and the WAA is an improvement over current methods for the posture tested. PMID:26145486
NASA Astrophysics Data System (ADS)
Rey, Michaël; Nikitin, Andrei V.; Tyuterev, Vladimir G.
2015-10-01
Accurate basis set convergence of first-principles predictions of rotationally resolved spectra at high energy range is a common challenging issue for variational methods. In this paper, a detailed convergence study for the methane spectra is presented both for vibrational and rotational degrees of freedom as well as for intensities. For this purpose, we use our previously reported nine-dimensional potential energy and dipole moment surfaces of the methane molecule [Nikitin et al. Chem Phys Lett 2011;501:179-86; 2013;565:5-11]. Vibration-rotation calculations were carried out using variational normal mode approach with a full account of the Td symmetry. The aim was to obtain accurate theoretical methane line lists for the wavenumber range beyond currently available spectra analyses. The focus of this paper is the complicated icosad range (6280-7900 cm-1) containing 20 bands and 134 sub-bands where over 90% of experimental lines still remain unassigned. We provide variational line lists converged to a "spectroscopic precision" for icosad transitions for T=80 K and T=296 K. The first one contains 70 300 lines and the second one 286 170 lines with the intensity cut-off 10-29cm-1 /(moleculecm-2) with Jmax=18. An average error in line positions of theoretical predictions up to J=15 is estimated as 0.2-0.5 cm-1 from the comparisons with currently analyzed bands. Ab initio line strength calculations give the integrated intensity 4.37 ×10-20cm-1 /(moleculecm-2) at T=80 K for the sum of all icosad bands. This is to be compared to the integrated intensity 4.36 ×10-20cm-1 /(moleculecm-2) of the experimental icosad line list recorded in Grenoble University [Campargue et al., J Mol Spectrosc 2013;291:16-22] using very sensitive laser techniques. The shapes of absorption bands are also in a good qualitative agreement with experimental spectra.
The Calculation of Theoretical Chromospheric Models and the Interpretation of the Solar Spectrum
NASA Technical Reports Server (NTRS)
Avrett, Eugene H.
1998-01-01
Since the early 1970s we have been developing the extensive computer programs needed to construct models of the solar atmosphere and to calculate detailed spectra for use in the interpretation of solar observations. This research involves two major related efforts: work by Avrett and Loeser on the Pandora computer program for non-LTE modeling of the solar atmosphere including a wide range of physical processes, and work by Rurucz on the detailed synthesis of the solar spectrum based on opacity data or over 58 million atomic and molecular lines. our goals are: to determine models of the various features observed on the Sun (sunspots, different components of quiet and active regions, and flares) by means of physically realistic models, and to calculate detailed spectra at all wavelengths that match observations of those features. These two goals are interrelated: discrepancies between calculated and observed spectra are used to determine improvements in the structure of the models, and in the detailed physical processes used in both the model calculations and the spectrum calculations. The atmospheric models obtained in this way provide not only the depth variation of various atmospheric parameters, but also a description of the internal physical processes that are responsible for non-radiative heating, and for solar activity in general.
Automated Routines for Calculating Whole-Stream Metabolism: Theoretical Background and User's Guide
Bales, Jerad D.; Nardi, Mark R.
2007-01-01
In order to standardize methods and facilitate rapid calculation and archival of stream-metabolism variables, the Stream Metabolism Program was developed to calculate gross primary production, net ecosystem production, respiration, and selected other variables from continuous measurements of dissolved-oxygen concentration, water temperature, and other user-supplied information. Methods for calculating metabolism from continuous measurements of dissolved-oxygen concentration and water temperature are fairly well known, but a standard set of procedures and computation software for all aspects of the calculations were not available previously. The Stream Metabolism Program addresses this deficiency with a stand-alone executable computer program written in Visual Basic.NET?, which runs in the Microsoft Windows? environment. All equations and assumptions used in the development of the software are documented in this report. Detailed guidance on application of the software is presented, along with a summary of the data required to use the software. Data from either a single station or paired (upstream, downstream) stations can be used with the software to calculate metabolism variables.
Theoretical calculation of the p-emitter length for snapback-free reverse-conducting IGBT
NASA Astrophysics Data System (ADS)
Liheng, Zhu; Xingbi, Chen
2014-06-01
A physically based equation for predicting required p-emitter length of a snapback-free reverse-conducting insulated gate bipolar transistor (RC-IGBT) with field-stop structure is proposed. The n-buffer resistances above the p-emitter region with anode geometries of linear strip, circular and annular type are calculated, and based on this, the minimum p-emitter lengths of those three geometries are given and verified by simulation. It is found that good agreement was achieved between the numerical calculation and simulation results. Moreover, the calculation results show that the annular case needs the shortest p-emitter length for RC-IGBT to be snapback-free.
Theoretical Calculation for the Ionization of Molecules by Short Strong Laser Pulses
Nagy, L.; Borbely, S.
2011-10-03
We have developed several calculation methods for the ionization of atoms and molecules by strong and ultrashort laser pulses, based on the numerical solution of the time dependent Schroedinger equation (TDSE) in the momentum space. We have performed calculations within the strong field approximation (Volkov) and using iterative and direct methods for solving the TDSE. The investigated molecules are H{sub 2}{sup +} and H{sub 2}O. In case of the ionization of diatomic molecules the interference effects in the ejected electron spectra due to the coherent addition of the waves associated to the electrons ejected from the vicinity of different nuclei were also analysed.
Vibrational spectra, structure, and theoretical calculations of 2-fluoro- and 3-fluoropyridine.
Boopalachandran, Praveenkumar; Laane, Jaan
2011-09-01
The infrared and Raman spectra of liquid and vapor-phase 2-fluoropyridine and 3-fluoropyridine have been recorded and assigned. Ab initio and DFT calculations were carried out to compute the molecular structures and to verify the vibrational assignments. The observed and calculated spectra agree extremely well. The ring bond distances of the fluoropyridines are very similar to those of pyridine except for a shortening of the C-N(F) bond in 2-fluoropyridine. The C-F bond stretching frequencies are similar to that in fluorobenzene reflecting the influence of the ring π bonding. PMID:21570895
Technology Transfer Automated Retrieval System (TEKTRAN)
Determination of absolute configuration (AC) is one of the most challenging features in the structure elucidation of chiral natural products, especially those with complex structures. With revolutionary advancements in the area of quantum chemical calculations of chiroptical spectroscopy over the pa...
NASA Astrophysics Data System (ADS)
Menapace, E.; Birattari, C.; Bonardi, M. L.; Groppi, F.; Morzenti, S.; Zona, C.
2005-05-01
The radionuclide production for biomedical applications has been brought up in the years, as a special nuclear application, at INFN LASA Laboratory, particularly in co-operation with the JRC-Ispra of EC. Mainly scientific aspects concerning radiation detection and the relevant instruments, the measurements of excitation functions of the involved nuclear reactions, the requested radiochemistry studies and further applications have been investigated. On the side of the nuclear data evaluations, based on nuclear model calculations and critically selected experimental data, the appropriate competence has been developed at ENEA Division for Advanced Physics Technologies. A series of high specific activity accelerator-produced radionuclides in no-carrier-added (NCA) form, for uses in metabolic radiotherapy and for PET radiodiagnostics, are investigated. In this work, last revised measurements and model calculations are reviewed for excitation functions of natZn(d,X)64Cu, 66Ga reactions, referring to irradiation experiments at K=38 variable energy Cyclotron of JRC-Ispra. Concerning the reaction data for producing 186gRe and 211At/211gPo (including significant emission spectra) and 210At, most recent and critically selected experimental results are considered and discussed in comparison with model calculations paying special care to pre-equilibrium effects estimate and to the appropriate overall parameterization. Model calculations are presented for 226Ra(p,2n)225Ac reaction, according to the working program of the ongoing IAEA CRP on the matter.
Menapace, E.; Birattari, C.; Bonardi, M.L.; Groppi, F.; Morzenti, S.; Zona, C.
2005-05-24
The radionuclide production for biomedical applications has been brought up in the years, as a special nuclear application, at INFN LASA Laboratory, particularly in co-operation with the JRC-Ispra of EC. Mainly scientific aspects concerning radiation detection and the relevant instruments, the measurements of excitation functions of the involved nuclear reactions, the requested radiochemistry studies and further applications have been investigated. On the side of the nuclear data evaluations, based on nuclear model calculations and critically selected experimental data, the appropriate competence has been developed at ENEA Division for Advanced Physics Technologies. A series of high specific activity accelerator-produced radionuclides in no-carrier-added (NCA) form, for uses in metabolic radiotherapy and for PET radiodiagnostics, are investigated. In this work, last revised measurements and model calculations are reviewed for excitation functions of natZn(d,X)64Cu, 66Ga reactions, referring to irradiation experiments at K=38 variable energy Cyclotron of JRC-Ispra. Concerning the reaction data for producing 186gRe and 211At/211gPo (including significant emission spectra) and 210At, most recent and critically selected experimental results are considered and discussed in comparison with model calculations paying special care to pre-equilibrium effects estimate and to the appropriate overall parameterization. Model calculations are presented for 226Ra(p,2n)225Ac reaction, according to the working program of the ongoing IAEA CRP on the matter.
Technology Transfer Automated Retrieval System (TEKTRAN)
Several flux-calculation (FC) schemes are available for determining soil-to-atmosphere emissions of nitrous oxide (N2O) and other trace gases using data from non-steady-state flux chambers. Recently developed methods claim to provide more accuracy in estimating the true pre-deployment flux (f0) comp...
Theoretical calculation of low-lying states of NaAr and NaXe
NASA Technical Reports Server (NTRS)
Laskowski, B. C.; Langhoff, S. R.; Stallcop, J. R.
1981-01-01
Potential curves as well as dipole moments and linking transition moments are calculated for the ground X 2 Sigma + and low lying excited A 2 Pi, B 2 Sigma +, C 2 Sigma +, (4) 2 Sigma +, (2) 2 Pi and (1) 2 Delta states of NaAr and NaXe. Calculations are performed using a self-consistent field plus configuration-interaction procedure with the core electrons replaced by an ab initio effective core potential. The potential curves obtained are found to be considerably less repulsive than the semiempirical curves of Pascale and Vandeplanque (1974) and to agree well with existing experimental data, although the binding energies of those states having potential minima due to van der Waals interactions are underestimated. Emission bands are also calculated for the X 2 Sigma + - C 2 Sigma + excimer transitions of NaAr and NaXe using the calculated transition moments and potential curves, and shown to agree well with experiment on the short-wavelength side of the maximum.
Calculations of band gaps in polyaniline from theoretical studies of oligomers
Kwon, O.; McKee, M.L.
2000-03-02
Geometries and band gaps of polyaniline oligomers up to decamer have been systematically calculated and analyzed using various computational techniques such as molecular mechanics, semiempirical, and ab initio methods. On the basis of fully optimized geometries of neural and charged forms of polyaniline oligomers, excitation energies are calculated at the semiempirical ZINDO (INDO/S) level and extrapolated to the band gap value of the infinite chain. Band gaps are also approximately by extrapolating the HOMO/LUMO difference calculated at the density functional level (B3LYP/6--31G*). The SINDO//AM1 band gaps in the reduced and oxidized form of polyaniline (4.3 and 2.7 eV) are in good agreement with experimental values (3.8 {+-} 2 and 1.8 {+-} 3 eV, respectively). The doped form of polyaniline (two positive charges per four aniline units) has been computed with a spin-unrestricted method (UAM1) and the band gap approximated from an extrapolation of the tetramer and octamer. The calculated band gap of 1.3 eV (UZINDO//UAM1) is in good agreement with experiment (1.5 eV). The influence of ring torsional angle and interchain interaction on the band gap of the polyaniline system are also discussed.
The three-fold theoretical basis of the Gravity Probe B gyro precession calculation
NASA Astrophysics Data System (ADS)
Adler, Ronald J.
2015-11-01
The Gravity Probe B (GP-B) experiment is complete and the results are in agreement with the predictions of general relativity (GR) for both the geodetic precession, 6.6 arcsec yr-1 to about 0.3%, and the Lense-Thirring precession, 39 marcsec to about 19%. This note is concerned with the theoretical basis for the predictions. The predictions depend on three elements of gravity theory, firstly that macroscopic gravity is described by a metric theory such as GR, secondly that the Lense-Thirring metric provides an approximate description of the gravitational field of the spinning Earth, and thirdly that the spin axis of a gyroscope is parallel displaced in spacetime, which gives its equation of motion. We look at each of these three elements to show how each is solidly based on previous experiments and well-tested theory. The agreement of GP-B with theory strengthens our belief that all three elements are correct and increases our confidence in applying GR to astrophysical phenomena. Conversely, if GP-B had not verified the predictions a major theoretical quandary would have occurred.
NASA Astrophysics Data System (ADS)
Latgé, A.; Porras-Montenegro, N.; de Dios-Leyva, M.; Oliveira, L. E.
1997-05-01
The acceptor-related photoluminescence of a GaAs-(Ga,Al)As superlattice, under the influence of a magnetic field applied parallel to the interfaces, is theoretically studied following a variational procedure within the effective-mass approximation. Electron and hole magnetic Landau levels and envelope wave functions were obtained by an expansion in terms of sine functions, whereas for the impurity levels the envelope functions were taken as products of sine and hydrogenic-like variational functions. Impurity binding energies and wave functions are obtained for acceptors at a general position in the superlattice and for different in-plane magnetic fields. Theoretical results corresponding to transitions from the conduction subband to states of acceptors (miniband-to-acceptor e-A0 transitions) at the edge and center positions of the GaAs quantum well compare well with available experimental data by Skromme et al. [Phys. Rev. Lett. 65, 2050 (1990)] on the magnetic-field dependence of the photoluminescence peak position of conduction miniband-to-acceptor transitions for different temperatures and values of the superlattice period.
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
Millimeterwave rotational spectrum and theoretical calculations of cis-propionic acid
NASA Astrophysics Data System (ADS)
Jaman, A. I.; Chakraborty, Shamik; Chakraborty, Rangana
2015-01-01
The millimeterwave rotational spectra of the cis conformer of propionic acid (C3H6O2) have been investigated in the ground vibrational state in the frequency range of 80.0-100.0 GHz. Many high J and K-1 (Jmax = 50, K-1 = 12) rotational lines have been assigned. A least-squares analysis of the measured and previously reported rotational transition frequencies resulted in the determination of an improved set of rotational and centrifugal distortion (CD) constants of the molecule. Detailed MP2 and DFT calculations were also carried out with various functional and basis sets to evaluate the spectroscopic constants, dipole moment, and various structural parameters of cis-propionic acid and compared with the corresponding experimental values. Potential energy surface has been calculated to identify other probable conformers in this molecule.
Theoretical Calculations of Refractive Properties for Hg3Te2Cl2 Crystals
NASA Astrophysics Data System (ADS)
Bokotey, O. V.
2016-05-01
This paper reviews the optical properties, such as refractive index, optical dielectric constant, and reflection coefficient of the Hg3Te2Cl2 crystals. The applications of the Hg3X2Y2 crystals as electronic, optical, and optoelectronic devices are very much determined by the nature and magnitude of these fundamental material properties. The origin of chemical bonding in the crystals is very important for definition of the physical and chemical properties. The main structural feature of the Hg3X2Y2 crystals is the presence of covalent pyramids [XHg3] and linear X-Hg-X groups. Optical properties are calculated according to the model proposed by Harrison. The refractive index in the spectral region far from the absorption edge is determined within the generalized single-oscillator model. The calculated results are found to be in good agreement with experimental data.
Theoretical Calculations of Refractive Properties for Hg3Te2Cl2 Crystals.
Bokotey, O V
2016-12-01
This paper reviews the optical properties, such as refractive index, optical dielectric constant, and reflection coefficient of the Hg3Te2Cl2 crystals. The applications of the Hg3X2Y2 crystals as electronic, optical, and optoelectronic devices are very much determined by the nature and magnitude of these fundamental material properties. The origin of chemical bonding in the crystals is very important for definition of the physical and chemical properties. The main structural feature of the Hg3X2Y2 crystals is the presence of covalent pyramids [XHg3] and linear X-Hg-X groups. Optical properties are calculated according to the model proposed by Harrison. The refractive index in the spectral region far from the absorption edge is determined within the generalized single-oscillator model. The calculated results are found to be in good agreement with experimental data. PMID:27184964
NASA Astrophysics Data System (ADS)
Ghazzali, Mohamed; Khattab, Sherine A. N.; Elnakady, Yasser A.; Al-Mekhlafi, Fahd A.; Al-Farhan, Khalid; El-Faham, Ayman
2013-08-01
A series of naphthyl and tolyl sulfonate ester were synthesized and characterized by H NMR. X-ray single crystal diffraction experiments established the molecular structure of three new sulfonate esters derivatives, and spectral data agree with these in solution. The observed hydrogen bonding is discussed on the basis of crystal structural analyses and DFT/MP2 geometry optimization quantum calculations. Antimicrobial activities were screened for selected compounds against three human cancer cell lines and Mosquito Culex pipiens larvae.
Shen, Ming; Roopchand, Rabia; Amoureux, Jean-Paul; Chen, Qun
2015-01-01
Quadrupolar echo NMR spectroscopy of solids often requires RF pulse excitation that covers spectral widths exceeding 100 kHz. In a recent work we found out that a four pulse, composite pulse COM-II ( 90180¯90135¯45 ), provided robust broadband excitation for deuterium quadrupolar echo spectroscopy. Moreover, when combined with an eight step phase cycle, spectral distortions arising from finite pulse widths were greatly supressed. In this paper we report on a theoretical analysis COM-II with 8-step phase cycle by average Hamiltonian theory. This treatment is combined with the fictitious spin-1 operator formalism, and the mechanism of the 8-step phase cycling that minimizes the spectral distortions is discussed. PMID:26681896
Silva, Mateus X; Galvão, Breno R L; Belchior, Jadson C
2014-05-21
Genetic algorithm is employed to survey an empirical potential energy surface for small Na(x)K(y) clusters with x + y ≤ 15, providing initial conditions for electronic structure methods. The minima of such empirical potential are assessed and corrected using high level ab initio methods such as CCSD(T), CR-CCSD(T)-L and MP2, and benchmark results are obtained for specific cases. The results are the first calculations for such small alloy clusters and may serve as a reference for further studies. The validity and choice of a proper functional and basis set for DFT calculations are then explored using the benchmark data, where it was found that the usual DFT approach may fail to provide the correct qualitative result for specific systems. The best general agreement to the benchmark calculations is achieved with def2-TZVPP basis set with SVWN5 functional, although the LANL2DZ basis set (with effective core potential) and SVWN5 functional provided the most cost-effective results. PMID:24691391
Theoretical calculations of ion acceleration in the vicinity of comet Giacobini-Zinner
NASA Technical Reports Server (NTRS)
Mckenzie, M. L.; Cravens, T. E.; Ye, G.
1994-01-01
Ionization of cometary neutral molecules produces ions which are picked up by the solar wind. The cometary ion pickup process for comet Giacobini-Zinner is studied in two ways: (1) with a test particle method in which trajectories are numerically calculated for several thousand ions whose initial locations were chosen randomly with probability proportional to the neutral density and (2) with a quasi-linear diffusion model. The cometary ion distribution function was calculated with the test particle model at several locations upstream of the bow shock and for several types of magnetic fluctuations (or waves). These waves were allowed to propagate in both directions along the magnetic field at the Alfven speed. Both pitch angle scattering and energy diffusion are evident in the derived ion distributions. The monochromatic waves result in less ion acceleration than turbulent fluctuations with about the same amplitude. The calculated ion distribution functions are in reasonable agreement with the distributions measured in the vicinity of comet Giacobini-Zinner in 1985 by particle detectors on the ICE spacecraft when the ratio of power in sunward propagating Alfven waves to the power in antisunward propagating waves is assumed to lie between about 20% and 50%. However, the quasi-linear diffusion model results agree best with the measured distribution functions when the power ratio is only about 3 percent.
Theoretical calculations on structural and electronic properties of BGaAsBi alloys
NASA Astrophysics Data System (ADS)
Aslan, Metin; Yalcin, Battal G.; Ustundag, Mehmet; Bagci, Sadik
2015-11-01
The structural and electronic properties of cubic B x Ga1- x As1- y Bi y alloys with bismuth (Bi) concentration of 0.0625, 0.125, 0.1875 and 0.25 are studied with various boron (B) compositions by means of density functional theory (DFT) within the Wu-Cohen (WC) exchange correlation potential based on generalized gradient approximation (GGA). For all studied alloy structures, we have implemented geometric optimization before the volume optimization calculations. The obtained equilibrium lattice constants and band gap of studied quaternary alloys are investigated for the first time in literature. While the lattice constant behavior changes linearly with boron concentration, increasing small amount of bismuth concentration alter the lattice constant nonlinearly. The present calculation shows that the band gap decreases with increasing bismuth concentration and direct band gap semiconductor alloy became an indirect band gap with increasing boron concentration. From the band offset calculation we have shown that increasing B and Bi concentration in host GaAs reduced the valance band offset in a heterostructure formed by GaAs and studied alloys.
Breiner, Joan; Gimeno, Benjamin S; Fenn, Mark
2007-01-01
Edaphic, foliar, and hydrologic forest nutrient status indicators from 15 mixed conifer forest stands in the Sierra Nevada, San Gabriel Mountains, and San Bernardino National Forest were used to estimate empirical or theoretical critical loads (CL) for nitrogen (N) as a nutrient. Soil acidification response to N deposition was also evaluated. Robust empirical relationships were found relating N deposition to plant N uptake (N in foliage), N fertility (litter C/N ratio), and soil acidification. However, no consistent empirical CL were obtained when the thresholds for parameters indicative of N excess from other types of ecosystems were used. Similarly, the highest theoretical CL for nutrient N calculated using the simple mass balance steady state model (estimates ranging from 1.4-8.8 kg N/ha/year) was approximately two times lower than the empirical observations. Further research is needed to derive the thresholds for indicators associated with the impairment of these mixed conifer forests exposed to chronic N deposition within a Mediterranean climate. Further development or parameterization of models for the calculation of theoretical critical loads suitable for these ecosystems will also be an important aspect of future critical loads research. PMID:17450298
Theoretical Model for the Calculation of Optical Properties of Gold Nanoparticles
NASA Astrophysics Data System (ADS)
Mendoza-García, A.; Romero-Depablos, A.; Ortega, M. A.; Paz, J. L.; Echevarría, L.
We have developed an analytical method to describe the optical properties of nanoparticles, whose results are in agreement with the observed experimental behavior according to the size of the nanoparticle under analysis. Our considerations to describe plasmonic absorption and dispersion are based on the combination of the two-level molecular system and the two-dimensional quantum box models. Employing the optical stochastic Bloch equations, we have determined the system's coherence, from which we have calculated expressions for the absorption coefficient and refractive index. The innovation of this methodology is that it allows us to take into account the solvent environment, which induce quantum effects not considered by classical treatments.
NASA Technical Reports Server (NTRS)
Avrett, E. H.
1985-01-01
Solar chromospheric models are described. The models included are based on the observed spectrum, and on the assumption of hydrostatic equilibrium. The calculations depend on realistic solutions of the radiative transfer and statistical equilibrium equations for optically thick lines and continua, and on including the effects of large numbers of lines throughout the spectrum. Although spectroheliograms show that the structure of the chromosphere is highly complex, one-dimensional models of particular features are reasonably successful in matching observed spectra. Such models were applied to the interpretation of chromospheric observations.
NASA Astrophysics Data System (ADS)
Arroudj, S.; Bouchouit, M.; Bouchouit, K.; Bouraiou, A.; Messaadia, L.; Kulyk, B.; Figa, V.; Bouacida, S.; Sofiani, Z.; Taboukhat, S.
2016-06-01
This paper explores the synthesis, structure characterization and optical properties of two new schiff bases. These compounds were obtained by condensation of o-tolidine with salicylaldehyde and cinnamaldehyde. The obtained ligands were characterized by UV, 1H and NMR. Their third-order NLO properties were measured using the third harmonic generation technique on thin films at 1064 nm. The electric dipole moment (μ), the polarizability (α) and the first hyperpolarizability (β) were calculated using the density functional B3LYP method with the lanl2dz basis set. For the results, the title compound shows nonzero β value revealing second order NLO behaviour.
NASA Astrophysics Data System (ADS)
Zhidong, Zhang; Zhiguang, Li; Jingli, Liu
The temperature-dependent pitch of a cholesteric phase is studied using both molecular field theory and the two-particle cluster theory. The interacting chiral molecules (as derived by van der Meer et al.) are placed at the sites of a three-dimensional, simple cubic lattice with orientations confined to two dimensions. The equilibrium pitch as functions of temperature is calculated and numerical results are compared with those predicted by Monte Carlo computer simulation. The two-particle cluster theory, taking into account short-range correlations between molecules, yields improved values compared with molecular field theory.
Munson, D.E.; Holcomb, D.J.; DeVries, K.L.; Brodsky, N.S.
1994-12-31
Cross-hole ultrasonic measurements were made in the immediate wall of the Air Intake Shaft of the Waste Isolation Pilot Plant facility. These measurements show that compressional wave speed markedly decreases at the shaft wall and then increases with radial distance from the shaft to eventually become that of solid or undamaged salt. This behavior is indicative of deformation damage or microfractures in the salt. These in situ data are compared to both laboratory measurements of wave speed as a function of volume dilatancy and to calculations based on the Multimechanism Deformation Coupled Fracture model, with reasonable agreement.
NASA Astrophysics Data System (ADS)
Kostko, Oleg; Troy, Tyler P.; Bandyopadhyay, Biswajit; Ahmed, Musahid
2015-03-01
Acetaldehyde, a probable human carcinogen and of environmental importance, upon solvation provides a test bed for understanding proton transfer pathways and catalytic mechanisms. In this study, we report on single photon vacuum ultraviolet photoionization of small acetaldehyde and acetaldehyde-water clusters. Appearance energies of protonated clusters are extracted from the experimental photoionization efficiency curves and compared to electronic structure calculations. The comparison of experimental data to computational results provides mechanistic insight into the fragmentation mechanisms of the observed mass spectra. Using deuterated water for isotopic tagging, we observe that proton transfer is mediated via acetaldehyde and not water in protonated acetaldehyde-water clusters.